subroutine SetSubgridLength(ixlo, ixhi, iylo, iyhi, izlo, izhi, $ nx,ny,nz,dx,dy,dz,len) IMPLICIT NONE integer ixlo, ixhi, iylo, iyhi, izlo, izhi DOUBLE PRECISION len(ixlo:ixhi,iylo:iyhi,izlo:izhi) INTEGER nx,ny,nz DOUBLE PRECISION dx,dy,dz integer i,j,k do k=izlo,izhi do j=iylo, iyhi do i=ixlo, ixhi len(i,j,k) = 0.0d0 enddo enddo enddo return end