#cls
text(In-Form Equivalent Of PLANT Case Z127
title
libref=727
  
  DISPLAY
  
  This In-Form case does what Z127 does, but without fortran coding
  or re-compilation.
  
  First, case z127 is loaded. Then NAMSAT = NONE is set, to 
  deactivate.

  It uses SOURCE In-Form statements.

  The Q1 contains PHOTON USE commands
  ENDDIS  
  PHOTON USE
  p



  gr ou z 1;vec z 1
  msg                  VELOCITY VECTORS
  msg Press Enter to continue
  pause
  red
  con c1 z 1 sh;in 30
  msg         NEAR GROUND CONCENTRATION OF 40 mcm AEROSOL
  msg Press Enter to continue
  pause;con off
  red
  con c2 z 1 sh;in 30
  msg         NEAR GROUND CONCENTRATION OF 120 mcm AEROSOL
  msg Press E to end
  ENDUSE

load(z127)
text(In-Form equivalent of PLANT case Z127
libref=727

  INFORM13BEGIN
    ** The next line de-activates PLANT **
namsat=none
CONG1=SKIP; VELG2=SKIP; CONG2=SKIP
  INFORM13END
  
     ** The following settings have been chosen as being
        preferable to those of z127, in respect of convergence
        and printout.                                         **
LSWEEP=22
RELAX(C2,FALSDT,1000.)
RELAX(WG1,FALSDT,0.01);RELAX(WG2,FALSDT,0.01)

     ** In-Form formulae which are equivalent to the
        PLANT formulae of z127                      **
  INFORM13BEGIN
PATCH(iCONG1,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP)
(SOURCE of C1 at iCONG1 is WG1/HEIGHT*(0-C1))

PATCH(iVELG2,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP)
(SOURCE of WG2 at iVELG2 is (COFA2+COFA3*WG2)*$
(AGRAV/(COFA2+COFA3*WG2)-WG2))

PATCH(iCONG2,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP)
(SOURCE of C2 at iCONG2 is WG2/HEIGHT*(0-C2))
  INFORM13END
DISTIL=T
EX(U1)=3.546E+00; EX(V1)=3.546E+00; EX(C1)=5.839E-03
EX(C2)=2.454E-03; EX(WG1)=4.844E-02; EX(WG2)=1.731E-01
STOP