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
This archive was generated by hypermail 2b30 : Thu Sep 12 2002 - 04:25:21 CDT