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
double precision sgs_len
sgs_len = 6.283185307/32
do k=izlo,izhi
do j=iylo, iyhi
do i=ixlo, ixhi
len(i,j,k) = sgs_len
enddo
enddo
enddo
return
end