construct_grid is does not handle no equidistant grids properly
Brought to you by:
tobson
Originally created by: wdobler@gmail.com
Originally owned by: wdobler@gmail.com
!FIXME: ! This dxmin/max concept is not compatible with the concept of a
spacially ! varying grid. !
dxmin = minval( (/dx,dy,dz,huge(dx)/), &
MASK=((/nxgrid,nygrid,nzgrid,2/) > 1) )
dxmax = maxval( (/dx,dy,dz,epsilon(dx)/), &
MASK=((/nxgrid,nygrid,nzgrid,2/) > 1) )