TALK=T;RUN(1,1)
display
Stirred reactor with a 2D population distribution, viz
NFLR reactednesses (or rather mixture fraction minus unburned-fuel
fraction) and NFLF mixture fractions.
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
DISPLAY
provided by l001.htm
Stirred reactor with a 1D population distribution, and
reactedness, ranging from zero to 1, as the population-
distinguishing attribute.
___________
It is supposed that two streams of fluid | |
enter a reactor which is sufficiently A ====> |
well-stirred for space-wise differences | stirred |
of conditions to be negligible, but not | \\\|/// C===>
sufficiently for micro-mixing to be | reactor |
complete. B ====> |
|___________|
The two streams have the same elemental
composition; but one may be more reacted than the other.
The flow is steady, and the total mass flow rate per unit volume
is 1 kg/s m**3 .
The chemical reaction rate is supposed to depend upon reactedness
in accordance with the formula:
rate (in kg/s m**3) = const * (1 - R) / R ** ratexp
where R = reactedness,
ratexp = a prescribed constant
const = conrea / integral [{(1 - R) / R ** ratexp} dR]
conrea = a second prescribed constant.
The micro-mixing rate (as defined in MFM theory) is taken as
conmix kg/s m**3 times the product of the mass fractions of
interacting fluids, in accordance with the "promiscuous Mendelian"
hypothesis, where conmix is again a prescribed constant.
The task is to compute:
(1) the fluid-population distribution in the reactor, and
(2) the population-average reactedness
(3) the root-mean-square of the reactedness fluctuations
These values depend upon 7 dimensionless parameters, namely:-
(a) the prescribed constants: RATEXP, CONREA & CONMIX;
(b) the relative proportions of the two mass flow rates,
namely FLOWA and FLOWB (= 1 - FLOWA);
(c) the reactednesses of the two streams, RA and RB;
and on the number of fluids in terms of which the population
distribution is to be described, NFLR.
In the following Q1, the values employed are as follows:
RATEXP CONREA CONMIX FLOWA RA RB NFLR
5.0 10.0 10.0 0.5 0.0 1.0 25
The opportunity is provided to change NFLR interactively in
the satellite, and CONREA and CONMIX during the EARTH run.
Cases 1-1 and 102 allow choice of FLOWA and RATEXP.
ENDDIS
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(2d Population, Transient, L113 )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 113
************************************************************
Group 2. Time dependence
STEADY = F
* Set overall time and no. of steps
TFIRST =0. ;TLAST =100.
FSTEP = 1 ;LSTEP = 100
TFRAC(1)=1.0E-02 ;TFRAC(2)=0.02
TFRAC(3)=0.03 ;TFRAC(4)=0.04
TFRAC(5)=0.05 ;TFRAC(6)=0.06
TFRAC(7)=0.07 ;TFRAC(8)=0.08
TFRAC(9)=0.09 ;TFRAC(10)=0.1
TFRAC(11)=0.11 ;TFRAC(12)=0.12
TFRAC(13)=0.13 ;TFRAC(14)=0.14
TFRAC(15)=0.15 ;TFRAC(16)=0.16
TFRAC(17)=0.17 ;TFRAC(18)=0.18
TFRAC(19)=0.19 ;TFRAC(20)=0.2
TFRAC(21)=0.21 ;TFRAC(22)=0.22
TFRAC(23)=0.23 ;TFRAC(24)=0.24
TFRAC(25)=0.25 ;TFRAC(26)=0.26
TFRAC(27)=0.27 ;TFRAC(28)=0.28
TFRAC(29)=0.29 ;TFRAC(30)=0.3
TFRAC(31)=0.31 ;TFRAC(32)=0.32
TFRAC(33)=0.33 ;TFRAC(34)=0.34
TFRAC(35)=0.35 ;TFRAC(36)=0.36
TFRAC(37)=0.37 ;TFRAC(38)=0.38
TFRAC(39)=0.39 ;TFRAC(40)=0.4
TFRAC(41)=0.41 ;TFRAC(42)=0.42
TFRAC(43)=0.43 ;TFRAC(44)=0.44
TFRAC(45)=0.45 ;TFRAC(46)=0.46
TFRAC(47)=0.47 ;TFRAC(48)=0.48
TFRAC(49)=0.49 ;TFRAC(50)=0.5
TFRAC(51)=0.51 ;TFRAC(52)=0.52
TFRAC(53)=0.53 ;TFRAC(54)=0.54
TFRAC(55)=0.55 ;TFRAC(56)=0.56
TFRAC(57)=0.57 ;TFRAC(58)=0.58
TFRAC(59)=0.59 ;TFRAC(60)=0.6
TFRAC(61)=0.61 ;TFRAC(62)=0.62
TFRAC(63)=0.63 ;TFRAC(64)=0.64
TFRAC(65)=0.65 ;TFRAC(66)=0.66
TFRAC(67)=0.67 ;TFRAC(68)=0.68
TFRAC(69)=0.69 ;TFRAC(70)=0.7
TFRAC(71)=0.71 ;TFRAC(72)=0.72
TFRAC(73)=0.73 ;TFRAC(74)=0.74
TFRAC(75)=0.75 ;TFRAC(76)=0.76
TFRAC(77)=0.77 ;TFRAC(78)=0.78
TFRAC(79)=0.79 ;TFRAC(80)=0.8
TFRAC(81)=0.81 ;TFRAC(82)=0.82
TFRAC(83)=0.83 ;TFRAC(84)=0.84
TFRAC(85)=0.85 ;TFRAC(86)=0.86
TFRAC(87)=0.87 ;TFRAC(88)=0.88
TFRAC(89)=0.89 ;TFRAC(90)=0.9
TFRAC(91)=0.91 ;TFRAC(92)=0.92
TFRAC(93)=0.93 ;TFRAC(94)=0.94
TFRAC(95)=0.95 ;TFRAC(96)=0.96
TFRAC(97)=0.97 ;TFRAC(98)=0.98
TFRAC(99)=0.99 ;TFRAC(100)=1.
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 1
XULAST =1.
XFRAC(1)=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(1)=P1 ;NAME(48)=F1
NAME(49)=F2 ;NAME(50)=F3
NAME(51)=F4 ;NAME(52)=F5
NAME(53)=F6 ;NAME(54)=F7
NAME(55)=F8 ;NAME(56)=F9
NAME(57)=F10 ;NAME(58)=F11
NAME(59)=F12 ;NAME(60)=F13
NAME(61)=F14 ;NAME(62)=F15
NAME(63)=F16 ;NAME(64)=F17
NAME(65)=F18 ;NAME(66)=F19
NAME(67)=F20 ;NAME(68)=F21
NAME(69)=F22 ;NAME(70)=F23
NAME(71)=F24 ;NAME(72)=F25
NAME(73)=F26 ;NAME(74)=F27
NAME(75)=F28 ;NAME(76)=F29
NAME(77)=F30 ;NAME(78)=F31
NAME(79)=F32 ;NAME(80)=F33
NAME(81)=F34 ;NAME(82)=F35
NAME(83)=F36 ;NAME(84)=F37
NAME(85)=F38 ;NAME(86)=F39
NAME(87)=F40 ;NAME(88)=F41
NAME(89)=F42 ;NAME(90)=F43
NAME(91)=F44 ;NAME(92)=F45
NAME(93)=F46 ;NAME(94)=F47
NAME(95)=F48 ;NAME(96)=F49
NAME(97)=F50 ;NAME(98)=F51
NAME(99)=F52 ;NAME(100)=F53
NAME(101)=F54 ;NAME(102)=F55
NAME(103)=F56 ;NAME(104)=F57
NAME(105)=F58 ;NAME(106)=F59
NAME(107)=F60 ;NAME(108)=F61
NAME(109)=F62 ;NAME(110)=F63
NAME(111)=F64 ;NAME(112)=F65
NAME(113)=F66 ;NAME(114)=F67
NAME(115)=F68 ;NAME(116)=F69
NAME(117)=F70 ;NAME(118)=F71
NAME(119)=F72 ;NAME(120)=F73
NAME(121)=F74 ;NAME(122)=F75
NAME(123)=F76 ;NAME(124)=F77
NAME(125)=F78 ;NAME(126)=F79
NAME(127)=F80 ;NAME(128)=F81
NAME(129)=F82 ;NAME(130)=F83
NAME(131)=F84 ;NAME(132)=F85
NAME(133)=F86 ;NAME(134)=F87
NAME(135)=F88 ;NAME(136)=F89
NAME(137)=F90 ;NAME(138)=F91
NAME(139)=F92 ;NAME(140)=F93
NAME(141)=F94 ;NAME(142)=F95
NAME(143)=F96 ;NAME(144)=F97
NAME(145)=F98 ;NAME(146)=F99
NAME(147)=F100 ;NAME(148)=AVER
NAME(149)=AVEF ;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(P1,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(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(F41,Y,Y,N,N,N,Y)
SOLUTN(F42,Y,Y,N,N,N,Y)
SOLUTN(F43,Y,Y,N,N,N,Y)
SOLUTN(F44,Y,Y,N,N,N,Y)
SOLUTN(F45,Y,Y,N,N,N,Y)
SOLUTN(F46,Y,Y,N,N,N,Y)
SOLUTN(F47,Y,Y,N,N,N,Y)
SOLUTN(F48,Y,Y,N,N,N,Y)
SOLUTN(F49,Y,Y,N,N,N,Y)
SOLUTN(F50,Y,Y,N,N,N,Y)
SOLUTN(F51,Y,Y,N,N,N,Y)
SOLUTN(F52,Y,Y,N,N,N,Y)
SOLUTN(F53,Y,Y,N,N,N,Y)
SOLUTN(F54,Y,Y,N,N,N,Y)
SOLUTN(F55,Y,Y,N,N,N,Y)
SOLUTN(F56,Y,Y,N,N,N,Y)
SOLUTN(F57,Y,Y,N,N,N,Y)
SOLUTN(F58,Y,Y,N,N,N,Y)
SOLUTN(F59,Y,Y,N,N,N,Y)
SOLUTN(F60,Y,Y,N,N,N,Y)
SOLUTN(F61,Y,Y,N,N,N,Y)
SOLUTN(F62,Y,Y,N,N,N,Y)
SOLUTN(F63,Y,Y,N,N,N,Y)
SOLUTN(F64,Y,Y,N,N,N,Y)
SOLUTN(F65,Y,Y,N,N,N,Y)
SOLUTN(F66,Y,Y,N,N,N,Y)
SOLUTN(F67,Y,Y,N,N,N,Y)
SOLUTN(F68,Y,Y,N,N,N,Y)
SOLUTN(F69,Y,Y,N,N,N,Y)
SOLUTN(F70,Y,Y,N,N,N,Y)
SOLUTN(F71,Y,Y,N,N,N,Y)
SOLUTN(F72,Y,Y,N,N,N,Y)
SOLUTN(F73,Y,Y,N,N,N,Y)
SOLUTN(F74,Y,Y,N,N,N,Y)
SOLUTN(F75,Y,Y,N,N,N,Y)
SOLUTN(F76,Y,Y,N,N,N,Y)
SOLUTN(F77,Y,Y,N,N,N,Y)
SOLUTN(F78,Y,Y,N,N,N,Y)
SOLUTN(F79,Y,Y,N,N,N,Y)
SOLUTN(F80,Y,Y,N,N,N,Y)
SOLUTN(F81,Y,Y,N,N,N,Y)
SOLUTN(F82,Y,Y,N,N,N,Y)
SOLUTN(F83,Y,Y,N,N,N,Y)
SOLUTN(F84,Y,Y,N,N,N,Y)
SOLUTN(F85,Y,Y,N,N,N,Y)
SOLUTN(F86,Y,Y,N,N,N,Y)
SOLUTN(F87,Y,Y,N,N,N,Y)
SOLUTN(F88,Y,Y,N,N,N,Y)
SOLUTN(F89,Y,Y,N,N,N,Y)
SOLUTN(F90,Y,Y,N,N,N,Y)
SOLUTN(F91,Y,Y,N,N,N,Y)
SOLUTN(F92,Y,Y,N,N,N,Y)
SOLUTN(F93,Y,Y,N,N,N,Y)
SOLUTN(F94,Y,Y,N,N,N,Y)
SOLUTN(F95,Y,Y,N,N,N,Y)
SOLUTN(F96,Y,Y,N,N,N,Y)
SOLUTN(F97,Y,Y,N,N,N,Y)
SOLUTN(F98,Y,Y,N,N,N,Y)
SOLUTN(F99,Y,Y,N,N,N,Y)
SOLUTN(F100,Y,Y,N,N,N,Y)
SOLUTN(AVER,Y,N,N,N,N,Y)
SOLUTN(AVEF,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(P1,Y,Y,Y,N,Y,Y)
TERMS(F1,N,Y,Y,Y,Y,Y)
TERMS(F2,N,Y,Y,Y,N,Y)
TERMS(F3,N,Y,Y,Y,Y,Y)
TERMS(F4,N,Y,Y,Y,N,Y)
TERMS(F5,N,Y,Y,Y,Y,Y)
TERMS(F6,N,Y,Y,Y,N,Y)
TERMS(F7,N,Y,Y,Y,Y,Y)
TERMS(F8,N,Y,Y,Y,N,Y)
TERMS(F9,N,Y,Y,Y,Y,Y)
TERMS(F10,N,Y,Y,Y,N,Y)
TERMS(F11,N,Y,Y,Y,Y,Y)
TERMS(F12,N,Y,Y,Y,N,Y)
TERMS(F13,N,Y,Y,Y,Y,Y)
TERMS(F14,N,Y,Y,Y,N,Y)
TERMS(F15,N,Y,Y,Y,Y,Y)
TERMS(F16,N,Y,Y,Y,N,Y)
TERMS(F17,N,Y,Y,Y,Y,Y)
TERMS(F18,N,Y,Y,Y,N,Y)
TERMS(F19,N,Y,Y,Y,Y,Y)
TERMS(F20,N,Y,Y,Y,N,Y)
TERMS(F21,N,Y,Y,Y,Y,Y)
TERMS(F22,N,Y,Y,Y,N,Y)
TERMS(F23,N,Y,Y,Y,Y,Y)
TERMS(F24,N,Y,Y,Y,N,Y)
TERMS(F25,N,Y,Y,Y,Y,Y)
TERMS(F26,N,Y,Y,Y,N,Y)
TERMS(F27,N,Y,Y,Y,Y,Y)
TERMS(F28,N,Y,Y,Y,N,Y)
TERMS(F29,N,Y,Y,Y,Y,Y)
TERMS(F30,N,Y,Y,Y,N,Y)
TERMS(F31,N,Y,Y,Y,Y,Y)
TERMS(F32,N,Y,Y,Y,N,Y)
TERMS(F33,N,Y,Y,Y,Y,Y)
TERMS(F34,N,Y,Y,Y,N,Y)
TERMS(F35,N,Y,Y,Y,Y,Y)
TERMS(F36,N,Y,Y,Y,N,Y)
TERMS(F37,N,Y,Y,Y,Y,Y)
TERMS(F38,N,Y,Y,Y,N,Y)
TERMS(F39,N,Y,Y,Y,Y,Y)
TERMS(F40,N,Y,Y,Y,N,Y)
TERMS(F41,N,Y,Y,Y,Y,Y)
TERMS(F42,N,Y,Y,Y,N,Y)
TERMS(F43,N,Y,Y,Y,Y,Y)
TERMS(F44,N,Y,Y,Y,N,Y)
TERMS(F45,N,Y,Y,Y,Y,Y)
TERMS(F46,N,Y,Y,Y,N,Y)
TERMS(F47,N,Y,Y,Y,Y,Y)
TERMS(F48,N,Y,Y,Y,N,Y)
TERMS(F49,N,Y,Y,Y,Y,Y)
TERMS(F50,N,Y,Y,Y,N,Y)
TERMS(F51,N,Y,Y,Y,Y,Y)
TERMS(F52,N,Y,Y,Y,N,Y)
TERMS(F53,N,Y,Y,Y,Y,Y)
TERMS(F54,N,Y,Y,Y,N,Y)
TERMS(F55,N,Y,Y,Y,Y,Y)
TERMS(F56,N,Y,Y,Y,N,Y)
TERMS(F57,N,Y,Y,Y,Y,Y)
TERMS(F58,N,Y,Y,Y,N,Y)
TERMS(F59,N,Y,Y,Y,Y,Y)
TERMS(F60,N,Y,Y,Y,N,Y)
TERMS(F61,N,Y,Y,Y,Y,Y)
TERMS(F62,N,Y,Y,Y,N,Y)
TERMS(F63,N,Y,Y,Y,Y,Y)
TERMS(F64,N,Y,Y,Y,N,Y)
TERMS(F65,N,Y,Y,Y,Y,Y)
TERMS(F66,N,Y,Y,Y,N,Y)
TERMS(F67,N,Y,Y,Y,Y,Y)
TERMS(F68,N,Y,Y,Y,N,Y)
TERMS(F69,N,Y,Y,Y,Y,Y)
TERMS(F70,N,Y,Y,Y,N,Y)
TERMS(F71,N,Y,Y,Y,Y,Y)
TERMS(F72,N,Y,Y,Y,N,Y)
TERMS(F73,N,Y,Y,Y,Y,Y)
TERMS(F74,N,Y,Y,Y,N,Y)
TERMS(F75,N,Y,Y,Y,Y,Y)
TERMS(F76,N,Y,Y,Y,N,Y)
TERMS(F77,N,Y,Y,Y,Y,Y)
TERMS(F78,N,Y,Y,Y,N,Y)
TERMS(F79,N,Y,Y,Y,Y,Y)
TERMS(F80,N,Y,Y,Y,N,Y)
TERMS(F81,N,Y,Y,Y,Y,Y)
TERMS(F82,N,Y,Y,Y,N,Y)
TERMS(F83,N,Y,Y,Y,Y,Y)
TERMS(F84,N,Y,Y,Y,N,Y)
TERMS(F85,N,Y,Y,Y,Y,Y)
TERMS(F86,N,Y,Y,Y,N,Y)
TERMS(F87,N,Y,Y,Y,Y,Y)
TERMS(F88,N,Y,Y,Y,N,Y)
TERMS(F89,N,Y,Y,Y,Y,Y)
TERMS(F90,N,Y,Y,Y,N,Y)
TERMS(F91,N,Y,Y,Y,Y,Y)
TERMS(F92,N,Y,Y,Y,N,Y)
TERMS(F93,N,Y,Y,Y,Y,Y)
TERMS(F94,N,Y,Y,Y,N,Y)
TERMS(F95,N,Y,Y,Y,Y,Y)
TERMS(F96,N,Y,Y,Y,N,Y)
TERMS(F97,N,Y,Y,Y,Y,Y)
TERMS(F98,N,Y,Y,Y,N,Y)
TERMS(F99,N,Y,Y,Y,Y,Y)
TERMS(F100,N,Y,Y,Y,N,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.0E-05 ;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.
PRNDTL(F41)=1. ;PRNDTL(F42)=1.
PRNDTL(F43)=1. ;PRNDTL(F44)=1.
PRNDTL(F45)=1. ;PRNDTL(F46)=1.
PRNDTL(F47)=1. ;PRNDTL(F48)=1.
PRNDTL(F49)=1. ;PRNDTL(F50)=1.
PRNDTL(F51)=1. ;PRNDTL(F52)=1.
PRNDTL(F53)=1. ;PRNDTL(F54)=1.
PRNDTL(F55)=1. ;PRNDTL(F56)=1.
PRNDTL(F57)=1. ;PRNDTL(F58)=1.
PRNDTL(F59)=1. ;PRNDTL(F60)=1.
PRNDTL(F61)=1. ;PRNDTL(F62)=1.
PRNDTL(F63)=1. ;PRNDTL(F64)=1.
PRNDTL(F65)=1. ;PRNDTL(F66)=1.
PRNDTL(F67)=1. ;PRNDTL(F68)=1.
PRNDTL(F69)=1. ;PRNDTL(F70)=1.
PRNDTL(F71)=1. ;PRNDTL(F72)=1.
PRNDTL(F73)=1. ;PRNDTL(F74)=1.
PRNDTL(F75)=1. ;PRNDTL(F76)=1.
PRNDTL(F77)=1. ;PRNDTL(F78)=1.
PRNDTL(F79)=1. ;PRNDTL(F80)=1.
PRNDTL(F81)=1. ;PRNDTL(F82)=1.
PRNDTL(F83)=1. ;PRNDTL(F84)=1.
PRNDTL(F85)=1. ;PRNDTL(F86)=1.
PRNDTL(F87)=1. ;PRNDTL(F88)=1.
PRNDTL(F89)=1. ;PRNDTL(F90)=1.
PRNDTL(F91)=1. ;PRNDTL(F92)=1.
PRNDTL(F93)=1. ;PRNDTL(F94)=1.
PRNDTL(F95)=1. ;PRNDTL(F96)=1.
PRNDTL(F97)=1. ;PRNDTL(F98)=1.
PRNDTL(F99)=1. ;PRNDTL(F100)=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.
PRT(F41)=1. ;PRT(F42)=1.
PRT(F43)=1. ;PRT(F44)=1.
PRT(F45)=1. ;PRT(F46)=1.
PRT(F47)=1. ;PRT(F48)=1.
PRT(F49)=1. ;PRT(F50)=1.
PRT(F51)=1. ;PRT(F52)=1.
PRT(F53)=1. ;PRT(F54)=1.
PRT(F55)=1. ;PRT(F56)=1.
PRT(F57)=1. ;PRT(F58)=1.
PRT(F59)=1. ;PRT(F60)=1.
PRT(F61)=1. ;PRT(F62)=1.
PRT(F63)=1. ;PRT(F64)=1.
PRT(F65)=1. ;PRT(F66)=1.
PRT(F67)=1. ;PRT(F68)=1.
PRT(F69)=1. ;PRT(F70)=1.
PRT(F71)=1. ;PRT(F72)=1.
PRT(F73)=1. ;PRT(F74)=1.
PRT(F75)=1. ;PRT(F76)=1.
PRT(F77)=1. ;PRT(F78)=1.
PRT(F79)=1. ;PRT(F80)=1.
PRT(F81)=1. ;PRT(F82)=1.
PRT(F83)=1. ;PRT(F84)=1.
PRT(F85)=1. ;PRT(F86)=1.
PRT(F87)=1. ;PRT(F88)=1.
PRT(F89)=1. ;PRT(F90)=1.
PRT(F91)=1. ;PRT(F92)=1.
PRT(F93)=1. ;PRT(F94)=1.
PRT(F95)=1. ;PRT(F96)=1.
PRT(F97)=1. ;PRT(F98)=1.
PRT(F99)=1. ;PRT(F100)=1.
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=1.0E-10 ;FIINIT(F1)=1.0E-10
FIINIT(F2)=1.0E-10 ;FIINIT(F3)=1.0E-10
FIINIT(F4)=1.0E-10 ;FIINIT(F5)=1.0E-10
FIINIT(F6)=1.0E-10 ;FIINIT(F7)=1.0E-10
FIINIT(F8)=1.0E-10 ;FIINIT(F9)=1.0E-10
FIINIT(F10)=1.0E-10 ;FIINIT(F11)=1.0E-10
FIINIT(F12)=1.0E-10 ;FIINIT(F13)=1.0E-10
FIINIT(F14)=1.0E-10 ;FIINIT(F15)=1.0E-10
FIINIT(F16)=1.0E-10 ;FIINIT(F17)=1.0E-10
FIINIT(F18)=1.0E-10 ;FIINIT(F19)=1.0E-10
FIINIT(F20)=1.0E-10 ;FIINIT(F21)=1.0E-10
FIINIT(F22)=1.0E-10 ;FIINIT(F23)=1.0E-10
FIINIT(F24)=1.0E-10 ;FIINIT(F25)=1.0E-10
FIINIT(F26)=1.0E-10 ;FIINIT(F27)=1.0E-10
FIINIT(F28)=1.0E-10 ;FIINIT(F29)=1.0E-10
FIINIT(F30)=1.0E-10 ;FIINIT(F31)=1.0E-10
FIINIT(F32)=1.0E-10 ;FIINIT(F33)=1.0E-10
FIINIT(F34)=1.0E-10 ;FIINIT(F35)=1.0E-10
FIINIT(F36)=1.0E-10 ;FIINIT(F37)=1.0E-10
FIINIT(F38)=1.0E-10 ;FIINIT(F39)=1.0E-10
FIINIT(F40)=1.0E-10 ;FIINIT(F41)=1.0E-10
FIINIT(F42)=1.0E-10 ;FIINIT(F43)=1.0E-10
FIINIT(F44)=1.0E-10 ;FIINIT(F45)=1.0E-10
FIINIT(F46)=1.0E-10 ;FIINIT(F47)=1.0E-10
FIINIT(F48)=1.0E-10 ;FIINIT(F49)=1.0E-10
FIINIT(F50)=1.0E-10 ;FIINIT(F51)=1.0E-10
FIINIT(F52)=1.0E-10 ;FIINIT(F53)=1.0E-10
FIINIT(F54)=1.0E-10 ;FIINIT(F55)=1.0E-10
FIINIT(F56)=1.0E-10 ;FIINIT(F57)=1.0E-10
FIINIT(F58)=1.0E-10 ;FIINIT(F59)=1.0E-10
FIINIT(F60)=1.0E-10 ;FIINIT(F61)=1.0E-10
FIINIT(F62)=1.0E-10 ;FIINIT(F63)=1.0E-10
FIINIT(F64)=1.0E-10 ;FIINIT(F65)=1.0E-10
FIINIT(F66)=1.0E-10 ;FIINIT(F67)=1.0E-10
FIINIT(F68)=1.0E-10 ;FIINIT(F69)=1.0E-10
FIINIT(F70)=1.0E-10 ;FIINIT(F71)=1.0E-10
FIINIT(F72)=1.0E-10 ;FIINIT(F73)=1.0E-10
FIINIT(F74)=1.0E-10 ;FIINIT(F75)=1.0E-10
FIINIT(F76)=1.0E-10 ;FIINIT(F77)=1.0E-10
FIINIT(F78)=1.0E-10 ;FIINIT(F79)=1.0E-10
FIINIT(F80)=1.0E-10 ;FIINIT(F81)=1.0E-10
FIINIT(F82)=1.0E-10 ;FIINIT(F83)=1.0E-10
FIINIT(F84)=1.0E-10 ;FIINIT(F85)=1.0E-10
FIINIT(F86)=1.0E-10 ;FIINIT(F87)=1.0E-10
FIINIT(F88)=1.0E-10 ;FIINIT(F89)=1.0E-10
FIINIT(F90)=1.0E-10 ;FIINIT(F91)=1.0E-10
FIINIT(F92)=1.0E-10 ;FIINIT(F93)=1.0E-10
FIINIT(F94)=1.0E-10 ;FIINIT(F95)=1.0E-10
FIINIT(F96)=1.0E-10 ;FIINIT(F97)=1.0E-10
FIINIT(F98)=1.0E-10 ;FIINIT(F99)=1.0E-10
FIINIT(F100)=1. ;FIINIT(AVER)=1.0E-10
FIINIT(AVEF)=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(INFLOWA ,CELL , 1, 1, 1, 1, 1, 1, 1, 100)
COVAL(INFLOWA ,P1 , FIXFLU ,0.5 )
COVAL(INFLOWA ,F1 ,0. ,1. )
PATCH(INFLOWB ,CELL , 1, 1, 1, 1, 1, 1, 1, 100)
COVAL(INFLOWB ,P1 , FIXFLU ,0.5 )
COVAL(INFLOWB ,F100,0. ,1. )
PATCH(OUTFLOW ,VOLUME, 1, 1, 1, 1, 1, 1, 1, 100)
COVAL(OUTFLOW ,P1 ,1. ,0. )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 25 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-20
************************************************************
Group 16. Terminate Iterations
LITER(P1)=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(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 ;LITER(F41)=20
LITER(F42)=20 ;LITER(F43)=20
LITER(F44)=20 ;LITER(F45)=20
LITER(F46)=20 ;LITER(F47)=20
LITER(F48)=20 ;LITER(F49)=20
LITER(F50)=20 ;LITER(F51)=20
LITER(F52)=20 ;LITER(F53)=20
LITER(F54)=20 ;LITER(F55)=20
LITER(F56)=20 ;LITER(F57)=20
LITER(F58)=20 ;LITER(F59)=20
LITER(F60)=20 ;LITER(F61)=20
LITER(F62)=20 ;LITER(F63)=20
LITER(F64)=20 ;LITER(F65)=20
LITER(F66)=20 ;LITER(F67)=20
LITER(F68)=20 ;LITER(F69)=20
LITER(F70)=20 ;LITER(F71)=20
LITER(F72)=20 ;LITER(F73)=20
LITER(F74)=20 ;LITER(F75)=20
LITER(F76)=20 ;LITER(F77)=20
LITER(F78)=20 ;LITER(F79)=20
LITER(F80)=20 ;LITER(F81)=20
LITER(F82)=20 ;LITER(F83)=20
LITER(F84)=20 ;LITER(F85)=20
LITER(F86)=20 ;LITER(F87)=20
LITER(F88)=20 ;LITER(F89)=20
LITER(F90)=20 ;LITER(F91)=20
LITER(F92)=20 ;LITER(F93)=20
LITER(F94)=20 ;LITER(F95)=20
LITER(F96)=20 ;LITER(F97)=20
LITER(F98)=20 ;LITER(F99)=20
LITER(F100)=20
ENDIT(P1)=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(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 ;ENDIT(F41)=1.0E-03
ENDIT(F42)=1.0E-03 ;ENDIT(F43)=1.0E-03
ENDIT(F44)=1.0E-03 ;ENDIT(F45)=1.0E-03
ENDIT(F46)=1.0E-03 ;ENDIT(F47)=1.0E-03
ENDIT(F48)=1.0E-03 ;ENDIT(F49)=1.0E-03
ENDIT(F50)=1.0E-03 ;ENDIT(F51)=1.0E-03
ENDIT(F52)=1.0E-03 ;ENDIT(F53)=1.0E-03
ENDIT(F54)=1.0E-03 ;ENDIT(F55)=1.0E-03
ENDIT(F56)=1.0E-03 ;ENDIT(F57)=1.0E-03
ENDIT(F58)=1.0E-03 ;ENDIT(F59)=1.0E-03
ENDIT(F60)=1.0E-03 ;ENDIT(F61)=1.0E-03
ENDIT(F62)=1.0E-03 ;ENDIT(F63)=1.0E-03
ENDIT(F64)=1.0E-03 ;ENDIT(F65)=1.0E-03
ENDIT(F66)=1.0E-03 ;ENDIT(F67)=1.0E-03
ENDIT(F68)=1.0E-03 ;ENDIT(F69)=1.0E-03
ENDIT(F70)=1.0E-03 ;ENDIT(F71)=1.0E-03
ENDIT(F72)=1.0E-03 ;ENDIT(F73)=1.0E-03
ENDIT(F74)=1.0E-03 ;ENDIT(F75)=1.0E-03
ENDIT(F76)=1.0E-03 ;ENDIT(F77)=1.0E-03
ENDIT(F78)=1.0E-03 ;ENDIT(F79)=1.0E-03
ENDIT(F80)=1.0E-03 ;ENDIT(F81)=1.0E-03
ENDIT(F82)=1.0E-03 ;ENDIT(F83)=1.0E-03
ENDIT(F84)=1.0E-03 ;ENDIT(F85)=1.0E-03
ENDIT(F86)=1.0E-03 ;ENDIT(F87)=1.0E-03
ENDIT(F88)=1.0E-03 ;ENDIT(F89)=1.0E-03
ENDIT(F90)=1.0E-03 ;ENDIT(F91)=1.0E-03
ENDIT(F92)=1.0E-03 ;ENDIT(F93)=1.0E-03
ENDIT(F94)=1.0E-03 ;ENDIT(F95)=1.0E-03
ENDIT(F96)=1.0E-03 ;ENDIT(F97)=1.0E-03
ENDIT(F98)=1.0E-03 ;ENDIT(F99)=1.0E-03
ENDIT(F100)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
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(F41,LINRLX,0.5)
RELAX(F42,LINRLX,0.5)
RELAX(F43,LINRLX,0.5)
RELAX(F44,LINRLX,0.5)
RELAX(F45,LINRLX,0.5)
RELAX(F46,LINRLX,0.5)
RELAX(F47,LINRLX,0.5)
RELAX(F48,LINRLX,0.5)
RELAX(F49,LINRLX,0.5)
RELAX(F50,LINRLX,0.5)
RELAX(F51,LINRLX,0.5)
RELAX(F52,LINRLX,0.5)
RELAX(F53,LINRLX,0.5)
RELAX(F54,LINRLX,0.5)
RELAX(F55,LINRLX,0.5)
RELAX(F56,LINRLX,0.5)
RELAX(F57,LINRLX,0.5)
RELAX(F58,LINRLX,0.5)
RELAX(F59,LINRLX,0.5)
RELAX(F60,LINRLX,0.5)
RELAX(F61,LINRLX,0.5)
RELAX(F62,LINRLX,0.5)
RELAX(F63,LINRLX,0.5)
RELAX(F64,LINRLX,0.5)
RELAX(F65,LINRLX,0.5)
RELAX(F66,LINRLX,0.5)
RELAX(F67,LINRLX,0.5)
RELAX(F68,LINRLX,0.5)
RELAX(F69,LINRLX,0.5)
RELAX(F70,LINRLX,0.5)
RELAX(F71,LINRLX,0.5)
RELAX(F72,LINRLX,0.5)
RELAX(F73,LINRLX,0.5)
RELAX(F74,LINRLX,0.5)
RELAX(F75,LINRLX,0.5)
RELAX(F76,LINRLX,0.5)
RELAX(F77,LINRLX,0.5)
RELAX(F78,LINRLX,0.5)
RELAX(F79,LINRLX,0.5)
RELAX(F80,LINRLX,0.5)
RELAX(F81,LINRLX,0.5)
RELAX(F82,LINRLX,0.5)
RELAX(F83,LINRLX,0.5)
RELAX(F84,LINRLX,0.5)
RELAX(F85,LINRLX,0.5)
RELAX(F86,LINRLX,0.5)
RELAX(F87,LINRLX,0.5)
RELAX(F88,LINRLX,0.5)
RELAX(F89,LINRLX,0.5)
RELAX(F90,LINRLX,0.5)
RELAX(F91,LINRLX,0.5)
RELAX(F92,LINRLX,0.5)
RELAX(F93,LINRLX,0.5)
RELAX(F94,LINRLX,0.5)
RELAX(F95,LINRLX,0.5)
RELAX(F96,LINRLX,0.5)
RELAX(F97,LINRLX,0.5)
RELAX(F98,LINRLX,0.5)
RELAX(F99,LINRLX,0.5)
RELAX(F100,LINRLX,0.5)
RELAX(AVER,LINRLX,1.)
RELAX(AVEF,LINRLX,1.)
RELAX(RATE,LINRLX,1.)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-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(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(F41)=1. ;VARMIN(F41)=0.
VARMAX(F42)=1. ;VARMIN(F42)=0.
VARMAX(F43)=1. ;VARMIN(F43)=0.
VARMAX(F44)=1. ;VARMIN(F44)=0.
VARMAX(F45)=1. ;VARMIN(F45)=0.
VARMAX(F46)=1. ;VARMIN(F46)=0.
VARMAX(F47)=1. ;VARMIN(F47)=0.
VARMAX(F48)=1. ;VARMIN(F48)=0.
VARMAX(F49)=1. ;VARMIN(F49)=0.
VARMAX(F50)=1. ;VARMIN(F50)=0.
VARMAX(F51)=1. ;VARMIN(F51)=0.
VARMAX(F52)=1. ;VARMIN(F52)=0.
VARMAX(F53)=1. ;VARMIN(F53)=0.
VARMAX(F54)=1. ;VARMIN(F54)=0.
VARMAX(F55)=1. ;VARMIN(F55)=0.
VARMAX(F56)=1. ;VARMIN(F56)=0.
VARMAX(F57)=1. ;VARMIN(F57)=0.
VARMAX(F58)=1. ;VARMIN(F58)=0.
VARMAX(F59)=1. ;VARMIN(F59)=0.
VARMAX(F60)=1. ;VARMIN(F60)=0.
VARMAX(F61)=1. ;VARMIN(F61)=0.
VARMAX(F62)=1. ;VARMIN(F62)=0.
VARMAX(F63)=1. ;VARMIN(F63)=0.
VARMAX(F64)=1. ;VARMIN(F64)=0.
VARMAX(F65)=1. ;VARMIN(F65)=0.
VARMAX(F66)=1. ;VARMIN(F66)=0.
VARMAX(F67)=1. ;VARMIN(F67)=0.
VARMAX(F68)=1. ;VARMIN(F68)=0.
VARMAX(F69)=1. ;VARMIN(F69)=0.
VARMAX(F70)=1. ;VARMIN(F70)=0.
VARMAX(F71)=1. ;VARMIN(F71)=0.
VARMAX(F72)=1. ;VARMIN(F72)=0.
VARMAX(F73)=1. ;VARMIN(F73)=0.
VARMAX(F74)=1. ;VARMIN(F74)=0.
VARMAX(F75)=1. ;VARMIN(F75)=0.
VARMAX(F76)=1. ;VARMIN(F76)=0.
VARMAX(F77)=1. ;VARMIN(F77)=0.
VARMAX(F78)=1. ;VARMIN(F78)=0.
VARMAX(F79)=1. ;VARMIN(F79)=0.
VARMAX(F80)=1. ;VARMIN(F80)=0.
VARMAX(F81)=1. ;VARMIN(F81)=0.
VARMAX(F82)=1. ;VARMIN(F82)=0.
VARMAX(F83)=1. ;VARMIN(F83)=0.
VARMAX(F84)=1. ;VARMIN(F84)=0.
VARMAX(F85)=1. ;VARMIN(F85)=0.
VARMAX(F86)=1. ;VARMIN(F86)=0.
VARMAX(F87)=1. ;VARMIN(F87)=0.
VARMAX(F88)=1. ;VARMIN(F88)=0.
VARMAX(F89)=1. ;VARMIN(F89)=0.
VARMAX(F90)=1. ;VARMIN(F90)=0.
VARMAX(F91)=1. ;VARMIN(F91)=0.
VARMAX(F92)=1. ;VARMIN(F92)=0.
VARMAX(F93)=1. ;VARMIN(F93)=0.
VARMAX(F94)=1. ;VARMIN(F94)=0.
VARMAX(F95)=1. ;VARMIN(F95)=0.
VARMAX(F96)=1. ;VARMIN(F96)=0.
VARMAX(F97)=1. ;VARMIN(F97)=0.
VARMAX(F98)=1. ;VARMIN(F98)=0.
VARMAX(F99)=1. ;VARMIN(F99)=0.
VARMAX(F100)=1. ;VARMIN(F100)=0.
VARMAX(AVER)=1. ;VARMIN(AVER)=0.
VARMAX(AVEF)=1. ;VARMIN(AVEF)=0.
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,MFM)
SPEDAT(SET,MFM,NFLUIDS,I,100)
SPEDAT(SET,MFM,NFLR,I,10)
SPEDAT(SET,MFM,NFLF,I,10)
SPEDAT(SET,MFM,CONREA,R,10.)
SPEDAT(SET,MFM,CONMIX,R,10.)
SPEDAT(SET,MFM,RATEXP,R,5.)
SPEDAT(SET,MFM,HIGHF,R,0.7)
SPEDAT(SET,MFM,LOWF,R,0.4)
SPEDAT(SET,MFM,VISCON,R,0.)
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=1. ;EX(F1)=0.09292
EX(F2)=0.01115 ;EX(F3)=1.955E-03
EX(F4)=1.797E-04 ;EX(F5)=9.862E-06
EX(F6)=2.526E-07 ;EX(F7)=0.
EX(F8)=0. ;EX(F9)=0.
EX(F10)=0. ;EX(F11)=0.
EX(F12)=0.01972 ;EX(F13)=0.02089
EX(F14)=0.01108 ;EX(F15)=3.299E-03
EX(F16)=3.901E-04 ;EX(F17)=1.864E-05
EX(F18)=0. ;EX(F19)=0.
EX(F20)=0. ;EX(F21)=0.
EX(F22)=0. ;EX(F23)=0.01117
EX(F24)=0.02233 ;EX(F25)=0.02312
EX(F26)=0.01431 ;EX(F27)=5.153E-03
EX(F28)=6.208E-04 ;EX(F29)=0.
EX(F30)=0. ;EX(F31)=0.
EX(F32)=0. ;EX(F33)=0.
EX(F34)=7.185E-03 ;EX(F35)=0.01668
EX(F36)=0.02581 ;EX(F37)=0.0262
EX(F38)=0.01866 ;EX(F39)=0.01015
EX(F40)=0. ;EX(F41)=0.
EX(F42)=0. ;EX(F43)=0.
EX(F44)=0. ;EX(F45)=8.458E-05
EX(F46)=1.026E-04 ;EX(F47)=1.217E-04
EX(F48)=1.505E-04 ;EX(F49)=2.286E-04
EX(F50)=0.1418 ;EX(F51)=0.
EX(F52)=0. ;EX(F53)=0.
EX(F54)=0. ;EX(F55)=0.
EX(F56)=3.179E-05 ;EX(F57)=4.554E-05
EX(F58)=6.335E-05 ;EX(F59)=1.197E-04
EX(F60)=0.1593 ;EX(F61)=0.
EX(F62)=0. ;EX(F63)=0.
EX(F64)=0. ;EX(F65)=0.
EX(F66)=0. ;EX(F67)=1.826E-05
EX(F68)=3.399E-05 ;EX(F69)=7.733E-05
EX(F70)=0.1226 ;EX(F71)=0.
EX(F72)=0. ;EX(F73)=0.
EX(F74)=0. ;EX(F75)=0.
EX(F76)=0. ;EX(F77)=0.
EX(F78)=6.367E-03 ;EX(F79)=0.01604
EX(F80)=0.058 ;EX(F81)=0.
EX(F82)=0. ;EX(F83)=0.
EX(F84)=0. ;EX(F85)=0.
EX(F86)=0. ;EX(F87)=0.
EX(F88)=0. ;EX(F89)=0.01153
EX(F90)=0.04496 ;EX(F91)=0.
EX(F92)=0. ;EX(F93)=0.
EX(F94)=0. ;EX(F95)=0.
EX(F96)=0. ;EX(F97)=0.
EX(F98)=0. ;EX(F99)=0.
EX(F100)=0.09541 ;EX(AVER)=0.7634
EX(AVEF)=0.5 ;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(P1,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(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(F41,Y,N,Y,Y,Y,Y)
OUTPUT(F42,Y,N,Y,Y,Y,Y)
OUTPUT(F43,Y,N,Y,Y,Y,Y)
OUTPUT(F44,Y,N,Y,Y,Y,Y)
OUTPUT(F45,Y,N,Y,Y,Y,Y)
OUTPUT(F46,Y,N,Y,Y,Y,Y)
OUTPUT(F47,Y,N,Y,Y,Y,Y)
OUTPUT(F48,Y,N,Y,Y,Y,Y)
OUTPUT(F49,Y,N,Y,Y,Y,Y)
OUTPUT(F50,Y,N,Y,Y,Y,Y)
OUTPUT(F51,Y,N,Y,Y,Y,Y)
OUTPUT(F52,Y,N,Y,Y,Y,Y)
OUTPUT(F53,Y,N,Y,Y,Y,Y)
OUTPUT(F54,Y,N,Y,Y,Y,Y)
OUTPUT(F55,Y,N,Y,Y,Y,Y)
OUTPUT(F56,Y,N,Y,Y,Y,Y)
OUTPUT(F57,Y,N,Y,Y,Y,Y)
OUTPUT(F58,Y,N,Y,Y,Y,Y)
OUTPUT(F59,Y,N,Y,Y,Y,Y)
OUTPUT(F60,Y,N,Y,Y,Y,Y)
OUTPUT(F61,Y,N,Y,Y,Y,Y)
OUTPUT(F62,Y,N,Y,Y,Y,Y)
OUTPUT(F63,Y,N,Y,Y,Y,Y)
OUTPUT(F64,Y,N,Y,Y,Y,Y)
OUTPUT(F65,Y,N,Y,Y,Y,Y)
OUTPUT(F66,Y,N,Y,Y,Y,Y)
OUTPUT(F67,Y,N,Y,Y,Y,Y)
OUTPUT(F68,Y,N,Y,Y,Y,Y)
OUTPUT(F69,Y,N,Y,Y,Y,Y)
OUTPUT(F70,Y,N,Y,Y,Y,Y)
OUTPUT(F71,Y,N,Y,Y,Y,Y)
OUTPUT(F72,Y,N,Y,Y,Y,Y)
OUTPUT(F73,Y,N,Y,Y,Y,Y)
OUTPUT(F74,Y,N,Y,Y,Y,Y)
OUTPUT(F75,Y,N,Y,Y,Y,Y)
OUTPUT(F76,Y,N,Y,Y,Y,Y)
OUTPUT(F77,Y,N,Y,Y,Y,Y)
OUTPUT(F78,Y,N,Y,Y,Y,Y)
OUTPUT(F79,Y,N,Y,Y,Y,Y)
OUTPUT(F80,Y,N,Y,Y,Y,Y)
OUTPUT(F81,Y,N,Y,Y,Y,Y)
OUTPUT(F82,Y,N,Y,Y,Y,Y)
OUTPUT(F83,Y,N,Y,Y,Y,Y)
OUTPUT(F84,Y,N,Y,Y,Y,Y)
OUTPUT(F85,Y,N,Y,Y,Y,Y)
OUTPUT(F86,Y,N,Y,Y,Y,Y)
OUTPUT(F87,Y,N,Y,Y,Y,Y)
OUTPUT(F88,Y,N,Y,Y,Y,Y)
OUTPUT(F89,Y,N,Y,Y,Y,Y)
OUTPUT(F90,Y,N,Y,Y,Y,Y)
OUTPUT(F91,Y,N,Y,Y,Y,Y)
OUTPUT(F92,Y,N,Y,Y,Y,Y)
OUTPUT(F93,Y,N,Y,Y,Y,Y)
OUTPUT(F94,Y,N,Y,Y,Y,Y)
OUTPUT(F95,Y,N,Y,Y,Y,Y)
OUTPUT(F96,Y,N,Y,Y,Y,Y)
OUTPUT(F97,Y,N,Y,Y,Y,Y)
OUTPUT(F98,Y,N,Y,Y,Y,Y)
OUTPUT(F99,Y,N,Y,Y,Y,Y)
OUTPUT(F100,Y,N,Y,Y,Y,Y)
OUTPUT(AVER,Y,N,Y,N,N,N)
OUTPUT(AVEF,Y,N,Y,N,N,N)
OUTPUT(RATE,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;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
NTPRIN = 100000 ;ISTPRF = 1 ;ISTPRL = 100000
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
No PATCHes yet used for this Group
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP