TALK=T;RUN(1,1)
Display
A turbulent combustible mixture flows steadily from left to right
towards and into a hot porous plug, which is represented by
fixing the mass fraction of the fully-reacted gas to unity there.
A distorted distance scale is used, so that the whole flame lies
between x=0 and x=1.
The RATE value is fixed arbitrarily, as is also the diffusivity
viscosity.
The case can be used for testing the influence of the fineness
of the population grid on the reaction-rate distribution.
ENDDIS
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(1D Flame Propgn; Conmix=1. Conrea=0.025 )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 4
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 40
XULAST =1.
XFRAC(1)=0.025 ;XFRAC(2)=0.05
XFRAC(3)=0.075 ;XFRAC(4)=0.1
XFRAC(5)=0.125 ;XFRAC(6)=0.15
XFRAC(7)=0.175 ;XFRAC(8)=0.2
XFRAC(9)=0.225 ;XFRAC(10)=0.25
XFRAC(11)=0.275 ;XFRAC(12)=0.3
XFRAC(13)=0.325 ;XFRAC(14)=0.35
XFRAC(15)=0.375 ;XFRAC(16)=0.4
XFRAC(17)=0.425 ;XFRAC(18)=0.45
XFRAC(19)=0.475 ;XFRAC(20)=0.5
XFRAC(21)=0.525 ;XFRAC(22)=0.55
XFRAC(23)=0.575 ;XFRAC(24)=0.6
XFRAC(25)=0.625 ;XFRAC(26)=0.65
XFRAC(27)=0.675 ;XFRAC(28)=0.7
XFRAC(29)=0.725 ;XFRAC(30)=0.75
XFRAC(31)=0.775 ;XFRAC(32)=0.8
XFRAC(33)=0.825 ;XFRAC(34)=0.85
XFRAC(35)=0.875 ;XFRAC(36)=0.9
XFRAC(37)=0.925 ;XFRAC(38)=0.95
XFRAC(39)=0.975 ;XFRAC(40)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 1
YVLAST =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(107)=F1 ;NAME(108)=F2
NAME(109)=F3 ;NAME(110)=F4
NAME(111)=F5 ;NAME(112)=F6
NAME(113)=F7 ;NAME(114)=F8
NAME(115)=F9 ;NAME(116)=F10
NAME(117)=F11 ;NAME(118)=F12
NAME(119)=F13 ;NAME(120)=F14
NAME(121)=F15 ;NAME(122)=F16
NAME(123)=F17 ;NAME(124)=F18
NAME(125)=F19 ;NAME(126)=F20
NAME(127)=F21 ;NAME(128)=F22
NAME(129)=F23 ;NAME(130)=F24
NAME(131)=F25 ;NAME(132)=F26
NAME(133)=F27 ;NAME(134)=F28
NAME(135)=F29 ;NAME(136)=F30
NAME(137)=F31 ;NAME(138)=F32
NAME(139)=F33 ;NAME(140)=F34
NAME(141)=F35 ;NAME(142)=F36
NAME(143)=F37 ;NAME(144)=F38
NAME(145)=F39 ;NAME(146)=F40
NAME(147)=MNSQ ;NAME(148)=AVER
NAME(149)=REAC ;NAME(150)=RATE
* Y in SOLUTN argument list denotes:
* 1-stored 2-solved 3-whole-field
* 4-point-by-point 5-explicit 6-harmonic averaging
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(F21,Y,Y,N,N,N,Y)
SOLUTN(F22,Y,Y,N,N,N,Y)
SOLUTN(F23,Y,Y,N,N,N,Y)
SOLUTN(F24,Y,Y,N,N,N,Y)
SOLUTN(F25,Y,Y,N,N,N,Y)
SOLUTN(F26,Y,Y,N,N,N,Y)
SOLUTN(F27,Y,Y,N,N,N,Y)
SOLUTN(F28,Y,Y,N,N,N,Y)
SOLUTN(F29,Y,Y,N,N,N,Y)
SOLUTN(F30,Y,Y,N,N,N,Y)
SOLUTN(F31,Y,Y,N,N,N,Y)
SOLUTN(F32,Y,Y,N,N,N,Y)
SOLUTN(F33,Y,Y,N,N,N,Y)
SOLUTN(F34,Y,Y,N,N,N,Y)
SOLUTN(F35,Y,Y,N,N,N,Y)
SOLUTN(F36,Y,Y,N,N,N,Y)
SOLUTN(F37,Y,Y,N,N,N,Y)
SOLUTN(F38,Y,Y,N,N,N,Y)
SOLUTN(F39,Y,Y,N,N,N,Y)
SOLUTN(F40,Y,Y,N,N,N,Y)
SOLUTN(MNSQ,Y,N,N,N,N,Y)
SOLUTN(AVER,Y,N,N,N,N,Y)
SOLUTN(REAC,Y,N,N,N,N,Y)
SOLUTN(RATE,Y,N,N,N,N,Y)
************************************************************
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(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(F21,N,Y,Y,Y,N,Y)
TERMS(F22,N,Y,Y,Y,Y,Y)
TERMS(F23,N,Y,Y,Y,N,Y)
TERMS(F24,N,Y,Y,Y,Y,Y)
TERMS(F25,N,Y,Y,Y,N,Y)
TERMS(F26,N,Y,Y,Y,Y,Y)
TERMS(F27,N,Y,Y,Y,N,Y)
TERMS(F28,N,Y,Y,Y,Y,Y)
TERMS(F29,N,Y,Y,Y,N,Y)
TERMS(F30,N,Y,Y,Y,Y,Y)
TERMS(F31,N,Y,Y,Y,N,Y)
TERMS(F32,N,Y,Y,Y,Y,Y)
TERMS(F33,N,Y,Y,Y,N,Y)
TERMS(F34,N,Y,Y,Y,Y,Y)
TERMS(F35,N,Y,Y,Y,N,Y)
TERMS(F36,N,Y,Y,Y,Y,Y)
TERMS(F37,N,Y,Y,Y,N,Y)
TERMS(F38,N,Y,Y,Y,Y,Y)
TERMS(F39,N,Y,Y,Y,N,Y)
TERMS(F40,N,Y,Y,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
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 =0.
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =0. ;DRH1DP =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
ENUL =1. ;ENUT =0.
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(F21)=1. ;PRNDTL(F22)=1.
PRNDTL(F23)=1. ;PRNDTL(F24)=1.
PRNDTL(F25)=1. ;PRNDTL(F26)=1.
PRNDTL(F27)=1. ;PRNDTL(F28)=1.
PRNDTL(F29)=1. ;PRNDTL(F30)=1.
PRNDTL(F31)=1. ;PRNDTL(F32)=1.
PRNDTL(F33)=1. ;PRNDTL(F34)=1.
PRNDTL(F35)=1. ;PRNDTL(F36)=1.
PRNDTL(F37)=1. ;PRNDTL(F38)=1.
PRNDTL(F39)=1. ;PRNDTL(F40)=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(F21)=1. ;PRT(F22)=1.
PRT(F23)=1. ;PRT(F24)=1.
PRT(F25)=1. ;PRT(F26)=1.
PRT(F27)=1. ;PRT(F28)=1.
PRT(F29)=1. ;PRT(F30)=1.
PRT(F31)=1. ;PRT(F32)=1.
PRT(F33)=1. ;PRT(F34)=1.
PRT(F35)=1. ;PRT(F36)=1.
PRT(F37)=1. ;PRT(F38)=1.
PRT(F39)=1. ;PRT(F40)=1.
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(F1)=0. ;FIINIT(F2)=0.
FIINIT(F3)=0. ;FIINIT(F4)=0.
FIINIT(F5)=0. ;FIINIT(F6)=0.
FIINIT(F7)=0. ;FIINIT(F8)=0.
FIINIT(F9)=0. ;FIINIT(F10)=0.
FIINIT(F11)=0. ;FIINIT(F12)=0.
FIINIT(F13)=0. ;FIINIT(F14)=0.
FIINIT(F15)=0. ;FIINIT(F16)=0.
FIINIT(F17)=0. ;FIINIT(F18)=0.
FIINIT(F19)=0. ;FIINIT(F20)=0.
FIINIT(F21)=0. ;FIINIT(F22)=0.
FIINIT(F23)=0. ;FIINIT(F24)=0.
FIINIT(F25)=0. ;FIINIT(F26)=0.
FIINIT(F27)=0. ;FIINIT(F28)=0.
FIINIT(F29)=0. ;FIINIT(F30)=0.
FIINIT(F31)=0. ;FIINIT(F32)=0.
FIINIT(F33)=0. ;FIINIT(F34)=0.
FIINIT(F35)=0. ;FIINIT(F36)=0.
FIINIT(F37)=0. ;FIINIT(F38)=0.
FIINIT(F39)=0. ;FIINIT(F40)=0.
FIINIT(MNSQ)=1.0E-10 ;FIINIT(AVER)=1.0E-10
FIINIT(REAC)=1.0E-10 ;FIINIT(RATE)=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(OUTFLOW ,VOLUME, 1, 1, 1, 1, 1, 1, 1, 1)
PATCH(COLD ,CELL , 1, 1, 1, 1, 1, 1, 1, 1)
COVAL(COLD ,F1 , FIXVAL ,1. )
COVAL(COLD ,F2 , FIXVAL ,0. )
COVAL(COLD ,F3 , FIXVAL ,0. )
COVAL(COLD ,F4 , FIXVAL ,0. )
COVAL(COLD ,F5 , FIXVAL ,0. )
COVAL(COLD ,F6 , FIXVAL ,0. )
COVAL(COLD ,F7 , FIXVAL ,0. )
COVAL(COLD ,F8 , FIXVAL ,0. )
COVAL(COLD ,F9 , FIXVAL ,0. )
COVAL(COLD ,F10 , FIXVAL ,0. )
COVAL(COLD ,F11 , FIXVAL ,0. )
COVAL(COLD ,F12 , FIXVAL ,0. )
COVAL(COLD ,F13 , FIXVAL ,0. )
COVAL(COLD ,F14 , FIXVAL ,0. )
COVAL(COLD ,F15 , FIXVAL ,0. )
COVAL(COLD ,F16 , FIXVAL ,0. )
COVAL(COLD ,F17 , FIXVAL ,0. )
COVAL(COLD ,F18 , FIXVAL ,0. )
COVAL(COLD ,F19 , FIXVAL ,0. )
COVAL(COLD ,F20 , FIXVAL ,0. )
COVAL(COLD ,F21 , FIXVAL ,0. )
COVAL(COLD ,F22 , FIXVAL ,0. )
COVAL(COLD ,F23 , FIXVAL ,0. )
COVAL(COLD ,F24 , FIXVAL ,0. )
COVAL(COLD ,F25 , FIXVAL ,0. )
COVAL(COLD ,F26 , FIXVAL ,0. )
COVAL(COLD ,F27 , FIXVAL ,0. )
COVAL(COLD ,F28 , FIXVAL ,0. )
COVAL(COLD ,F29 , FIXVAL ,0. )
COVAL(COLD ,F30 , FIXVAL ,0. )
COVAL(COLD ,F31 , FIXVAL ,0. )
COVAL(COLD ,F32 , FIXVAL ,0. )
COVAL(COLD ,F33 , FIXVAL ,0. )
COVAL(COLD ,F34 , FIXVAL ,0. )
COVAL(COLD ,F35 , FIXVAL ,0. )
COVAL(COLD ,F36 , FIXVAL ,0. )
COVAL(COLD ,F37 , FIXVAL ,0. )
COVAL(COLD ,F38 , FIXVAL ,0. )
COVAL(COLD ,F39 , FIXVAL ,0. )
COVAL(COLD ,F40 , FIXVAL ,0. )
PATCH(HOT ,CELL , 40, 40, 1, 1, 1, 1, 1, 1)
COVAL(HOT ,F1 , FIXVAL ,0. )
COVAL(HOT ,F2 , FIXVAL ,0. )
COVAL(HOT ,F3 , FIXVAL ,0. )
COVAL(HOT ,F4 , FIXVAL ,0. )
COVAL(HOT ,F5 , FIXVAL ,0. )
COVAL(HOT ,F6 , FIXVAL ,0. )
COVAL(HOT ,F7 , FIXVAL ,0. )
COVAL(HOT ,F8 , FIXVAL ,0. )
COVAL(HOT ,F9 , FIXVAL ,0. )
COVAL(HOT ,F10 , FIXVAL ,0. )
COVAL(HOT ,F11 , FIXVAL ,0. )
COVAL(HOT ,F12 , FIXVAL ,0. )
COVAL(HOT ,F13 , FIXVAL ,0. )
COVAL(HOT ,F14 , FIXVAL ,0. )
COVAL(HOT ,F15 , FIXVAL ,0. )
COVAL(HOT ,F16 , FIXVAL ,0. )
COVAL(HOT ,F17 , FIXVAL ,0. )
COVAL(HOT ,F18 , FIXVAL ,0. )
COVAL(HOT ,F19 , FIXVAL ,0. )
COVAL(HOT ,F20 , FIXVAL ,0. )
COVAL(HOT ,F21 , FIXVAL ,0. )
COVAL(HOT ,F22 , FIXVAL ,0. )
COVAL(HOT ,F23 , FIXVAL ,0. )
COVAL(HOT ,F24 , FIXVAL ,0. )
COVAL(HOT ,F25 , FIXVAL ,0. )
COVAL(HOT ,F26 , FIXVAL ,0. )
COVAL(HOT ,F27 , FIXVAL ,0. )
COVAL(HOT ,F28 , FIXVAL ,0. )
COVAL(HOT ,F29 , FIXVAL ,0. )
COVAL(HOT ,F30 , FIXVAL ,0. )
COVAL(HOT ,F31 , FIXVAL ,0. )
COVAL(HOT ,F32 , FIXVAL ,0. )
COVAL(HOT ,F33 , FIXVAL ,0. )
COVAL(HOT ,F34 , FIXVAL ,0. )
COVAL(HOT ,F35 , FIXVAL ,0. )
COVAL(HOT ,F36 , FIXVAL ,0. )
COVAL(HOT ,F37 , FIXVAL ,0. )
COVAL(HOT ,F38 , FIXVAL ,0. )
COVAL(HOT ,F39 , FIXVAL ,0. )
COVAL(HOT ,F40 ,1000. ,1. )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 5000 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-20
************************************************************
Group 16. Terminate Iterations
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(F21)=20 ;LITER(F22)=20
LITER(F23)=20 ;LITER(F24)=20
LITER(F25)=20 ;LITER(F26)=20
LITER(F27)=20 ;LITER(F28)=20
LITER(F29)=20 ;LITER(F30)=20
LITER(F31)=20 ;LITER(F32)=20
LITER(F33)=20 ;LITER(F34)=20
LITER(F35)=20 ;LITER(F36)=20
LITER(F37)=20 ;LITER(F38)=20
LITER(F39)=20 ;LITER(F40)=20
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(F21)=1.0E-03 ;ENDIT(F22)=1.0E-03
ENDIT(F23)=1.0E-03 ;ENDIT(F24)=1.0E-03
ENDIT(F25)=1.0E-03 ;ENDIT(F26)=1.0E-03
ENDIT(F27)=1.0E-03 ;ENDIT(F28)=1.0E-03
ENDIT(F29)=1.0E-03 ;ENDIT(F30)=1.0E-03
ENDIT(F31)=1.0E-03 ;ENDIT(F32)=1.0E-03
ENDIT(F33)=1.0E-03 ;ENDIT(F34)=1.0E-03
ENDIT(F35)=1.0E-03 ;ENDIT(F36)=1.0E-03
ENDIT(F37)=1.0E-03 ;ENDIT(F38)=1.0E-03
ENDIT(F39)=1.0E-03 ;ENDIT(F40)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(F1,LINRLX,0.5)
RELAX(F2,LINRLX,0.5)
RELAX(F3,LINRLX,0.5)
RELAX(F4,LINRLX,0.5)
RELAX(F5,LINRLX,0.5)
RELAX(F6,LINRLX,0.5)
RELAX(F7,LINRLX,0.5)
RELAX(F8,LINRLX,0.5)
RELAX(F9,LINRLX,0.5)
RELAX(F10,LINRLX,0.5)
RELAX(F11,LINRLX,0.5)
RELAX(F12,LINRLX,0.5)
RELAX(F13,LINRLX,0.5)
RELAX(F14,LINRLX,0.5)
RELAX(F15,LINRLX,0.5)
RELAX(F16,LINRLX,0.5)
RELAX(F17,LINRLX,0.5)
RELAX(F18,LINRLX,0.5)
RELAX(F19,LINRLX,0.5)
RELAX(F20,LINRLX,0.5)
RELAX(F21,LINRLX,0.5)
RELAX(F22,LINRLX,0.5)
RELAX(F23,LINRLX,0.5)
RELAX(F24,LINRLX,0.5)
RELAX(F25,LINRLX,0.5)
RELAX(F26,LINRLX,0.5)
RELAX(F27,LINRLX,0.5)
RELAX(F28,LINRLX,0.5)
RELAX(F29,LINRLX,0.5)
RELAX(F30,LINRLX,0.5)
RELAX(F31,LINRLX,0.5)
RELAX(F32,LINRLX,0.5)
RELAX(F33,LINRLX,0.5)
RELAX(F34,LINRLX,0.5)
RELAX(F35,LINRLX,0.5)
RELAX(F36,LINRLX,0.5)
RELAX(F37,LINRLX,0.5)
RELAX(F38,LINRLX,0.5)
RELAX(F39,LINRLX,0.5)
RELAX(F40,LINRLX,0.5)
RELAX(MNSQ,LINRLX,1.)
RELAX(AVER,LINRLX,1.)
RELAX(REAC,LINRLX,1.)
RELAX(RATE,LINRLX,1.)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
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(F21)=1. ;VARMIN(F21)=0.
VARMAX(F22)=1. ;VARMIN(F22)=0.
VARMAX(F23)=1. ;VARMIN(F23)=0.
VARMAX(F24)=1. ;VARMIN(F24)=0.
VARMAX(F25)=1. ;VARMIN(F25)=0.
VARMAX(F26)=1. ;VARMIN(F26)=0.
VARMAX(F27)=1. ;VARMIN(F27)=0.
VARMAX(F28)=1. ;VARMIN(F28)=0.
VARMAX(F29)=1. ;VARMIN(F29)=0.
VARMAX(F30)=1. ;VARMIN(F30)=0.
VARMAX(F31)=1. ;VARMIN(F31)=0.
VARMAX(F32)=1. ;VARMIN(F32)=0.
VARMAX(F33)=1. ;VARMIN(F33)=0.
VARMAX(F34)=1. ;VARMIN(F34)=0.
VARMAX(F35)=1. ;VARMIN(F35)=0.
VARMAX(F36)=1. ;VARMIN(F36)=0.
VARMAX(F37)=1. ;VARMIN(F37)=0.
VARMAX(F38)=1. ;VARMIN(F38)=0.
VARMAX(F39)=1. ;VARMIN(F39)=0.
VARMAX(F40)=1. ;VARMIN(F40)=0.
VARMAX(MNSQ)=1. ;VARMIN(MNSQ)=0.
VARMAX(AVER)=1. ;VARMIN(AVER)=0.
VARMAX(REAC)=1.0E+10 ;VARMIN(REAC)=-1.0E+10
VARMAX(RATE)=1.0E+10 ;VARMIN(RATE)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
PARSOL = F
ISG62 = 1
SPEDAT(SET,MFM,MFMMOD,C,XFL1D)
SPEDAT(SET,MFM,NFLUIDS,I,40)
SPEDAT(SET,MFM,NFLR,I,40)
SPEDAT(SET,MFM,NFLF,I,1)
SPEDAT(SET,MFM,CONREA,R,0.025)
SPEDAT(SET,MFM,CONMIX,R,1.)
SPEDAT(SET,MFM,RATEXP,R,5.)
SPEDAT(SET,MFM,HIGHF,R,1.)
SPEDAT(SET,MFM,LOWF,R,0.)
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(F1)=0.4197 ;EX(F2)=0.01603
EX(F3)=7.67E-03 ;EX(F4)=6.074E-03
EX(F5)=5.199E-03 ;EX(F6)=4.627E-03
EX(F7)=4.217E-03 ;EX(F8)=3.905E-03
EX(F9)=3.656E-03 ;EX(F10)=3.452E-03
EX(F11)=3.281E-03 ;EX(F12)=3.133E-03
EX(F13)=3.003E-03 ;EX(F14)=2.888E-03
EX(F15)=2.785E-03 ;EX(F16)=2.692E-03
EX(F17)=2.606E-03 ;EX(F18)=2.529E-03
EX(F19)=2.457E-03 ;EX(F20)=2.392E-03
EX(F21)=2.333E-03 ;EX(F22)=2.279E-03
EX(F23)=2.231E-03 ;EX(F24)=2.188E-03
EX(F25)=2.152E-03 ;EX(F26)=2.122E-03
EX(F27)=2.098E-03 ;EX(F28)=2.082E-03
EX(F29)=2.075E-03 ;EX(F30)=2.077E-03
EX(F31)=2.09E-03 ;EX(F32)=2.117E-03
EX(F33)=2.162E-03 ;EX(F34)=2.229E-03
EX(F35)=2.328E-03 ;EX(F36)=2.475E-03
EX(F37)=2.701E-03 ;EX(F38)=3.078E-03
EX(F39)=4.043E-03 ;EX(F40)=0.4528
EX(MNSQ)=0. ;EX(AVER)=0.
EX(REAC)=2.263E-03 ;EX(RATE)=1.
************************************************************
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(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(F21,Y,N,Y,Y,Y,Y)
OUTPUT(F22,Y,N,Y,Y,Y,Y)
OUTPUT(F23,Y,N,Y,Y,Y,Y)
OUTPUT(F24,Y,N,Y,Y,Y,Y)
OUTPUT(F25,Y,N,Y,Y,Y,Y)
OUTPUT(F26,Y,N,Y,Y,Y,Y)
OUTPUT(F27,Y,N,Y,Y,Y,Y)
OUTPUT(F28,Y,N,Y,Y,Y,Y)
OUTPUT(F29,Y,N,Y,Y,Y,Y)
OUTPUT(F30,Y,N,Y,Y,Y,Y)
OUTPUT(F31,Y,N,Y,Y,Y,Y)
OUTPUT(F32,Y,N,Y,Y,Y,Y)
OUTPUT(F33,Y,N,Y,Y,Y,Y)
OUTPUT(F34,Y,N,Y,Y,Y,Y)
OUTPUT(F35,Y,N,Y,Y,Y,Y)
OUTPUT(F36,Y,N,Y,Y,Y,Y)
OUTPUT(F37,Y,N,Y,Y,Y,Y)
OUTPUT(F38,Y,N,Y,Y,Y,Y)
OUTPUT(F39,Y,N,Y,Y,Y,Y)
OUTPUT(F40,Y,N,Y,Y,Y,Y)
OUTPUT(MNSQ,Y,N,Y,N,N,N)
OUTPUT(AVER,Y,N,Y,N,N,N)
OUTPUT(REAC,Y,N,Y,N,N,N)
OUTPUT(RATE,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 20 ;IYMON = 1 ;IZMON = 1
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
PATCH(AVES ,PROFIL, 1, 40, 1, 1, 1, 1, 1, 1)
PLOT(AVES ,REAC,0. ,0. )
PATCH(EXT ,PROFIL, 1, 40, 1, 1, 1, 1, 1, 1)
PLOT(EXT ,F1 ,0. ,1. )
PLOT(EXT ,F40 ,0. ,0. )
PATCH(INT ,PROFIL, 1, 40, 1, 1, 1, 1, 1, 1)
PLOT(INT ,F2 ,0. ,0.03 )
PLOT(INT ,F11 ,0. ,0.03 )
PLOT(INT ,F20 ,0. ,0.03 )
PLOT(INT ,F29 ,0. ,0.03 )
PLOT(INT ,F38 ,0. ,0.03 )
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP