#cls
text(In-Form Equivalent Of PLANT Case Z126
title
libref=726
DISPLAY
This In-Form case does what Z126 does, but without fortran coding
or re-compilation.
First, case z126 is loaded. Then NAMSAT = NONE is set, to
deactivate.
It uses SOURCE In-Form statements.
The Q1 contains PHOTON USE commands
ENDDIS
load(z126)
text(In-Form equivalent of PLANT case Z126
libref=726
INFORM13BEGIN
** The next line de-activates PLANT **
namsat=none
SVEL=SKIP; INT1=SKIP; INT2=SKIP
LSWEEP=300
** In-Form formulae which are equivalent
to the PLANT formulae of z126 **
PATCH(ISVEL,CELL,1,NX,1,NY,1,NZ,1,LSTEP)
(SOURCE of SLIP at ISVEL is SQRT((V1-V2)^2+(W1-W2)^2) with FIXV)
PATCH(IINT1,PHASEM,1,NX,1,NY,1,NZ,1,LSTEP)
(SOURCE of H1 at IINT1 is 0.05*sqrt(SLIP)*(H2-H1))
PATCH(IINT2,PHASEM,1,NX,1,NY,1,NZ,1,LSTEP)
(SOURCE of H2 at IINT2 is 0.05*sqrt(SLIP)*(H1-H2))
INFORM13END
DISTIL=T
EX(P1 )=1.361E+02; EX(V1 )=2.559E-02; EX(V2 )=2.480E-01
EX(W1 )=1.022E+00; EX(W2 )=1.044E+00; EX(R1 )=9.717E-01
EX(R2 )=2.828E-02; EX(H1 )=3.536E-01; EX(H2 )=6.799E-01
EX(CMDO)=8.340E+00; EX(W2CR)=8.668E-01; EX(V2CR)=4.656E-01
EX(SLIP)=2.195E-01; EX(WCRT)=8.614E-01; EX(VCRT)=3.030E-01
PHOTON USE
p
msg SOIL DRIFT POLLUTION IN A RIVER
gr ou x 1
msg Velocity vectors for the water
vec x 1 sh
msg Press Enter to continue
pause
vec off;red
set vec comp - v2cr w2cr
msg Velocity vectors for the soil particles
vec x 1 sh
msg Press Enter to continue
pause
vec off;red
msg Slip velocity
con slip x 1 fil;0.01
msg Press Enter to continue
pause
con off;red
msg Volume fraction of soil particles
con r2 x 1 fil;0.01
msg Press Enter to continue
pause
con off;red
msg Sedimentation mass transfer rate
con cmdo x 1 fil;0.01
msg Press Enter to continue
pause
con off;red
msg Water temperature
con h1 x 1 fil;0.01
msg Press Enter to continue
pause
con off;red
msg Soil temperature
con h2 x 1 fil;0.01
msg Press e to END
ENDUSE
STOP