TALK=T;RUN(1,1)
  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
   PLANTBEGIN
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
         ** Ask PLANT to introduce first time step, 5 seconds.
       DT=5.
      REGION(1,1,1,1,1,1,1,1)
    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)
    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
    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)
    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)
    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
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(TIME STEP CALCULATIONS                  )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 619
 ************************************************************
  Group 2. Time dependence
 STEADY = F
    * Set overall time and no. of steps
 TFIRST =0. ;TLAST = GRND
 FSTEP = 1 ;LSTEP = 2
 TFRAC(1)=1.
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 2
 XULAST =2.
 XFRAC(1)=0.5 ;XFRAC(2)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 2
 YVLAST =4.
 YFRAC(1)=0.5 ;YFRAC(2)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 2
 ZWLAST =6.
 ZFRAC(1)=0.5 ;ZFRAC(2)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(3)=U1 ;NAME(5)=V1
 NAME(7)=W1 ;NAME(14)=H1
    * Y in SOLUTN argument list denotes:
    * 1-stored 2-solved 3-whole-field
    * 4-point-by-point 5-explicit 6-harmonic averaging 
 SOLUTN(U1,Y,N,N,N,N,Y)
 SOLUTN(V1,Y,N,N,N,N,Y)
 SOLUTN(W1,Y,N,N,N,N,Y)
 SOLUTN(H1,Y,Y,N,N,N,Y)
 ************************************************************
  Group 8. Terms & Devices
    * Y in TERMS argument list denotes:
    * 1-built-in source 2-convection 3-diffusion 4-transient
    * 5-first phase variable 6-interphase transport         
 TERMS(H1,Y,Y,Y,Y,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 HUNIT =1.
 ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
 ************************************************************
  Group 9. Properties used if PRPS is not
  stored, and where PRPS = -1.0 if it is!
 RHO1 =1. ;TMP1 =0. ;EL1 =0.
 TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
 DVO1DT =0. ;DRH1DP =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 ENUL =1.0E-05 ;ENUT =0.
 PRNDTL(H1)=1.
 PRT(H1)=1.
 CP1 =1. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(U1)=1. ;FIINIT(V1)=0.5
 FIINIT(W1)=0.1 ;FIINIT(H1)=1.0E-10
 
 PATCH(MAXVEL ,INIVAL, 1, 1, 2, 2, 2, 2, 1, 1)
 INIT(MAXVEL ,U1 ,0. ,10. )
 INIADD = F
 FSWEEP = 1
 NAMFI =CHAM
 ************************************************************
  Group 12. Patchwise adjustment of terms
  Patches for this group are printed with those
  for Group 13.
  Their names begin either with GP12 or &
 ************************************************************
  Group 13. Boundary & Special Sources
   No PATCHes yet used for this Group
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 2 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-05
 ************************************************************
  Group 16. Terminate Iterations
 LITER(H1)=20
 ENDIT(H1)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(U1,LINRLX,1.)
 RELAX(V1,LINRLX,1.)
 RELAX(W1,LINRLX,1.)
 RELAX(H1,FALSDT,1.0E+09)
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
 VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
 VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
 VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 NAMSAT =MOSG
 PARSOL = F
 ISG62 = 1
 SPEDAT(SET,GXMONI,TRANSIENT,L,F)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 RG( 1) =1.0E+20
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(U1)=3.25 ;EX(V1)=0.5
 EX(W1)=0.1 ;EX(H1)=5.445E-04
 ************************************************************
  Group 21. Print-out of Variables
 INIFLD = F ;SUBWGR = F
    * Y in OUTPUT argument list denotes:
    * 1-field 2-correction-eq. monitor 3-selective dumping      
    * 4-whole-field residual 5-spot-value table 6-residual table
 OUTPUT(U1,Y,N,Y,N,N,N)
 OUTPUT(V1,Y,N,Y,N,N,N)
 OUTPUT(W1,Y,N,Y,N,N,N)
 OUTPUT(H1,Y,N,Y,Y,Y,Y)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 1 ;IYMON = 1 ;IZMON = 1
 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = 1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NTPRIN = 100000 ;ISTPRF = 1 ;ISTPRL = 100000
 NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
 NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
 NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
 XZPR = F ;YZPR = F
 IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
 ISWPRF = 1 ;ISWPRL = 100000
 ITABL = 3 ;IPROF = 1
 ABSIZ =0.5 ;ORSIZ =0.4
 NTZPRF = 1 ;NCOLPF = 50
 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
   No PATCHes yet used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP