>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
  DISPLAY
    PLANT is used for specifying two time-step  sizes.  
    
    The first is set equal to 5 sec.
    The second time step is computed by dividing the smallest cell 
    size by the largest velocity.
 
    The solved-for variable is H1; but, since the actual solution is  
    of not of interest, no initial or boundary conditions are
    introduced.
 
    The non-uniform grid and velocity field are  initialised
    for expected size of the second time step to be 0.1 sec.
    
    PLANT information :
     * Data input groups used:  2, 19
     * Ground groups planted :  1, 2, 19-3
     * Headings used  : SCTS??
     * Functions used : SUM
     * Commands used  : IF, REGION, TEXT.
  ENDDIS
 
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
 
    GROUP 1. Run title and other preliminaries
TEXT(TIME STEP CALCULATIONS
    GROUP 2. Transience; time-step specification
STEADY=F
LSTEP=2
TLAST=GRND
  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
    The above  three  lines will instruct EARTH to perform 2
    time steps,  the size of which being determined  by  the
    settings PLANTed in group 2.

   PLANTBEGIN
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
         ** Ask PLANT to introduce first time step, 5 seconds.
       DT=5.
      REGION(1,1,1,1,1,1,1,1)
  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
    REGION command  is  used to control this simple setting.
    Its first six arguments limit DO loop to perform  single
    action over first time step defined by the last argument
    pair.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
  ** Choose the smallest cell size at the start of each iz-slab
     for the last sweep of the first time step.
RG(1)=GREAT
              DXG2D - cell size in X-direction,
              DYG2D - cell size in Y-direction,
              DZWNZ - cell size in Z-direction,
       RG(1)=AMIN1(RG(1),DXG2D,DYG2D,DZWNZ)
      IF(ISTEP.EQ.1.AND.ISWEEP.EQ.LSWEEP)
  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
    Command IF restricts the whole-domain-default extents of
    PLANT statement.  The arguments of IF  command  instruct
    PLANT   to  introduce  the  logical  conditions  of  the
    bracketed expression.  The execution  of  the  statement
    returns   the   auxilliary  variable,  RG(1),  with  the
    smallest cell size over the whole domain.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
      
      ** Output of smallest cell size using summation in one
         cell at 1st time step, last sweep  for checking.
       SIZMIN=SUM(RG(1))
      TEXT(Smallest cell size)
      REGION(1,1,1,1,1,1,1,1) /ISWEEP.EQ.LSWEEP
  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
    The function SUM is used in above three  lines  just  to
    dump into globcalc file the value of RG(1) headed by the
    character arguments of TEXT command.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
         ** Choose the largest velocity value just before the
            second time step calculation
       RG(2)=AMAX1(RG(2),U1,V1,W1)
      IF(ISTEP.EQ.2)
         ** Output the largest velocity using summation in one
            cell at 2nd time step for checking.
       VELMAX=SUM(RG(2))
      TEXT(Largest velocity)
      REGION(1,1,1,1,1,1,2,2)
  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
    Auxiliary variable,  RG(2), will be the largest velocity
    value  as  result  of the above statements at the second
    time moment.  It will be dumped  in  globcalc  file  and
    appropriately headed.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
       DT=SIZMIN/VELMAX
      REGION(1,1,1,1,1,1,2,2)
  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
    The final stage of time step settings:  second time step
    size, DT, is set as division of smallest size by largest
    velocity.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
   PLANTEND

    ** The simple non-uniform grid with smallest cell size of 1 m.
    GROUP 3. X-direction grid specification
GRDPWR(X,2,2.,1.)
    GROUP 4. Y-direction grid specification
GRDPWR(Y,2,4.,1.)
    GROUP 5. Z-direction grid specification
GRDPWR(Z,2,6.,1.)
    GROUP 7. Variables stored, solved & named
SOLVE(H1)
STORE(U1,V1,W1)
    GROUP 11. Initialization of variable or porosity fields
    ** The simple non-uniform velocity field with maximum
       value, 10.0m/s, at east-south-high corner of the domain.
INIADD=F
FIINIT(U1)=1.;FIINIT(V1)=0.5;FIINIT(W1)=0.1
PATCH(MAXVEL,INIVAL,NX-1,NX-1,NY,NY,NZ,NZ,1,1)
COVAL(MAXVEL,U1,0.0,10.)
    GROUP 15. Termination of sweeps
LSWEEP=2
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
    GROUP 19. Data communicated by satellite to GROUND
NAMSAT=MOSG
dmpstk=t
DISTIL=T
EX(U1)=3.250E+00; EX(V1)=5.000E-01
EX(W1)=1.000E-01; EX(H1)=5.445E-04
 LIBREF=619
STOP