TALK=T;RUN(1,1)
DISPLAY
A hydrocarbon fuel enters steadily a duct of uniform cross-
section steadily, through the porous wall of which is supplied a
steady stream of air, at a low rate for the first 3/4 of the
length, and at an increased rate over the last 1/4 .
The air and fuel are supposed to react in a mixed-is-burned
manner; but smoke is presumed to be generated at a rate which is
proportional to the concentration of unburned fuel times the
absolute temperature raised to the power SMOEXP.
The multi-fuel model of turbulence is used, but with
arbitrarily chosen values of the micro-mixing rate, CONMIX * RATE.
The situation is discussed in the
report by DB Spalding, dated 1998
ENDDIS
PHOTON USE
AUTOPLOT
file; phi 5
cl; da 1; prod; da 1; oxid; da 1; fuel; scale; cola 1; col3 2
colf 3
msg mass fractions of prod - yellow, oxid - blue, free fuel -red
msg
msg flow is from left to right
pause
cl; da 1; smos; da 1; smom; f11; scale; cola 1; col3 2;
msg mass fractions of smoke according to:
msg
msg single-fluid - yellow; multi-fluid - blue
pause
cl; da 1; f1; da 1; f10; da 1; f20; scale; cola 1; col3 2; colf 3
msg mass fractions of f1 - yellow, f10 - blue, f20 - red
pause
cl; da 1; f1; da 1; f3; da 1; f5; da 1; f7; da 1; f9; da 1; f11;
da 1; f13; da 1; f15; da 1; f17; da 1; f19; da 1; f20; scale;
cola 1; col1 2; col1 3; col1 4; col1 5; col3 6; col1 7;
col1 8; col1 9; col1 10; colf 11
msg mass fractions of some intermediate fluids also
pause
cl; da 1; mnsq; cola 1
msg root-mean-square fluctuations
pause
cl; da 1; tmp1; cola 1 ; da 1; tmpm; col3 2
msg temperatures: single fluid - yellow; multi-fluid - blue
pause
cl; da 1; fuel; cola 1; da 1; frfm; col3 2;
msg free fuel from single- (yellow) and multi-fluid models (blue)
pause
enduse
DISPLAY
Notes:
(1) four macros are provided in order reduce the repetiveness of
MFM q1s.
(2) they are named mfm, mfm1, mfm2 and mfm3, and are normally
called in that order at the top of the q1 file, which then
appears as:
#mfm
settings of nflr, nflf, etc
#mfm1
settings of
#mfm2
#mfm3
(3) the character variable mfm is declared and set in the
always-loaded core-library macro 014.htm, which, in its turn,
declares and sets the character variables mfm1, mfm2, mfm3.
(4) the values are:
mfm = $L004
mfm1 = $L003
mfm2 = $L002
mfm3 = $L001
ENDDIS
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(MFM; 1dsmoke; conmix=5. )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 4
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 50
XULAST =1.
XFRAC(1)=0.02 ;XFRAC(2)=0.04
XFRAC(3)=0.06 ;XFRAC(4)=0.08
XFRAC(5)=0.1 ;XFRAC(6)=0.12
XFRAC(7)=0.14 ;XFRAC(8)=0.16
XFRAC(9)=0.18 ;XFRAC(10)=0.2
XFRAC(11)=0.22 ;XFRAC(12)=0.24
XFRAC(13)=0.26 ;XFRAC(14)=0.28
XFRAC(15)=0.3 ;XFRAC(16)=0.32
XFRAC(17)=0.34 ;XFRAC(18)=0.36
XFRAC(19)=0.38 ;XFRAC(20)=0.4
XFRAC(21)=0.42 ;XFRAC(22)=0.44
XFRAC(23)=0.46 ;XFRAC(24)=0.48
XFRAC(25)=0.5 ;XFRAC(26)=0.52
XFRAC(27)=0.54 ;XFRAC(28)=0.56
XFRAC(29)=0.58 ;XFRAC(30)=0.6
XFRAC(31)=0.62 ;XFRAC(32)=0.64
XFRAC(33)=0.66 ;XFRAC(34)=0.68
XFRAC(35)=0.7 ;XFRAC(36)=0.72
XFRAC(37)=0.74 ;XFRAC(38)=0.76
XFRAC(39)=0.78 ;XFRAC(40)=0.8
XFRAC(41)=0.82 ;XFRAC(42)=0.84
XFRAC(43)=0.86 ;XFRAC(44)=0.88
XFRAC(45)=0.9 ;XFRAC(46)=0.92
XFRAC(47)=0.94 ;XFRAC(48)=0.96
XFRAC(49)=0.98 ;XFRAC(50)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 1
YVLAST =0.1
YFRAC(1)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 1
ZWLAST =1.
ZFRAC(1)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(3)=U1
NAME(14)=H1 ;NAME(115)=PRDM
NAME(116)=OXDM ;NAME(117)=FRFM
NAME(118)=TMPM ;NAME(119)=SMOS
NAME(120)=SMOM ;NAME(121)=F1
NAME(122)=F2 ;NAME(123)=F3
NAME(124)=F4 ;NAME(125)=F5
NAME(126)=F6 ;NAME(127)=F7
NAME(128)=F8 ;NAME(129)=F9
NAME(130)=F10 ;NAME(131)=F11
NAME(132)=F12 ;NAME(133)=F13
NAME(134)=F14 ;NAME(135)=F15
NAME(136)=F16 ;NAME(137)=F17
NAME(138)=F18 ;NAME(139)=F19
NAME(140)=F20 ;NAME(141)=MNSQ
NAME(142)=AVEF ;NAME(143)=SPH1
NAME(144)=RHO1 ;NAME(145)=TMP1
NAME(146)=PROD ;NAME(147)=OXID
NAME(148)=FUEL ;NAME(149)=RATE
NAME(150)=MIXF
* 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,N,N,N,Y)
SOLUTN(U1,Y,Y,N,N,N,Y)
SOLUTN(H1,Y,Y,N,N,N,Y)
SOLUTN(PRDM,Y,N,N,N,N,Y)
SOLUTN(OXDM,Y,N,N,N,N,Y)
SOLUTN(FRFM,Y,N,N,N,N,Y)
SOLUTN(TMPM,Y,N,N,N,N,Y)
SOLUTN(SMOS,Y,Y,N,N,N,Y)
SOLUTN(SMOM,Y,Y,N,N,N,Y)
SOLUTN(F1,Y,Y,N,N,N,Y)
SOLUTN(F2,Y,Y,N,N,N,Y)
SOLUTN(F3,Y,Y,N,N,N,Y)
SOLUTN(F4,Y,Y,N,N,N,Y)
SOLUTN(F5,Y,Y,N,N,N,Y)
SOLUTN(F6,Y,Y,N,N,N,Y)
SOLUTN(F7,Y,Y,N,N,N,Y)
SOLUTN(F8,Y,Y,N,N,N,Y)
SOLUTN(F9,Y,Y,N,N,N,Y)
SOLUTN(F10,Y,Y,N,N,N,Y)
SOLUTN(F11,Y,Y,N,N,N,Y)
SOLUTN(F12,Y,Y,N,N,N,Y)
SOLUTN(F13,Y,Y,N,N,N,Y)
SOLUTN(F14,Y,Y,N,N,N,Y)
SOLUTN(F15,Y,Y,N,N,N,Y)
SOLUTN(F16,Y,Y,N,N,N,Y)
SOLUTN(F17,Y,Y,N,N,N,Y)
SOLUTN(F18,Y,Y,N,N,N,Y)
SOLUTN(F19,Y,Y,N,N,N,Y)
SOLUTN(F20,Y,Y,N,N,N,Y)
SOLUTN(MNSQ,Y,N,N,N,N,Y)
SOLUTN(AVEF,Y,N,N,N,N,Y)
SOLUTN(SPH1,Y,N,N,N,N,Y)
SOLUTN(RHO1,Y,N,N,N,N,Y)
SOLUTN(TMP1,Y,N,N,N,N,Y)
SOLUTN(PROD,Y,N,N,N,N,Y)
SOLUTN(OXID,Y,N,N,N,N,Y)
SOLUTN(FUEL,Y,N,N,N,N,Y)
SOLUTN(RATE,Y,N,N,N,N,Y)
SOLUTN(MIXF,Y,Y,N,N,N,Y)
DEN1 = 144
TEMP1 = 145
************************************************************
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(U1,Y,Y,Y,Y,Y,Y)
TERMS(H1,N,Y,Y,N,Y,N)
TERMS(SMOS,N,Y,Y,Y,N,Y)
TERMS(SMOM,N,Y,Y,Y,Y,Y)
TERMS(F1,N,Y,Y,Y,N,Y)
TERMS(F2,N,Y,Y,Y,Y,Y)
TERMS(F3,N,Y,Y,Y,N,Y)
TERMS(F4,N,Y,Y,Y,Y,Y)
TERMS(F5,N,Y,Y,Y,N,Y)
TERMS(F6,N,Y,Y,Y,Y,Y)
TERMS(F7,N,Y,Y,Y,N,Y)
TERMS(F8,N,Y,Y,Y,Y,Y)
TERMS(F9,N,Y,Y,Y,N,Y)
TERMS(F10,N,Y,Y,Y,Y,Y)
TERMS(F11,N,Y,Y,Y,N,Y)
TERMS(F12,N,Y,Y,Y,Y,Y)
TERMS(F13,N,Y,Y,Y,N,Y)
TERMS(F14,N,Y,Y,Y,Y,Y)
TERMS(F15,N,Y,Y,Y,N,Y)
TERMS(F16,N,Y,Y,Y,Y,Y)
TERMS(F17,N,Y,Y,Y,N,Y)
TERMS(F18,N,Y,Y,Y,Y,Y)
TERMS(F19,N,Y,Y,Y,N,Y)
TERMS(F20,N,Y,Y,Y,Y,Y)
TERMS(MIXF,N,Y,Y,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
NEWRH1 = T
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
DENPCO = T
************************************************************
Group 9. Properties used if PRPS is not
stored, and where PRPS = -1.0 if it is!
RHO1 = GRND6 ;TMP1 = GRND7 ;EL1 =0.
TSURR =0. ;TEMP0 =0. ;PRESS0 =1.0E+05
DVO1DT =0. ;DRH1DP =0.
RHO1A =16. ;RHO1B =29. ;RHO1C =28.
TMP1A =1500. ;TMP1B =1500. ;TMP1C =1500.
TMP2A =0.054825 ;TMP2B =4.9E+07 ;TMP2C =1.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
ENUL =1.0E-05 ;ENUT =0.
PRNDTL(U1)=1. ;PRNDTL(H1)=1.
PRNDTL(SMOS)=1. ;PRNDTL(SMOM)=1.
PRNDTL(F1)=1. ;PRNDTL(F2)=1.
PRNDTL(F3)=1. ;PRNDTL(F4)=1.
PRNDTL(F5)=1. ;PRNDTL(F6)=1.
PRNDTL(F7)=1. ;PRNDTL(F8)=1.
PRNDTL(F9)=1. ;PRNDTL(F10)=1.
PRNDTL(F11)=1. ;PRNDTL(F12)=1.
PRNDTL(F13)=1. ;PRNDTL(F14)=1.
PRNDTL(F15)=1. ;PRNDTL(F16)=1.
PRNDTL(F17)=1. ;PRNDTL(F18)=1.
PRNDTL(F19)=1. ;PRNDTL(F20)=1.
PRNDTL(MIXF)=1.
PRT(U1)=1. ;PRT(H1)=1.
PRT(SMOS)=1. ;PRT(SMOM)=1.
PRT(F1)=1. ;PRT(F2)=1.
PRT(F3)=1. ;PRT(F4)=1.
PRT(F5)=1. ;PRT(F6)=1.
PRT(F7)=1. ;PRT(F8)=1.
PRT(F9)=1. ;PRT(F10)=1.
PRT(F11)=1. ;PRT(F12)=1.
PRT(F13)=1. ;PRT(F14)=1.
PRT(F15)=1. ;PRT(F16)=1.
PRT(F17)=1. ;PRT(F18)=1.
PRT(F19)=1. ;PRT(F20)=1.
PRT(MIXF)=1.
CP1 = GRND10 ;CP2 =1.
CP1A =1500. ;CP1B =1500. ;CP1C =1500. ;CP1D =0.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=1.0E-10 ;FIINIT(U1)=2.
FIINIT(H1)=1.1595E+06 ;FIINIT(PRDM)=1.0E-10
FIINIT(OXDM)=1.0E-10 ;FIINIT(FRFM)=1.0E-10
FIINIT(TMPM)=1.0E-10 ;FIINIT(SMOS)=0.
FIINIT(SMOM)=0. ;FIINIT(F1)=0.05
FIINIT(F2)=0.05 ;FIINIT(F3)=0.05
FIINIT(F4)=0.05 ;FIINIT(F5)=0.05
FIINIT(F6)=0.05 ;FIINIT(F7)=0.05
FIINIT(F8)=0.05 ;FIINIT(F9)=0.05
FIINIT(F10)=0.05 ;FIINIT(F11)=0.05
FIINIT(F12)=0.05 ;FIINIT(F13)=0.05
FIINIT(F14)=0.05 ;FIINIT(F15)=0.05
FIINIT(F16)=0.05 ;FIINIT(F17)=0.05
FIINIT(F18)=0.05 ;FIINIT(F19)=0.05
FIINIT(F20)=0.05 ;FIINIT(MNSQ)=1.0E-10
FIINIT(AVEF)=1.0E-10 ;FIINIT(SPH1)=1.0E-10
FIINIT(RHO1)=0.353169 ;FIINIT(TMP1)=1.0E-10
FIINIT(PROD)=1.0E-10 ;FIINIT(OXID)=1.0E-10
FIINIT(FUEL)=1. ;FIINIT(RATE)=100.
FIINIT(MIXF)=1.
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(FUIN ,WEST , 1, 1, 1, 1, 1, 1, 1, 1)
COVAL(FUIN ,P1 , FIXFLU ,0.389704 )
COVAL(FUIN ,U1 ,0. ,2. )
COVAL(FUIN ,H1 ,0. ,5.01595E+07 )
COVAL(FUIN ,F20 ,0. ,1. )
COVAL(FUIN ,MIXF,0. ,1. )
PATCH(OUTLET ,EAST , 50, 50, 1, 1, 1, 1, 1, 1)
COVAL(OUTLET ,P1 ,1. ,0. )
COVAL(OUTLET ,SMOS,0. , SAME )
COVAL(OUTLET ,SMOM,0. , SAME )
COVAL(OUTLET ,F1 ,0. , SAME )
COVAL(OUTLET ,F2 ,0. , SAME )
COVAL(OUTLET ,F3 ,0. , SAME )
COVAL(OUTLET ,F4 ,0. , SAME )
COVAL(OUTLET ,F5 ,0. , SAME )
COVAL(OUTLET ,F6 ,0. , SAME )
COVAL(OUTLET ,F7 ,0. , SAME )
COVAL(OUTLET ,F8 ,0. , SAME )
COVAL(OUTLET ,F9 ,0. , SAME )
COVAL(OUTLET ,F10 ,0. , SAME )
COVAL(OUTLET ,F11 ,0. , SAME )
COVAL(OUTLET ,F12 ,0. , SAME )
COVAL(OUTLET ,F13 ,0. , SAME )
COVAL(OUTLET ,F14 ,0. , SAME )
COVAL(OUTLET ,F15 ,0. , SAME )
COVAL(OUTLET ,F16 ,0. , SAME )
COVAL(OUTLET ,F17 ,0. , SAME )
COVAL(OUTLET ,F18 ,0. , SAME )
COVAL(OUTLET ,F19 ,0. , SAME )
COVAL(OUTLET ,F20 ,0. , SAME )
PATCH(AIRIN ,NORTH , 1, 50, 1, 1, 1, 1, 1, 1)
COVAL(AIRIN ,P1 , FIXFLU ,0.988874 )
COVAL(AIRIN ,H1 ,0. ,1.1595E+06 )
COVAL(AIRIN ,F1 ,0. ,1. )
PATCH(AIRIN2 ,NORTH , 37, 50, 1, 1, 1, 1, 1, 1)
COVAL(AIRIN2 ,P1 , FIXFLU ,28.253557 )
COVAL(AIRIN2 ,H1 ,0. ,1.1595E+06 )
COVAL(AIRIN2 ,F1 ,0. ,1. )
PATCH(SMOKE ,PHASEM, 1, 50, 1, 1, 1, 1, 1, 1)
COVAL(SMOKE ,SMOS, FIXFLU , GRND )
COVAL(SMOKE ,SMOM, FIXFLU , GRND )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 500 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-03
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(U1)=20
LITER(H1)=20 ;LITER(SMOS)=20
LITER(SMOM)=20 ;LITER(F1)=20
LITER(F2)=20 ;LITER(F3)=20
LITER(F4)=20 ;LITER(F5)=20
LITER(F6)=20 ;LITER(F7)=20
LITER(F8)=20 ;LITER(F9)=20
LITER(F10)=20 ;LITER(F11)=20
LITER(F12)=20 ;LITER(F13)=20
LITER(F14)=20 ;LITER(F15)=20
LITER(F16)=20 ;LITER(F17)=20
LITER(F18)=20 ;LITER(F19)=20
LITER(F20)=20 ;LITER(MIXF)=20
ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
ENDIT(H1)=1.0E-03 ;ENDIT(SMOS)=1.0E-03
ENDIT(SMOM)=1.0E-03 ;ENDIT(F1)=1.0E-03
ENDIT(F2)=1.0E-03 ;ENDIT(F3)=1.0E-03
ENDIT(F4)=1.0E-03 ;ENDIT(F5)=1.0E-03
ENDIT(F6)=1.0E-03 ;ENDIT(F7)=1.0E-03
ENDIT(F8)=1.0E-03 ;ENDIT(F9)=1.0E-03
ENDIT(F10)=1.0E-03 ;ENDIT(F11)=1.0E-03
ENDIT(F12)=1.0E-03 ;ENDIT(F13)=1.0E-03
ENDIT(F14)=1.0E-03 ;ENDIT(F15)=1.0E-03
ENDIT(F16)=1.0E-03 ;ENDIT(F17)=1.0E-03
ENDIT(F18)=1.0E-03 ;ENDIT(F19)=1.0E-03
ENDIT(F20)=1.0E-03 ;ENDIT(MIXF)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,0.5)
RELAX(U1,FALSDT,1.)
RELAX(H1,FALSDT,10.)
RELAX(PRDM,LINRLX,1.)
RELAX(OXDM,LINRLX,1.)
RELAX(FRFM,LINRLX,1.)
RELAX(TMPM,LINRLX,1.)
RELAX(SMOS,FALSDT,1.0E+09)
RELAX(SMOM,FALSDT,1.0E+09)
RELAX(F1,LINRLX,0.25)
RELAX(F2,LINRLX,0.25)
RELAX(F3,LINRLX,0.25)
RELAX(F4,LINRLX,0.25)
RELAX(F5,LINRLX,0.25)
RELAX(F6,LINRLX,0.25)
RELAX(F7,LINRLX,0.25)
RELAX(F8,LINRLX,0.25)
RELAX(F9,LINRLX,0.25)
RELAX(F10,LINRLX,0.25)
RELAX(F11,LINRLX,0.25)
RELAX(F12,LINRLX,0.25)
RELAX(F13,LINRLX,0.25)
RELAX(F14,LINRLX,0.25)
RELAX(F15,LINRLX,0.25)
RELAX(F16,LINRLX,0.25)
RELAX(F17,LINRLX,0.25)
RELAX(F18,LINRLX,0.25)
RELAX(F19,LINRLX,0.25)
RELAX(F20,LINRLX,0.25)
RELAX(MNSQ,LINRLX,1.)
RELAX(AVEF,LINRLX,1.)
RELAX(SPH1,LINRLX,1.)
RELAX(RHO1,LINRLX,0.1)
RELAX(TMP1,LINRLX,1.)
RELAX(PROD,LINRLX,1.)
RELAX(OXID,LINRLX,1.)
RELAX(FUEL,LINRLX,1.)
RELAX(RATE,LINRLX,1.)
RELAX(MIXF,FALSDT,1.0E+09)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10
VARMAX(PRDM)=1.0E+10 ;VARMIN(PRDM)=-1.0E+10
VARMAX(OXDM)=1.0E+10 ;VARMIN(OXDM)=-1.0E+10
VARMAX(FRFM)=1.0E+10 ;VARMIN(FRFM)=-1.0E+10
VARMAX(TMPM)=1.0E+10 ;VARMIN(TMPM)=-1.0E+10
VARMAX(SMOS)=1.0E+10 ;VARMIN(SMOS)=-1.0E+10
VARMAX(SMOM)=1.0E+10 ;VARMIN(SMOM)=-1.0E+10
VARMAX(F1)=1. ;VARMIN(F1)=0.
VARMAX(F2)=1. ;VARMIN(F2)=0.
VARMAX(F3)=1. ;VARMIN(F3)=0.
VARMAX(F4)=1. ;VARMIN(F4)=0.
VARMAX(F5)=1. ;VARMIN(F5)=0.
VARMAX(F6)=1. ;VARMIN(F6)=0.
VARMAX(F7)=1. ;VARMIN(F7)=0.
VARMAX(F8)=1. ;VARMIN(F8)=0.
VARMAX(F9)=1. ;VARMIN(F9)=0.
VARMAX(F10)=1. ;VARMIN(F10)=0.
VARMAX(F11)=1. ;VARMIN(F11)=0.
VARMAX(F12)=1. ;VARMIN(F12)=0.
VARMAX(F13)=1. ;VARMIN(F13)=0.
VARMAX(F14)=1. ;VARMIN(F14)=0.
VARMAX(F15)=1. ;VARMIN(F15)=0.
VARMAX(F16)=1. ;VARMIN(F16)=0.
VARMAX(F17)=1. ;VARMIN(F17)=0.
VARMAX(F18)=1. ;VARMIN(F18)=0.
VARMAX(F19)=1. ;VARMIN(F19)=0.
VARMAX(F20)=1. ;VARMIN(F20)=0.
VARMAX(MNSQ)=1. ;VARMIN(MNSQ)=0.
VARMAX(AVEF)=1. ;VARMIN(AVEF)=0.
VARMAX(SPH1)=1.0E+10 ;VARMIN(SPH1)=-1.0E+10
VARMAX(RHO1)=0.706339 ;VARMIN(RHO1)=0.035317
VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=773.
VARMAX(PROD)=1.0E+10 ;VARMIN(PROD)=-1.0E+10
VARMAX(OXID)=1.0E+10 ;VARMIN(OXID)=-1.0E+10
VARMAX(FUEL)=1.0E+10 ;VARMIN(FUEL)=-1.0E+10
VARMAX(RATE)=1.0E+10 ;VARMIN(RATE)=-1.0E+10
VARMAX(MIXF)=1.0E+10 ;VARMIN(MIXF)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
PARSOL = F
ISG62 = 1
SPEDAT(SET,MFM,FFUIN,R,1.)
SPEDAT(SET,MFM,TAIRIN,R,773.)
SPEDAT(SET,MFM,TFUIN,R,773.)
SPEDAT(SET,MFM,TSTOIC,R,2606.314941)
SPEDAT(SET,MFM,FACST,R,0.054825)
SPEDAT(SET,MFM,GASCON,R,8314.)
SPEDAT(SET,MFM,SPVLMX,R,100.)
SPEDAT(SET,MFM,SPVLMN,R,1.0E-02)
SPEDAT(SET,MFM,SMOCON,R,1.)
SPEDAT(SET,MFM,SMOEXP,R,7.)
SPEDAT(SET,MFM,MFMMOD,C,SCRS)
SPEDAT(SET,MFM,NFLUIDS,I,20)
SPEDAT(SET,MFM,NFLR,I,1)
SPEDAT(SET,MFM,NFLF,I,20)
SPEDAT(SET,MFM,CONREA,R,0.)
SPEDAT(SET,MFM,CONMIX,R,5.)
SPEDAT(SET,MFM,VISCON,R,0.)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=3.786E+04 ;EX(U1)=77.07
EX(H1)=9.408E+06 ;EX(PRDM)=0.3831
EX(OXDM)=0. ;EX(FRFM)=0.1397
EX(TMPM)=1460. ;EX(SMOS)=3.547E-04
EX(SMOM)=1.579E-04 ;EX(F1)=0.
EX(F2)=0.1142 ;EX(F3)=0.07434
EX(F4)=0.05543 ;EX(F5)=0.043
EX(F6)=0.03415 ;EX(F7)=0.02767
EX(F8)=0.02282 ;EX(F9)=0.01912
EX(F10)=0.01626 ;EX(F11)=0.01401
EX(F12)=0.01224 ;EX(F13)=0.01083
EX(F14)=9.705E-03 ;EX(F15)=8.832E-03
EX(F16)=8.191E-03 ;EX(F17)=7.813E-03
EX(F18)=7.901E-03 ;EX(F19)=0.01022
EX(F20)=0.03056 ;EX(MNSQ)=0.1773
EX(AVEF)=0.1683 ;EX(SPH1)=1500.
EX(RHO1)=0.2311 ;EX(TMP1)=1997.
EX(PROD)=0.6835 ;EX(OXID)=0.1856
EX(FUEL)=0.1309 ;EX(RATE)=100.
EX(MIXF)=0.1683
************************************************************
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(U1,Y,N,Y,Y,Y,Y)
OUTPUT(H1,Y,N,Y,Y,Y,Y)
OUTPUT(PRDM,Y,N,Y,N,N,N)
OUTPUT(OXDM,Y,N,Y,N,N,N)
OUTPUT(FRFM,Y,N,Y,N,N,N)
OUTPUT(TMPM,Y,N,Y,N,N,N)
OUTPUT(SMOS,Y,N,Y,Y,Y,Y)
OUTPUT(SMOM,Y,N,Y,Y,Y,Y)
OUTPUT(F1,Y,N,Y,Y,Y,Y)
OUTPUT(F2,Y,N,Y,Y,Y,Y)
OUTPUT(F3,Y,N,Y,Y,Y,Y)
OUTPUT(F4,Y,N,Y,Y,Y,Y)
OUTPUT(F5,Y,N,Y,Y,Y,Y)
OUTPUT(F6,Y,N,Y,Y,Y,Y)
OUTPUT(F7,Y,N,Y,Y,Y,Y)
OUTPUT(F8,Y,N,Y,Y,Y,Y)
OUTPUT(F9,Y,N,Y,Y,Y,Y)
OUTPUT(F10,Y,N,Y,Y,Y,Y)
OUTPUT(F11,Y,N,Y,Y,Y,Y)
OUTPUT(F12,Y,N,Y,Y,Y,Y)
OUTPUT(F13,Y,N,Y,Y,Y,Y)
OUTPUT(F14,Y,N,Y,Y,Y,Y)
OUTPUT(F15,Y,N,Y,Y,Y,Y)
OUTPUT(F16,Y,N,Y,Y,Y,Y)
OUTPUT(F17,Y,N,Y,Y,Y,Y)
OUTPUT(F18,Y,N,Y,Y,Y,Y)
OUTPUT(F19,Y,N,Y,Y,Y,Y)
OUTPUT(F20,Y,N,Y,Y,Y,Y)
OUTPUT(MNSQ,Y,N,Y,N,N,N)
OUTPUT(AVEF,Y,N,Y,N,N,N)
OUTPUT(SPH1,Y,N,Y,N,N,N)
OUTPUT(RHO1,Y,N,Y,N,N,N)
OUTPUT(TMP1,Y,N,Y,N,N,N)
OUTPUT(PROD,Y,N,Y,N,N,N)
OUTPUT(OXID,Y,N,Y,N,N,N)
OUTPUT(FUEL,Y,N,Y,N,N,N)
OUTPUT(RATE,Y,N,Y,N,N,N)
OUTPUT(MIXF,Y,N,Y,Y,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 37 ;IYMON = 1 ;IZMON = 1
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 1000 ;NUMCLS = 5
NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
IPLTF = 2 ;IPLTL = 1000 ;NPLT = 5
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.2
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 40
PATCH(LONGPLOT,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
PLOT(LONGPLOT,PROD,0. ,1. )
PLOT(LONGPLOT,OXID,0. ,1. )
PLOT(LONGPLOT,FUEL,0. ,1. )
PLOT(LONGPLOT,MIXF,0. ,1. )
PATCH(LONGPLT2,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
PLOT(LONGPLT2,P1 ,0. ,0. )
PLOT(LONGPLT2,U1 ,0. ,0. )
PLOT(LONGPLT2,TMP1,0. ,0. )
PATCH(FPROF1 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
PLOT(FPROF1 ,F1 ,0. ,0. )
PATCH(FPROF2 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
PLOT(FPROF2 ,F2 ,0. ,0. )
PATCH(FPROF3 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
PLOT(FPROF3 ,F3 ,0. ,0. )
PATCH(FPROF4 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
PLOT(FPROF4 ,F4 ,0. ,0. )
PATCH(FPROF5 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
PLOT(FPROF5 ,F5 ,0. ,0. )
PATCH(FPROF6 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
PLOT(FPROF6 ,F6 ,0. ,0. )
PATCH(FPROF7 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
PLOT(FPROF7 ,F7 ,0. ,0. )
PATCH(FPROF8 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
PLOT(FPROF8 ,F8 ,0. ,0. )
PATCH(FPROF9 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
PLOT(FPROF9 ,F9 ,0. ,0. )
PATCH(FPROF10 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
PLOT(FPROF10 ,F10 ,0. ,0. )
PATCH(FPROF11 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
PLOT(FPROF11 ,F11 ,0. ,0. )
PATCH(FPROF12 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
PLOT(FPROF12 ,F12 ,0. ,0. )
PATCH(FPROF13 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
PLOT(FPROF13 ,F13 ,0. ,0. )
PATCH(FPROF14 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
PLOT(FPROF14 ,F14 ,0. ,0. )
PATCH(FPROF15 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
PLOT(FPROF15 ,F15 ,0. ,0. )
PATCH(FPROF16 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
PLOT(FPROF16 ,F16 ,0. ,0. )
PATCH(FPROF17 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
PLOT(FPROF17 ,F17 ,0. ,0. )
PATCH(FPROF18 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
PLOT(FPROF18 ,F18 ,0. ,0. )
PATCH(FPROF19 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
PLOT(FPROF19 ,F19 ,0. ,0. )
PATCH(FPROF20 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
PLOT(FPROF20 ,F20 ,0. ,0. )
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP