From ff9920f375f6956f43ee4a00fd0a63fd03f5aecb Mon Sep 17 00:00:00 2001
From: Pedro Costa
Date: Tue, 29 Aug 2023 14:29:14 +0200
Subject: [PATCH] Fixed OOB issue in `out1d` when averaging along `x` or `y`.
(#90)
---
src/output.f90 | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/src/output.f90 b/src/output.f90
index 5558141d..27b05f5a 100644
--- a/src/output.f90
+++ b/src/output.f90
@@ -64,7 +64,8 @@ subroutine out1d(fname,ng,lo,hi,idir,l,dl,z_g,dz,p)
integer , intent(in), dimension(3) :: ng,lo,hi
integer , intent(in) :: idir
real(rp), intent(in), dimension(3) :: l,dl
- real(rp), intent(in), dimension(0:) :: z_g,dz
+ real(rp), intent(in), dimension(0: ) :: z_g
+ real(rp), intent(in), dimension(0:lo(3)-1) :: dz
real(rp), intent(in), dimension(lo(1)-1:,lo(2)-1:,lo(3)-1:) :: p
real(rp), allocatable, dimension(:) :: p1d
integer :: i,j,k