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