TALK=T;RUN(1,1)
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
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
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 )
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 )
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
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)
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)
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)
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+F$
8+F10+F12+F14+F16)*F2
COVAL(MIX,F2, FIXFLU,GRND)
* Fluid 3
VAL=2.*RG(1)*EPKE*(F2*F4+F1*F5)- RG(1)*EPKE*(F$
1+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)*E$
PKE*(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)- R$
G(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+F1$
2+F14+F16)*F6
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
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+F1$
0+F12+F14+F16)*F8
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+ F$
5*F13+F4*F14+F3*F15+F2*F16+F1*F17)- RG(1)*EPKE*(F1+F3+F5+$
F7+F17+F11+F13+F15)*F9
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
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+F$
5+F7+F9+F17+F13+F15)*F11
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
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
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)- R$
G(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)*EPK$
E*(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 )
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+$
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
4./16.*F13+ 3./16.*F14+ 2./16.*F15+$
1./16.*F16+ 0./16.*F17
IF(ISWEEP.EQ.LSWEEP)
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)
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(C$
AV-12./16)*F5 + ABS(CAV-11./16.)*F6 + ABS(CAV-10./16)$
*F7 + ABS(CAV- 9./16.)*F8 +$
ABS(CAV-12./16)*F5 + ABS(CAV-11./16.)*F6 +$
ABS(CAV-10./16)*F7 + ABS(CAV- 9./16.)*F8 + ABS(C$
AV- 8./16)*F9 + ABS(CAV- 7./16.)*F10+ ABS(CAV- 6./16)$
*F11+ ABS(CAV- 5./16.)*F12+$
ABS(CAV- 8./16)*F9 + ABS(CAV- 7./16.)*F10+$
ABS(CAV- 6./16)*F11+ ABS(CAV- 5./16.)*F12+ ABS(C$
AV- 4./16)*F13+ ABS(CAV- 3./16.)*F14+ ABS(CAV- 2./16)$
*F15+ ABS(CAV- 1./16.)*F16+$
ABS(CAV- 4./16)*F13+ ABS(CAV- 3./16.)*F14+$
ABS(CAV- 2./16)*F15+ ABS(CAV- 1./16.)*F16+ ABS(C$
AV- 0./16)*F17
ABS(CAV- 0./16)*F17
IF(ISWEEP.EQ.LSWEEP)
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)
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
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
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
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(17-fluid turbulence model )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 781
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = F
NX = 1
XULAST =0.1
XFRAC(1)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 15
YVLAST =0.15
YFRAC(1)=0.066667 ;YFRAC(2)=0.133333
YFRAC(3)=0.2 ;YFRAC(4)=0.266667
YFRAC(5)=0.333333 ;YFRAC(6)=0.4
YFRAC(7)=0.466667 ;YFRAC(8)=0.533333
YFRAC(9)=0.6 ;YFRAC(10)=0.666667
YFRAC(11)=0.733333 ;YFRAC(12)=0.8
YFRAC(13)=0.866667 ;YFRAC(14)=0.933333
YFRAC(15)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 20
ZWLAST =3.
ZFRAC(1)=0.05 ;ZFRAC(2)=0.1
ZFRAC(3)=0.15 ;ZFRAC(4)=0.2
ZFRAC(5)=0.25 ;ZFRAC(6)=0.3
ZFRAC(7)=0.35 ;ZFRAC(8)=0.4
ZFRAC(9)=0.45 ;ZFRAC(10)=0.5
ZFRAC(11)=0.55 ;ZFRAC(12)=0.6
ZFRAC(13)=0.65 ;ZFRAC(14)=0.7
ZFRAC(15)=0.75 ;ZFRAC(16)=0.8
ZFRAC(17)=0.85 ;ZFRAC(18)=0.9
ZFRAC(19)=0.95 ;ZFRAC(20)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(5)=V1
NAME(7)=W1 ;NAME(12)=KE
NAME(13)=EP ;NAME(14)=H1
NAME(117)=FPD ;NAME(118)=GF
NAME(119)=GAV ;NAME(120)=MAS
NAME(121)=CAV ;NAME(122)=F17
NAME(123)=F16 ;NAME(124)=F15
NAME(125)=F14 ;NAME(126)=F13
NAME(127)=F12 ;NAME(128)=F11
NAME(129)=F10 ;NAME(130)=F9
NAME(131)=F8 ;NAME(132)=F7
NAME(133)=F6 ;NAME(134)=F5
NAME(135)=F4 ;NAME(136)=F3
NAME(137)=F2 ;NAME(138)=F1
NAME(139)=GGF ;NAME(140)=GG
NAME(141)=DFYN ;NAME(142)=DFZH
NAME(143)=DFY ;NAME(144)=DFZ
NAME(145)=GENG ;NAME(146)=EPKE
NAME(147)=GEN1 ;NAME(148)=LEN1
NAME(149)=ENUT ;NAME(150)=G
* Y in SOLUTN argument list denotes:
* 1-stored 2-solved 3-whole-field
* 4-point-by-point 5-explicit 6-harmonic averaging
SOLUTN(P1,Y,Y,Y,N,N,N)
SOLUTN(V1,Y,Y,N,N,N,Y)
SOLUTN(W1,Y,Y,N,N,N,Y)
SOLUTN(KE,Y,Y,N,N,N,N)
SOLUTN(EP,Y,Y,N,N,N,N)
SOLUTN(H1,Y,Y,N,N,N,Y)
SOLUTN(FPD,Y,N,N,N,N,Y)
SOLUTN(GF,Y,N,N,N,N,Y)
SOLUTN(GAV,Y,N,N,N,N,Y)
SOLUTN(MAS,Y,N,N,N,N,Y)
SOLUTN(CAV,Y,N,N,N,N,Y)
SOLUTN(F17,Y,Y,N,N,N,Y)
SOLUTN(F16,Y,Y,N,N,N,Y)
SOLUTN(F15,Y,Y,N,N,N,Y)
SOLUTN(F14,Y,Y,N,N,N,Y)
SOLUTN(F13,Y,Y,N,N,N,Y)
SOLUTN(F12,Y,Y,N,N,N,Y)
SOLUTN(F11,Y,Y,N,N,N,Y)
SOLUTN(F10,Y,Y,N,N,N,Y)
SOLUTN(F9,Y,Y,N,N,N,Y)
SOLUTN(F8,Y,Y,N,N,N,Y)
SOLUTN(F7,Y,Y,N,N,N,Y)
SOLUTN(F6,Y,Y,N,N,N,Y)
SOLUTN(F5,Y,Y,N,N,N,Y)
SOLUTN(F4,Y,Y,N,N,N,Y)
SOLUTN(F3,Y,Y,N,N,N,Y)
SOLUTN(F2,Y,Y,N,N,N,Y)
SOLUTN(F1,Y,Y,N,N,N,Y)
SOLUTN(GGF,Y,N,N,N,N,Y)
SOLUTN(GG,Y,N,N,N,N,Y)
SOLUTN(DFYN,Y,N,N,N,N,Y)
SOLUTN(DFZH,Y,N,N,N,N,Y)
SOLUTN(DFY,Y,N,N,N,N,Y)
SOLUTN(DFZ,Y,N,N,N,N,Y)
SOLUTN(GENG,Y,N,N,N,N,Y)
SOLUTN(EPKE,Y,N,N,N,N,Y)
SOLUTN(GEN1,Y,N,N,N,N,Y)
SOLUTN(LEN1,Y,N,N,N,N,Y)
SOLUTN(ENUT,Y,N,N,N,N,Y)
SOLUTN(G,Y,Y,N,N,N,Y)
VIST = 149
LEN1 = 148
************************************************************
Group 8. Terms & Devices
* Y in TERMS argument list denotes:
* 1-built-in source 2-convection 3-diffusion 4-transient
* 5-first phase variable 6-interphase transport
TERMS(P1,Y,Y,Y,N,Y,Y)
TERMS(V1,Y,Y,Y,Y,Y,Y)
TERMS(W1,Y,Y,Y,Y,Y,Y)
TERMS(KE,N,Y,Y,Y,Y,N)
TERMS(EP,N,Y,Y,Y,Y,N)
TERMS(H1,N,Y,Y,Y,Y,Y)
TERMS(F17,N,Y,Y,Y,Y,Y)
TERMS(F16,N,Y,Y,Y,Y,Y)
TERMS(F15,N,Y,Y,Y,Y,Y)
TERMS(F14,N,Y,Y,Y,Y,Y)
TERMS(F13,N,Y,Y,Y,Y,Y)
TERMS(F12,N,Y,Y,Y,Y,Y)
TERMS(F11,N,Y,Y,Y,Y,Y)
TERMS(F10,N,Y,Y,Y,Y,Y)
TERMS(F9,N,Y,Y,Y,Y,Y)
TERMS(F8,N,Y,Y,Y,Y,Y)
TERMS(F7,N,Y,Y,Y,Y,Y)
TERMS(F6,N,Y,Y,Y,Y,Y)
TERMS(F5,N,Y,Y,Y,Y,Y)
TERMS(F4,N,Y,Y,Y,Y,Y)
TERMS(F3,N,Y,Y,Y,Y,Y)
TERMS(F2,N,Y,Y,Y,Y,Y)
TERMS(F1,N,Y,Y,Y,Y,Y)
TERMS(G,N,Y,Y,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
NEWENT = T
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
************************************************************
Group 9. Properties used if PRPS is not
stored, and where PRPS = -1.0 if it is!
RHO1 =1. ;TMP1 =0. ;EL1 = GRND4
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =0. ;DRH1DP =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
EL1A =0. ;EL1B =0. ;EL1C =0.
ENUL =3.0E-06 ;ENUT = GRND3
ENUTA =0. ;ENUTB =0. ;ENUTC =0.
IENUTA = 0
PRNDTL(V1)=1. ;PRNDTL(W1)=1.
PRNDTL(KE)=1. ;PRNDTL(EP)=1.
PRNDTL(H1)=0.71 ;PRNDTL(F17)=0.71
PRNDTL(F16)=0.71 ;PRNDTL(F15)=0.71
PRNDTL(F14)=0.71 ;PRNDTL(F13)=0.71
PRNDTL(F12)=0.71 ;PRNDTL(F11)=0.71
PRNDTL(F10)=0.71 ;PRNDTL(F9)=0.71
PRNDTL(F8)=0.71 ;PRNDTL(F7)=0.71
PRNDTL(F6)=0.71 ;PRNDTL(F5)=0.71
PRNDTL(F4)=0.71 ;PRNDTL(F3)=0.71
PRNDTL(F2)=0.71 ;PRNDTL(F1)=0.71
PRNDTL(G)=0.7
PRT(V1)=1. ;PRT(W1)=1.
PRT(KE)=1. ;PRT(EP)=1.314
PRT(H1)=0.86 ;PRT(F17)=0.86
PRT(F16)=0.86 ;PRT(F15)=0.86
PRT(F14)=0.86 ;PRT(F13)=0.86
PRT(F12)=0.86 ;PRT(F11)=0.86
PRT(F10)=0.86 ;PRT(F9)=0.86
PRT(F8)=0.86 ;PRT(F7)=0.86
PRT(F6)=0.86 ;PRT(F5)=0.86
PRT(F4)=0.86 ;PRT(F3)=0.86
PRT(F2)=0.86 ;PRT(F1)=0.86
PRT(G)=0.7
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=1.3E-04 ;FIINIT(V1)=1.0E-10
FIINIT(W1)=6. ;FIINIT(KE)=0.234
FIINIT(EP)=15.149659 ;FIINIT(H1)=1.
FIINIT(FPD)=0. ;FIINIT(GF)=0.
FIINIT(GAV)=1.0E-10 ;FIINIT(MAS)=1.0E-10
FIINIT(CAV)=1.0E-10 ;FIINIT(F17)=1.0E-10
FIINIT(F16)=1.0E-10 ;FIINIT(F15)=1.0E-10
FIINIT(F14)=1.0E-10 ;FIINIT(F13)=1.0E-10
FIINIT(F12)=1.0E-10 ;FIINIT(F11)=1.0E-10
FIINIT(F10)=1.0E-10 ;FIINIT(F9)=1.0E-10
FIINIT(F8)=1.0E-10 ;FIINIT(F7)=1.0E-10
FIINIT(F6)=1.0E-10 ;FIINIT(F5)=1.0E-10
FIINIT(F4)=1.0E-10 ;FIINIT(F3)=1.0E-10
FIINIT(F2)=1.0E-10 ;FIINIT(F1)=0.
FIINIT(GGF)=0. ;FIINIT(GG)=1.0E-10
FIINIT(DFYN)=0. ;FIINIT(DFZH)=0.
FIINIT(DFY)=0. ;FIINIT(DFZ)=0.
FIINIT(GENG)=1.0E-10 ;FIINIT(EPKE)=1.0E-10
FIINIT(GEN1)=1.0E-10 ;FIINIT(LEN1)=0.015
FIINIT(ENUT)=0.015 ;FIINIT(G)=1.0E-10
No PATCHes yet used for this Group
INIADD = F
FSWEEP = 1
NAMFI =CHAM
************************************************************
Group 12. Patchwise adjustment of terms
Patches for this group are printed with those
for Group 13.
Their names begin either with GP12 or &
************************************************************
Group 13. Boundary & Special Sources
PATCH(KESOURCE,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1)
COVAL(KESOURCE,KE , GRND4 , GRND4 )
COVAL(KESOURCE,EP , GRND4 , GRND4 )
PATCH(IN1 ,LOW , 1, 1, 1, 3, 1, 1, 1, 1)
COVAL(IN1 ,P1 , FIXFLU ,10. )
COVAL(IN1 ,V1 ,0. ,0. )
COVAL(IN1 ,W1 ,0. ,10. )
COVAL(IN1 ,KE ,0. ,0.45 )
COVAL(IN1 ,EP ,0. ,30.058847 )
COVAL(IN1 ,H1 ,0. ,1. )
COVAL(IN1 ,F17 ,0. ,0. )
COVAL(IN1 ,F16 ,0. ,0. )
COVAL(IN1 ,F15 ,0. ,0. )
COVAL(IN1 ,F14 ,0. ,0. )
COVAL(IN1 ,F13 ,0. ,0. )
COVAL(IN1 ,F12 ,0. ,0. )
COVAL(IN1 ,F11 ,0. ,0. )
COVAL(IN1 ,F10 ,0. ,0. )
COVAL(IN1 ,F9 ,0. ,0. )
COVAL(IN1 ,F8 ,0. ,0. )
COVAL(IN1 ,F7 ,0. ,0. )
COVAL(IN1 ,F6 ,0. ,0. )
COVAL(IN1 ,F5 ,0. ,0. )
COVAL(IN1 ,F4 ,0. ,0. )
COVAL(IN1 ,F3 ,0. ,0. )
COVAL(IN1 ,F2 ,0. ,0. )
COVAL(IN1 ,F1 ,0. ,1. )
COVAL(IN1 ,G ,0. ,0. )
PATCH(IN2 ,LOW , 1, 1, 4, 15, 1, 1, 1, 1)
COVAL(IN2 ,P1 , FIXFLU ,2. )
COVAL(IN2 ,V1 ,0. ,0. )
COVAL(IN2 ,W1 ,0. ,2. )
COVAL(IN2 ,KE ,0. ,0.018 )
COVAL(IN2 ,EP ,0. ,0.240471 )
COVAL(IN2 ,H1 ,0. ,0. )
COVAL(IN2 ,F17 ,0. ,1. )
COVAL(IN2 ,F16 ,0. ,0. )
COVAL(IN2 ,F15 ,0. ,0. )
COVAL(IN2 ,F14 ,0. ,0. )
COVAL(IN2 ,F13 ,0. ,0. )
COVAL(IN2 ,F12 ,0. ,0. )
COVAL(IN2 ,F11 ,0. ,0. )
COVAL(IN2 ,F10 ,0. ,0. )
COVAL(IN2 ,F9 ,0. ,0. )
COVAL(IN2 ,F8 ,0. ,0. )
COVAL(IN2 ,F7 ,0. ,0. )
COVAL(IN2 ,F6 ,0. ,0. )
COVAL(IN2 ,F5 ,0. ,0. )
COVAL(IN2 ,F4 ,0. ,0. )
COVAL(IN2 ,F3 ,0. ,0. )
COVAL(IN2 ,F2 ,0. ,0. )
COVAL(IN2 ,F1 ,0. ,0. )
COVAL(IN2 ,G ,0. ,0. )
PATCH(OUTLET ,HIGH , 1, 1, 1, 15, 20, 20, 1, 1)
COVAL(OUTLET ,P1 ,1.0E+05 ,0. )
COVAL(OUTLET ,V1 ,0. ,0. )
COVAL(OUTLET ,W1 ,0. ,0. )
COVAL(OUTLET ,KE ,0. ,0. )
COVAL(OUTLET ,EP ,0. ,0. )
PATCH(WFNN ,NWALL , 1, 1, 15, 15, 1, 20, 1, 1)
COVAL(WFNN ,W1 , GRND2 ,0. )
COVAL(WFNN ,KE , GRND2 , GRND2 )
COVAL(WFNN ,EP , GRND2 , GRND2 )
PATCH(ISORG ,VOLUME, 1, 1, 1, 15, 1, 20, 1, 1)
COVAL(ISORG ,G ,In-Form:source - see Grp 19)
PATCH(IWG ,VOLUME, 1, 1, 15, 15, 1, 20, 1, 1)
COVAL(IWG ,G ,In-Form:source - see Grp 19)
PATCH(PAT1 ,CELL , 1, 1, 1, 15, 1, 19, 1, 1)
PATCH(PAT2 ,CELL , 1, 1, 1, 14, 1, 20, 1, 1)
PATCH(PAT3 ,CELL , 1, 1, 1, 15, 20, 20, 1, 1)
PATCH(PAT4 ,CELL , 1, 1, 15, 15, 1, 20, 1, 1)
PATCH(IMIX ,PHASEM, 1, 1, 1, 15, 1, 20, 1, 1)
COVAL(IMIX ,F17 ,In-Form:source - see Grp 19)
COVAL(IMIX ,F16 ,In-Form:source - see Grp 19)
COVAL(IMIX ,F15 ,In-Form:source - see Grp 19)
COVAL(IMIX ,F14 ,In-Form:source - see Grp 19)
COVAL(IMIX ,F13 ,In-Form:source - see Grp 19)
COVAL(IMIX ,F12 ,In-Form:source - see Grp 19)
COVAL(IMIX ,F11 ,In-Form:source - see Grp 19)
COVAL(IMIX ,F10 ,In-Form:source - see Grp 19)
COVAL(IMIX ,F9 ,In-Form:source - see Grp 19)
COVAL(IMIX ,F8 ,In-Form:source - see Grp 19)
COVAL(IMIX ,F7 ,In-Form:source - see Grp 19)
COVAL(IMIX ,F6 ,In-Form:source - see Grp 19)
COVAL(IMIX ,F5 ,In-Form:source - see Grp 19)
COVAL(IMIX ,F4 ,In-Form:source - see Grp 19)
COVAL(IMIX ,F3 ,In-Form:source - see Grp 19)
COVAL(IMIX ,F2 ,In-Form:source - see Grp 19)
COVAL(IMIX ,F1 ,In-Form:source - see Grp 19)
PATCH(PATCH1 ,CELL , 1, 1, 1, 15, 1, 1, 1, 1)
PATCH(PATCH2 ,CELL , 1, 1, 1, 15, 2, 2, 1, 1)
PATCH(PATCH3 ,CELL , 1, 1, 1, 15, 3, 3, 1, 1)
PATCH(PATCH4 ,CELL , 1, 1, 1, 15, 4, 4, 1, 1)
PATCH(PATCH5 ,CELL , 1, 1, 1, 15, 5, 5, 1, 1)
PATCH(PATCH6 ,CELL , 1, 1, 1, 15, 6, 6, 1, 1)
PATCH(PATCH7 ,CELL , 1, 1, 1, 15, 7, 7, 1, 1)
PATCH(PATCH8 ,CELL , 1, 1, 1, 15, 8, 8, 1, 1)
PATCH(PATCH9 ,CELL , 1, 1, 1, 15, 9, 9, 1, 1)
PATCH(PATCH10 ,CELL , 1, 1, 1, 15, 10, 10, 1, 1)
PATCH(PATCH11 ,CELL , 1, 1, 1, 15, 11, 11, 1, 1)
PATCH(PATCH12 ,CELL , 1, 1, 1, 15, 12, 12, 1, 1)
PATCH(PATCH13 ,CELL , 1, 1, 1, 15, 13, 13, 1, 1)
PATCH(PATCH14 ,CELL , 1, 1, 1, 15, 14, 14, 1, 1)
PATCH(PATCH15 ,CELL , 1, 1, 1, 15, 15, 15, 1, 1)
PATCH(PATCH16 ,CELL , 1, 1, 1, 15, 16, 16, 1, 1)
PATCH(PATCH17 ,CELL , 1, 1, 1, 15, 17, 17, 1, 1)
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 250 ;ISWC1 = 1
LITHYD = 10 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-02
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(V1)=10
LITER(W1)=10 ;LITER(KE)=20
LITER(EP)=20 ;LITER(H1)=20
LITER(F17)=20 ;LITER(F16)=20
LITER(F15)=20 ;LITER(F14)=20
LITER(F13)=20 ;LITER(F12)=20
LITER(F11)=20 ;LITER(F10)=20
LITER(F9)=20 ;LITER(F8)=20
LITER(F7)=20 ;LITER(F6)=20
LITER(F5)=20 ;LITER(F4)=20
LITER(F3)=20 ;LITER(F2)=20
LITER(F1)=20 ;LITER(G)=20
ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03
ENDIT(W1)=1.0E-03 ;ENDIT(KE)=1.0E-03
ENDIT(EP)=1.0E-03 ;ENDIT(H1)=1.0E-03
ENDIT(F17)=1.0E-03 ;ENDIT(F16)=1.0E-03
ENDIT(F15)=1.0E-03 ;ENDIT(F14)=1.0E-03
ENDIT(F13)=1.0E-03 ;ENDIT(F12)=1.0E-03
ENDIT(F11)=1.0E-03 ;ENDIT(F10)=1.0E-03
ENDIT(F9)=1.0E-03 ;ENDIT(F8)=1.0E-03
ENDIT(F7)=1.0E-03 ;ENDIT(F6)=1.0E-03
ENDIT(F5)=1.0E-03 ;ENDIT(F4)=1.0E-03
ENDIT(F3)=1.0E-03 ;ENDIT(F2)=1.0E-03
ENDIT(F1)=1.0E-03 ;ENDIT(G)=1.0E-03
************************************************************
Group 17. Relaxation
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(H1,FALSDT,1.0E+09)
RELAX(FPD,LINRLX,1.)
RELAX(GF,LINRLX,1.)
RELAX(GAV,LINRLX,1.)
RELAX(MAS,LINRLX,1.)
RELAX(CAV,LINRLX,1.)
RELAX(F17,LINRLX,0.15)
RELAX(F16,LINRLX,0.15)
RELAX(F15,LINRLX,0.15)
RELAX(F14,LINRLX,0.15)
RELAX(F13,LINRLX,0.15)
RELAX(F12,LINRLX,0.15)
RELAX(F11,LINRLX,0.15)
RELAX(F10,LINRLX,0.15)
RELAX(F9,LINRLX,0.15)
RELAX(F8,LINRLX,0.15)
RELAX(F7,LINRLX,0.15)
RELAX(F6,LINRLX,0.15)
RELAX(F5,LINRLX,0.15)
RELAX(F4,LINRLX,0.15)
RELAX(F3,LINRLX,0.15)
RELAX(F2,LINRLX,0.15)
RELAX(F1,LINRLX,0.15)
RELAX(GGF,LINRLX,1.)
RELAX(GG,LINRLX,1.)
RELAX(DFYN,LINRLX,1.)
RELAX(DFZH,LINRLX,1.)
RELAX(DFY,LINRLX,1.)
RELAX(DFZ,LINRLX,1.)
RELAX(GENG,LINRLX,1.)
RELAX(EPKE,LINRLX,1.)
RELAX(GEN1,LINRLX,1.)
RELAX(LEN1,LINRLX,1.)
RELAX(ENUT,LINRLX,1.)
RELAX(G,FALSDT,0.025)
KELIN = 3
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10
VARMAX(FPD)=1.0E+10 ;VARMIN(FPD)=-1.0E+10
VARMAX(GF)=1.0E+10 ;VARMIN(GF)=-1.0E+10
VARMAX(GAV)=1.0E+10 ;VARMIN(GAV)=-1.0E+10
VARMAX(MAS)=1.0E+10 ;VARMIN(MAS)=-1.0E+10
VARMAX(CAV)=1.0E+10 ;VARMIN(CAV)=-1.0E+10
VARMAX(F17)=1. ;VARMIN(F17)=0.
VARMAX(F16)=1. ;VARMIN(F16)=0.
VARMAX(F15)=1. ;VARMIN(F15)=0.
VARMAX(F14)=1. ;VARMIN(F14)=0.
VARMAX(F13)=1. ;VARMIN(F13)=0.
VARMAX(F12)=1. ;VARMIN(F12)=0.
VARMAX(F11)=1. ;VARMIN(F11)=0.
VARMAX(F10)=1. ;VARMIN(F10)=0.
VARMAX(F9)=1. ;VARMIN(F9)=0.
VARMAX(F8)=1. ;VARMIN(F8)=0.
VARMAX(F7)=1. ;VARMIN(F7)=0.
VARMAX(F6)=1. ;VARMIN(F6)=0.
VARMAX(F5)=1. ;VARMIN(F5)=0.
VARMAX(F4)=1. ;VARMIN(F4)=0.
VARMAX(F3)=1. ;VARMIN(F3)=0.
VARMAX(F2)=1. ;VARMIN(F2)=0.
VARMAX(F1)=1. ;VARMIN(F1)=0.
VARMAX(GGF)=1.0E+10 ;VARMIN(GGF)=-1.0E+10
VARMAX(GG)=1.0E+10 ;VARMIN(GG)=-1.0E+10
VARMAX(DFYN)=1.0E+10 ;VARMIN(DFYN)=-1.0E+10
VARMAX(DFZH)=1.0E+10 ;VARMIN(DFZH)=-1.0E+10
VARMAX(DFY)=1.0E+10 ;VARMIN(DFY)=-1.0E+10
VARMAX(DFZ)=1.0E+10 ;VARMIN(DFZ)=-1.0E+10
VARMAX(GENG)=1.0E+10 ;VARMIN(GENG)=-1.0E+10
VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
VARMAX(GEN1)=1.0E+10 ;VARMIN(GEN1)=-1.0E+10
VARMAX(LEN1)=1.0E+10 ;VARMIN(LEN1)=-1.0E+10
VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
VARMAX(G)=1.0E+10 ;VARMIN(G)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
GENK = T
PARSOL = F
ISG62 = 1
SPEDAT(SET,SOURCE,G!ISORG,C,=2.0*1.*EPKE*(GENG/(2.0*1.*EPKE+TINY)$)
SPEDAT(SET,SOURCE,G!ISORG,C,-G))
SPEDAT(SET,SOURCE,G!IWG,C,=GENG/(2.0*1.*EPKE+TINY)!FIXVAL)
SPEDAT(SET,MAKE,RES,C,=1.E-10)
SPEDAT(SET,STORED,RES,C,=SUM(VOL*(GENG-2.*1.*EPKE*G)/(15*20))!ZSL$)
SPEDAT(SET,STORED,RES,C,FIN)
SPEDAT(SET,STORED,RESREF(G),C,=RES!ZSLFIN)
SPEDAT(SET,STORED,DFZ!PAT1,C,=((H1[&&+1]-H1)/DZG)^2!ZSLFIN)
SPEDAT(SET,STORED,DFY!PAT2,C,=((H1[&+1&]-H1)/DYG)^2!ZSLFIN)
SPEDAT(SET,STORED,DFZH!PAT3,C,=((H1-H1[&&-1])/DZG[&&-1])^2!ZSLFIN)
SPEDAT(SET,STORED,DFYN!PAT4,C,=((H1-H1[&-1&])/DYG[&-1&])^2!ZSLFIN)
SPEDAT(SET,STORED,GENG,C,=2.8*1.*ENUT*(DFZ+DFY+DFZH+DFYN)!ZSLFIN)
SPEDAT(SET,STORED,GG,C,=SQRT(G)!IF(ISWEEP.EQ.250)!ZSLFIN)
SPEDAT(SET,STORED,GGF,C,=GG/(H1+1.0E-20)!IF(ISWEEP.EQ.250)!ZSLFIN)
SPEDAT(SET,SOURCE,F1!IMIX,C,=5.*EPKE*(F3+F5+F7+F9+F11+F13+F15+F17$)
SPEDAT(SET,SOURCE,F1!IMIX,C,)*(0.-F1)!LINE)
SPEDAT(SET,SOURCE,F2!IMIX,C,=2.*5.*EPKE*(F1*F3)-5.*EPKE*(F4+F6+F8$)
SPEDAT(SET,SOURCE,F2!IMIX,C,+F10+F12+F14+F16)*F2)
SPEDAT(SET,SOURCE,F3!IMIX,C,=2.*5.*EPKE*(F2*F4+F1*F5)-5.*EPKE*(F1$)
SPEDAT(SET,SOURCE,F3!IMIX,C,+F17+F5+F7+F9+F11+F13+F15)*F3)
SPEDAT(SET,SOURCE,F4!IMIX,C,=2.*5.*EPKE*(F3*F5+F2*F6+F1*F7)-5.*EP$)
SPEDAT(SET,SOURCE,F4!IMIX,C,KE*(F2+F6+F8+F10+F12+F14+F16)*F4)
SPEDAT(SET,SOURCE,F5!IMIX,C,=2.*5.*EPKE*(F4*F6+F3*F7+F2*F8+F1*F9)$)
SPEDAT(SET,SOURCE,F5!IMIX,C,-5.*EPKE*(F1+F3+F17+F7+F9+F11+F13+F15$)
SPEDAT(SET,SOURCE,F5!IMIX,C,)*F5)
SPEDAT(SET,SOURCE,F6!IMIX,C,=2.*5.*EPKE*(F5*F7+F4*F8+F3*F9+F2*F10$)
SPEDAT(SET,SOURCE,F6!IMIX,C,+F1*F11)-5.*EPKE*(F2+F4+F8+F10+F12+F1$)
SPEDAT(SET,SOURCE,F6!IMIX,C,4+F16)*F6)
SPEDAT(SET,SOURCE,F7!IMIX,C,=2.*5.*EPKE*(F6*F8+F5*F9+F4*F10+F3*F1$)
SPEDAT(SET,SOURCE,F7!IMIX,C,1+F2*F12+F1*F13)-5.*EPKE*(F1+F3+F5+F1$)
SPEDAT(SET,SOURCE,F7!IMIX,C,7+F9+F11+F13+F15)*F7)
SPEDAT(SET,SOURCE,F8!IMIX,C,=2.*5.*EPKE*(F7*F9+F6*F10+F5*F11+F4*F$)
SPEDAT(SET,SOURCE,F8!IMIX,C,12+F3*F13+F2*F14+F1*F15)-5.*EPKE*(F2+$)
SPEDAT(SET,SOURCE,F8!IMIX,C,F4+F6+F10+F12+F14+F16)*F8)
SPEDAT(SET,SOURCE,F9!IMIX,C,=2.*5.*EPKE*(F8*F10+F7*F11+F6*F12+F5*$)
SPEDAT(SET,SOURCE,F9!IMIX,C,F13+F4*F14+F3*F15+F2*F16+F1*F17)-5.*E$)
SPEDAT(SET,SOURCE,F9!IMIX,C,PKE*(F1+F3+F5+F7+F17+F11+F13+F15)*F9)
SPEDAT(SET,SOURCE,F10!IMIX,C,=2.*5.*EPKE*(F9*F11+F8*F12+F7*F13+F6$)
SPEDAT(SET,SOURCE,F10!IMIX,C,*F14+F5*F15+F4*F16+F3*F17)-5.*EPKE*($)
SPEDAT(SET,SOURCE,F10!IMIX,C,F2+F4+F6+F8+F12+F14+F16)*F10)
SPEDAT(SET,SOURCE,F11!IMIX,C,=2.*5.*EPKE*(F10*F12+F9*F13+F8*F14+F$)
SPEDAT(SET,SOURCE,F11!IMIX,C,7*F15+F6*F16+F5*F17)-5.*EPKE*(F1+F3+$)
SPEDAT(SET,SOURCE,F11!IMIX,C,F5+F7+F9+F17+F13+F15)*F11)
SPEDAT(SET,SOURCE,F12!IMIX,C,=2.*5.*EPKE*(F11*F13+F10*F14+F9*F15+$)
SPEDAT(SET,SOURCE,F12!IMIX,C,F8*F16+F7*F17)-5.*EPKE*(F2+F4+F6+F8+$)
SPEDAT(SET,SOURCE,F12!IMIX,C,F10+F14+F16)*F12)
SPEDAT(SET,SOURCE,F13!IMIX,C,=2.*5.*EPKE*(F12*F14+F11*F15+F10*F16$)
SPEDAT(SET,SOURCE,F13!IMIX,C,+F9*F17)-5.*EPKE*(F1+F3+F5+F7+F9+F11$)
SPEDAT(SET,SOURCE,F13!IMIX,C,+F15+F17)*F13)
SPEDAT(SET,SOURCE,F14!IMIX,C,=2.*5.*EPKE*(F13*F15+F12*F16+F11*F17$)
SPEDAT(SET,SOURCE,F14!IMIX,C,)-5.*EPKE*(F2+F4+F6+F8+F10+F12+F16)*$)
SPEDAT(SET,SOURCE,F14!IMIX,C,F14)
SPEDAT(SET,SOURCE,F15!IMIX,C,=2.*5.*EPKE*(F14*F16+F13*F17)-5.*EPK$)
SPEDAT(SET,SOURCE,F15!IMIX,C,E*(F1+F3+F5+F7+F9+F11+F13+F17)*F15)
SPEDAT(SET,SOURCE,F16!IMIX,C,=2.*5.*EPKE*(F15*F17)-5.*EPKE*(F2+F4$)
SPEDAT(SET,SOURCE,F16!IMIX,C,+F6+F8+F10+F12+F14)*F16)
SPEDAT(SET,SOURCE,F17!IMIX,C,=5.*EPKE*(F1+F3+F5+F7+F9+F11+F13+F15$)
SPEDAT(SET,SOURCE,F17!IMIX,C,)*(0.-F17)!LINE)
SPEDAT(SET,STORED,CAV,C,=16./16.*F1+15./16.*F2+14./16.*F3+13./16.$)
SPEDAT(SET,STORED,CAV,C,*F4+12./16.*F5+11./16.*F6+10./16.*F7+9./1$)
SPEDAT(SET,STORED,CAV,C,6.*F8+8./16.*F9!IF(ISWEEP.EQ.250)!ZSLFIN)
SPEDAT(SET,STORED,CAV,C,=CAV+7./16.*F10+6./16.*F11+5./16.*F12+4./$)
SPEDAT(SET,STORED,CAV,C,16.*F13+3./16.*F14+2./16.*F15+1./16.*F16+$)
SPEDAT(SET,STORED,CAV,C,0./16.*F17!IF(ISWEEP.EQ.250)!ZSLFIN)
SPEDAT(SET,STORED,MAS,C,=F1+F2+F3+F4+F5+F6+F7+F8+F9+F10+F11+F12+F$)
SPEDAT(SET,STORED,MAS,C,13+F14+F15+F16+F17!IF(ISWEEP.EQ.250)!ZSLF$)
SPEDAT(SET,STORED,MAS,C,IN)
SPEDAT(SET,STORED,GAV,C,=ABS(CAV-16./16)*F1+ABS(CAV-15./16.)*F2+A$)
SPEDAT(SET,STORED,GAV,C,BS(CAV-14./16)*F3+ABS(CAV-13./16.)*F4+ABS$)
SPEDAT(SET,STORED,GAV,C,(CAV-12./16)*F5+ABS(CAV-11./16.)*F6!ZSLFI$)
SPEDAT(SET,STORED,GAV,C,N)
SPEDAT(SET,STORED,GAV,C,=GAV+ABS(CAV-10./16)*F7+ABS(CAV-9./16.)*F$)
SPEDAT(SET,STORED,GAV,C,8+ABS(CAV-8./16)*F9+ABS(CAV-7./16.)*F10+A$)
SPEDAT(SET,STORED,GAV,C,BS(CAV-6./16)*F11+ABS(CAV-5./16.)*F12!ZSL$)
SPEDAT(SET,STORED,GAV,C,FIN)
SPEDAT(SET,STORED,GAV,C,=GAV+ABS(CAV-4./16)*F13+ABS(CAV-3./16.)*F$)
SPEDAT(SET,STORED,GAV,C,14+ABS(CAV-2./16)*F15+ABS(CAV-1./16.)*F16$)
SPEDAT(SET,STORED,GAV,C,+ABS(CAV-0./16)*F17!ZSLFIN)
SPEDAT(SET,STORED,GF,C,=GAV/(CAV+1.0E-20)!IF(ISWEEP.EQ.250)!ZSLFIN)
SPEDAT(SET,STORED,FPD!PATCH1,C,=F1[1&4&4]*MAX(ABS(F1[1&4&4]-YV)/($)
SPEDAT(SET,STORED,FPD!PATCH1,C,F1[1&4&4]-YV)&0.)!IF(ISWEEP.EQ.250$)
SPEDAT(SET,STORED,FPD!PATCH1,C,)!ZSLFIN)
SPEDAT(SET,STORED,FPD!PATCH2,C,=F2[1&4&4]*MAX(ABS(F2[1&4&4]-YV)/($)
SPEDAT(SET,STORED,FPD!PATCH2,C,F2[1&4&4]-YV)&0.0)!IF(ISWEEP.EQ.25$)
SPEDAT(SET,STORED,FPD!PATCH2,C,0)!ZSLFIN)
SPEDAT(SET,STORED,FPD!PATCH3,C,=F3[1&4&4]*MAX(ABS(F3[1&4&4]-YV)/($)
SPEDAT(SET,STORED,FPD!PATCH3,C,F3[1&4&4]-YV)&0.0)!IF(ISWEEP.EQ.25$)
SPEDAT(SET,STORED,FPD!PATCH3,C,0)!ZSLFIN)
SPEDAT(SET,STORED,FPD!PATCH4,C,=F4[1&4&4]*MAX(ABS(F4[1&4&4]-YV)/($)
SPEDAT(SET,STORED,FPD!PATCH4,C,F4[1&4&4]-YV)&0.0)!IF(ISWEEP.EQ.25$)
SPEDAT(SET,STORED,FPD!PATCH4,C,0)!ZSLFIN)
SPEDAT(SET,STORED,FPD!PATCH5,C,=F5[1&4&4]*MAX(ABS(F5[1&4&4]-YV)/($)
SPEDAT(SET,STORED,FPD!PATCH5,C,F5[1&4&4]-YV)&0.0)!IF(ISWEEP.EQ.25$)
SPEDAT(SET,STORED,FPD!PATCH5,C,0)!ZSLFIN)
SPEDAT(SET,STORED,FPD!PATCH6,C,=F6[1&4&4]*MAX(ABS(F6[1&4&4]-YV)/($)
SPEDAT(SET,STORED,FPD!PATCH6,C,F6[1&4&4]-YV)&0.0)!IF(ISWEEP.EQ.25$)
SPEDAT(SET,STORED,FPD!PATCH6,C,0)!ZSLFIN)
SPEDAT(SET,STORED,FPD!PATCH7,C,=F7[1&4&4]*MAX(ABS(F7[1&4&4]-YV)/($)
SPEDAT(SET,STORED,FPD!PATCH7,C,F7[1&4&4]-YV)&0.0)!IF(ISWEEP.EQ.25$)
SPEDAT(SET,STORED,FPD!PATCH7,C,0)!ZSLFIN)
SPEDAT(SET,STORED,FPD!PATCH8,C,=F8[1&4&4]*MAX(ABS(F8[1&4&4]-YV)/($)
SPEDAT(SET,STORED,FPD!PATCH8,C,F8[1&4&4]-YV)&0.0)!IF(ISWEEP.EQ.25$)
SPEDAT(SET,STORED,FPD!PATCH8,C,0)!ZSLFIN)
SPEDAT(SET,STORED,FPD!PATCH9,C,=F9[1&4&4]*MAX(ABS(F9[1&4&4]-YV)/($)
SPEDAT(SET,STORED,FPD!PATCH9,C,F9[1&4&4]-YV)&0.0)!IF(ISWEEP.EQ.25$)
SPEDAT(SET,STORED,FPD!PATCH9,C,0)!ZSLFIN)
SPEDAT(SET,STORED,FPD!PATCH10,C,=F10[1&4&4]*MAX(ABS(F10[1&4&4]-YV$)
SPEDAT(SET,STORED,FPD!PATCH10,C,)/(F10[1&4&4]-YV)&0.0)!IF(ISWEEP.$)
SPEDAT(SET,STORED,FPD!PATCH10,C,EQ.250)!ZSLFIN)
SPEDAT(SET,STORED,FPD!PATCH11,C,=F11[1&4&4]*MAX(ABS(F11[1&4&4]-YV$)
SPEDAT(SET,STORED,FPD!PATCH11,C,)/(F11[1&4&4]-YV)&0.0)!IF(ISWEEP.$)
SPEDAT(SET,STORED,FPD!PATCH11,C,EQ.250)!ZSLFIN)
SPEDAT(SET,STORED,FPD!PATCH12,C,=F12[1&4&4]*MAX(ABS(F12[1&4&4]-YV$)
SPEDAT(SET,STORED,FPD!PATCH12,C,)/(F12[1&4&4]-YV)&0.0)!IF(ISWEEP.$)
SPEDAT(SET,STORED,FPD!PATCH12,C,EQ.250)!ZSLFIN)
SPEDAT(SET,STORED,FPD!PATCH13,C,=F13[1&4&4]*MAX(ABS(F13[1&4&4]-YV$)
SPEDAT(SET,STORED,FPD!PATCH13,C,)/(F13[1&4&4]-YV)&0.0)!IF(ISWEEP.$)
SPEDAT(SET,STORED,FPD!PATCH13,C,EQ.250)!ZSLFIN)
SPEDAT(SET,STORED,FPD!PATCH14,C,=F14[1&4&4]*MAX(ABS(F14[1&4&4]-YV$)
SPEDAT(SET,STORED,FPD!PATCH14,C,)/(F14[1&4&4]-YV)&0.0)!IF(ISWEEP.$)
SPEDAT(SET,STORED,FPD!PATCH14,C,EQ.250)!ZSLFIN)
SPEDAT(SET,STORED,FPD!PATCH15,C,=F15[1&4&4]*MAX(ABS(F15[1&4&4]-YV$)
SPEDAT(SET,STORED,FPD!PATCH15,C,)/(F15[1&4&4]-YV)&0.0)!IF(ISWEEP.$)
SPEDAT(SET,STORED,FPD!PATCH15,C,EQ.250)!ZSLFIN)
SPEDAT(SET,STORED,FPD!PATCH16,C,=F16[1&4&4]*MAX(ABS(F16[1&4&4]-YV$)
SPEDAT(SET,STORED,FPD!PATCH16,C,)/(F16[1&4&4]-YV)&0.0)!IF(ISWEEP.$)
SPEDAT(SET,STORED,FPD!PATCH16,C,EQ.250)!ZSLFIN)
SPEDAT(SET,STORED,FPD!PATCH17,C,=F17[1&4&4]*MAX(ABS(F17[1&4&4]-YV$)
SPEDAT(SET,STORED,FPD!PATCH17,C,)/(F17[1&4&4]-YV)&0.0)!IF(ISWEEP.$)
SPEDAT(SET,STORED,FPD!PATCH17,C,EQ.250)!ZSLFIN)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
IG( 1) = 4
IG( 2) = 4
RG( 1) =5.
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=0.5613 ;EX(V1)=0.02532
EX(W1)=2.84 ;EX(KE)=0.4442
EX(EP)=6.611 ;EX(H1)=0.1939
EX(FPD)=0.02724 ;EX(GF)=0.4936
EX(GAV)=0.05976 ;EX(MAS)=1.
EX(CAV)=0.1939 ;EX(F17)=0.1402
EX(F16)=0.07228 ;EX(F15)=0.2095
EX(F14)=0.3191 ;EX(F13)=0.1075
EX(F12)=0.04737 ;EX(F11)=0.0269
EX(F10)=0.01746 ;EX(F9)=0.01188
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)=0.01499
EX(GGF)=2.83E+13 ;EX(GG)=0.05768
EX(DFYN)=0.3957 ;EX(DFZH)=3.238E-06
EX(DFY)=17.65 ;EX(DFZ)=0.04226
EX(GENG)=0.2216 ;EX(EPKE)=8.359
EX(GEN1)=1724. ;EX(LEN1)=0.01567
EX(ENUT)=4.963E-03 ;EX(G)=5.591E-03
************************************************************
Group 21. Print-out of Variables
INIFLD = F ;SUBWGR = F
* Y in OUTPUT argument list denotes:
* 1-field 2-correction-eq. monitor 3-selective dumping
* 4-whole-field residual 5-spot-value table 6-residual table
OUTPUT(P1,Y,N,Y,Y,Y,Y)
OUTPUT(V1,Y,N,Y,Y,Y,Y)
OUTPUT(W1,Y,N,Y,Y,Y,Y)
OUTPUT(KE,Y,Y,Y,Y,Y,Y)
OUTPUT(EP,Y,N,Y,Y,Y,Y)
OUTPUT(H1,Y,Y,Y,Y,Y,Y)
OUTPUT(FPD,Y,N,Y,N,N,N)
OUTPUT(GF,Y,N,Y,N,N,N)
OUTPUT(GAV,Y,N,Y,N,N,N)
OUTPUT(MAS,Y,N,Y,N,N,N)
OUTPUT(CAV,Y,N,Y,N,N,N)
OUTPUT(F17,Y,N,Y,Y,Y,Y)
OUTPUT(F16,Y,N,Y,Y,Y,Y)
OUTPUT(F15,Y,N,Y,Y,Y,Y)
OUTPUT(F14,Y,N,Y,Y,Y,Y)
OUTPUT(F13,Y,N,Y,Y,Y,Y)
OUTPUT(F12,Y,N,Y,Y,Y,Y)
OUTPUT(F11,Y,N,Y,Y,Y,Y)
OUTPUT(F10,Y,N,Y,Y,Y,Y)
OUTPUT(F9,Y,N,Y,Y,Y,Y)
OUTPUT(F8,Y,N,Y,Y,Y,Y)
OUTPUT(F7,Y,N,Y,Y,Y,Y)
OUTPUT(F6,Y,N,Y,Y,Y,Y)
OUTPUT(F5,Y,N,Y,Y,Y,Y)
OUTPUT(F4,Y,N,Y,Y,Y,Y)
OUTPUT(F3,Y,N,Y,Y,Y,Y)
OUTPUT(F2,Y,N,Y,Y,Y,Y)
OUTPUT(F1,Y,N,Y,Y,Y,Y)
OUTPUT(GGF,Y,N,Y,N,N,N)
OUTPUT(GG,Y,N,Y,N,N,N)
OUTPUT(DFYN,Y,N,Y,N,N,N)
OUTPUT(DFZH,Y,N,Y,N,N,N)
OUTPUT(DFY,Y,N,Y,N,N,N)
OUTPUT(DFZ,Y,N,Y,N,N,N)
OUTPUT(GENG,Y,N,Y,N,N,N)
OUTPUT(EPKE,Y,N,Y,N,N,N)
OUTPUT(GEN1,Y,N,Y,N,N,N)
OUTPUT(LEN1,Y,N,Y,N,N,N)
OUTPUT(ENUT,Y,N,Y,N,N,N)
OUTPUT(G,Y,N,Y,Y,Y,Y)
WALPRN = T
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 14 ;IZMON = 19
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NYPRIN = 1 ;IYPRF = 1 ;IYPRL = 30
NZPRIN = 1 ;IZPRF = 1 ;IZPRL = 10000
XZPR = F ;YZPR = F
IPLTF = 1 ;IPLTL = -1 ;NPLT = 1
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.2
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
PATCH(PROF1 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
PLOT(PROF1 ,F1 ,0. ,0. )
PATCH(PROF2 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
PLOT(PROF2 ,F2 ,0. ,0. )
PATCH(PROF3 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
PLOT(PROF3 ,F3 ,0. ,0. )
PATCH(PROF4 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
PLOT(PROF4 ,F4 ,0. ,0. )
PATCH(PROF5 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
PLOT(PROF5 ,F5 ,0. ,0. )
PATCH(PROF6 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
PLOT(PROF6 ,F6 ,0. ,0. )
PATCH(PROF7 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
PLOT(PROF7 ,F7 ,0. ,0. )
PATCH(PROF8 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
PLOT(PROF8 ,F8 ,0. ,0. )
PATCH(PROF9 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
PLOT(PROF9 ,F9 ,0. ,0. )
PATCH(PROF10 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
PLOT(PROF10 ,F10 ,0. ,0. )
PATCH(PROF11 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
PLOT(PROF11 ,F11 ,0. ,0. )
PATCH(PROF12 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
PLOT(PROF12 ,F12 ,0. ,0. )
PATCH(PROF13 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
PLOT(PROF13 ,F13 ,0. ,0. )
PATCH(PROF14 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
PLOT(PROF14 ,F14 ,0. ,0. )
PATCH(PROF15 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
PLOT(PROF15 ,F15 ,0. ,0. )
PATCH(PROF16 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
PLOT(PROF16 ,F16 ,0. ,0. )
PATCH(PROF17 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
PLOT(PROF17 ,F17 ,0. ,0. )
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP