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

  It uses the In-Form statements:
      STORED and SOURCE.

  The Q1 contains PHOTON USE commands
  ENDDIS  
  PHOTON USE
  p



  gr ou z 1
  con fliq z 1 fi;0.001
  vec z 1
  msg Velocity vectors and liquid volume fraction.
  ENDUSE

load(z106) 
text(In-Form equivalent of PLANT case Z106
libref=706
  INFORM13BEGIN
    ** The next line de-activates PLANT **
namsat=none
MUSHY=SKIP
  INFORM13END

     ** The following settings have been chosen as being
        preferable to those of z106, in respect of convergence
        and printout.                                         **
SELREF=F
RELAX(U1,FALSDT,2.E-4); RELAX(V1,FALSDT,2.E-4)
  
     ** In-Form formulae which are equivalent to the
        PLANT formulae of z106                      **
  INFORM7BEGIN
(STORED VAR FLIQ IS MAX(1.e-5&MIN(1.&(TEMP-:RG(1):)/(:RG(2):$
-:RG(1):))))
  INFORM7END
  INFORM13BEGIN
PATCH (iMUSHY,VOLUME,1,NX,1,NY,1,NZ,1,1)
(SOURCE OF U1 AT iMUSHY IS (0-U1)*:RG(3):*(1-FLIQ)/FLIQ)
(SOURCE OF V1 AT iMUSHY IS (0-V1)*:RG(3):*(1-FLIQ)/FLIQ)
  INFORM13END
DISTIL=T
EX(P1)=2.877E+01; EX(U1)=8.874E-02; EX(V1)=4.652E-02
EX(FLIQ)=6.456E-01; EX(TEMP)=2.736E+02; EX(VISL)=1.000E-03
STOP