TALK=T;RUN(1,1)
  DISPLAY
  Stirred reactor, with 1D (reactedness population), transient.

  The is initially full of reacted fluid, measured by F:NFLUIDS: .

  Completely-unreacted fluid, measured by F1, enters in the A
  stream; no fluid at all enters in the B stream.

  The final steady state is of interest; and, in particular,
  whether this consists almost completely of low-numbered fluids,
  F1, F2 etc, in which case the reactor is said to have been
  'extinguished', or whether significant mass fractions of
  higher-numbered fluids are to be found.

  The final value of the variable AVER  (i.e. average reactedness)
  is the simplest indicator. When extinction occurs, it tends
  asymptotically to zero; otherwise it maintains a finite value of
  between (say) 35 and 100 %.

  Its precise value depends upon the 'physical' constants CONMIX,
  CONREA, and RATEXP; and also on the 'numerical-grid-related'
  quantities, LSTEP and NFLR.
  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(AWSR transient, conmix=100. conrea=1.   )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 4
 ************************************************************
  Group 2. Time dependence
 STEADY = F
    * Set overall time and no. of steps
 TFIRST =0. ;TLAST =1.
 FSTEP = 1 ;LSTEP = 50
 TFRAC(1)=0.02 ;TFRAC(2)=0.04
 TFRAC(3)=0.06 ;TFRAC(4)=0.08
 TFRAC(5)=0.1 ;TFRAC(6)=0.12
 TFRAC(7)=0.14 ;TFRAC(8)=0.16
 TFRAC(9)=0.18 ;TFRAC(10)=0.2
 TFRAC(11)=0.22 ;TFRAC(12)=0.24
 TFRAC(13)=0.26 ;TFRAC(14)=0.28
 TFRAC(15)=0.3 ;TFRAC(16)=0.32
 TFRAC(17)=0.34 ;TFRAC(18)=0.36
 TFRAC(19)=0.38 ;TFRAC(20)=0.4
 TFRAC(21)=0.42 ;TFRAC(22)=0.44
 TFRAC(23)=0.46 ;TFRAC(24)=0.48
 TFRAC(25)=0.5 ;TFRAC(26)=0.52
 TFRAC(27)=0.54 ;TFRAC(28)=0.56
 TFRAC(29)=0.58 ;TFRAC(30)=0.6
 TFRAC(31)=0.62 ;TFRAC(32)=0.64
 TFRAC(33)=0.66 ;TFRAC(34)=0.68
 TFRAC(35)=0.7 ;TFRAC(36)=0.72
 TFRAC(37)=0.74 ;TFRAC(38)=0.76
 TFRAC(39)=0.78 ;TFRAC(40)=0.8
 TFRAC(41)=0.82 ;TFRAC(42)=0.84
 TFRAC(43)=0.86 ;TFRAC(44)=0.88
 TFRAC(45)=0.9 ;TFRAC(46)=0.92
 TFRAC(47)=0.94 ;TFRAC(48)=0.96
 TFRAC(49)=0.98 ;TFRAC(50)=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(47)=RCTD
 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)=MNSQ ;NAME(149)=AVER
 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(RCTD,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(MNSQ,Y,N,N,N,N,Y)
 SOLUTN(AVER,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(RCTD,N,Y,Y,Y,N,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(RCTD)=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(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(RCTD)=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(RCTD)=1.
 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(MNSQ)=1.0E-10 ;FIINIT(AVER)=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, 50)
 COVAL(INFLOWA ,P1 , FIXFLU ,1. )
 COVAL(INFLOWA ,F1 ,0. ,1. )
 
 PATCH(OUTFLOW ,VOLUME, 1, 1, 1, 1, 1, 1, 1, 50)
 COVAL(OUTFLOW ,P1 ,1. ,0. )
 
 PATCH(INFLOWB ,CELL , 1, 1, 1, 1, 1, 1, 1, 50)
 COVAL(INFLOWB ,F100,0. ,1. )
 
 PATCH(WHOLE ,VOLUME, 1, 1, 1, 1, 1, 1, 1, 50)
 COVAL(WHOLE ,RCTD,In-Form:source - see Grp 19)
 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(RCTD)=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(RCTD)=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(RCTD,LINRLX,0.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(MNSQ,LINRLX,1.)
 RELAX(AVER,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(RCTD)=1. ;VARMIN(RCTD)=0.
 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(MNSQ)=1. ;VARMIN(MNSQ)=0.
 VARMAX(AVER)=1. ;VARMIN(AVER)=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,100)
 SPEDAT(SET,MFM,NFLF,I,1)
 SPEDAT(SET,MFM,CONREA,R,1.)
 SPEDAT(SET,MFM,CONMIX,R,100.)
 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,TRANSIENT,L,F)
 SPEDAT(SET,SOURCE,RCTD!WHOLE,C,=42.*(1-RCTD)*RCTD^5.!LINE)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=1. ;EX(RCTD)=0.9735
 EX(F1)=0.02001 ;EX(F2)=5.298E-04
 EX(F3)=5.134E-04 ;EX(F4)=5.163E-04
 EX(F5)=5.226E-04 ;EX(F6)=5.305E-04
 EX(F7)=5.397E-04 ;EX(F8)=5.498E-04
 EX(F9)=5.607E-04 ;EX(F10)=5.723E-04
 EX(F11)=5.845E-04 ;EX(F12)=5.973E-04
 EX(F13)=6.107E-04 ;EX(F14)=6.247E-04
 EX(F15)=6.392E-04 ;EX(F16)=6.544E-04
 EX(F17)=6.701E-04 ;EX(F18)=6.865E-04
 EX(F19)=7.034E-04 ;EX(F20)=7.21E-04
 EX(F21)=7.391E-04 ;EX(F22)=7.579E-04
 EX(F23)=7.774E-04 ;EX(F24)=7.975E-04
 EX(F25)=8.183E-04 ;EX(F26)=8.398E-04
 EX(F27)=8.62E-04 ;EX(F28)=8.849E-04
 EX(F29)=9.085E-04 ;EX(F30)=9.329E-04
 EX(F31)=9.58E-04 ;EX(F32)=9.84E-04
 EX(F33)=1.011E-03 ;EX(F34)=1.038E-03
 EX(F35)=1.067E-03 ;EX(F36)=1.096E-03
 EX(F37)=1.126E-03 ;EX(F38)=1.157E-03
 EX(F39)=1.189E-03 ;EX(F40)=1.222E-03
 EX(F41)=1.256E-03 ;EX(F42)=1.292E-03
 EX(F43)=1.328E-03 ;EX(F44)=1.365E-03
 EX(F45)=1.404E-03 ;EX(F46)=1.444E-03
 EX(F47)=1.485E-03 ;EX(F48)=1.528E-03
 EX(F49)=1.572E-03 ;EX(F50)=1.618E-03
 EX(F51)=1.666E-03 ;EX(F52)=1.715E-03
 EX(F53)=1.767E-03 ;EX(F54)=1.82E-03
 EX(F55)=1.876E-03 ;EX(F56)=1.934E-03
 EX(F57)=1.995E-03 ;EX(F58)=2.059E-03
 EX(F59)=2.126E-03 ;EX(F60)=2.196E-03
 EX(F61)=2.269E-03 ;EX(F62)=2.347E-03
 EX(F63)=2.429E-03 ;EX(F64)=2.516E-03
 EX(F65)=2.608E-03 ;EX(F66)=2.706E-03
 EX(F67)=2.81E-03 ;EX(F68)=2.921E-03
 EX(F69)=3.04E-03 ;EX(F70)=3.167E-03
 EX(F71)=3.305E-03 ;EX(F72)=3.453E-03
 EX(F73)=3.613E-03 ;EX(F74)=3.788E-03
 EX(F75)=3.978E-03 ;EX(F76)=4.187E-03
 EX(F77)=4.416E-03 ;EX(F78)=4.669E-03
 EX(F79)=4.95E-03 ;EX(F80)=5.263E-03
 EX(F81)=5.614E-03 ;EX(F82)=6.011E-03
 EX(F83)=6.462E-03 ;EX(F84)=6.978E-03
 EX(F85)=7.575E-03 ;EX(F86)=8.272E-03
 EX(F87)=9.093E-03 ;EX(F88)=0.01007
 EX(F89)=0.01126 ;EX(F90)=0.01273
 EX(F91)=0.01458 ;EX(F92)=0.01696
 EX(F93)=0.02012 ;EX(F94)=0.02451
 EX(F95)=0.0309 ;EX(F96)=0.04093
 EX(F97)=0.05854 ;EX(F98)=0.09591
 EX(F99)=0.1887 ;EX(F100)=0.2724
 EX(MNSQ)=0.2086 ;EX(AVER)=0.8947
 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(RCTD,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(MNSQ,Y,N,Y,N,N,N)
 OUTPUT(AVER,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
 
 PATCH(PROFILES,PROFIL, 1, 1, 1, 1, 1, 1, 1, 50)
 PLOT(PROFILES,RCTD,0. ,1. )
 PLOT(PROFILES,F94 ,0. ,0.1 )
 PLOT(PROFILES,F96 ,0. ,0.1 )
 PLOT(PROFILES,F98 ,0. ,0.1 )
 PLOT(PROFILES,F100,0. ,1. )
 PLOT(PROFILES,AVER,0. ,1. )
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP