#cls
text(In-Form Equivalent Of PLANT Case Z125
title
libref=725
DISPLAY
This In-Form case does what Z125 does, but without fortran coding
or re-compilation.
First, case z125 is loaded. Then NAMSAT = NONE is set, to
deactivate.
It uses PROPERTY In-Form statement.
The Q1 contains PHOTON USE commands
ENDDIS
PHOTON USE
P
rot z ang 180;rot x ang 90
gr ou x 1
set vec ref 2.
vec x 1 sh
msg( LN2 velocity distribution
msg( Press ENTER to continue
pause;
vec off;red
set vec comp - v2 w2
vec x 1 sh
msg( N2 velocity distribution
msg( Press ENTER to continue
pause;
vec off;red
con r2 x 1 fil;.0001
msg( N2 volume fraction contours
msg( Press E to finish
ENDUSE
load(z125)
text(In-Form equivalent of PLANT case Z125
libref=725
INFORM1BEGIN
** The next line de-activates PLANT **
namsat=none
INFORM1END
** The following settings have been chosen as being
preferable to those of z125, in respect of convergence
and printout. **
LSWEEP=750
RELAX(P1,LINRLX,0.95);RELAX(R1,LINRLX,0.95);RELAX(R2,LINRLX,0.95)
RELAX(V1,FALSDT,0.0001);RELAX(W1,FALSDT,0.0001)
RELAX(V2,FALSDT,0.0001);RELAX(W2,FALSDT,0.0001)
** In-Form property formulae which are equivalent to the
PLANT formulae of z125 **
INFORM9BEGIN
(PROPERTY CFIPS is 500.*R2*MASS1)
INFORM9END
DISTIL=T
EX(P1)=7.487E+02; EX(V1)=1.854E-01; EX(V2)=2.891E-01
EX(W1)=5.604E-01; EX(W2)=5.735E-01; EX(R1)=7.816E-01
EX(R2)=2.184E-01; EX(VPOR)=8.521E-01
store(prps)
EX(PRPS)=3.014E+01
STOP