#cls
text(In-Form Equivalent Of PLANT Case Z131
title
libref=731
  
  DISPLAY
  
  This In-Form case does what Z131 does, but without fortran coding
  or re-compilation.
  
  First, case z131 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
  SET VEC REF 0.005;ve Z 1 sh
  msg                    Velocity distribution:
  msg Press  to continue
  pause;clear
  gr ou z 1
  msg                   Temperature distribution :
  con h1 z 1 fil;0.0004
  msg Press e to END
  ENDUSE

load(z131) 
text(In-Form equivalent of PLANT case Z131
libref=731

  INFORM13BEGIN
    ** The next line de-activates PLANT **
namsat=none
ERG2=SKIP; ERG4=SKIP
  
     ** In-Form formulae which are equivalent to the
        PLANT formulae of z131                      **
PATCH(POROS,INIVAL,1,NX,1,NY,1,NZ,1,1)
COVAL(POROS,EPOR,0.,POROS)
COVAL(POROS,NPOR,0.,POROS)
COVAL(POROS,VPOR,0.,POROS)

PATCH(iERG2,PHASEM,1,NX,1,NY,1,NZ,1,1000)
(SOURCE of U1 at iERG2 is :RG(2):*ABS(U1)*(0-U1))

PATCH(iERG4,PHASEM,1,NX,1,NY,1,NZ,1,1000)
(SOURCE of V1 at iERG4 is :RG(2):*ABS(V1)*(0-V1))
  INFORM13END
DISTIL=T
EX(P1)=8.714E-03; EX(U1)=7.371E-04; EX(V1)=7.447E-04
EX(H1)=5.021E+02; EX(NPOR)=5.000E-01; EX(VPOR)=4.644E-01
EX(EPOR)=5.000E-01; EX(RHO1)=6.643E-01
STORE(PRPS)
EX(PRPS)=1.501E+01
STOP