#369
grdpwr(t,lstep,tlast,0.5) ! make time grid of power-law type
grdpwr(x,nx,-xulast,1.1) ! make x-grid of geometric-progression
! type
clr
display
This case is like 369; but both the x-wise and t-wise grids
are strongly non-uniform.
The motion is nevertheless imparted to the fluid, because
whenever the box finds itself in a cell which is larger than
itself, it "expands" so as to fill the cell.
This expansion is performed automatically within 'EARTH'.
ENDDIS
distil=t
store(rsu1,aep1)
ex(p1)=1.147E+00; ex(rsu1)=2.200E-07; ex(aep1)=2.095E-01
libref=374
stop