User: rhaas Date: 2014/04/29 07:52 PM
Modified: /trunk/src/ NullInterp_MaskInit.F90
Log: Fix check for guard cells
before it was using unitialized memory.
Fix provided by Bela Szilagyi
File Changes:
Directory: /trunk/src/ ======================
File [modified]: NullInterp_MaskInit.F90 Delta lines: +3 -3 =================================================================== --- trunk/src/NullInterp_MaskInit.F90 2014-04-25 19:58:03 UTC (rev 16) +++ trunk/src/NullInterp_MaskInit.F90 2014-04-30 00:52:32 UTC (rev 17) @@ -99,9 +99,9 @@ ! for guard points near the patch boundaries
if (((lbnd(1).eq.0).and.any(EG_mask(1:N_ang_stencil_size,:).ne.0)) .or. & - ((ubnd(1).eq.gsh(1)-1).and.any(guard_mask(lsh(1)-N_ang_stencil_size+1:,:).ne.0)) .or. & - ((lbnd(2).eq.0).and.any(guard_mask(:,1:N_ang_stencil_size).ne.0)) .or. & - ((ubnd(2).eq.gsh(2)-1).and.any(guard_mask(:,lsh(2)-N_ang_stencil_size+1:).ne.0))) & + ((ubnd(1).eq.gsh(1)-1).and.any(EG_mask(lsh(1)-N_ang_stencil_size+1:,:).ne.0)) .or. & + ((lbnd(2).eq.0).and.any(EG_mask(:,1:N_ang_stencil_size).ne.0)) .or. & + ((ubnd(2).eq.gsh(2)-1).and.any(EG_mask(:,lsh(2)-N_ang_stencil_size+1:).ne.0))) & call CCTK_WARN(0, "mask setup error");
! mark guard points