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
  DISPLAY
  This case extends case
  230
  by treating the temperatures and smoke concentrations as
  continuously-varying attributes.

  Their names are T1C, T2C etc and S1C, S2C etc
    
  The scalar variable RADFACT is introduced in order that the
  effects of radiative heat losses associated with smoke content
  may be qualitatively simulated.

  ENDDIS
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(MFM; 1dsmoke; conmix=5.                 )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 230
 ************************************************************
  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(75)=T1C
 NAME(76)=T2C ;NAME(77)=T3C
 NAME(78)=T4C ;NAME(79)=T5C
 NAME(80)=T6C ;NAME(81)=T7C
 NAME(82)=T8C ;NAME(83)=T9C
 NAME(84)=T10C ;NAME(85)=T11C
 NAME(86)=T12C ;NAME(87)=T13C
 NAME(88)=T14C ;NAME(89)=T15C
 NAME(90)=T16C ;NAME(91)=T17C
 NAME(92)=T18C ;NAME(93)=T19C
 NAME(94)=T20C ;NAME(95)=S1C
 NAME(96)=S2C ;NAME(97)=S3C
 NAME(98)=S4C ;NAME(99)=S5C
 NAME(100)=S6C ;NAME(101)=S7C
 NAME(102)=S8C ;NAME(103)=S9C
 NAME(104)=S10C ;NAME(105)=S11C
 NAME(106)=S12C ;NAME(107)=S13C
 NAME(108)=S14C ;NAME(109)=S15C
 NAME(110)=S16C ;NAME(111)=S17C
 NAME(112)=S18C ;NAME(113)=S19C
 NAME(114)=S20C ;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(T1C,Y,Y,N,N,N,Y)
 SOLUTN(T2C,Y,Y,N,N,N,Y)
 SOLUTN(T3C,Y,Y,N,N,N,Y)
 SOLUTN(T4C,Y,Y,N,N,N,Y)
 SOLUTN(T5C,Y,Y,N,N,N,Y)
 SOLUTN(T6C,Y,Y,N,N,N,Y)
 SOLUTN(T7C,Y,Y,N,N,N,Y)
 SOLUTN(T8C,Y,Y,N,N,N,Y)
 SOLUTN(T9C,Y,Y,N,N,N,Y)
 SOLUTN(T10C,Y,Y,N,N,N,Y)
 SOLUTN(T11C,Y,Y,N,N,N,Y)
 SOLUTN(T12C,Y,Y,N,N,N,Y)
 SOLUTN(T13C,Y,Y,N,N,N,Y)
 SOLUTN(T14C,Y,Y,N,N,N,Y)
 SOLUTN(T15C,Y,Y,N,N,N,Y)
 SOLUTN(T16C,Y,Y,N,N,N,Y)
 SOLUTN(T17C,Y,Y,N,N,N,Y)
 SOLUTN(T18C,Y,Y,N,N,N,Y)
 SOLUTN(T19C,Y,Y,N,N,N,Y)
 SOLUTN(T20C,Y,Y,N,N,N,Y)
 SOLUTN(S1C,Y,Y,N,N,N,Y)
 SOLUTN(S2C,Y,Y,N,N,N,Y)
 SOLUTN(S3C,Y,Y,N,N,N,Y)
 SOLUTN(S4C,Y,Y,N,N,N,Y)
 SOLUTN(S5C,Y,Y,N,N,N,Y)
 SOLUTN(S6C,Y,Y,N,N,N,Y)
 SOLUTN(S7C,Y,Y,N,N,N,Y)
 SOLUTN(S8C,Y,Y,N,N,N,Y)
 SOLUTN(S9C,Y,Y,N,N,N,Y)
 SOLUTN(S10C,Y,Y,N,N,N,Y)
 SOLUTN(S11C,Y,Y,N,N,N,Y)
 SOLUTN(S12C,Y,Y,N,N,N,Y)
 SOLUTN(S13C,Y,Y,N,N,N,Y)
 SOLUTN(S14C,Y,Y,N,N,N,Y)
 SOLUTN(S15C,Y,Y,N,N,N,Y)
 SOLUTN(S16C,Y,Y,N,N,N,Y)
 SOLUTN(S17C,Y,Y,N,N,N,Y)
 SOLUTN(S18C,Y,Y,N,N,N,Y)
 SOLUTN(S19C,Y,Y,N,N,N,Y)
 SOLUTN(S20C,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(T1C,N,Y,Y,Y,N,Y)
 TERMS(T2C,N,Y,Y,Y,Y,Y)
 TERMS(T3C,N,Y,Y,Y,N,Y)
 TERMS(T4C,N,Y,Y,Y,Y,Y)
 TERMS(T5C,N,Y,Y,Y,N,Y)
 TERMS(T6C,N,Y,Y,Y,Y,Y)
 TERMS(T7C,N,Y,Y,Y,N,Y)
 TERMS(T8C,N,Y,Y,Y,Y,Y)
 TERMS(T9C,N,Y,Y,Y,N,Y)
 TERMS(T10C,N,Y,Y,Y,Y,Y)
 TERMS(T11C,N,Y,Y,Y,N,Y)
 TERMS(T12C,N,Y,Y,Y,Y,Y)
 TERMS(T13C,N,Y,Y,Y,N,Y)
 TERMS(T14C,N,Y,Y,Y,Y,Y)
 TERMS(T15C,N,Y,Y,Y,N,Y)
 TERMS(T16C,N,Y,Y,Y,Y,Y)
 TERMS(T17C,N,Y,Y,Y,N,Y)
 TERMS(T18C,N,Y,Y,Y,Y,Y)
 TERMS(T19C,N,Y,Y,Y,N,Y)
 TERMS(T20C,N,Y,Y,Y,Y,Y)
 TERMS(S1C,N,Y,Y,Y,N,Y)
 TERMS(S2C,N,Y,Y,Y,Y,Y)
 TERMS(S3C,N,Y,Y,Y,N,Y)
 TERMS(S4C,N,Y,Y,Y,Y,Y)
 TERMS(S5C,N,Y,Y,Y,N,Y)
 TERMS(S6C,N,Y,Y,Y,Y,Y)
 TERMS(S7C,N,Y,Y,Y,N,Y)
 TERMS(S8C,N,Y,Y,Y,Y,Y)
 TERMS(S9C,N,Y,Y,Y,N,Y)
 TERMS(S10C,N,Y,Y,Y,Y,Y)
 TERMS(S11C,N,Y,Y,Y,N,Y)
 TERMS(S12C,N,Y,Y,Y,Y,Y)
 TERMS(S13C,N,Y,Y,Y,N,Y)
 TERMS(S14C,N,Y,Y,Y,Y,Y)
 TERMS(S15C,N,Y,Y,Y,N,Y)
 TERMS(S16C,N,Y,Y,Y,Y,Y)
 TERMS(S17C,N,Y,Y,Y,N,Y)
 TERMS(S18C,N,Y,Y,Y,Y,Y)
 TERMS(S19C,N,Y,Y,Y,N,Y)
 TERMS(S20C,N,Y,Y,Y,Y,Y)
 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(T1C)=1. ;PRNDTL(T2C)=1.
 PRNDTL(T3C)=1. ;PRNDTL(T4C)=1.
 PRNDTL(T5C)=1. ;PRNDTL(T6C)=1.
 PRNDTL(T7C)=1. ;PRNDTL(T8C)=1.
 PRNDTL(T9C)=1. ;PRNDTL(T10C)=1.
 PRNDTL(T11C)=1. ;PRNDTL(T12C)=1.
 PRNDTL(T13C)=1. ;PRNDTL(T14C)=1.
 PRNDTL(T15C)=1. ;PRNDTL(T16C)=1.
 PRNDTL(T17C)=1. ;PRNDTL(T18C)=1.
 PRNDTL(T19C)=1. ;PRNDTL(T20C)=1.
 PRNDTL(S1C)=1. ;PRNDTL(S2C)=1.
 PRNDTL(S3C)=1. ;PRNDTL(S4C)=1.
 PRNDTL(S5C)=1. ;PRNDTL(S6C)=1.
 PRNDTL(S7C)=1. ;PRNDTL(S8C)=1.
 PRNDTL(S9C)=1. ;PRNDTL(S10C)=1.
 PRNDTL(S11C)=1. ;PRNDTL(S12C)=1.
 PRNDTL(S13C)=1. ;PRNDTL(S14C)=1.
 PRNDTL(S15C)=1. ;PRNDTL(S16C)=1.
 PRNDTL(S17C)=1. ;PRNDTL(S18C)=1.
 PRNDTL(S19C)=1. ;PRNDTL(S20C)=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(T1C)=1. ;PRT(T2C)=1.
 PRT(T3C)=1. ;PRT(T4C)=1.
 PRT(T5C)=1. ;PRT(T6C)=1.
 PRT(T7C)=1. ;PRT(T8C)=1.
 PRT(T9C)=1. ;PRT(T10C)=1.
 PRT(T11C)=1. ;PRT(T12C)=1.
 PRT(T13C)=1. ;PRT(T14C)=1.
 PRT(T15C)=1. ;PRT(T16C)=1.
 PRT(T17C)=1. ;PRT(T18C)=1.
 PRT(T19C)=1. ;PRT(T20C)=1.
 PRT(S1C)=1. ;PRT(S2C)=1.
 PRT(S3C)=1. ;PRT(S4C)=1.
 PRT(S5C)=1. ;PRT(S6C)=1.
 PRT(S7C)=1. ;PRT(S8C)=1.
 PRT(S9C)=1. ;PRT(S10C)=1.
 PRT(S11C)=1. ;PRT(S12C)=1.
 PRT(S13C)=1. ;PRT(S14C)=1.
 PRT(S15C)=1. ;PRT(S16C)=1.
 PRT(S17C)=1. ;PRT(S18C)=1.
 PRT(S19C)=1. ;PRT(S20C)=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(T1C)=1.0E-10
 FIINIT(T2C)=1.0E-10 ;FIINIT(T3C)=1.0E-10
 FIINIT(T4C)=1.0E-10 ;FIINIT(T5C)=1.0E-10
 FIINIT(T6C)=1.0E-10 ;FIINIT(T7C)=1.0E-10
 FIINIT(T8C)=1.0E-10 ;FIINIT(T9C)=1.0E-10
 FIINIT(T10C)=1.0E-10 ;FIINIT(T11C)=1.0E-10
 FIINIT(T12C)=1.0E-10 ;FIINIT(T13C)=1.0E-10
 FIINIT(T14C)=1.0E-10 ;FIINIT(T15C)=1.0E-10
 FIINIT(T16C)=1.0E-10 ;FIINIT(T17C)=1.0E-10
 FIINIT(T18C)=1.0E-10 ;FIINIT(T19C)=1.0E-10
 FIINIT(T20C)=1.0E-10 ;FIINIT(S1C)=1.0E-10
 FIINIT(S2C)=1.0E-10 ;FIINIT(S3C)=1.0E-10
 FIINIT(S4C)=1.0E-10 ;FIINIT(S5C)=1.0E-10
 FIINIT(S6C)=1.0E-10 ;FIINIT(S7C)=1.0E-10
 FIINIT(S8C)=1.0E-10 ;FIINIT(S9C)=1.0E-10
 FIINIT(S10C)=1.0E-10 ;FIINIT(S11C)=1.0E-10
 FIINIT(S12C)=1.0E-10 ;FIINIT(S13C)=1.0E-10
 FIINIT(S14C)=1.0E-10 ;FIINIT(S15C)=1.0E-10
 FIINIT(S16C)=1.0E-10 ;FIINIT(S17C)=1.0E-10
 FIINIT(S18C)=1.0E-10 ;FIINIT(S19C)=1.0E-10
 FIINIT(S20C)=1.0E-10 ;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 ,S1C ,0. ,0. )
 COVAL(FUIN ,S2C ,0. ,0. )
 COVAL(FUIN ,S3C ,0. ,0. )
 COVAL(FUIN ,S4C ,0. ,0. )
 COVAL(FUIN ,S5C ,0. ,0. )
 COVAL(FUIN ,S6C ,0. ,0. )
 COVAL(FUIN ,S7C ,0. ,0. )
 COVAL(FUIN ,S8C ,0. ,0. )
 COVAL(FUIN ,S9C ,0. ,0. )
 COVAL(FUIN ,S10C,0. ,0. )
 COVAL(FUIN ,S11C,0. ,0. )
 COVAL(FUIN ,S12C,0. ,0. )
 COVAL(FUIN ,S13C,0. ,0. )
 COVAL(FUIN ,S14C,0. ,0. )
 COVAL(FUIN ,S15C,0. ,0. )
 COVAL(FUIN ,S16C,0. ,0. )
 COVAL(FUIN ,S17C,0. ,0. )
 COVAL(FUIN ,S18C,0. ,0. )
 COVAL(FUIN ,S19C,0. ,0. )
 COVAL(FUIN ,S20C,0. ,0. )
 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 ,S1C ,0. ,0. )
 COVAL(AIRIN ,S2C ,0. ,0. )
 COVAL(AIRIN ,S3C ,0. ,0. )
 COVAL(AIRIN ,S4C ,0. ,0. )
 COVAL(AIRIN ,S5C ,0. ,0. )
 COVAL(AIRIN ,S6C ,0. ,0. )
 COVAL(AIRIN ,S7C ,0. ,0. )
 COVAL(AIRIN ,S8C ,0. ,0. )
 COVAL(AIRIN ,S9C ,0. ,0. )
 COVAL(AIRIN ,S10C,0. ,0. )
 COVAL(AIRIN ,S11C,0. ,0. )
 COVAL(AIRIN ,S12C,0. ,0. )
 COVAL(AIRIN ,S13C,0. ,0. )
 COVAL(AIRIN ,S14C,0. ,0. )
 COVAL(AIRIN ,S15C,0. ,0. )
 COVAL(AIRIN ,S16C,0. ,0. )
 COVAL(AIRIN ,S17C,0. ,0. )
 COVAL(AIRIN ,S18C,0. ,0. )
 COVAL(AIRIN ,S19C,0. ,0. )
 COVAL(AIRIN ,S20C,0. ,0. )
 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 ,S1C ,0. ,0. )
 COVAL(AIRIN2 ,S2C ,0. ,0. )
 COVAL(AIRIN2 ,S3C ,0. ,0. )
 COVAL(AIRIN2 ,S4C ,0. ,0. )
 COVAL(AIRIN2 ,S5C ,0. ,0. )
 COVAL(AIRIN2 ,S6C ,0. ,0. )
 COVAL(AIRIN2 ,S7C ,0. ,0. )
 COVAL(AIRIN2 ,S8C ,0. ,0. )
 COVAL(AIRIN2 ,S9C ,0. ,0. )
 COVAL(AIRIN2 ,S10C,0. ,0. )
 COVAL(AIRIN2 ,S11C,0. ,0. )
 COVAL(AIRIN2 ,S12C,0. ,0. )
 COVAL(AIRIN2 ,S13C,0. ,0. )
 COVAL(AIRIN2 ,S14C,0. ,0. )
 COVAL(AIRIN2 ,S15C,0. ,0. )
 COVAL(AIRIN2 ,S16C,0. ,0. )
 COVAL(AIRIN2 ,S17C,0. ,0. )
 COVAL(AIRIN2 ,S18C,0. ,0. )
 COVAL(AIRIN2 ,S19C,0. ,0. )
 COVAL(AIRIN2 ,S20C,0. ,0. )
 COVAL(AIRIN2 ,F1 ,0. ,1. )
 
 PATCH(SMOKE ,PHASEM, 1, 50, 1, 1, 1, 1, 1, 1)
 COVAL(SMOKE ,T1C , GRND , GRND )
 COVAL(SMOKE ,T2C , GRND , GRND )
 COVAL(SMOKE ,T3C , GRND , GRND )
 COVAL(SMOKE ,T4C , GRND , GRND )
 COVAL(SMOKE ,T5C , GRND , GRND )
 COVAL(SMOKE ,T6C , GRND , GRND )
 COVAL(SMOKE ,T7C , GRND , GRND )
 COVAL(SMOKE ,T8C , GRND , GRND )
 COVAL(SMOKE ,T9C , GRND , GRND )
 COVAL(SMOKE ,T10C, GRND , GRND )
 COVAL(SMOKE ,T11C, GRND , GRND )
 COVAL(SMOKE ,T12C, GRND , GRND )
 COVAL(SMOKE ,T13C, GRND , GRND )
 COVAL(SMOKE ,T14C, GRND , GRND )
 COVAL(SMOKE ,T15C, GRND , GRND )
 COVAL(SMOKE ,T16C, GRND , GRND )
 COVAL(SMOKE ,T17C, GRND , GRND )
 COVAL(SMOKE ,T18C, GRND , GRND )
 COVAL(SMOKE ,T19C, GRND , GRND )
 COVAL(SMOKE ,T20C, GRND , GRND )
 COVAL(SMOKE ,S1C , GRND , GRND )
 COVAL(SMOKE ,S2C , GRND , GRND )
 COVAL(SMOKE ,S3C , GRND , GRND )
 COVAL(SMOKE ,S4C , GRND , GRND )
 COVAL(SMOKE ,S5C , GRND , GRND )
 COVAL(SMOKE ,S6C , GRND , GRND )
 COVAL(SMOKE ,S7C , GRND , GRND )
 COVAL(SMOKE ,S8C , GRND , GRND )
 COVAL(SMOKE ,S9C , GRND , GRND )
 COVAL(SMOKE ,S10C, GRND , GRND )
 COVAL(SMOKE ,S11C, GRND , GRND )
 COVAL(SMOKE ,S12C, GRND , GRND )
 COVAL(SMOKE ,S13C, GRND , GRND )
 COVAL(SMOKE ,S14C, GRND , GRND )
 COVAL(SMOKE ,S15C, GRND , GRND )
 COVAL(SMOKE ,S16C, GRND , GRND )
 COVAL(SMOKE ,S17C, GRND , GRND )
 COVAL(SMOKE ,S18C, GRND , GRND )
 COVAL(SMOKE ,S19C, GRND , GRND )
 COVAL(SMOKE ,S20C, GRND , GRND )
 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(T1C)=20
 LITER(T2C)=20 ;LITER(T3C)=20
 LITER(T4C)=20 ;LITER(T5C)=20
 LITER(T6C)=20 ;LITER(T7C)=20
 LITER(T8C)=20 ;LITER(T9C)=20
 LITER(T10C)=20 ;LITER(T11C)=20
 LITER(T12C)=20 ;LITER(T13C)=20
 LITER(T14C)=20 ;LITER(T15C)=20
 LITER(T16C)=20 ;LITER(T17C)=20
 LITER(T18C)=20 ;LITER(T19C)=20
 LITER(T20C)=20 ;LITER(S1C)=20
 LITER(S2C)=20 ;LITER(S3C)=20
 LITER(S4C)=20 ;LITER(S5C)=20
 LITER(S6C)=20 ;LITER(S7C)=20
 LITER(S8C)=20 ;LITER(S9C)=20
 LITER(S10C)=20 ;LITER(S11C)=20
 LITER(S12C)=20 ;LITER(S13C)=20
 LITER(S14C)=20 ;LITER(S15C)=20
 LITER(S16C)=20 ;LITER(S17C)=20
 LITER(S18C)=20 ;LITER(S19C)=20
 LITER(S20C)=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(T1C)=1.0E-03
 ENDIT(T2C)=1.0E-03 ;ENDIT(T3C)=1.0E-03
 ENDIT(T4C)=1.0E-03 ;ENDIT(T5C)=1.0E-03
 ENDIT(T6C)=1.0E-03 ;ENDIT(T7C)=1.0E-03
 ENDIT(T8C)=1.0E-03 ;ENDIT(T9C)=1.0E-03
 ENDIT(T10C)=1.0E-03 ;ENDIT(T11C)=1.0E-03
 ENDIT(T12C)=1.0E-03 ;ENDIT(T13C)=1.0E-03
 ENDIT(T14C)=1.0E-03 ;ENDIT(T15C)=1.0E-03
 ENDIT(T16C)=1.0E-03 ;ENDIT(T17C)=1.0E-03
 ENDIT(T18C)=1.0E-03 ;ENDIT(T19C)=1.0E-03
 ENDIT(T20C)=1.0E-03 ;ENDIT(S1C)=1.0E-03
 ENDIT(S2C)=1.0E-03 ;ENDIT(S3C)=1.0E-03
 ENDIT(S4C)=1.0E-03 ;ENDIT(S5C)=1.0E-03
 ENDIT(S6C)=1.0E-03 ;ENDIT(S7C)=1.0E-03
 ENDIT(S8C)=1.0E-03 ;ENDIT(S9C)=1.0E-03
 ENDIT(S10C)=1.0E-03 ;ENDIT(S11C)=1.0E-03
 ENDIT(S12C)=1.0E-03 ;ENDIT(S13C)=1.0E-03
 ENDIT(S14C)=1.0E-03 ;ENDIT(S15C)=1.0E-03
 ENDIT(S16C)=1.0E-03 ;ENDIT(S17C)=1.0E-03
 ENDIT(S18C)=1.0E-03 ;ENDIT(S19C)=1.0E-03
 ENDIT(S20C)=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(T1C,FALSDT,1.0E+09)
 RELAX(T2C,FALSDT,1.0E+09)
 RELAX(T3C,FALSDT,1.0E+09)
 RELAX(T4C,FALSDT,1.0E+09)
 RELAX(T5C,FALSDT,1.0E+09)
 RELAX(T6C,FALSDT,1.0E+09)
 RELAX(T7C,FALSDT,1.0E+09)
 RELAX(T8C,FALSDT,1.0E+09)
 RELAX(T9C,FALSDT,1.0E+09)
 RELAX(T10C,FALSDT,1.0E+09)
 RELAX(T11C,FALSDT,1.0E+09)
 RELAX(T12C,FALSDT,1.0E+09)
 RELAX(T13C,FALSDT,1.0E+09)
 RELAX(T14C,FALSDT,1.0E+09)
 RELAX(T15C,FALSDT,1.0E+09)
 RELAX(T16C,FALSDT,1.0E+09)
 RELAX(T17C,FALSDT,1.0E+09)
 RELAX(T18C,FALSDT,1.0E+09)
 RELAX(T19C,FALSDT,1.0E+09)
 RELAX(T20C,FALSDT,1.0E+09)
 RELAX(S1C,FALSDT,1.0E+09)
 RELAX(S2C,FALSDT,1.0E+09)
 RELAX(S3C,FALSDT,1.0E+09)
 RELAX(S4C,FALSDT,1.0E+09)
 RELAX(S5C,FALSDT,1.0E+09)
 RELAX(S6C,FALSDT,1.0E+09)
 RELAX(S7C,FALSDT,1.0E+09)
 RELAX(S8C,FALSDT,1.0E+09)
 RELAX(S9C,FALSDT,1.0E+09)
 RELAX(S10C,FALSDT,1.0E+09)
 RELAX(S11C,FALSDT,1.0E+09)
 RELAX(S12C,FALSDT,1.0E+09)
 RELAX(S13C,FALSDT,1.0E+09)
 RELAX(S14C,FALSDT,1.0E+09)
 RELAX(S15C,FALSDT,1.0E+09)
 RELAX(S16C,FALSDT,1.0E+09)
 RELAX(S17C,FALSDT,1.0E+09)
 RELAX(S18C,FALSDT,1.0E+09)
 RELAX(S19C,FALSDT,1.0E+09)
 RELAX(S20C,FALSDT,1.0E+09)
 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(T1C)=1.0E+10 ;VARMIN(T1C)=-1.0E+10
 VARMAX(T2C)=1.0E+10 ;VARMIN(T2C)=-1.0E+10
 VARMAX(T3C)=1.0E+10 ;VARMIN(T3C)=-1.0E+10
 VARMAX(T4C)=1.0E+10 ;VARMIN(T4C)=-1.0E+10
 VARMAX(T5C)=1.0E+10 ;VARMIN(T5C)=-1.0E+10
 VARMAX(T6C)=1.0E+10 ;VARMIN(T6C)=-1.0E+10
 VARMAX(T7C)=1.0E+10 ;VARMIN(T7C)=-1.0E+10
 VARMAX(T8C)=1.0E+10 ;VARMIN(T8C)=-1.0E+10
 VARMAX(T9C)=1.0E+10 ;VARMIN(T9C)=-1.0E+10
 VARMAX(T10C)=1.0E+10 ;VARMIN(T10C)=-1.0E+10
 VARMAX(T11C)=1.0E+10 ;VARMIN(T11C)=-1.0E+10
 VARMAX(T12C)=1.0E+10 ;VARMIN(T12C)=-1.0E+10
 VARMAX(T13C)=1.0E+10 ;VARMIN(T13C)=-1.0E+10
 VARMAX(T14C)=1.0E+10 ;VARMIN(T14C)=-1.0E+10
 VARMAX(T15C)=1.0E+10 ;VARMIN(T15C)=-1.0E+10
 VARMAX(T16C)=1.0E+10 ;VARMIN(T16C)=-1.0E+10
 VARMAX(T17C)=1.0E+10 ;VARMIN(T17C)=-1.0E+10
 VARMAX(T18C)=1.0E+10 ;VARMIN(T18C)=-1.0E+10
 VARMAX(T19C)=1.0E+10 ;VARMIN(T19C)=-1.0E+10
 VARMAX(T20C)=1.0E+10 ;VARMIN(T20C)=-1.0E+10
 VARMAX(S1C)=1.0E+10 ;VARMIN(S1C)=-1.0E+10
 VARMAX(S2C)=1.0E+10 ;VARMIN(S2C)=-1.0E+10
 VARMAX(S3C)=1.0E+10 ;VARMIN(S3C)=-1.0E+10
 VARMAX(S4C)=1.0E+10 ;VARMIN(S4C)=-1.0E+10
 VARMAX(S5C)=1.0E+10 ;VARMIN(S5C)=-1.0E+10
 VARMAX(S6C)=1.0E+10 ;VARMIN(S6C)=-1.0E+10
 VARMAX(S7C)=1.0E+10 ;VARMIN(S7C)=-1.0E+10
 VARMAX(S8C)=1.0E+10 ;VARMIN(S8C)=-1.0E+10
 VARMAX(S9C)=1.0E+10 ;VARMIN(S9C)=-1.0E+10
 VARMAX(S10C)=1.0E+10 ;VARMIN(S10C)=-1.0E+10
 VARMAX(S11C)=1.0E+10 ;VARMIN(S11C)=-1.0E+10
 VARMAX(S12C)=1.0E+10 ;VARMIN(S12C)=-1.0E+10
 VARMAX(S13C)=1.0E+10 ;VARMIN(S13C)=-1.0E+10
 VARMAX(S14C)=1.0E+10 ;VARMIN(S14C)=-1.0E+10
 VARMAX(S15C)=1.0E+10 ;VARMIN(S15C)=-1.0E+10
 VARMAX(S16C)=1.0E+10 ;VARMIN(S16C)=-1.0E+10
 VARMAX(S17C)=1.0E+10 ;VARMIN(S17C)=-1.0E+10
 VARMAX(S18C)=1.0E+10 ;VARMIN(S18C)=-1.0E+10
 VARMAX(S19C)=1.0E+10 ;VARMIN(S19C)=-1.0E+10
 VARMAX(S20C)=1.0E+10 ;VARMIN(S20C)=-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.)
 SPEDAT(SET,MFM,RADFACT,R,0.5)
 ************************************************************
  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(T1C)=695.700012
 EX(T2C)=2280. ;EX(T3C)=2258.
 EX(T4C)=2166. ;EX(T5C)=2074.
 EX(T6C)=1982. ;EX(T7C)=1890.
 EX(T8C)=1798. ;EX(T9C)=1706.
 EX(T10C)=1614. ;EX(T11C)=1523.
 EX(T12C)=1431. ;EX(T13C)=1339.
 EX(T14C)=1247. ;EX(T15C)=1155.
 EX(T16C)=1063. ;EX(T17C)=971.299988
 EX(T18C)=879.5 ;EX(T19C)=787.599976
 EX(T20C)=695.700012 ;EX(S1C)=0.
 EX(S2C)=1.805E-05 ;EX(S3C)=1.53E-04
 EX(S4C)=2.297E-04 ;EX(S5C)=2.581E-04
 EX(S6C)=2.542E-04 ;EX(S7C)=2.309E-04
 EX(S8C)=1.978E-04 ;EX(S9C)=1.616E-04
 EX(S10C)=1.267E-04 ;EX(S11C)=9.567E-05
 EX(S12C)=6.961E-05 ;EX(S13C)=4.881E-05
 EX(S14C)=3.297E-05 ;EX(S15C)=2.142E-05
 EX(S16C)=1.337E-05 ;EX(S17C)=8.002E-06
 EX(S18C)=4.586E-06 ;EX(S19C)=2.451E-06
 EX(S20C)=8.065E-07 ;EX(PRDM)=0.3831
 EX(OXDM)=0. ;EX(FRFM)=0.1397
 EX(TMPM)=1314. ;EX(SMOS)=1.696E-04
 EX(SMOM)=7.553E-05 ;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(T1C,Y,N,Y,Y,Y,Y)
 OUTPUT(T2C,Y,N,Y,Y,Y,Y)
 OUTPUT(T3C,Y,N,Y,Y,Y,Y)
 OUTPUT(T4C,Y,N,Y,Y,Y,Y)
 OUTPUT(T5C,Y,N,Y,Y,Y,Y)
 OUTPUT(T6C,Y,N,Y,Y,Y,Y)
 OUTPUT(T7C,Y,N,Y,Y,Y,Y)
 OUTPUT(T8C,Y,N,Y,Y,Y,Y)
 OUTPUT(T9C,Y,N,Y,Y,Y,Y)
 OUTPUT(T10C,Y,N,Y,Y,Y,Y)
 OUTPUT(T11C,Y,N,Y,Y,Y,Y)
 OUTPUT(T12C,Y,N,Y,Y,Y,Y)
 OUTPUT(T13C,Y,N,Y,Y,Y,Y)
 OUTPUT(T14C,Y,N,Y,Y,Y,Y)
 OUTPUT(T15C,Y,N,Y,Y,Y,Y)
 OUTPUT(T16C,Y,N,Y,Y,Y,Y)
 OUTPUT(T17C,Y,N,Y,Y,Y,Y)
 OUTPUT(T18C,Y,N,Y,Y,Y,Y)
 OUTPUT(T19C,Y,N,Y,Y,Y,Y)
 OUTPUT(T20C,Y,N,Y,Y,Y,Y)
 OUTPUT(S1C,Y,N,Y,Y,Y,Y)
 OUTPUT(S2C,Y,N,Y,Y,Y,Y)
 OUTPUT(S3C,Y,N,Y,Y,Y,Y)
 OUTPUT(S4C,Y,N,Y,Y,Y,Y)
 OUTPUT(S5C,Y,N,Y,Y,Y,Y)
 OUTPUT(S6C,Y,N,Y,Y,Y,Y)
 OUTPUT(S7C,Y,N,Y,Y,Y,Y)
 OUTPUT(S8C,Y,N,Y,Y,Y,Y)
 OUTPUT(S9C,Y,N,Y,Y,Y,Y)
 OUTPUT(S10C,Y,N,Y,Y,Y,Y)
 OUTPUT(S11C,Y,N,Y,Y,Y,Y)
 OUTPUT(S12C,Y,N,Y,Y,Y,Y)
 OUTPUT(S13C,Y,N,Y,Y,Y,Y)
 OUTPUT(S14C,Y,N,Y,Y,Y,Y)
 OUTPUT(S15C,Y,N,Y,Y,Y,Y)
 OUTPUT(S16C,Y,N,Y,Y,Y,Y)
 OUTPUT(S17C,Y,N,Y,Y,Y,Y)
 OUTPUT(S18C,Y,N,Y,Y,Y,Y)
 OUTPUT(S19C,Y,N,Y,Y,Y,Y)
 OUTPUT(S20C,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. )
 
 PATCH(TPROF1 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(TPROF1 ,T1C ,0. ,0. )
 
 PATCH(TPROF2 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(TPROF2 ,T2C ,0. ,0. )
 
 PATCH(TPROF3 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(TPROF3 ,T3C ,0. ,0. )
 
 PATCH(TPROF4 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(TPROF4 ,T4C ,0. ,0. )
 
 PATCH(TPROF5 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(TPROF5 ,T5C ,0. ,0. )
 
 PATCH(TPROF6 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(TPROF6 ,T6C ,0. ,0. )
 
 PATCH(TPROF7 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(TPROF7 ,T7C ,0. ,0. )
 
 PATCH(TPROF8 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(TPROF8 ,T8C ,0. ,0. )
 
 PATCH(TPROF9 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(TPROF9 ,T9C ,0. ,0. )
 
 PATCH(TPROF10 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(TPROF10 ,T10C,0. ,0. )
 
 PATCH(TPROF11 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(TPROF11 ,T11C,0. ,0. )
 
 PATCH(TPROF12 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(TPROF12 ,T12C,0. ,0. )
 
 PATCH(TPROF13 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(TPROF13 ,T13C,0. ,0. )
 
 PATCH(TPROF14 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(TPROF14 ,T14C,0. ,0. )
 
 PATCH(TPROF15 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(TPROF15 ,T15C,0. ,0. )
 
 PATCH(TPROF16 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(TPROF16 ,T16C,0. ,0. )
 
 PATCH(TPROF17 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(TPROF17 ,T17C,0. ,0. )
 
 PATCH(TPROF18 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(TPROF18 ,T18C,0. ,0. )
 
 PATCH(TPROF19 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(TPROF19 ,T19C,0. ,0. )
 
 PATCH(TPROF20 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(TPROF20 ,T20C,0. ,0. )
 
 PATCH(SPROF1 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(SPROF1 ,S1C ,0. ,0. )
 
 PATCH(SPROF2 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(SPROF2 ,S2C ,0. ,0. )
 
 PATCH(SPROF3 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(SPROF3 ,S3C ,0. ,0. )
 
 PATCH(SPROF4 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(SPROF4 ,S4C ,0. ,0. )
 
 PATCH(SPROF5 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(SPROF5 ,S5C ,0. ,0. )
 
 PATCH(SPROF6 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(SPROF6 ,S6C ,0. ,0. )
 
 PATCH(SPROF7 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(SPROF7 ,S7C ,0. ,0. )
 
 PATCH(SPROF8 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(SPROF8 ,S8C ,0. ,0. )
 
 PATCH(SPROF9 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(SPROF9 ,S9C ,0. ,0. )
 
 PATCH(SPROF10 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(SPROF10 ,S10C,0. ,0. )
 
 PATCH(SPROF11 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(SPROF11 ,S11C,0. ,0. )
 
 PATCH(SPROF12 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(SPROF12 ,S12C,0. ,0. )
 
 PATCH(SPROF13 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(SPROF13 ,S13C,0. ,0. )
 
 PATCH(SPROF14 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(SPROF14 ,S14C,0. ,0. )
 
 PATCH(SPROF15 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(SPROF15 ,S15C,0. ,0. )
 
 PATCH(SPROF16 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(SPROF16 ,S16C,0. ,0. )
 
 PATCH(SPROF17 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(SPROF17 ,S17C,0. ,0. )
 
 PATCH(SPROF18 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(SPROF18 ,S18C,0. ,0. )
 
 PATCH(SPROF19 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(SPROF19 ,S19C,0. ,0. )
 
 PATCH(SPROF20 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1)
 PLOT(SPROF20 ,S20C,0. ,0. )
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP