TALK=T;RUN(1,1)
  Display
  A turbulent combustible mixture flows steadily from left to right
  towards and into a hot porous plug, which is represented by
  fixing the mass fraction of the fully-reacted gas to unity there.

  A distorted distance scale is used, so that the whole flame lies
  between x=0 and x=1.

  The RATE value is fixed arbitrarily, as is also the diffusivity
  viscosity.

  The case can be used for testing the influence of the fineness
  of the population grid on the reaction-rate distribution.

  ENDDIS
  DISPLAY
  Notes:

  (1) four macros are provided in order reduce the repetiveness of
      MFM q1s.

  (2) they are named mfm, mfm1, mfm2 and mfm3, and are normally
      called in that order at the top of the q1 file, which then
      appears  as:

      #mfm
      settings of nflr, nflf, etc
      #mfm1
      settings of
      #mfm2
      #mfm3

  (3) the character variable mfm is declared and set in the
      always-loaded core-library macro 014.htm, which, in its turn,
      declares and sets the character variables mfm1, mfm2, mfm3.

  (4) the values are:
      mfm  = $L004
      mfm1 = $L003
      mfm2 = $L002
      mfm3 = $L001
  ENDDIS
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(1D Flame Propgn; Conmix=1. Conrea=0.025 )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 4
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 40
 XULAST =1.
 XFRAC(1)=0.025 ;XFRAC(2)=0.05
 XFRAC(3)=0.075 ;XFRAC(4)=0.1
 XFRAC(5)=0.125 ;XFRAC(6)=0.15
 XFRAC(7)=0.175 ;XFRAC(8)=0.2
 XFRAC(9)=0.225 ;XFRAC(10)=0.25
 XFRAC(11)=0.275 ;XFRAC(12)=0.3
 XFRAC(13)=0.325 ;XFRAC(14)=0.35
 XFRAC(15)=0.375 ;XFRAC(16)=0.4
 XFRAC(17)=0.425 ;XFRAC(18)=0.45
 XFRAC(19)=0.475 ;XFRAC(20)=0.5
 XFRAC(21)=0.525 ;XFRAC(22)=0.55
 XFRAC(23)=0.575 ;XFRAC(24)=0.6
 XFRAC(25)=0.625 ;XFRAC(26)=0.65
 XFRAC(27)=0.675 ;XFRAC(28)=0.7
 XFRAC(29)=0.725 ;XFRAC(30)=0.75
 XFRAC(31)=0.775 ;XFRAC(32)=0.8
 XFRAC(33)=0.825 ;XFRAC(34)=0.85
 XFRAC(35)=0.875 ;XFRAC(36)=0.9
 XFRAC(37)=0.925 ;XFRAC(38)=0.95
 XFRAC(39)=0.975 ;XFRAC(40)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 1
 YVLAST =1.
 YFRAC(1)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 1
 ZWLAST =1.
 ZFRAC(1)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(107)=F1 ;NAME(108)=F2
 NAME(109)=F3 ;NAME(110)=F4
 NAME(111)=F5 ;NAME(112)=F6
 NAME(113)=F7 ;NAME(114)=F8
 NAME(115)=F9 ;NAME(116)=F10
 NAME(117)=F11 ;NAME(118)=F12
 NAME(119)=F13 ;NAME(120)=F14
 NAME(121)=F15 ;NAME(122)=F16
 NAME(123)=F17 ;NAME(124)=F18
 NAME(125)=F19 ;NAME(126)=F20
 NAME(127)=F21 ;NAME(128)=F22
 NAME(129)=F23 ;NAME(130)=F24
 NAME(131)=F25 ;NAME(132)=F26
 NAME(133)=F27 ;NAME(134)=F28
 NAME(135)=F29 ;NAME(136)=F30
 NAME(137)=F31 ;NAME(138)=F32
 NAME(139)=F33 ;NAME(140)=F34
 NAME(141)=F35 ;NAME(142)=F36
 NAME(143)=F37 ;NAME(144)=F38
 NAME(145)=F39 ;NAME(146)=F40
 NAME(147)=MNSQ ;NAME(148)=AVER
 NAME(149)=REAC ;NAME(150)=RATE
    * Y in SOLUTN argument list denotes:
    * 1-stored 2-solved 3-whole-field
    * 4-point-by-point 5-explicit 6-harmonic averaging 
 SOLUTN(F1,Y,Y,N,N,N,Y)
 SOLUTN(F2,Y,Y,N,N,N,Y)
 SOLUTN(F3,Y,Y,N,N,N,Y)
 SOLUTN(F4,Y,Y,N,N,N,Y)
 SOLUTN(F5,Y,Y,N,N,N,Y)
 SOLUTN(F6,Y,Y,N,N,N,Y)
 SOLUTN(F7,Y,Y,N,N,N,Y)
 SOLUTN(F8,Y,Y,N,N,N,Y)
 SOLUTN(F9,Y,Y,N,N,N,Y)
 SOLUTN(F10,Y,Y,N,N,N,Y)
 SOLUTN(F11,Y,Y,N,N,N,Y)
 SOLUTN(F12,Y,Y,N,N,N,Y)
 SOLUTN(F13,Y,Y,N,N,N,Y)
 SOLUTN(F14,Y,Y,N,N,N,Y)
 SOLUTN(F15,Y,Y,N,N,N,Y)
 SOLUTN(F16,Y,Y,N,N,N,Y)
 SOLUTN(F17,Y,Y,N,N,N,Y)
 SOLUTN(F18,Y,Y,N,N,N,Y)
 SOLUTN(F19,Y,Y,N,N,N,Y)
 SOLUTN(F20,Y,Y,N,N,N,Y)
 SOLUTN(F21,Y,Y,N,N,N,Y)
 SOLUTN(F22,Y,Y,N,N,N,Y)
 SOLUTN(F23,Y,Y,N,N,N,Y)
 SOLUTN(F24,Y,Y,N,N,N,Y)
 SOLUTN(F25,Y,Y,N,N,N,Y)
 SOLUTN(F26,Y,Y,N,N,N,Y)
 SOLUTN(F27,Y,Y,N,N,N,Y)
 SOLUTN(F28,Y,Y,N,N,N,Y)
 SOLUTN(F29,Y,Y,N,N,N,Y)
 SOLUTN(F30,Y,Y,N,N,N,Y)
 SOLUTN(F31,Y,Y,N,N,N,Y)
 SOLUTN(F32,Y,Y,N,N,N,Y)
 SOLUTN(F33,Y,Y,N,N,N,Y)
 SOLUTN(F34,Y,Y,N,N,N,Y)
 SOLUTN(F35,Y,Y,N,N,N,Y)
 SOLUTN(F36,Y,Y,N,N,N,Y)
 SOLUTN(F37,Y,Y,N,N,N,Y)
 SOLUTN(F38,Y,Y,N,N,N,Y)
 SOLUTN(F39,Y,Y,N,N,N,Y)
 SOLUTN(F40,Y,Y,N,N,N,Y)
 SOLUTN(MNSQ,Y,N,N,N,N,Y)
 SOLUTN(AVER,Y,N,N,N,N,Y)
 SOLUTN(REAC,Y,N,N,N,N,Y)
 SOLUTN(RATE,Y,N,N,N,N,Y)
 ************************************************************
  Group 8. Terms & Devices
    * Y in TERMS argument list denotes:
    * 1-built-in source 2-convection 3-diffusion 4-transient
    * 5-first phase variable 6-interphase transport         
 TERMS(F1,N,Y,Y,Y,N,Y)
 TERMS(F2,N,Y,Y,Y,Y,Y)
 TERMS(F3,N,Y,Y,Y,N,Y)
 TERMS(F4,N,Y,Y,Y,Y,Y)
 TERMS(F5,N,Y,Y,Y,N,Y)
 TERMS(F6,N,Y,Y,Y,Y,Y)
 TERMS(F7,N,Y,Y,Y,N,Y)
 TERMS(F8,N,Y,Y,Y,Y,Y)
 TERMS(F9,N,Y,Y,Y,N,Y)
 TERMS(F10,N,Y,Y,Y,Y,Y)
 TERMS(F11,N,Y,Y,Y,N,Y)
 TERMS(F12,N,Y,Y,Y,Y,Y)
 TERMS(F13,N,Y,Y,Y,N,Y)
 TERMS(F14,N,Y,Y,Y,Y,Y)
 TERMS(F15,N,Y,Y,Y,N,Y)
 TERMS(F16,N,Y,Y,Y,Y,Y)
 TERMS(F17,N,Y,Y,Y,N,Y)
 TERMS(F18,N,Y,Y,Y,Y,Y)
 TERMS(F19,N,Y,Y,Y,N,Y)
 TERMS(F20,N,Y,Y,Y,Y,Y)
 TERMS(F21,N,Y,Y,Y,N,Y)
 TERMS(F22,N,Y,Y,Y,Y,Y)
 TERMS(F23,N,Y,Y,Y,N,Y)
 TERMS(F24,N,Y,Y,Y,Y,Y)
 TERMS(F25,N,Y,Y,Y,N,Y)
 TERMS(F26,N,Y,Y,Y,Y,Y)
 TERMS(F27,N,Y,Y,Y,N,Y)
 TERMS(F28,N,Y,Y,Y,Y,Y)
 TERMS(F29,N,Y,Y,Y,N,Y)
 TERMS(F30,N,Y,Y,Y,Y,Y)
 TERMS(F31,N,Y,Y,Y,N,Y)
 TERMS(F32,N,Y,Y,Y,Y,Y)
 TERMS(F33,N,Y,Y,Y,N,Y)
 TERMS(F34,N,Y,Y,Y,Y,Y)
 TERMS(F35,N,Y,Y,Y,N,Y)
 TERMS(F36,N,Y,Y,Y,Y,Y)
 TERMS(F37,N,Y,Y,Y,N,Y)
 TERMS(F38,N,Y,Y,Y,Y,Y)
 TERMS(F39,N,Y,Y,Y,N,Y)
 TERMS(F40,N,Y,Y,Y,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
 ************************************************************
  Group 9. Properties used if PRPS is not
  stored, and where PRPS = -1.0 if it is!
 RHO1 =1. ;TMP1 =0.
 EL1 =0.
 TSURR =0. ;TEMP0 =0.
 PRESS0 =0.
 DVO1DT =0. ;DRH1DP =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 ENUL =1. ;ENUT =0.
 PRNDTL(F1)=1. ;PRNDTL(F2)=1.
 PRNDTL(F3)=1. ;PRNDTL(F4)=1.
 PRNDTL(F5)=1. ;PRNDTL(F6)=1.
 PRNDTL(F7)=1. ;PRNDTL(F8)=1.
 PRNDTL(F9)=1. ;PRNDTL(F10)=1.
 PRNDTL(F11)=1. ;PRNDTL(F12)=1.
 PRNDTL(F13)=1. ;PRNDTL(F14)=1.
 PRNDTL(F15)=1. ;PRNDTL(F16)=1.
 PRNDTL(F17)=1. ;PRNDTL(F18)=1.
 PRNDTL(F19)=1. ;PRNDTL(F20)=1.
 PRNDTL(F21)=1. ;PRNDTL(F22)=1.
 PRNDTL(F23)=1. ;PRNDTL(F24)=1.
 PRNDTL(F25)=1. ;PRNDTL(F26)=1.
 PRNDTL(F27)=1. ;PRNDTL(F28)=1.
 PRNDTL(F29)=1. ;PRNDTL(F30)=1.
 PRNDTL(F31)=1. ;PRNDTL(F32)=1.
 PRNDTL(F33)=1. ;PRNDTL(F34)=1.
 PRNDTL(F35)=1. ;PRNDTL(F36)=1.
 PRNDTL(F37)=1. ;PRNDTL(F38)=1.
 PRNDTL(F39)=1. ;PRNDTL(F40)=1.
 PRT(F1)=1. ;PRT(F2)=1.
 PRT(F3)=1. ;PRT(F4)=1.
 PRT(F5)=1. ;PRT(F6)=1.
 PRT(F7)=1. ;PRT(F8)=1.
 PRT(F9)=1. ;PRT(F10)=1.
 PRT(F11)=1. ;PRT(F12)=1.
 PRT(F13)=1. ;PRT(F14)=1.
 PRT(F15)=1. ;PRT(F16)=1.
 PRT(F17)=1. ;PRT(F18)=1.
 PRT(F19)=1. ;PRT(F20)=1.
 PRT(F21)=1. ;PRT(F22)=1.
 PRT(F23)=1. ;PRT(F24)=1.
 PRT(F25)=1. ;PRT(F26)=1.
 PRT(F27)=1. ;PRT(F28)=1.
 PRT(F29)=1. ;PRT(F30)=1.
 PRT(F31)=1. ;PRT(F32)=1.
 PRT(F33)=1. ;PRT(F34)=1.
 PRT(F35)=1. ;PRT(F36)=1.
 PRT(F37)=1. ;PRT(F38)=1.
 PRT(F39)=1. ;PRT(F40)=1.
 CP1 =1. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(F1)=0. ;FIINIT(F2)=0.
 FIINIT(F3)=0. ;FIINIT(F4)=0.
 FIINIT(F5)=0. ;FIINIT(F6)=0.
 FIINIT(F7)=0. ;FIINIT(F8)=0.
 FIINIT(F9)=0. ;FIINIT(F10)=0.
 FIINIT(F11)=0. ;FIINIT(F12)=0.
 FIINIT(F13)=0. ;FIINIT(F14)=0.
 FIINIT(F15)=0. ;FIINIT(F16)=0.
 FIINIT(F17)=0. ;FIINIT(F18)=0.
 FIINIT(F19)=0. ;FIINIT(F20)=0.
 FIINIT(F21)=0. ;FIINIT(F22)=0.
 FIINIT(F23)=0. ;FIINIT(F24)=0.
 FIINIT(F25)=0. ;FIINIT(F26)=0.
 FIINIT(F27)=0. ;FIINIT(F28)=0.
 FIINIT(F29)=0. ;FIINIT(F30)=0.
 FIINIT(F31)=0. ;FIINIT(F32)=0.
 FIINIT(F33)=0. ;FIINIT(F34)=0.
 FIINIT(F35)=0. ;FIINIT(F36)=0.
 FIINIT(F37)=0. ;FIINIT(F38)=0.
 FIINIT(F39)=0. ;FIINIT(F40)=0.
 FIINIT(MNSQ)=1.0E-10 ;FIINIT(AVER)=1.0E-10
 FIINIT(REAC)=1.0E-10 ;FIINIT(RATE)=1.
   No PATCHes yet used for this Group
 INIADD = F
 FSWEEP = 1
 NAMFI =CHAM
 ************************************************************
  Group 12. Patchwise adjustment of terms
  Patches for this group are printed with those
  for Group 13.
  Their names begin either with GP12 or &
 ************************************************************
  Group 13. Boundary & Special Sources
 
 PATCH(OUTFLOW ,VOLUME, 1, 1, 1, 1, 1, 1, 1, 1)
 
 PATCH(COLD ,CELL , 1, 1, 1, 1, 1, 1, 1, 1)
 COVAL(COLD ,F1 , FIXVAL ,1. )
 COVAL(COLD ,F2 , FIXVAL ,0. )
 COVAL(COLD ,F3 , FIXVAL ,0. )
 COVAL(COLD ,F4 , FIXVAL ,0. )
 COVAL(COLD ,F5 , FIXVAL ,0. )
 COVAL(COLD ,F6 , FIXVAL ,0. )
 COVAL(COLD ,F7 , FIXVAL ,0. )
 COVAL(COLD ,F8 , FIXVAL ,0. )
 COVAL(COLD ,F9 , FIXVAL ,0. )
 COVAL(COLD ,F10 , FIXVAL ,0. )
 COVAL(COLD ,F11 , FIXVAL ,0. )
 COVAL(COLD ,F12 , FIXVAL ,0. )
 COVAL(COLD ,F13 , FIXVAL ,0. )
 COVAL(COLD ,F14 , FIXVAL ,0. )
 COVAL(COLD ,F15 , FIXVAL ,0. )
 COVAL(COLD ,F16 , FIXVAL ,0. )
 COVAL(COLD ,F17 , FIXVAL ,0. )
 COVAL(COLD ,F18 , FIXVAL ,0. )
 COVAL(COLD ,F19 , FIXVAL ,0. )
 COVAL(COLD ,F20 , FIXVAL ,0. )
 COVAL(COLD ,F21 , FIXVAL ,0. )
 COVAL(COLD ,F22 , FIXVAL ,0. )
 COVAL(COLD ,F23 , FIXVAL ,0. )
 COVAL(COLD ,F24 , FIXVAL ,0. )
 COVAL(COLD ,F25 , FIXVAL ,0. )
 COVAL(COLD ,F26 , FIXVAL ,0. )
 COVAL(COLD ,F27 , FIXVAL ,0. )
 COVAL(COLD ,F28 , FIXVAL ,0. )
 COVAL(COLD ,F29 , FIXVAL ,0. )
 COVAL(COLD ,F30 , FIXVAL ,0. )
 COVAL(COLD ,F31 , FIXVAL ,0. )
 COVAL(COLD ,F32 , FIXVAL ,0. )
 COVAL(COLD ,F33 , FIXVAL ,0. )
 COVAL(COLD ,F34 , FIXVAL ,0. )
 COVAL(COLD ,F35 , FIXVAL ,0. )
 COVAL(COLD ,F36 , FIXVAL ,0. )
 COVAL(COLD ,F37 , FIXVAL ,0. )
 COVAL(COLD ,F38 , FIXVAL ,0. )
 COVAL(COLD ,F39 , FIXVAL ,0. )
 COVAL(COLD ,F40 , FIXVAL ,0. )
 
 PATCH(HOT ,CELL , 40, 40, 1, 1, 1, 1, 1, 1)
 COVAL(HOT ,F1 , FIXVAL ,0. )
 COVAL(HOT ,F2 , FIXVAL ,0. )
 COVAL(HOT ,F3 , FIXVAL ,0. )
 COVAL(HOT ,F4 , FIXVAL ,0. )
 COVAL(HOT ,F5 , FIXVAL ,0. )
 COVAL(HOT ,F6 , FIXVAL ,0. )
 COVAL(HOT ,F7 , FIXVAL ,0. )
 COVAL(HOT ,F8 , FIXVAL ,0. )
 COVAL(HOT ,F9 , FIXVAL ,0. )
 COVAL(HOT ,F10 , FIXVAL ,0. )
 COVAL(HOT ,F11 , FIXVAL ,0. )
 COVAL(HOT ,F12 , FIXVAL ,0. )
 COVAL(HOT ,F13 , FIXVAL ,0. )
 COVAL(HOT ,F14 , FIXVAL ,0. )
 COVAL(HOT ,F15 , FIXVAL ,0. )
 COVAL(HOT ,F16 , FIXVAL ,0. )
 COVAL(HOT ,F17 , FIXVAL ,0. )
 COVAL(HOT ,F18 , FIXVAL ,0. )
 COVAL(HOT ,F19 , FIXVAL ,0. )
 COVAL(HOT ,F20 , FIXVAL ,0. )
 COVAL(HOT ,F21 , FIXVAL ,0. )
 COVAL(HOT ,F22 , FIXVAL ,0. )
 COVAL(HOT ,F23 , FIXVAL ,0. )
 COVAL(HOT ,F24 , FIXVAL ,0. )
 COVAL(HOT ,F25 , FIXVAL ,0. )
 COVAL(HOT ,F26 , FIXVAL ,0. )
 COVAL(HOT ,F27 , FIXVAL ,0. )
 COVAL(HOT ,F28 , FIXVAL ,0. )
 COVAL(HOT ,F29 , FIXVAL ,0. )
 COVAL(HOT ,F30 , FIXVAL ,0. )
 COVAL(HOT ,F31 , FIXVAL ,0. )
 COVAL(HOT ,F32 , FIXVAL ,0. )
 COVAL(HOT ,F33 , FIXVAL ,0. )
 COVAL(HOT ,F34 , FIXVAL ,0. )
 COVAL(HOT ,F35 , FIXVAL ,0. )
 COVAL(HOT ,F36 , FIXVAL ,0. )
 COVAL(HOT ,F37 , FIXVAL ,0. )
 COVAL(HOT ,F38 , FIXVAL ,0. )
 COVAL(HOT ,F39 , FIXVAL ,0. )
 COVAL(HOT ,F40 ,1000. ,1. )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 5000 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-20
 ************************************************************
  Group 16. Terminate Iterations
 LITER(F1)=20 ;LITER(F2)=20
 LITER(F3)=20 ;LITER(F4)=20
 LITER(F5)=20 ;LITER(F6)=20
 LITER(F7)=20 ;LITER(F8)=20
 LITER(F9)=20 ;LITER(F10)=20
 LITER(F11)=20 ;LITER(F12)=20
 LITER(F13)=20 ;LITER(F14)=20
 LITER(F15)=20 ;LITER(F16)=20
 LITER(F17)=20 ;LITER(F18)=20
 LITER(F19)=20 ;LITER(F20)=20
 LITER(F21)=20 ;LITER(F22)=20
 LITER(F23)=20 ;LITER(F24)=20
 LITER(F25)=20 ;LITER(F26)=20
 LITER(F27)=20 ;LITER(F28)=20
 LITER(F29)=20 ;LITER(F30)=20
 LITER(F31)=20 ;LITER(F32)=20
 LITER(F33)=20 ;LITER(F34)=20
 LITER(F35)=20 ;LITER(F36)=20
 LITER(F37)=20 ;LITER(F38)=20
 LITER(F39)=20 ;LITER(F40)=20
 ENDIT(F1)=1.0E-03 ;ENDIT(F2)=1.0E-03
 ENDIT(F3)=1.0E-03 ;ENDIT(F4)=1.0E-03
 ENDIT(F5)=1.0E-03 ;ENDIT(F6)=1.0E-03
 ENDIT(F7)=1.0E-03 ;ENDIT(F8)=1.0E-03
 ENDIT(F9)=1.0E-03 ;ENDIT(F10)=1.0E-03
 ENDIT(F11)=1.0E-03 ;ENDIT(F12)=1.0E-03
 ENDIT(F13)=1.0E-03 ;ENDIT(F14)=1.0E-03
 ENDIT(F15)=1.0E-03 ;ENDIT(F16)=1.0E-03
 ENDIT(F17)=1.0E-03 ;ENDIT(F18)=1.0E-03
 ENDIT(F19)=1.0E-03 ;ENDIT(F20)=1.0E-03
 ENDIT(F21)=1.0E-03 ;ENDIT(F22)=1.0E-03
 ENDIT(F23)=1.0E-03 ;ENDIT(F24)=1.0E-03
 ENDIT(F25)=1.0E-03 ;ENDIT(F26)=1.0E-03
 ENDIT(F27)=1.0E-03 ;ENDIT(F28)=1.0E-03
 ENDIT(F29)=1.0E-03 ;ENDIT(F30)=1.0E-03
 ENDIT(F31)=1.0E-03 ;ENDIT(F32)=1.0E-03
 ENDIT(F33)=1.0E-03 ;ENDIT(F34)=1.0E-03
 ENDIT(F35)=1.0E-03 ;ENDIT(F36)=1.0E-03
 ENDIT(F37)=1.0E-03 ;ENDIT(F38)=1.0E-03
 ENDIT(F39)=1.0E-03 ;ENDIT(F40)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(F1,LINRLX,0.5)
 RELAX(F2,LINRLX,0.5)
 RELAX(F3,LINRLX,0.5)
 RELAX(F4,LINRLX,0.5)
 RELAX(F5,LINRLX,0.5)
 RELAX(F6,LINRLX,0.5)
 RELAX(F7,LINRLX,0.5)
 RELAX(F8,LINRLX,0.5)
 RELAX(F9,LINRLX,0.5)
 RELAX(F10,LINRLX,0.5)
 RELAX(F11,LINRLX,0.5)
 RELAX(F12,LINRLX,0.5)
 RELAX(F13,LINRLX,0.5)
 RELAX(F14,LINRLX,0.5)
 RELAX(F15,LINRLX,0.5)
 RELAX(F16,LINRLX,0.5)
 RELAX(F17,LINRLX,0.5)
 RELAX(F18,LINRLX,0.5)
 RELAX(F19,LINRLX,0.5)
 RELAX(F20,LINRLX,0.5)
 RELAX(F21,LINRLX,0.5)
 RELAX(F22,LINRLX,0.5)
 RELAX(F23,LINRLX,0.5)
 RELAX(F24,LINRLX,0.5)
 RELAX(F25,LINRLX,0.5)
 RELAX(F26,LINRLX,0.5)
 RELAX(F27,LINRLX,0.5)
 RELAX(F28,LINRLX,0.5)
 RELAX(F29,LINRLX,0.5)
 RELAX(F30,LINRLX,0.5)
 RELAX(F31,LINRLX,0.5)
 RELAX(F32,LINRLX,0.5)
 RELAX(F33,LINRLX,0.5)
 RELAX(F34,LINRLX,0.5)
 RELAX(F35,LINRLX,0.5)
 RELAX(F36,LINRLX,0.5)
 RELAX(F37,LINRLX,0.5)
 RELAX(F38,LINRLX,0.5)
 RELAX(F39,LINRLX,0.5)
 RELAX(F40,LINRLX,0.5)
 RELAX(MNSQ,LINRLX,1.)
 RELAX(AVER,LINRLX,1.)
 RELAX(REAC,LINRLX,1.)
 RELAX(RATE,LINRLX,1.)
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(F1)=1. ;VARMIN(F1)=0.
 VARMAX(F2)=1. ;VARMIN(F2)=0.
 VARMAX(F3)=1. ;VARMIN(F3)=0.
 VARMAX(F4)=1. ;VARMIN(F4)=0.
 VARMAX(F5)=1. ;VARMIN(F5)=0.
 VARMAX(F6)=1. ;VARMIN(F6)=0.
 VARMAX(F7)=1. ;VARMIN(F7)=0.
 VARMAX(F8)=1. ;VARMIN(F8)=0.
 VARMAX(F9)=1. ;VARMIN(F9)=0.
 VARMAX(F10)=1. ;VARMIN(F10)=0.
 VARMAX(F11)=1. ;VARMIN(F11)=0.
 VARMAX(F12)=1. ;VARMIN(F12)=0.
 VARMAX(F13)=1. ;VARMIN(F13)=0.
 VARMAX(F14)=1. ;VARMIN(F14)=0.
 VARMAX(F15)=1. ;VARMIN(F15)=0.
 VARMAX(F16)=1. ;VARMIN(F16)=0.
 VARMAX(F17)=1. ;VARMIN(F17)=0.
 VARMAX(F18)=1. ;VARMIN(F18)=0.
 VARMAX(F19)=1. ;VARMIN(F19)=0.
 VARMAX(F20)=1. ;VARMIN(F20)=0.
 VARMAX(F21)=1. ;VARMIN(F21)=0.
 VARMAX(F22)=1. ;VARMIN(F22)=0.
 VARMAX(F23)=1. ;VARMIN(F23)=0.
 VARMAX(F24)=1. ;VARMIN(F24)=0.
 VARMAX(F25)=1. ;VARMIN(F25)=0.
 VARMAX(F26)=1. ;VARMIN(F26)=0.
 VARMAX(F27)=1. ;VARMIN(F27)=0.
 VARMAX(F28)=1. ;VARMIN(F28)=0.
 VARMAX(F29)=1. ;VARMIN(F29)=0.
 VARMAX(F30)=1. ;VARMIN(F30)=0.
 VARMAX(F31)=1. ;VARMIN(F31)=0.
 VARMAX(F32)=1. ;VARMIN(F32)=0.
 VARMAX(F33)=1. ;VARMIN(F33)=0.
 VARMAX(F34)=1. ;VARMIN(F34)=0.
 VARMAX(F35)=1. ;VARMIN(F35)=0.
 VARMAX(F36)=1. ;VARMIN(F36)=0.
 VARMAX(F37)=1. ;VARMIN(F37)=0.
 VARMAX(F38)=1. ;VARMIN(F38)=0.
 VARMAX(F39)=1. ;VARMIN(F39)=0.
 VARMAX(F40)=1. ;VARMIN(F40)=0.
 VARMAX(MNSQ)=1. ;VARMIN(MNSQ)=0.
 VARMAX(AVER)=1. ;VARMIN(AVER)=0.
 VARMAX(REAC)=1.0E+10 ;VARMIN(REAC)=-1.0E+10
 VARMAX(RATE)=1.0E+10 ;VARMIN(RATE)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 PARSOL = F
 ISG62 = 1
 SPEDAT(SET,MFM,MFMMOD,C,XFL1D)
 SPEDAT(SET,MFM,NFLUIDS,I,40)
 SPEDAT(SET,MFM,NFLR,I,40)
 SPEDAT(SET,MFM,NFLF,I,1)
 SPEDAT(SET,MFM,CONREA,R,0.025)
 SPEDAT(SET,MFM,CONMIX,R,1.)
 SPEDAT(SET,MFM,RATEXP,R,5.)
 SPEDAT(SET,MFM,HIGHF,R,1.)
 SPEDAT(SET,MFM,LOWF,R,0.)
 SPEDAT(SET,MFM,VISCON,R,0.)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(F1)=0.4197 ;EX(F2)=0.01603
 EX(F3)=7.67E-03 ;EX(F4)=6.074E-03
 EX(F5)=5.199E-03 ;EX(F6)=4.627E-03
 EX(F7)=4.217E-03 ;EX(F8)=3.905E-03
 EX(F9)=3.656E-03 ;EX(F10)=3.452E-03
 EX(F11)=3.281E-03 ;EX(F12)=3.133E-03
 EX(F13)=3.003E-03 ;EX(F14)=2.888E-03
 EX(F15)=2.785E-03 ;EX(F16)=2.692E-03
 EX(F17)=2.606E-03 ;EX(F18)=2.529E-03
 EX(F19)=2.457E-03 ;EX(F20)=2.392E-03
 EX(F21)=2.333E-03 ;EX(F22)=2.279E-03
 EX(F23)=2.231E-03 ;EX(F24)=2.188E-03
 EX(F25)=2.152E-03 ;EX(F26)=2.122E-03
 EX(F27)=2.098E-03 ;EX(F28)=2.082E-03
 EX(F29)=2.075E-03 ;EX(F30)=2.077E-03
 EX(F31)=2.09E-03 ;EX(F32)=2.117E-03
 EX(F33)=2.162E-03 ;EX(F34)=2.229E-03
 EX(F35)=2.328E-03 ;EX(F36)=2.475E-03
 EX(F37)=2.701E-03 ;EX(F38)=3.078E-03
 EX(F39)=4.043E-03 ;EX(F40)=0.4528
 EX(MNSQ)=0. ;EX(AVER)=0.
 EX(REAC)=2.263E-03 ;EX(RATE)=1.
 ************************************************************
  Group 21. Print-out of Variables
 INIFLD = F ;SUBWGR = F
    * Y in OUTPUT argument list denotes:
    * 1-field 2-correction-eq. monitor 3-selective dumping      
    * 4-whole-field residual 5-spot-value table 6-residual table
 OUTPUT(F1,Y,N,Y,Y,Y,Y)
 OUTPUT(F2,Y,N,Y,Y,Y,Y)
 OUTPUT(F3,Y,N,Y,Y,Y,Y)
 OUTPUT(F4,Y,N,Y,Y,Y,Y)
 OUTPUT(F5,Y,N,Y,Y,Y,Y)
 OUTPUT(F6,Y,N,Y,Y,Y,Y)
 OUTPUT(F7,Y,N,Y,Y,Y,Y)
 OUTPUT(F8,Y,N,Y,Y,Y,Y)
 OUTPUT(F9,Y,N,Y,Y,Y,Y)
 OUTPUT(F10,Y,N,Y,Y,Y,Y)
 OUTPUT(F11,Y,N,Y,Y,Y,Y)
 OUTPUT(F12,Y,N,Y,Y,Y,Y)
 OUTPUT(F13,Y,N,Y,Y,Y,Y)
 OUTPUT(F14,Y,N,Y,Y,Y,Y)
 OUTPUT(F15,Y,N,Y,Y,Y,Y)
 OUTPUT(F16,Y,N,Y,Y,Y,Y)
 OUTPUT(F17,Y,N,Y,Y,Y,Y)
 OUTPUT(F18,Y,N,Y,Y,Y,Y)
 OUTPUT(F19,Y,N,Y,Y,Y,Y)
 OUTPUT(F20,Y,N,Y,Y,Y,Y)
 OUTPUT(F21,Y,N,Y,Y,Y,Y)
 OUTPUT(F22,Y,N,Y,Y,Y,Y)
 OUTPUT(F23,Y,N,Y,Y,Y,Y)
 OUTPUT(F24,Y,N,Y,Y,Y,Y)
 OUTPUT(F25,Y,N,Y,Y,Y,Y)
 OUTPUT(F26,Y,N,Y,Y,Y,Y)
 OUTPUT(F27,Y,N,Y,Y,Y,Y)
 OUTPUT(F28,Y,N,Y,Y,Y,Y)
 OUTPUT(F29,Y,N,Y,Y,Y,Y)
 OUTPUT(F30,Y,N,Y,Y,Y,Y)
 OUTPUT(F31,Y,N,Y,Y,Y,Y)
 OUTPUT(F32,Y,N,Y,Y,Y,Y)
 OUTPUT(F33,Y,N,Y,Y,Y,Y)
 OUTPUT(F34,Y,N,Y,Y,Y,Y)
 OUTPUT(F35,Y,N,Y,Y,Y,Y)
 OUTPUT(F36,Y,N,Y,Y,Y,Y)
 OUTPUT(F37,Y,N,Y,Y,Y,Y)
 OUTPUT(F38,Y,N,Y,Y,Y,Y)
 OUTPUT(F39,Y,N,Y,Y,Y,Y)
 OUTPUT(F40,Y,N,Y,Y,Y,Y)
 OUTPUT(MNSQ,Y,N,Y,N,N,N)
 OUTPUT(AVER,Y,N,Y,N,N,N)
 OUTPUT(REAC,Y,N,Y,N,N,N)
 OUTPUT(RATE,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 20 ;IYMON = 1 ;IZMON = 1
 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
 IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
 ISWPRF = 1 ;ISWPRL = 100000
 ITABL = 3 ;IPROF = 1
 ABSIZ =0.5 ;ORSIZ =0.4
 NTZPRF = 1 ;NCOLPF = 50
 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
 
 PATCH(AVES ,PROFIL, 1, 40, 1, 1, 1, 1, 1, 1)
 PLOT(AVES ,REAC,0. ,0. )
 
 PATCH(EXT ,PROFIL, 1, 40, 1, 1, 1, 1, 1, 1)
 PLOT(EXT ,F1 ,0. ,1. )
 PLOT(EXT ,F40 ,0. ,0. )
 
 PATCH(INT ,PROFIL, 1, 40, 1, 1, 1, 1, 1, 1)
 PLOT(INT ,F2 ,0. ,0.03 )
 PLOT(INT ,F11 ,0. ,0.03 )
 PLOT(INT ,F20 ,0. ,0.03 )
 PLOT(INT ,F29 ,0. ,0.03 )
 PLOT(INT ,F38 ,0. ,0.03 )
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP