TALK=F;RUN( 1, 1)

TEXT(Library case Y603: ADAPTIVE Z-DIRECTION STEP SIZE

  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    PLANT information :
     * Data input groups used:  5
     * Ground groups planted :  5
     * Headings used  : SCZS??
     * Functions used : None
     * Commands used  : IF
  <<<<<<<<<<<<<<<<<<<<<<  Comment ends <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
 LOAD(150)
    GROUP 5. Z-direction grid specification
NAMSAT=MOSG
AZDZ=GRND
   DZ=0.1*YVLAST
  IF(IZ.LE.2)
   DZ=DZL*(1.-0.02*TEMP[,1,-1]/:TJET:)
  IF(IZ.GT.2)
  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
   The above  two  PLANT blocks overwrite the correspondong actions in
   library PIL.

   They make the step size,  DZ,  a fixed fraction of the  calculation
   domain width,  YVLAST,  for the steps less or equal 2.  If the step
   number is greater than 2,  the  each  next  step  size  is  made  a
   fraction of the thickness for previous z-slab, DZL.

   The fraction  factor  is  supposed  to  be adaptive to the upstream
   temperature at the first cell in Y-direction.  The  consequence  is
   that the step size gradually reduces towards the domain end.
  <<<<<<<<<<<<<<<<<<<<<<  Comment ends <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
  PHOTON USE
  p;  parphi;;;

  msg        PLANE JET
  msg
  msg        ADAPTIVE Z-DIRECTION STEP SIZE
  msg        Planted non-uniform grid  and vectors
  msg            -
  msg Press  to continue
  gr  x 1;vec x 1 sh
  msg Press e to END
  enduse
STOP