PHOTON USE
P;phi; 10 1;;
msg( Single fluid central-jet concentration contours
con h1 x 1 fil;.001
pause
con cl; red
msg( Averaged 17-fluid contours
con cav x 1 fil;.001
pause
con cl; red
msg( Concentration fluctuations by transport equation
con gg x 1 fil;.001
pause
con cl; red
msg( Averaged 17-fluid concentration fluctuation
con gav x 1 fil;.001
msg
msg Hit Enter for FPD hystogram
msg
pause
p
20 1
con fpd x 1 fil;10.
msg Hit Enter to continue
ENDUSE
>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
DISPLAY
In this case, PLANT is used to introduce a turbulence
model for concentration fluctuations which employs no
conservation equations for statistical properties of
the fluctuations. It is a variant of Multi-Fluid
concept of Brian Spalding as employed by Sergei Zhubrin.
The 17-fluid model is considered here to simulate the
turbulent mixing resulting from the admission of two
separate, isothermal coaxial jets of different
composition into a concentric duct as depicted
diagrammatically below.
ENDDIS
Entrance /////////////// Duct wall /////////////////Exit
-> :::::::::::: ->
surrounding -> :::::::::: -->
jet ------->:::::Mixing layer spreading to wall--->
" -------> ---->
central jet ------->-.-.-.-.- Symmetry axis .-.-.-.-.-.----->
In parallel with above model, the conservation equation
for the transport of the square of fluctuating
concentration component is also PLANTed for comparative
purposes.
At the final stage of computations the results are
processed to get the Fluid Population Distribution
histogram which may be viewed by PHOTON. The use file
for the latter is supplied.
The example employs the number of PLANT features:
- specific sources introduction,
- reference residuals calculations,
- intervention in calculation,
- processing the results and
- special print-out preparations,
PLANT information :
* Data input groups used: 13, 19
* Ground groups planted : 13, 19-6
* Headings used : SORC??, SC06??
* Functions used : SUM
* Commands used : IF, REGION
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
GROUP 1. Run title and other preliminaries
TEXT(CONFINED JET FLOW: 17 FLUID MFM
-------------------------------
REAL(HIN,GMIXL,CLEN,WIDTH,WIN1,WIN2,REYNO,WD2)
REAL(TKEIN1,EPIN1,TKEIN2,EPIN2)
INTEGER(IYJ);IYJ=3
REYNO=1.E6;WIDTH=0.3;HIN=1.
WIN1=10.;WIN2=2.0
GROUP 3. X-direction grid specification
CARTES=F;XULAST=0.1
GROUP 4. Y-direction grid specification
NY=15;WD2=0.5*WIDTH;GRDPWR(Y,NY,WD2,1.0)
GROUP 5. Z-direction grid specification
NZ=20;CLEN=20.*WD2;GRDPWR(Z,NZ,CLEN,1.0)
GROUP 7. Variables stored, solved & named
** H1 - single fluid concentration variable;
G - square of concentration fluctuation;
SOLVE(P1,W1,V1,H1,G)
** GENG - production for the square of concentration fluctuation
STORE(ENUT,LEN1,GEN1,EPKE,GENG)
SOLUTN(P1,Y,Y,Y,N,N,N)
TURMOD(KEMODL)
GROUP 8. Terms (in differential equations) & devices
TERMS(H1,N,Y,Y,Y,Y,Y)
TERMS( G,N,Y,Y,Y,Y,Y)
GROUP 9. Properties of the medium (or media)
RHO1=1.0;ENUL=WIN1*WIDTH/REYNO
PRT(H1)= 0.86;PRNDTL(H1)= 0.71
PRT(G) = 0.7 ;PRNDTL(G) = 0.7
GROUP 11. Initialization of variable or porosity fields
FIINIT(W1)=0.5*(WIN1+WIN2);FIINIT(H1)=HIN;FIINIT(LEN1)=0.1*YVLAST
FIINIT(ENUT)=0.01*WIN1*YVLAST
** TKEIN = 0.25*WIN1*WIN1*FRIC where FRIC=0.018 AT REYNO=1.E5
TKEIN1=0.25*WIN1*WIN1*0.018
TKEIN2=0.25*WIN2*WIN2*0.018
FIINIT(KE)=0.5*(TKEIN1+TKEIN2)
** EPIN = 0.1643*KIN**1.5/LMIX where LMIX=0.045*WIDTH
GMIXL=0.011*WD2
EPIN2=TKEIN2**1.5/GMIXL*0.1643
EPIN1=TKEIN1**1.5/GMIXL*0.1643
FIINIT(EP)=0.5*(EPIN1+EPIN2)
FIINIT(P1)=1.3E-4
GROUP 13. Boundary conditions and special sources
** Inlet Boundaries
INLET(IN1,LOW,1,1,1,IYJ,1,1,1,1)
VALUE(IN1,P1 , WIN1)
VALUE(IN1,W1 , WIN1)
VALUE(IN1,H1 , 1.0)
VALUE(IN1,KE , TKEIN1)
VALUE(IN1,EP , EPIN1)
VALUE(IN1,G , 0.0)
INLET(IN2,LOW,1,1,IYJ+1,NY,1,1,1,1)
VALUE(IN2,P1, WIN2)
VALUE(IN2,W1, WIN2)
VALUE(IN2,H1, 0.0)
VALUE(IN2,KE, TKEIN2)
VALUE(IN2,EP, EPIN2)
VALUE(IN2,G , 0.0)
**Outlet boundary
PATCH(OUTLET,HIGH,1,NX,1,NY,NZ,NZ,1,1)
COVAL(OUTLET,P1,1.0e05,0.0)
COVAL(OUTLET,W1,ONLYMS,0.0);COVAL(OUTLET,V1,ONLYMS,0.0)
COVAL(OUTLET,KE,ONLYMS,0.0);COVAL(OUTLET,EP,ONLYMS,0.0)
**North-Wall boundary
WALL (WFNN,NORTH,1,NX,NY,NY,1,NZ,1,1)
GROUP 15. Termination of sweeps
LSWEEP=250
RESFAC=0.01
GROUP 16. Termination of iterations
LITHYD=10
GROUP 17. Under-relaxation devices
KELIN=3
RELAX(P1,LINRLX,0.25)
RELAX(V1,FALSDT,0.025);RELAX(W1,FALSDT,0.025)
RELAX(KE,FALSDT,0.025);RELAX(EP,FALSDT,0.025)
RELAX(G,FALSDT ,0.025)
GROUP 19. Data communicated by SATELLITE to GROUND
GROUP 21. Print-out of variables
WALPRN=T;OUTPUT(KE,Y,Y,Y,Y,Y,Y);OUTPUT(H1,Y,Y,Y,Y,Y,Y)
GROUP 22. Monitor print-out
IZMON=NZ-1;IYMON=NY-1;UWATCH=T
GROUP 23. Field print-out and plot control
NPLT=1;NZPRIN=1
NYPRIN=1;IYPRF=1;IYPRL=30
TSTSWP=-1
NAMSAT=MOSG
PLANTBEGIN
**Source term for g
PATCH(SORG,VOLUME,1,NX,1,NY,1,NZ,1,1)
CO=2.0*:RHO1:*EPKE
VAL=GENG/(2.0*:RHO1:*EPKE+TINY)
COVAL(SORG,G , GRND ,GRND )
>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
Above statements contain the formulae for combined
source/sink term for the production and dissipation of
the concentration fluctuations.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
PATCH(WG,VOLUME,1,NX,NY,NY,1,NZ,1,1)
VAL=GENG/(2.0*:RHO1:*EPKE+TINY)
COVAL(WG,G , FIXVAL ,GRND )
>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
Production of G is made equal to its dissipation at
North-Wall boundary by FIXVALing its value to the
production rate divided by twice product of density and
EPKE. The latter is built-in variable standing for
EP/KE.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
** Provide the re-calculation of reference residuals for G
RES=SUM(VOL*(GENG-2.*:RHO1:*EPKE*G)/(NY*NZ))
RESREF(G)=RES
>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
By above two statements the reference residuals for G is
calculated at the and of z-slab as a sum of its
generation rate per cell.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
**Calculation of GENG
* Auxiliary variables
STORE(DFZ,DFY,DFZH,DFYN)
FIINIT(DFZ) =0.0;FIINIT(DFY) =0.0
FIINIT(DFZH)=0.0;FIINIT(DFYN)=0.0
DFZ=((H1[,,+1]-H1)/DZGNZ)**2
REGION(1,NX,1,NY,1,NZ-1,1,1)
DFY=((H1[,+1,]-H1)/DYG2D)**2
REGION(1,1,1,NY-1,1,NZ,1,1)
DFZH=((H1-H1[,,-1])/DZGNZ[,,-1])**2
REGION(1,NX,1,NY,NZ,NZ,1,1)
DFYN=((H1-H1[,-1,])/DYG2D[,-1,])**2
REGION(1,NX,NY,NY,1,NZ,1,1)
>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
The above statements calculate the square concentration
derivatives separately for internal and near domain
bounadary regions.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
GENG=2.8*:RHO1:*ENUT*(DFZ+DFY+DFZH+DFYN)
>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
The sum of the radial and longitudinal derivatives is
multiplied by density and turbulent viscosity times 2.8
to get the generation term.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
**Output calculations
GG - concentration fluctuation;
GGF - concentration fluctuation normalised by local
concentration of central jet fluid.
STORE(GG,GGF)
FIINIT(GGF)=0.0
GG=SQRT(G)
IF(ISWEEP.EQ.LSWEEP)
GGF=GG/(H1+TINY)
IF(ISWEEP.EQ.LSWEEP)
>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
At the end of z-slab for the last sweep the
concentration fluctuation is calculated and normalized
by the local average concentration.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
******************** MFTM section ***************
** Number of fluids in population
INTEGER(NFLUIDS)
NFLUIDS=17
** Micro-mixing constant
REAL(MMC)
MMC = 5.0 ; RG(1) = MMC
** Solve for fluid mass fractions F1, F2, ..., F17
DO II=1,NFLUIDS
SOLVE(F:II:)
TERMS(F:II:,N,Y,y,y,y,y)
PRT(F:II:)= 0.86;PRNDTL(F:II:)= 0.71
RELAX(F:II:,linrlx,0.15)
VARMIN(F:II:)=0.0;VARMAX(F:II:)=1.0
PATCH(PROF:II:,PROFIL,1,1,1,1,1,20,1,1)
PLOT(PROF:II:,F:II:,0.000E+00, 0.000E+00)
ENDDO
ABSIZ=0.5; ORSIZ=0.2
** Fluid population boundary conditions
INLET(IN1,LOW,1,1,1,3,1,1,1,1)
INLET(IN2,LOW,1,1,4,NY,1,1,1,1)
DO II=1,NFLUIDS
VALUE(IN1,F:II:,0.0)
VALUE(IN2,F:II:,0.0)
ENDDO
VALUE(IN1,F1 , 1.0); fiinit(f1)=0
VALUE(IN2,F:NFLUIDS:, 1.0)
** Coupling/splitting rates
PATCH(MIX,PHASEM,1,NX,1,NY,1,NZ,1,1)
* Fluid 1
CO=RG(1)*EPKE*(F3+F5+F7+F9+F11+F13+F15+F17)
COVAL(MIX,F1 , GRND ,0.0 )
* Fluid 2
VAL=2.*RG(1)*EPKE*(F1*F3)-$
RG(1)*EPKE*(F4+F6+F8+F10+F12+F14+F16)*F2
COVAL(MIX,F2, FIXFLU,GRND)
* Fluid 3
VAL=2.*RG(1)*EPKE*(F2*F4+F1*F5)-$
RG(1)*EPKE*(F1+F17+F5+F7+F9+F11+F13+F15)*F3
COVAL(MIX,F3, FIXFLU,GRND)
* Fluid 4
VAL=2.*RG(1)*EPKE*(F3*F5+F2*F6+F1*F7)-$
RG(1)*EPKE*(F2+F6+F8+F10+F12+F14+F16)*F4
COVAL(MIX,F4, FIXFLU,GRND)
* Fluid 5
VAL=2.*RG(1)*EPKE*(F4*F6+F3*F7+F2*F8+F1*F9)-$
RG(1)*EPKE*(F1+F3+F17+F7+F9+F11+F13+F15)*F5
COVAL(MIX,F5, FIXFLU,GRND)
* Fluid 6
VAL=2.*RG(1)*EPKE*(F5*F7+F4*F8+F3*F9+$
F2*F10+F1*F11)-$
RG(1)*EPKE*(F2+F4+F8+F10+F12+F14+F16)*F6
COVAL(MIX,F6, FIXFLU,GRND)
* Fluid 7
VAL=2.*RG(1)*EPKE*(F6*F8+F5*F9+F4*F10+$
F3*F11+F2*F12+F1*F13)-$
RG(1)*EPKE*(F1+F3+F5+F17+F9+F11+F13+F15)*F7
COVAL(MIX,F7, FIXFLU,GRND)
* Fluid 8
VAL=2.*RG(1)*EPKE*(F7*F9+F6*F10+F5*F11+$
F4*F12+F3*F13+F2*F14+F1*F15)-$
RG(1)*EPKE*(F2+F4+F6+F10+F12+F14+F16)*F8
COVAL(MIX,F8, FIXFLU,GRND)
* Fluid 9
VAL=2.*RG(1)*EPKE*(F8*F10+F7*F11+F6*F12+$
F5*F13+F4*F14+F3*F15+F2*F16+F1*F17)-$
RG(1)*EPKE*(F1+F3+F5+F7+F17+F11+F13+F15)*F9
COVAL(MIX,F9, FIXFLU,GRND)
* Fluid 10
VAL=2.*RG(1)*EPKE*(F9*F11+F8*F12+F7*F13+$
F6*F14+F5*F15+F4*F16+F3*F17)-$
RG(1)*EPKE*(F2+F4+F6+F8+F12+F14+F16)*F10
COVAL(MIX,F10, FIXFLU,GRND)
* Fluid 11
VAL=2.*RG(1)*EPKE*(F10*F12+F9*F13+F8*F14+$
F7*F15+F6*F16+F5*F17)-$
RG(1)*EPKE*(F1+F3+F5+F7+F9+F17+F13+F15)*F11
COVAL(MIX,F11, FIXFLU,GRND)
* Fluid 12
VAL=2.*RG(1)*EPKE*(F11*F13+F10*F14+F9*F15+$
F8*F16+F7*F17)-$
RG(1)*EPKE*(F2+F4+F6+F8+F10+F14+F16)*F12
COVAL(MIX,F12, FIXFLU,GRND)
* Fluid 13
VAL=2.*RG(1)*EPKE*(F12*F14+F11*F15+$
F10*F16+F9*F17)-$
RG(1)*EPKE*(F1+F3+F5+F7+F9+F11+F15+F17)*F13
COVAL(MIX,F13, FIXFLU,GRND)
* Fluid 14
VAL=2.*RG(1)*EPKE*(F13*F15+F12*F16+F11*F17)-$
RG(1)*EPKE*(F2+F4+F6+F8+F10+F12+F16)*F14
COVAL(MIX,F14, FIXFLU,GRND)
* Fluid 15
VAL=2.*RG(1)*EPKE*(F14*F16+F13*F17)-$
RG(1)*EPKE*(F1+F3+F5+F7+F9+F11+F13+F17)*F15
COVAL(MIX,F15, FIXFLU,GRND)
* Fluid 16
VAL=2.*RG(1)*EPKE*(F15*F17)-$
RG(1)*EPKE*(F2+F4+F6+F8+F10+F12+F14)*F16
COVAL(MIX,F16, FIXFLU,GRND)
* Fluid 17
CO=RG(1)*EPKE*(F1+F3+F5+F7+F9+F11+F13+F15)
COVAL(MIX,F17 , GRND ,0.0 )
>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
The above source/sink terms in the fluid-mass-fraction
equations are shared according to a coupling/splitting
scheme derived from Spalding concept.
The scheme hypotheses is that the coupling may only take
place between those parent fluids which would produce
the appropriate offsprings inheriting the ATTRIBUTES of
either parent in EQUAL proportion.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
** Output calculations
STORE(CAV,MAS,GAV,GF)
FIINIT(GF)=0.0
CAV=16./16.*F1 + 15./16.*F2 + 14./16.*F3 +$
13./16.*F4 + 12./16.*F5 + 11./16.*F6 +$
10./16.*F7 + 9./16.*F8 + 8./16.*F9 +$
7./16.*F10+ 6./16.*F11+ 5./16.*F12+$
4./16.*F13+ 3./16.*F14+ 2./16.*F15+$
1./16.*F16+ 0./16.*F17
IF(ISWEEP.EQ.LSWEEP)
>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
At the end of the iz-slab for the last sweep, CAV,
averaged concentration of central jet fluid, is
calculated from the individual fluid mass-fractionsand
their arrributes;
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
MAS=F1+F2+F3+F4+F5+F6+F7+F8+F9+F10+F11+F12+$
F13+F14+F15+F16+F17
IF(ISWEEP.EQ.LSWEEP)
>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
At the end of the iz-slab for the last sweep, MAS, sum of
fluid mass fractions, is calculated to check its equality
to unity;
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
GAV=ABS(CAV-16./16)*F1 + ABS(CAV-15./16.)*F2 +$
ABS(CAV-14./16)*F3 + ABS(CAV-13./16.)*F4 +$
ABS(CAV-12./16)*F5 + ABS(CAV-11./16.)*F6 +$
ABS(CAV-10./16)*F7 + ABS(CAV- 9./16.)*F8 +$
ABS(CAV- 8./16)*F9 + ABS(CAV- 7./16.)*F10+$
ABS(CAV- 6./16)*F11+ ABS(CAV- 5./16.)*F12+$
ABS(CAV- 4./16)*F13+ ABS(CAV- 3./16.)*F14+$
ABS(CAV- 2./16)*F15+ ABS(CAV- 1./16.)*F16+$
ABS(CAV- 0./16)*F17
IF(ISWEEP.EQ.LSWEEP)
>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
At the end of the iz-slab for the last sweep, GAV,
averaged concentration fluctuation, is calculated as the
sum of local deviations of averaged concentrations from
the individual concentration attributes;
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
GF =GAV/(CAV+TINY)
IF(ISWEEP.EQ.LSWEEP)
>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
At the end of the iz-slab for the last sweep, GF,
averaged concentration fluctuation normalised by local
averaged concentration of central jet fluid, is
calculated.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
** Output data processing for plotting PDF
* Specify the cell in question: IY=IG(1), IZ=IG(2)
IG(1)=4; IG(2)=4
STORE(FPD);FIINIT(FPD)=0.0
FPD=F1[1,IG(1),IG(2)]*AMAX1(ABS(F1[1,IG(1),IG(2)]$
-YV2D)/(F1[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
REGION(1,1,1,NY,1,1) /ISWEEP.EQ.LSWEEP
FPD=F2[1,IG(1),IG(2)]*AMAX1(ABS(F2[1,IG(1),IG(2)]$
-YV2D)/(F2[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
REGION(1,1,1,NY,2,2) /ISWEEP.EQ.LSWEEP
FPD=F3[1,IG(1),IG(2)]*AMAX1(ABS(F3[1,IG(1),IG(2)]$
-YV2D)/(F3[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
REGION(1,1,1,NY,3,3) /ISWEEP.EQ.LSWEEP
FPD=F4[1,IG(1),IG(2)]*AMAX1(ABS(F4[1,IG(1),IG(2)]$
-YV2D)/(F4[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
REGION(1,1,1,NY,4,4) /ISWEEP.EQ.LSWEEP
FPD=F5[1,IG(1),IG(2)]*AMAX1(ABS(F5[1,IG(1),IG(2)]$
-YV2D)/(F5[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
REGION(1,1,1,NY,5,5) /ISWEEP.EQ.LSWEEP
FPD=F6[1,IG(1),IG(2)]*AMAX1(ABS(F6[1,IG(1),IG(2)]$
-YV2D)/(F6[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
REGION(1,1,1,NY,6,6) /ISWEEP.EQ.LSWEEP
FPD=F7[1,IG(1),IG(2)]*AMAX1(ABS(F7[1,IG(1),IG(2)]$
-YV2D)/(F7[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
REGION(1,1,1,NY,7,7) /ISWEEP.EQ.LSWEEP
FPD=F8[1,IG(1),IG(2)]*AMAX1(ABS(F8[1,IG(1),IG(2)]$
-YV2D)/(F8[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
REGION(1,1,1,NY,8,8) /ISWEEP.EQ.LSWEEP
FPD=F9[1,IG(1),IG(2)]*AMAX1(ABS(F9[1,IG(1),IG(2)]$
-YV2D)/(F9[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
REGION(1,1,1,NY,9,9) /ISWEEP.EQ.LSWEEP
FPD=F10[1,IG(1),IG(2)]*AMAX1(ABS(F10[1,IG(1),IG(2)]$
-YV2D)/(F10[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
REGION(1,1,1,NY,10,10) /ISWEEP.EQ.LSWEEP
FPD=F11[1,IG(1),IG(2)]*AMAX1(ABS(F11[1,IG(1),IG(2)]$
-YV2D)/(F11[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
REGION(1,1,1,NY,11,11) /ISWEEP.EQ.LSWEEP
FPD=F12[1,IG(1),IG(2)]*AMAX1(ABS(F12[1,IG(1),IG(2)]$
-YV2D)/(F12[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
REGION(1,1,1,NY,12,12) /ISWEEP.EQ.LSWEEP
FPD=F13[1,IG(1),IG(2)]*AMAX1(ABS(F13[1,IG(1),IG(2)]$
-YV2D)/(F13[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
REGION(1,1,1,NY,13,13) /ISWEEP.EQ.LSWEEP
FPD=F14[1,IG(1),IG(2)]*AMAX1(ABS(F14[1,IG(1),IG(2)]$
-YV2D)/(F14[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
REGION(1,1,1,NY,14,14) /ISWEEP.EQ.LSWEEP
FPD=F15[1,IG(1),IG(2)]*AMAX1(ABS(F15[1,IG(1),IG(2)]$
-YV2D)/(F15[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
REGION(1,1,1,NY,15,15) /ISWEEP.EQ.LSWEEP
FPD=F16[1,IG(1),IG(2)]*AMAX1(ABS(F16[1,IG(1),IG(2)]$
-YV2D)/(F16[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
REGION(1,1,1,NY,16,16) /ISWEEP.EQ.LSWEEP
FPD=F17[1,IG(1),IG(2)]*AMAX1(ABS(F17[1,IG(1),IG(2)]$
-YV2D)/(F17[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
REGION(1,1,1,NY,17,17) /ISWEEP.EQ.LSWEEP
>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
The above operations are made at the end of the iz-slab
for the last sweep to fill each IY-column of the domain
by Fi value.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
PLANTEND
dmpstk=t
DISTIL=T
EX(P1)=5.613E-01; EX(V1)=2.532E-02; EX(W1)=2.840E+00
EX(KE)=4.442E-01; EX(EP)=6.611E+00; EX(H1)=1.939E-01
EX(FPD)=2.724E-02; EX(GF)=4.936E-01; EX(GAV)=5.976E-02
EX(MAS)=1.000E+00; EX(CAV)=1.939E-01; EX(F17)=1.402E-01
EX(F16)=7.228E-02; EX(F15)=2.095E-01; EX(F14)=3.191E-01
EX(F13)=1.075E-01; EX(F12)=4.737E-02; EX(F11)=2.690E-02
EX(F10)=1.746E-02; EX(F9)=1.188E-02; EX(F8)=8.836E-03
EX(F7)=6.076E-03; EX(F6)=5.269E-03; EX(F5)=3.765E-03
EX(F4)=3.483E-03; EX(F3)=2.355E-03; EX(F2)=3.138E-03
EX(F1)=1.499E-02; EX(GGF)=5.634E+13; EX(GG)=5.768E-02
EX(DFYN)=3.957E-01; EX(DFZH)=3.238E-06; EX(DFY)=1.765E+01
EX(DFZ)=4.226E-02; EX(GENG)=2.216E-01; EX(EPKE)=8.359E+00
EX(GEN1)=1.724E+03; EX(LEN1)=1.567E-02
EX(ENUT)=4.963E-03; EX(G)=5.590E-03
LIBREF=616
STOP