TALK=T;RUN(1,1)
  DISPLAY
  A hydrocarbon fuel enters steadily a duct of uniform cross-
  section steadily, through the porous wall of which is supplied a
  steady stream of air, at a low rate for the first 3/4 of the
  length, and at an increased rate over the last 1/4 .

  The air and fuel are supposed to react in a mixed-is-burned
  manner; but smoke is presumed to be generated at a rate which is
  proportional to the concentration of unburned fuel times the
  absolute temperature raised to the power SMOEXP.

  The multi-fuel model of turbulence is used, but with
  arbitrarily chosen values of the micro-mixing rate, CONMIX * RATE.

  The situation is discussed in the
  
  report by DB Spalding, dated 1998




  ENDDIS
  PHOTON USE
  AUTOPLOT
  file; phi 5

  cl; da 1; prod; da 1; oxid; da 1; fuel; scale; cola 1; col3 2
  colf 3
  msg  mass fractions of prod - yellow, oxid - blue, free fuel -red
  msg
  msg  flow is from left to right
  pause

  cl; da 1; smos; da 1; smom;  f11; scale; cola 1; col3 2;
  msg  mass fractions of smoke according to:
  msg
  msg  single-fluid - yellow;                   multi-fluid - blue
  pause

  cl; da 1; f1; da 1; f10; da 1; f20; scale; cola 1; col3 2; colf 3
  msg  mass fractions of f1 - yellow, f10 - blue, f20 - red
  pause

  cl; da 1; f1; da 1; f3; da 1; f5; da 1; f7; da 1; f9; da 1; f11;
      da 1; f13; da 1; f15; da 1; f17; da 1; f19; da 1; f20; scale;
      cola 1; col1 2; col1 3; col1 4; col1 5; col3 6; col1 7;
      col1 8; col1 9; col1 10; colf 11
  msg mass fractions of some intermediate fluids also
  pause

  cl; da 1; mnsq; cola 1
  msg root-mean-square fluctuations
  pause

  cl; da 1; tmp1; cola 1 ; da 1; tmpm; col3 2
  msg temperatures: single fluid - yellow; multi-fluid - blue
  pause

  cl; da 1; fuel; cola 1; da 1; frfm; col3 2;
  msg free fuel from single- (yellow) and multi-fluid models (blue)
  pause
  enduse
  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(MFM; 1dsmoke; conmix=5.                 )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 4
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 50
 XULAST =1.
 XFRAC(1)=0.02 ;XFRAC(2)=0.04
 XFRAC(3)=0.06 ;XFRAC(4)=0.08
 XFRAC(5)=0.1 ;XFRAC(6)=0.12
 XFRAC(7)=0.14 ;XFRAC(8)=0.16
 XFRAC(9)=0.18 ;XFRAC(10)=0.2
 XFRAC(11)=0.22 ;XFRAC(12)=0.24
 XFRAC(13)=0.26 ;XFRAC(14)=0.28
 XFRAC(15)=0.3 ;XFRAC(16)=0.32
 XFRAC(17)=0.34 ;XFRAC(18)=0.36
 XFRAC(19)=0.38 ;XFRAC(20)=0.4
 XFRAC(21)=0.42 ;XFRAC(22)=0.44
 XFRAC(23)=0.46 ;XFRAC(24)=0.48
 XFRAC(25)=0.5 ;XFRAC(26)=0.52
 XFRAC(27)=0.54 ;XFRAC(28)=0.56
 XFRAC(29)=0.58 ;XFRAC(30)=0.6
 XFRAC(31)=0.62 ;XFRAC(32)=0.64
 XFRAC(33)=0.66 ;XFRAC(34)=0.68
 XFRAC(35)=0.7 ;XFRAC(36)=0.72
 XFRAC(37)=0.74 ;XFRAC(38)=0.76
 XFRAC(39)=0.78 ;XFRAC(40)=0.8
 XFRAC(41)=0.82 ;XFRAC(42)=0.84
 XFRAC(43)=0.86 ;XFRAC(44)=0.88
 XFRAC(45)=0.9 ;XFRAC(46)=0.92
 XFRAC(47)=0.94 ;XFRAC(48)=0.96
 XFRAC(49)=0.98 ;XFRAC(50)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 1
 YVLAST =0.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(3)=U1
 NAME(14)=H1 ;NAME(115)=PRDM
 NAME(116)=OXDM ;NAME(117)=FRFM
 NAME(118)=TMPM ;NAME(119)=SMOS
 NAME(120)=SMOM ;NAME(121)=F1
 NAME(122)=F2 ;NAME(123)=F3
 NAME(124)=F4 ;NAME(125)=F5
 NAME(126)=F6 ;NAME(127)=F7
 NAME(128)=F8 ;NAME(129)=F9
 NAME(130)=F10 ;NAME(131)=F11
 NAME(132)=F12 ;NAME(133)=F13
 NAME(134)=F14 ;NAME(135)=F15
 NAME(136)=F16 ;NAME(137)=F17
 NAME(138)=F18 ;NAME(139)=F19
 NAME(140)=F20 ;NAME(141)=MNSQ
 NAME(142)=AVEF ;NAME(143)=SPH1
 NAME(144)=RHO1 ;NAME(145)=TMP1
 NAME(146)=PROD ;NAME(147)=OXID
 NAME(148)=FUEL ;NAME(149)=RATE
 NAME(150)=MIXF
    * 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(U1,Y,Y,N,N,N,Y)
 SOLUTN(H1,Y,Y,N,N,N,Y)
 SOLUTN(PRDM,Y,N,N,N,N,Y)
 SOLUTN(OXDM,Y,N,N,N,N,Y)
 SOLUTN(FRFM,Y,N,N,N,N,Y)
 SOLUTN(TMPM,Y,N,N,N,N,Y)
 SOLUTN(SMOS,Y,Y,N,N,N,Y)
 SOLUTN(SMOM,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(MNSQ,Y,N,N,N,N,Y)
 SOLUTN(AVEF,Y,N,N,N,N,Y)
 SOLUTN(SPH1,Y,N,N,N,N,Y)
 SOLUTN(RHO1,Y,N,N,N,N,Y)
 SOLUTN(TMP1,Y,N,N,N,N,Y)
 SOLUTN(PROD,Y,N,N,N,N,Y)
 SOLUTN(OXID,Y,N,N,N,N,Y)
 SOLUTN(FUEL,Y,N,N,N,N,Y)
 SOLUTN(RATE,Y,N,N,N,N,Y)
 SOLUTN(MIXF,Y,Y,N,N,N,Y)
 DEN1 = 144
 TEMP1 = 145
 ************************************************************
  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(U1,Y,Y,Y,Y,Y,Y)
 TERMS(H1,N,Y,Y,N,Y,N)
 TERMS(SMOS,N,Y,Y,Y,N,Y)
 TERMS(SMOM,N,Y,Y,Y,Y,Y)
 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(MIXF,N,Y,Y,Y,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWRH1 = T
 ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
 DENPCO = T
 ************************************************************
  Group 9. Properties used if PRPS is not
  stored, and where PRPS = -1.0 if it is!
 RHO1 = GRND6 ;TMP1 = GRND7
 EL1 =0.
 TSURR =0. ;TEMP0 =0.
 PRESS0 =1.0E+05
 DVO1DT =0. ;DRH1DP =0.
 RHO1A =16. ;RHO1B =29.
 RHO1C =28.
 TMP1A =1500. ;TMP1B =1500.
 TMP1C =1500.
 TMP2A =0.054825 ;TMP2B =4.9E+07
 TMP2C =1.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 ENUL =1.0E-05 ;ENUT =0.
 PRNDTL(U1)=1. ;PRNDTL(H1)=1.
 PRNDTL(SMOS)=1. ;PRNDTL(SMOM)=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(MIXF)=1.
 PRT(U1)=1. ;PRT(H1)=1.
 PRT(SMOS)=1. ;PRT(SMOM)=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(MIXF)=1.
 CP1 = GRND10 ;CP2 =1.
 CP1A =1500. ;CP1B =1500.
 CP1C =1500. ;CP1D =0.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1.0E-10 ;FIINIT(U1)=2.
 FIINIT(H1)=1.1595E+06 ;FIINIT(PRDM)=1.0E-10
 FIINIT(OXDM)=1.0E-10 ;FIINIT(FRFM)=1.0E-10
 FIINIT(TMPM)=1.0E-10 ;FIINIT(SMOS)=0.
 FIINIT(SMOM)=0. ;FIINIT(F1)=0.05
 FIINIT(F2)=0.05 ;FIINIT(F3)=0.05
 FIINIT(F4)=0.05 ;FIINIT(F5)=0.05
 FIINIT(F6)=0.05 ;FIINIT(F7)=0.05
 FIINIT(F8)=0.05 ;FIINIT(F9)=0.05
 FIINIT(F10)=0.05 ;FIINIT(F11)=0.05
 FIINIT(F12)=0.05 ;FIINIT(F13)=0.05
 FIINIT(F14)=0.05 ;FIINIT(F15)=0.05
 FIINIT(F16)=0.05 ;FIINIT(F17)=0.05
 FIINIT(F18)=0.05 ;FIINIT(F19)=0.05
 FIINIT(F20)=0.05 ;FIINIT(MNSQ)=1.0E-10
 FIINIT(AVEF)=1.0E-10 ;FIINIT(SPH1)=1.0E-10
 FIINIT(RHO1)=0.353169 ;FIINIT(TMP1)=1.0E-10
 FIINIT(PROD)=1.0E-10 ;FIINIT(OXID)=1.0E-10
 FIINIT(FUEL)=1. ;FIINIT(RATE)=100.
 FIINIT(MIXF)=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(FUIN ,WEST , 1, 1, 1, 1, 1, 1, 1, 1)
 COVAL(FUIN ,P1 , FIXFLU ,0.389704 )
 COVAL(FUIN ,U1 ,0. ,2. )
 COVAL(FUIN ,H1 ,0. ,5.01595E+07 )
 COVAL(FUIN ,F20 ,0. ,1. )
 COVAL(FUIN ,MIXF,0. ,1. )
 
 PATCH(OUTLET ,EAST , 50, 50, 1, 1, 1, 1, 1, 1)
 COVAL(OUTLET ,P1 ,1. ,0. )
 COVAL(OUTLET ,SMOS,0. , SAME )
 COVAL(OUTLET ,SMOM,0. , SAME )
 COVAL(OUTLET ,F1 ,0. , SAME )
 COVAL(OUTLET ,F2 ,0. , SAME )
 COVAL(OUTLET ,F3 ,0. , SAME )
 COVAL(OUTLET ,F4 ,0. , SAME )
 COVAL(OUTLET ,F5 ,0. , SAME )
 COVAL(OUTLET ,F6 ,0. , SAME )
 COVAL(OUTLET ,F7 ,0. , SAME )
 COVAL(OUTLET ,F8 ,0. , SAME )
 COVAL(OUTLET ,F9 ,0. , SAME )
 COVAL(OUTLET ,F10 ,0. , SAME )
 COVAL(OUTLET ,F11 ,0. , SAME )
 COVAL(OUTLET ,F12 ,0. , SAME )
 COVAL(OUTLET ,F13 ,0. , SAME )
 COVAL(OUTLET ,F14 ,0. , SAME )
 COVAL(OUTLET ,F15 ,0. , SAME )
 COVAL(OUTLET ,F16 ,0. , SAME )
 COVAL(OUTLET ,F17 ,0. , SAME )
 COVAL(OUTLET ,F18 ,0. , SAME )
 COVAL(OUTLET ,F19 ,0. , SAME )
 COVAL(OUTLET ,F20 ,0. , SAME )
 
 PATCH(AIRIN ,NORTH , 1, 50, 1, 1, 1, 1, 1, 1)
 COVAL(AIRIN ,P1 , FIXFLU ,0.988874 )
 COVAL(AIRIN ,H1 ,0. ,1.1595E+06 )
 COVAL(AIRIN ,F1 ,0. ,1. )
 
 PATCH(AIRIN2 ,NORTH , 37, 50, 1, 1, 1, 1, 1, 1)
 COVAL(AIRIN2 ,P1 , FIXFLU ,28.253557 )
 COVAL(AIRIN2 ,H1 ,0. ,1.1595E+06 )
 COVAL(AIRIN2 ,F1 ,0. ,1. )
 
 PATCH(SMOKE ,PHASEM, 1, 50, 1, 1, 1, 1, 1, 1)
 COVAL(SMOKE ,SMOS, FIXFLU , GRND )
 COVAL(SMOKE ,SMOM, FIXFLU , GRND )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 500 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-03
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(U1)=20
 LITER(H1)=20 ;LITER(SMOS)=20
 LITER(SMOM)=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(MIXF)=20
 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
 ENDIT(H1)=1.0E-03 ;ENDIT(SMOS)=1.0E-03
 ENDIT(SMOM)=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(MIXF)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,0.5)
 RELAX(U1,FALSDT,1.)
 RELAX(H1,FALSDT,10.)
 RELAX(PRDM,LINRLX,1.)
 RELAX(OXDM,LINRLX,1.)
 RELAX(FRFM,LINRLX,1.)
 RELAX(TMPM,LINRLX,1.)
 RELAX(SMOS,FALSDT,1.0E+09)
 RELAX(SMOM,FALSDT,1.0E+09)
 RELAX(F1,LINRLX,0.25)
 RELAX(F2,LINRLX,0.25)
 RELAX(F3,LINRLX,0.25)
 RELAX(F4,LINRLX,0.25)
 RELAX(F5,LINRLX,0.25)
 RELAX(F6,LINRLX,0.25)
 RELAX(F7,LINRLX,0.25)
 RELAX(F8,LINRLX,0.25)
 RELAX(F9,LINRLX,0.25)
 RELAX(F10,LINRLX,0.25)
 RELAX(F11,LINRLX,0.25)
 RELAX(F12,LINRLX,0.25)
 RELAX(F13,LINRLX,0.25)
 RELAX(F14,LINRLX,0.25)
 RELAX(F15,LINRLX,0.25)
 RELAX(F16,LINRLX,0.25)
 RELAX(F17,LINRLX,0.25)
 RELAX(F18,LINRLX,0.25)
 RELAX(F19,LINRLX,0.25)
 RELAX(F20,LINRLX,0.25)
 RELAX(MNSQ,LINRLX,1.)
 RELAX(AVEF,LINRLX,1.)
 RELAX(SPH1,LINRLX,1.)
 RELAX(RHO1,LINRLX,0.1)
 RELAX(TMP1,LINRLX,1.)
 RELAX(PROD,LINRLX,1.)
 RELAX(OXID,LINRLX,1.)
 RELAX(FUEL,LINRLX,1.)
 RELAX(RATE,LINRLX,1.)
 RELAX(MIXF,FALSDT,1.0E+09)
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
 VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
 VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10
 VARMAX(PRDM)=1.0E+10 ;VARMIN(PRDM)=-1.0E+10
 VARMAX(OXDM)=1.0E+10 ;VARMIN(OXDM)=-1.0E+10
 VARMAX(FRFM)=1.0E+10 ;VARMIN(FRFM)=-1.0E+10
 VARMAX(TMPM)=1.0E+10 ;VARMIN(TMPM)=-1.0E+10
 VARMAX(SMOS)=1.0E+10 ;VARMIN(SMOS)=-1.0E+10
 VARMAX(SMOM)=1.0E+10 ;VARMIN(SMOM)=-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(MNSQ)=1. ;VARMIN(MNSQ)=0.
 VARMAX(AVEF)=1. ;VARMIN(AVEF)=0.
 VARMAX(SPH1)=1.0E+10 ;VARMIN(SPH1)=-1.0E+10
 VARMAX(RHO1)=0.706339 ;VARMIN(RHO1)=0.035317
 VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=773.
 VARMAX(PROD)=1.0E+10 ;VARMIN(PROD)=-1.0E+10
 VARMAX(OXID)=1.0E+10 ;VARMIN(OXID)=-1.0E+10
 VARMAX(FUEL)=1.0E+10 ;VARMIN(FUEL)=-1.0E+10
 VARMAX(RATE)=1.0E+10 ;VARMIN(RATE)=-1.0E+10
 VARMAX(MIXF)=1.0E+10 ;VARMIN(MIXF)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 PARSOL = F
 ISG62 = 1
 SPEDAT(SET,MFM,FFUIN,R,1.)
 SPEDAT(SET,MFM,TAIRIN,R,773.)
 SPEDAT(SET,MFM,TFUIN,R,773.)
 SPEDAT(SET,MFM,TSTOIC,R,2606.314941)
 SPEDAT(SET,MFM,FACST,R,0.054825)
 SPEDAT(SET,MFM,GASCON,R,8314.)
 SPEDAT(SET,MFM,SPVLMX,R,100.)
 SPEDAT(SET,MFM,SPVLMN,R,1.0E-02)
 SPEDAT(SET,MFM,SMOCON,R,1.)
 SPEDAT(SET,MFM,SMOEXP,R,7.)
 SPEDAT(SET,MFM,MFMMOD,C,SCRS)
 SPEDAT(SET,MFM,NFLUIDS,I,20)
 SPEDAT(SET,MFM,NFLR,I,1)
 SPEDAT(SET,MFM,NFLF,I,20)
 SPEDAT(SET,MFM,CONREA,R,0.)
 SPEDAT(SET,MFM,CONMIX,R,5.)
 SPEDAT(SET,MFM,VISCON,R,0.)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=3.786E+04 ;EX(U1)=77.07
 EX(H1)=9.408E+06 ;EX(PRDM)=0.3831
 EX(OXDM)=0. ;EX(FRFM)=0.1397
 EX(TMPM)=1460. ;EX(SMOS)=3.547E-04
 EX(SMOM)=1.579E-04 ;EX(F1)=0.
 EX(F2)=0.1142 ;EX(F3)=0.07434
 EX(F4)=0.05543 ;EX(F5)=0.043
 EX(F6)=0.03415 ;EX(F7)=0.02767
 EX(F8)=0.02282 ;EX(F9)=0.01912
 EX(F10)=0.01626 ;EX(F11)=0.01401
 EX(F12)=0.01224 ;EX(F13)=0.01083
 EX(F14)=9.705E-03 ;EX(F15)=8.832E-03
 EX(F16)=8.191E-03 ;EX(F17)=7.813E-03
 EX(F18)=7.901E-03 ;EX(F19)=0.01022
 EX(F20)=0.03056 ;EX(MNSQ)=0.1773
 EX(AVEF)=0.1683 ;EX(SPH1)=1500.
 EX(RHO1)=0.2311 ;EX(TMP1)=1997.
 EX(PROD)=0.6835 ;EX(OXID)=0.1856
 EX(FUEL)=0.1309 ;EX(RATE)=100.
 EX(MIXF)=0.1683
 ************************************************************
  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(U1,Y,N,Y,Y,Y,Y)
 OUTPUT(H1,Y,N,Y,Y,Y,Y)
 OUTPUT(PRDM,Y,N,Y,N,N,N)
 OUTPUT(OXDM,Y,N,Y,N,N,N)
 OUTPUT(FRFM,Y,N,Y,N,N,N)
 OUTPUT(TMPM,Y,N,Y,N,N,N)
 OUTPUT(SMOS,Y,N,Y,Y,Y,Y)
 OUTPUT(SMOM,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(MNSQ,Y,N,Y,N,N,N)
 OUTPUT(AVEF,Y,N,Y,N,N,N)
 OUTPUT(SPH1,Y,N,Y,N,N,N)
 OUTPUT(RHO1,Y,N,Y,N,N,N)
 OUTPUT(TMP1,Y,N,Y,N,N,N)
 OUTPUT(PROD,Y,N,Y,N,N,N)
 OUTPUT(OXID,Y,N,Y,N,N,N)
 OUTPUT(FUEL,Y,N,Y,N,N,N)
 OUTPUT(RATE,Y,N,Y,N,N,N)
 OUTPUT(MIXF,Y,N,Y,Y,Y,Y)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 37 ;IYMON = 1 ;IZMON = 1
 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 1000 ;NUMCLS = 5
 NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
 IPLTF = 2 ;IPLTL = 1000 ;NPLT = 5
 ISWPRF = 1 ;ISWPRL = 100000
 ITABL = 3 ;IPROF = 1
 ABSIZ =0.5 ;ORSIZ =0.2
 NTZPRF = 1 ;NCOLPF = 50
 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 40
 
 PATCH(LONGPLOT,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(LONGPLOT,PROD,0. ,1. )
 PLOT(LONGPLOT,OXID,0. ,1. )
 PLOT(LONGPLOT,FUEL,0. ,1. )
 PLOT(LONGPLOT,MIXF,0. ,1. )
 
 PATCH(LONGPLT2,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(LONGPLT2,P1 ,0. ,0. )
 PLOT(LONGPLT2,U1 ,0. ,0. )
 PLOT(LONGPLT2,TMP1,0. ,0. )
 
 PATCH(FPROF1 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(FPROF1 ,F1 ,0. ,0. )
 
 PATCH(FPROF2 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(FPROF2 ,F2 ,0. ,0. )
 
 PATCH(FPROF3 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(FPROF3 ,F3 ,0. ,0. )
 
 PATCH(FPROF4 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(FPROF4 ,F4 ,0. ,0. )
 
 PATCH(FPROF5 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(FPROF5 ,F5 ,0. ,0. )
 
 PATCH(FPROF6 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(FPROF6 ,F6 ,0. ,0. )
 
 PATCH(FPROF7 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(FPROF7 ,F7 ,0. ,0. )
 
 PATCH(FPROF8 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(FPROF8 ,F8 ,0. ,0. )
 
 PATCH(FPROF9 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(FPROF9 ,F9 ,0. ,0. )
 
 PATCH(FPROF10 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(FPROF10 ,F10 ,0. ,0. )
 
 PATCH(FPROF11 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(FPROF11 ,F11 ,0. ,0. )
 
 PATCH(FPROF12 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(FPROF12 ,F12 ,0. ,0. )
 
 PATCH(FPROF13 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(FPROF13 ,F13 ,0. ,0. )
 
 PATCH(FPROF14 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(FPROF14 ,F14 ,0. ,0. )
 
 PATCH(FPROF15 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(FPROF15 ,F15 ,0. ,0. )
 
 PATCH(FPROF16 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(FPROF16 ,F16 ,0. ,0. )
 
 PATCH(FPROF17 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(FPROF17 ,F17 ,0. ,0. )
 
 PATCH(FPROF18 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(FPROF18 ,F18 ,0. ,0. )
 
 PATCH(FPROF19 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(FPROF19 ,F19 ,0. ,0. )
 
 PATCH(FPROF20 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(FPROF20 ,F20 ,0. ,0. )
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP