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

  The Q1 uses the In-Form statements:
      PROPERTY, INITIAL and SOURCE.

  The Q1 contains PHOTON USE commands
  ENDDIS  
  PHOTON USE
  p



  up z
  g ou ce x 1
  g ou ce x 1 y 1 5 z 9 12
  ve x 1 sh
  msg                    Velocity distribution:
  msg Press  to continue
  pause;ve off;red
  msg                   Temperature distribution :
  con temp x 1 fil;0.001
  g ou ce x 1 y 1 5 z 9 12 col 4
  msg Press e to END
  ENDUSE

load(z124) 
text(In-Form equivalent of PLANT case Z124
libref=724

  INFORM13BEGIN
    ** The next line de-activates PLANT **
namsat=none
TMP=SKIP; BODY=SKIP
PRNDTL(TEMP)=-GRND2
LSWEEP=500
RELAX(P1,LINRLX,0.5)
RELAX(V1,FALSDT,1.e-05);RELAX(W1,FALSDT,1.e-05)
RELAX(TEMP,FALSDT,1.e-05)

  INFORM13END
  
     ** In-Form formulae which are equivalent to the
        PLANT formulae of z124                      **
  INFORM9BEGIN
enul=1 ! cancel call to GROUND
(PROPERTY var VISL is VISC)

(PROPERTY var PRNDTL(TEMP) is DIFS)
  INFORM9END

  INFORM11BEGIN
PATCH(iTMP,INIVAL,1,1,1,NY,1,NZ,1,1)
(INITIAL of TEMP at iTMP is ZG/0.3)
  INFORM11END
     Buoyancy forces
  INFORM13BEGIN
PATCH(iBODY,VOLUME,1,1,1,NY,1,NZ,1,1)
(SOURCE of W1 at iBODY is 100.*5.4*TEMP)
  INFORM13END
DISTIL=T
EX(P1)=2.953E+02; EX(V1)=5.011E-05; EX(W1)=8.162E-05
EX(DIFS)=6.090E+01; EX(VISC)=1.000E+09; EX(TEMP)=5.000E-01
SOLUTN(P1,P,P,Y,P,P,P)
SOLUTN(TEMP,P,P,Y,P,P,P)
STOP