[FLASH-BUGS] Bug in MarkRefLib.F90

[Sebastian Heinz]

Dear FLASH team,

There is a small bug in MarkRefLib.F90:

lines 185 through 194 should read

!=====================================================================

        if (ndim .eq. 3) then
           if ( (coords(3)-zc) > 0 ) then
              call dBaseGetCoords(znl, zCoord, block_no, face_zcoords)
              zcell = nguard+1
           else if ( (coords(3)-zc) < 0 ) then
              call dBaseGetCoords(znr, zCoord, block_no, face_zcoords)
              zcell = nguard+nzb
           else 
              call dBaseGetCoords(znr, zCoord, block_no, face_zcoords)
              zcell = nguard+(nzb/2)
           endif
        endif

!=====================================================================



Best regards,

Sebastian