TALK=T;RUN(1,1)
  PHOTON USE
  p;;;;

  view x
  msg compare
  msg fluid-average values: avef, tmpm, frfm, prdm, oxdm, smom
  msg with corresponding single-fluid EBU values
  msg                    mixf, tmp1, fuel, prod, oxid, smos
  pause
  msg x-plane mixture fraction avef according to MFM
  con avef x 1 fi;0.001;pause;con off;red
  msg x-plane mixture fraction mixf according to EBU
  msg This is, as it should be, similar to avef
  con mixf x 1 fi;0.001;pause;con off;red

  msg x-plane                  tmpm according to MFM
  con tmpm x 1 fi;0.001;pause;con off;red
  msg x-plane                  tmp1 according to EBU
  con tmp1 x 1 fi;0.001;pause;con off;red

  msg x-plane                  frfm according to MFM
  con frfm x 1 fi;0.001;pause;con off;red
  msg x-plane                  fuel according to EBU
  con fuel x 1 fi;0.001;pause;con off;red

  msg x-plane                  prdm according to MFM
  con prdm x 1 fi;0.001;pause;con off;red
  msg x-plane                  prod according to EBU
  con prod x 1 fi;0.001;pause;con off;red

  msg x-plane                  oxdm according to MFM
  con oxdm x 1 fi;0.001;pause;con off;red
  msg x-plane                  oxid according to EBU
  con oxidd x 1 fi;0.001;pause;con off;red

  msg x-plane smoke concentration according to eddy-break-up model
  con smos x 1 fi;0.001;pause;con off;red
  msg x-plane smoke concentration according to multi-fluid model
  con smom x 1 fi;0.001;pause;con off;red

  msg x-plane MFM fluid-1 concentration
  con f1 x 1 fi;0.001;pause;con off;red
  msg x-plane MFM fluid-11 concentration
  con f11 x 1 fi;0.001;pause;con off;red
  msg x-plane MFM fluid-6 concentration
  con f6 x 1 fi;0.001
  ENDUSE
  PHOTON USE
  p;;;;

  view z
  msg exit smoke concentration according to single-fluid model
  con smos z m fi;0.001
  pause;con off;red
  msg exit smoke concentration according to multi-fluid model
  con smom z m fi;0.001;pause;con off;red;view x

  msg x-plane smoke concentration according to single-fluid model
  con smos x 1 fi;0.001;pause;con off;red
  msg x-plane smoke concentration according to multi-fluid model
  con smom x 1 fi;0.001;pause;con off;red
  msg x-plane mixture fraction avef according to multi-fluid model
  con avef x 1 fi;0.001;pause;con off;red
  msg x-plane mixture fraction mixf according to single-fluid model
  con mixf x 1 fi;0.001;pause;con off;red
  msg x-plane fluid-1 concentration according to multi-fluid model
  con f1 x 1 fi;0.001;pause;con off;red
  msg x-plane fluid-11 concentration according to multi-fluid model
  con f11 x 1 fi;0.001;pause;con off;red
  msg x-plane fluid-6 concentration according to multi-fluid model
  con f6 x 1 fi;0.001
  ENDUSE
  DISPLAY
  This simulation shows how the multi-fluid model of turbulence
  predicts a smoke-production pattern in a gas-turbine combustor
  which is different from, and probably more realistic than, that
  predicted by a conventional single-fluid model.

  It is explained and discussed extensively in a report which is
  published on CHAM's website www.cham.co.uk (Click on short-cuts,
  then MFM then 1998b).

  The geometry is that of library case 492, which is loaded first.
  Note the $ in #$492; it ensures that the above PHOTON USE
  commands remain in the Q1 file, and therefore take priority over
  those which appear in case 493.

  The injected fuel-air composition, which is stoichiometric in
  case 492, is made twice stoichiometric here.

  To reduce computation time, it is best to answer N to the question
  about radiation.
  ENDDIS

  PHOTON USE
  ext;;;;

  gr ou z 1; use patgeo
  msg boundary condition patches. Press RETURN
  pause
  msg temperature contours. Press RETURN
  con tmp1 y m fi;0.001;con tmp1 x m fi;0.001;con tmp1 z m fi;0.001
  gr ou y m;gr ou z m
  msg press RETURN for view x
  pause;  con off;  view x;  con tmp1 x m fi;0.001
  con tmp1 y m fi;0.001;  gr ou y m;  gr ou x m;  gr ou z 1
  msg This is view x. Press RETURN for view y
  pause;  con off;  view y;  con tmp1 y m fi;0.001
  msg This is view y. Press RETURN for view z
  pause;  con off; view z;  con tmp1 z m fi;0.001;  gr ou z m
  msg This is view z. Press RETURN for view x and velocity vectors
  pause;  con off;  view x
  msg velocity vectors. Type menu for menu and further possibilities
  msg Press e to end. Otherwise enter photon-readable commands
  vec x m sh;  gr ou y m;  gr ou x m;  gr ou z 1
  ENDUSE
  DISPLAY
  READQ1_BEGIN
  The following text is provided as an example of what, by use of
  the readq1 command, can be transmitted to EARTH for writing near
  the top of the RESULT file. Such text should leave columns 1 and 2
  blank; and it should not extend beyond the 68th column.

    ****************************************************************
    * This library case dates from the earliest days of PHOENICS,  *
    * when Professor WU Chung-Hua ('turbomachinery Wu', who had    *
    * returned to China from the USA) visited CHAM in 1982.        *
    *                                                              *
    * The configuration of the combustion chamber, and its being   *
    * supplied with premixed fuel vapour and air, was proposed by  *
    * Professor Wu's accompanying assistant.                       *
    ****************************************************************
  READQ1_END
  The shape of the combustion chamber is as shown.
  Pre-mixed fuel vapour and air enter near the axis on left.
  Secondary and dilution air enter through holes in outer wall.
  A 36-degree sector is simulated.

                 secondary ox. inlets   dilution inlet
           ____________ 1,2 _______________     ________
  blocked /
  region /
        /
       /                                             outlet
       |
       |_____   _
  fuel-ox. inlet|            Symmetry axis
       -- - -- -| -- - -- -  -- - -- - -- - -- - -- - -- -

  The flow is turbulent; the Simple Chemical Reaction Scheme is
  used; and the reaction-rate is physically controlled by means of
  the Eddy-Breakup Model.
  ENDDIS
  DISPLAY
  Notes:

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

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

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

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

  (4) the values are:
      mfm  = $L004
      mfm1 = $L003
      mfm2 = $L002
      mfm3 = $L001
  ENDDIS
  display
  This 2D version of the 3D combustor facilitates exploration of the
  influences of the model variables, and of the boundary conditions.

  Some of these are explored in the following cases, 401, etc.
  enddis
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(2D Combustor With Smoke Generation      )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 400
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = F
 NX = 1
 XULAST =0.628319
   Method of pairs used for grid setting.
 XFRAC(1)=-6. ;XFRAC(2)=0.166666
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 10
 YVLAST =1.
   Method of pairs used for grid setting.
 YFRAC(1)=-5. ;YFRAC(2)=6.5E-03
 YFRAC(3)=1. ;YFRAC(4)=5.0E-03
 YFRAC(5)=1. ;YFRAC(6)=8.0E-03
 YFRAC(7)=3. ;YFRAC(8)=6.5E-03
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 13
 ZWLAST =1.
   Method of pairs used for grid setting.
 ZFRAC(1)=-2. ;ZFRAC(3)=1.
 ZFRAC(5)=1. ;ZFRAC(7)=1.
 ZFRAC(9)=1. ;ZFRAC(11)=7.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(1)=P1 ;NAME(5)=V1
 NAME(7)=W1 ;NAME(12)=KE
 NAME(13)=EP ;NAME(14)=H1
 NAME(115)=PRPS ;NAME(116)=PRDM
 NAME(117)=OXDM ;NAME(118)=FRFM
 NAME(119)=TMPM ;NAME(120)=SMOS
 NAME(121)=SMOM ;NAME(122)=F1
 NAME(123)=F2 ;NAME(124)=F3
 NAME(125)=F4 ;NAME(126)=F5
 NAME(127)=F6 ;NAME(128)=F7
 NAME(129)=F8 ;NAME(130)=F9
 NAME(131)=F10 ;NAME(132)=F11
 NAME(133)=MNSQ ;NAME(134)=AVEF
 NAME(135)=MIXL ;NAME(136)=RATE
 NAME(137)=DEGF ;NAME(138)=RRAT
 NAME(139)=EPKE ;NAME(140)=HPOR
 NAME(141)=NPOR ;NAME(142)=EPOR
 NAME(143)=VPOR ;NAME(144)=ENUT
 NAME(145)=RHO1 ;NAME(146)=TMP1
 NAME(147)=PROD ;NAME(148)=OXID
 NAME(149)=FUEL ;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,Y,N,N,N)
 SOLUTN(V1,Y,Y,N,N,N,Y)
 SOLUTN(W1,Y,Y,N,N,N,Y)
 SOLUTN(KE,Y,Y,N,N,N,N)
 SOLUTN(EP,Y,Y,N,N,N,N)
 SOLUTN(H1,Y,Y,N,N,N,Y)
 SOLUTN(PRPS,Y,N,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(MNSQ,Y,N,N,N,N,Y)
 SOLUTN(AVEF,Y,N,N,N,N,Y)
 SOLUTN(MIXL,Y,N,N,N,N,N)
 SOLUTN(RATE,Y,N,N,N,N,Y)
 SOLUTN(DEGF,Y,N,N,N,N,N)
 SOLUTN(RRAT,Y,N,N,N,N,N)
 SOLUTN(EPKE,Y,N,N,N,N,N)
 SOLUTN(HPOR,Y,N,N,N,N,N)
 SOLUTN(NPOR,Y,N,N,N,N,N)
 SOLUTN(EPOR,Y,N,N,N,N,N)
 SOLUTN(VPOR,Y,N,N,N,N,N)
 SOLUTN(ENUT,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(MIXF,Y,Y,N,N,N,Y)
 DEN1 = 145
 VIST = 144
 EPOR = 142 ;HPOR = 140 ;NPOR = 141 ;VPOR = 143
 TEMP1 = 146
 PRPS = 115
 ************************************************************
  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(V1,Y,Y,Y,Y,Y,Y)
 TERMS(W1,Y,Y,Y,Y,Y,Y)
 TERMS(KE,N,Y,Y,Y,Y,N)
 TERMS(EP,N,Y,Y,Y,Y,N)
 TERMS(H1,N,Y,Y,N,Y,N)
 TERMS(SMOS,N,Y,Y,Y,Y,Y)
 TERMS(SMOM,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(MIXF,N,Y,Y,Y,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWRH1 = T
 NEWENT = 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 = GRND4
 TSURR =0. ;TEMP0 =0.
 PRESS0 =8.0E+05
 DVO1DT =0. ;DRH1DP =0.
 RHO1A =16. ;RHO1B =29.
 RHO1C =28.
 TMP1A =0. ;TMP1B =0.
 TMP1C =0.
 TMP2A =0.054825 ;TMP2B =4.9E+07
 TMP2C =0.054825
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 EL1A =0. ;EL1B =0.
 EL1C =0.
 ENUL =9.983361E-06 ;ENUT = GRND3
 ENUTA =0. ;ENUTB =0.
 ENUTC =0.
 IENUTA = 0
 PRNDTL(V1)=1. ;PRNDTL(W1)=1.
 PRNDTL(KE)=1. ;PRNDTL(EP)=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(MIXF)=1.
 PRT(V1)=1. ;PRT(W1)=1.
 PRT(KE)=1. ;PRT(EP)=1.314
 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(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(V1)=1.0E-10
 FIINIT(W1)=10. ;FIINIT(KE)=11.25
 FIINIT(EP)=353.255707 ;FIINIT(H1)=1.1595E+06
 FIINIT(PRPS)=-1. ;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.090909
 FIINIT(F2)=0.090909 ;FIINIT(F3)=0.090909
 FIINIT(F4)=0.090909 ;FIINIT(F5)=0.090909
 FIINIT(F6)=0.090909 ;FIINIT(F7)=0.090909
 FIINIT(F8)=0.090909 ;FIINIT(F9)=0.090909
 FIINIT(F10)=0.090909 ;FIINIT(F11)=0.090909
 FIINIT(MNSQ)=1.0E-10 ;FIINIT(AVEF)=1.0E-10
 FIINIT(MIXL)=1.0E-10 ;FIINIT(RATE)=31.400507
 FIINIT(DEGF)=1.0E-10 ;FIINIT(RRAT)=1.0E-10
 FIINIT(EPKE)=1.0E-10 ;FIINIT(HPOR)=1.
 FIINIT(NPOR)=1. ;FIINIT(EPOR)=1.
 FIINIT(VPOR)=1. ;FIINIT(ENUT)=1.0E-10
 FIINIT(RHO1)=3.606 ;FIINIT(TMP1)=773.
 FIINIT(PROD)=1.0E-10 ;FIINIT(OXID)=1.0E-10
 FIINIT(FUEL)=0.109649 ;FIINIT(MIXF)=0.109649
 
 PATCH(CMP1 ,INIVAL, 1, 6, 1, 2, 1, 3, 1, 1)
 INIT(CMP1 ,PRPS,0. ,199. )
 
 PATCH(CMP2 ,INIVAL, 1, 1, 7, 7, 1, 1, 1, 1)
 
 PATCH(CMP3 ,INIVAL, 1, 6, 7, 7, 1, 1, 1, 1)
 INIT(CMP3 ,VPOR,0. ,0.709 )
 
 PATCH(CMP4 ,INIVAL, 1, 6, 7, 7, 1, 1, 1, 1)
 INIT(CMP4 ,NPOR,0. ,0.418 )
 
 PATCH(CMP5 ,INIVAL, 1, 6, 7, 7, 1, 1, 1, 1)
 INIT(CMP5 ,HPOR,0. ,1. )
 
 PATCH(CMP6 ,INIVAL, 1, 1, 8, 8, 1, 1, 1, 1)
 
 PATCH(CMP7 ,INIVAL, 1, 6, 8, 8, 1, 1, 1, 1)
 INIT(CMP7 ,VPOR,0. ,0.185 )
 
 PATCH(CMP8 ,INIVAL, 1, 6, 8, 8, 1, 1, 1, 1)
 INIT(CMP8 ,V1 ,0. ,0. )
 INIT(CMP8 ,NPOR,0. ,0. )
 
 PATCH(CMP9 ,INIVAL, 1, 6, 8, 8, 1, 1, 1, 1)
 INIT(CMP9 ,HPOR,0. ,0.877 )
 
 PATCH(CMP10 ,INIVAL, 1, 1, 9, 10, 1, 1, 1, 1)
 
 PATCH(CMP11 ,INIVAL, 1, 6, 9, 10, 1, 1, 1, 1)
 INIT(CMP11 ,PRPS,0. ,199. )
 
 PATCH(CMP12 ,INIVAL, 1, 6, 9, 9, 1, 1, 1, 1)
 INIT(CMP12 ,V1 ,0. ,0. )
 INIT(CMP12 ,NPOR,0. ,0. )
 
 PATCH(CMP13 ,INIVAL, 1, 6, 9, 10, 1, 1, 1, 1)
 INIT(CMP13 ,W1 ,0. ,0. )
 INIT(CMP13 ,HPOR,0. ,0. )
 
 PATCH(CMP14 ,INIVAL, 1, 1, 8, 8, 2, 2, 1, 1)
 
 PATCH(CMP15 ,INIVAL, 1, 6, 8, 8, 2, 2, 1, 1)
 INIT(CMP15 ,VPOR,0. ,0.997 )
 
 PATCH(CMP16 ,INIVAL, 1, 6, 8, 8, 2, 2, 1, 1)
 INIT(CMP16 ,NPOR,0. ,0.946 )
 
 PATCH(CMP17 ,INIVAL, 1, 6, 8, 8, 2, 2, 1, 1)
 INIT(CMP17 ,HPOR,0. ,1. )
 
 PATCH(CMP18 ,INIVAL, 1, 1, 9, 9, 2, 2, 1, 1)
 
 PATCH(CMP19 ,INIVAL, 1, 6, 9, 9, 2, 2, 1, 1)
 INIT(CMP19 ,VPOR,0. ,0.706 )
 
 PATCH(CMP20 ,INIVAL, 1, 6, 9, 9, 2, 2, 1, 1)
 INIT(CMP20 ,NPOR,0. ,0.473 )
 
 PATCH(CMP21 ,INIVAL, 1, 6, 9, 9, 2, 2, 1, 1)
 INIT(CMP21 ,HPOR,0. ,1. )
 
 PATCH(CMP22 ,INIVAL, 1, 1, 10, 10, 2, 2, 1, 1)
 
 PATCH(CMP23 ,INIVAL, 1, 6, 10, 10, 2, 2, 1, 1)
 INIT(CMP23 ,VPOR,0. ,0.236 )
 
 PATCH(CMP24 ,INIVAL, 1, 6, 10, 10, 2, 2, 1, 1)
 
 PATCH(CMP25 ,INIVAL, 1, 6, 10, 10, 2, 2, 1, 1)
 INIT(CMP25 ,HPOR,0. ,1. )
 
 PATCH(WALL4 ,INIVAL, 1, 6, 2, 2, 1, 3, 1, 1)
 INIT(WALL4 ,V1 ,0. ,0. )
 INIT(WALL4 ,NPOR,0. ,0. )
 
 PATCH(WALL6 ,INIVAL, 1, 6, 1, 2, 3, 3, 1, 1)
 INIT(WALL6 ,W1 ,0. ,0. )
 INIT(WALL6 ,HPOR,0. ,0. )
 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(KESOURCE,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1)
 COVAL(KESOURCE,KE , GRND4 , GRND4 )
 COVAL(KESOURCE,EP , GRND4 , GRND4 )
 
 PATCH(FOIN ,SOUTH , 1, 6, 3, 3, 3, 3, 1, 1)
 COVAL(FOIN ,P1 , FIXFLU ,180.300003 )
 COVAL(FOIN ,V1 ,0. ,50. )
 COVAL(FOIN ,W1 ,0. ,0. )
 COVAL(FOIN ,KE ,0. ,0.0625 )
 COVAL(FOIN ,EP ,0. ,3.672038 )
 COVAL(FOIN ,H1 ,0. ,4.001728E+06 )
 COVAL(FOIN ,F11 ,0. ,1. )
 COVAL(FOIN ,FUEL,0. ,0.058005 )
 COVAL(FOIN ,MIXF,0. ,0.109649 )
 
 PATCH(SOIN1 ,CELL , 1, 1, 10, 10, 7, 7, 1, 1)
 COVAL(SOIN1 ,P1 , FIXFLU ,8.157E-03 )
 COVAL(SOIN1 ,V1 ,0. ,-40. )
 COVAL(SOIN1 ,W1 ,0. ,0. )
 COVAL(SOIN1 ,KE ,0. ,0.123 )
 COVAL(SOIN1 ,EP ,0. ,8.429416 )
 COVAL(SOIN1 ,H1 ,0. ,1.1595E+06 )
 COVAL(SOIN1 ,F1 ,0. ,1. )
 COVAL(SOIN1 ,FUEL,0. ,0. )
 COVAL(SOIN1 ,MIXF,0. ,0. )
 
 PATCH(SOIN2 ,CELL , 4, 4, 10, 10, 7, 7, 1, 1)
 COVAL(SOIN2 ,P1 , FIXFLU ,8.157E-03 )
 COVAL(SOIN2 ,V1 ,0. ,-40. )
 COVAL(SOIN2 ,W1 ,0. ,0. )
 COVAL(SOIN2 ,KE ,0. ,0.123 )
 COVAL(SOIN2 ,EP ,0. ,8.429416 )
 COVAL(SOIN2 ,H1 ,0. ,1.1595E+06 )
 COVAL(SOIN2 ,F1 ,0. ,1. )
 COVAL(SOIN2 ,FUEL,0. ,0. )
 COVAL(SOIN2 ,MIXF,0. ,0. )
 
 PATCH(DILUIN ,CELL , 1, 1, 10, 10, 10, 10, 1, 1)
 COVAL(DILUIN ,P1 , FIXFLU ,0.0145 )
 COVAL(DILUIN ,V1 ,0. ,-40. )
 COVAL(DILUIN ,W1 ,0. ,0. )
 COVAL(DILUIN ,KE ,0. ,0.123 )
 COVAL(DILUIN ,EP ,0. ,6.299674 )
 COVAL(DILUIN ,H1 ,0. ,1.1595E+06 )
 COVAL(DILUIN ,F1 ,0. ,1. )
 COVAL(DILUIN ,FUEL,0. ,0. )
 COVAL(DILUIN ,MIXF,0. ,0. )
 
 PATCH(OUTLET ,HIGH , 1, 6, 1, 10, 13, 13, 1, 1)
 COVAL(OUTLET ,P1 ,1000. ,0. )
 COVAL(OUTLET ,V1 ,0. ,0. )
 COVAL(OUTLET ,W1 ,0. ,0. )
 COVAL(OUTLET ,KE ,0. , SAME )
 COVAL(OUTLET ,EP ,0. , SAME )
 COVAL(OUTLET ,H1 ,0. , SAME )
 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 ,FUEL,0. , SAME )
 COVAL(OUTLET ,MIXF,0. , SAME )
 
 PATCH(WALL1 ,NWALL , 1, 6, 10, 10, 3, 6, 1, 1)
 COVAL(WALL1 ,W1 , GRND2 ,0. )
 COVAL(WALL1 ,KE , GRND2 , GRND2 )
 COVAL(WALL1 ,EP , GRND2 , GRND2 )
 
 PATCH(WALL2 ,NWALL , 1, 6, 10, 10, 8, 9, 1, 1)
 COVAL(WALL2 ,W1 , GRND2 ,0. )
 COVAL(WALL2 ,KE , GRND2 , GRND2 )
 COVAL(WALL2 ,EP , GRND2 , GRND2 )
 
 PATCH(WALL3 ,NWALL , 1, 6, 10, 10, 11, 13, 1, 1)
 COVAL(WALL3 ,W1 , GRND2 ,0. )
 COVAL(WALL3 ,KE , GRND2 , GRND2 )
 COVAL(WALL3 ,EP , GRND2 , GRND2 )
 
 PATCH(WALL4-NW,NWALL , 1, 6, 2, 2, 1, 3, 1, 1)
 COVAL(WALL4-NW,W1 , GRND2 ,0. )
 COVAL(WALL4-NW,KE , GRND2 , GRND2 )
 COVAL(WALL4-NW,EP , GRND2 , GRND2 )
 
 PATCH(WALL4-SW,SWALL , 1, 6, 3, 3, 1, 3, 1, 1)
 COVAL(WALL4-SW,W1 , GRND2 ,0. )
 COVAL(WALL4-SW,KE , GRND2 , GRND2 )
 COVAL(WALL4-SW,EP , GRND2 , GRND2 )
 
 PATCH(WALL5 ,LWALL , 1, 6, 3, 6, 1, 1, 1, 1)
 COVAL(WALL5 ,V1 , GRND2 ,0. )
 COVAL(WALL5 ,KE , GRND2 , GRND2 )
 COVAL(WALL5 ,EP , GRND2 , GRND2 )
 
 PATCH(WALL6-HW,HWALL , 1, 6, 1, 2, 3, 3, 1, 1)
 COVAL(WALL6-HW,V1 , GRND2 ,0. )
 COVAL(WALL6-HW,KE , GRND2 , GRND2 )
 COVAL(WALL6-HW,EP , GRND2 , GRND2 )
 
 PATCH(WALL6-LW,LWALL , 1, 6, 1, 2, 4, 4, 1, 1)
 COVAL(WALL6-LW,V1 , GRND2 ,0. )
 COVAL(WALL6-LW,KE , GRND2 , GRND2 )
 COVAL(WALL6-LW,EP , GRND2 , GRND2 )
 
 PATCH(WHOLE ,PHASEM, 1, 6, 1, 10, 1, 13, 1, 1)
 COVAL(WHOLE ,FUEL,In-Form:source - see Grp 19)
 
 PATCH(SMOKE ,PHASEM, 1, 6, 1, 10, 1, 13, 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
 SARAH =1.
 SELREF = T
 RESFAC =1.0E-03
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(V1)=20
 LITER(W1)=20 ;LITER(KE)=20
 LITER(EP)=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(MIXF)=20
 ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03
 ENDIT(W1)=1.0E-03 ;ENDIT(KE)=1.0E-03
 ENDIT(EP)=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(MIXF)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(V1,FALSDT,0.02)
 RELAX(W1,FALSDT,0.02)
 RELAX(KE,LINRLX,0.5)
 RELAX(EP,LINRLX,0.5)
 RELAX(H1,FALSDT,10.)
 RELAX(PRPS,LINRLX,1.)
 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,1.)
 RELAX(F2,LINRLX,1.)
 RELAX(F3,LINRLX,1.)
 RELAX(F4,LINRLX,1.)
 RELAX(F5,LINRLX,1.)
 RELAX(F6,LINRLX,1.)
 RELAX(F7,LINRLX,1.)
 RELAX(F8,LINRLX,1.)
 RELAX(F9,LINRLX,1.)
 RELAX(F10,LINRLX,1.)
 RELAX(F11,LINRLX,1.)
 RELAX(MNSQ,LINRLX,1.)
 RELAX(AVEF,LINRLX,1.)
 RELAX(MIXL,LINRLX,1.)
 RELAX(RATE,LINRLX,1.)
 RELAX(DEGF,LINRLX,1.)
 RELAX(RRAT,LINRLX,1.)
 RELAX(EPKE,LINRLX,1.)
 RELAX(ENUT,LINRLX,0.1)
 RELAX(RHO1,LINRLX,1.)
 RELAX(TMP1,LINRLX,1.)
 RELAX(PROD,LINRLX,1.)
 RELAX(OXID,LINRLX,1.)
 RELAX(FUEL,LINRLX,1.)
 RELAX(MIXF,FALSDT,1.0E+09)
 KELIN = 3
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
 VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
 VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
 VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
 VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
 VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10
 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-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(MNSQ)=1. ;VARMIN(MNSQ)=0.
 VARMAX(AVEF)=1. ;VARMIN(AVEF)=0.
 VARMAX(MIXL)=1.0E+10 ;VARMIN(MIXL)=-1.0E+10
 VARMAX(RATE)=1.0E+10 ;VARMIN(RATE)=-1.0E+10
 VARMAX(DEGF)=1.0E+10 ;VARMIN(DEGF)=-1.0E+10
 VARMAX(RRAT)=1.0E+10 ;VARMIN(RRAT)=-1.0E+10
 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
 VARMAX(HPOR)=1.0E+10 ;VARMIN(HPOR)=-1.0E+10
 VARMAX(NPOR)=1.0E+10 ;VARMIN(NPOR)=-1.0E+10
 VARMAX(EPOR)=1.0E+10 ;VARMIN(EPOR)=-1.0E+10
 VARMAX(VPOR)=1.0E+10 ;VARMIN(VPOR)=-1.0E+10
 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
 VARMAX(RHO1)=1.0E+10 ;VARMIN(RHO1)=-1.0E+10
 VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=773.
 VARMAX(PROD)=1.0E+10 ;VARMIN(PROD)=0.
 VARMAX(OXID)=1.0E+10 ;VARMIN(OXID)=0.
 VARMAX(FUEL)=1.0E+10 ;VARMIN(FUEL)=0.
 VARMAX(MIXF)=1.0E+10 ;VARMIN(MIXF)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 GENK = T
 READQ1 = T
 PARSOL = F
 ISG62 = 1
 CHSOA =0.054825 ;CHSOB =1.
 SPEDAT(SET,PROPERTY,TMP1,C,=MAX(773.&(H1-4.9E+07*FUEL)/(1500.*FUE$)
 SPEDAT(SET,PROPERTY,TMP1,C,L+1500.*OXID+1500.*PROD))!IMAT<100)
 SPEDAT(SET,PROPERTY,RHO1,C,=(P1+8.0E+05)*1.0/(FUEL/16.+OXID/29.+P$)
 SPEDAT(SET,PROPERTY,RHO1,C,ROD/28.)/(TMP1*8313.4))
 SPEDAT(SET,STORED,PROD,C,=(MIXF-FUEL)*(1+17.24)!ZSLFIN)
 SPEDAT(SET,STORED,OXID,C,=1-FUEL-PROD!ZSLFIN)
 SPEDAT(SET,SOURCE,FUEL!WHOLE,C,=COVAL(1.*EPKE&0.0))
 SPEDAT(SET,STORED,RRAT,C,=1.*EPKE*FUEL!ZSLFIN)
 SPEDAT(SET,LONGNAME,TMP1,C,absolute_temperature_of_the_gas_Kelvin)
 SPEDAT(SET,STORED,DEGF,C,=(9./5.)*(TMP1-273)!ZSLFIN)
 SPEDAT(SET,LONGNAME,DEGF,C,degrees_Fahrenheit)
 SPEDAT(SET,LONGNAME,RRAT,C,rate_of_consumption_of_fuel_kg/m^3sec)
 SPEDAT(SET,STORED,MIXL,C,=LEN1!ZSLFIN)
 SPEDAT(SET,STORED,RATE,C,=EPKE!ZSLFIN)
 SPEDAT(SET,MFM,FFUIN,R,0.109649)
 SPEDAT(SET,MFM,TAIRIN,R,773.)
 SPEDAT(SET,MFM,TFUIN,R,773.)
 SPEDAT(SET,MFM,TSTOIC,R,2106.90918)
 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,11)
 SPEDAT(SET,MFM,NFLR,I,1)
 SPEDAT(SET,MFM,NFLF,I,11)
 SPEDAT(SET,MFM,CONREA,R,0.)
 SPEDAT(SET,MFM,CONMIX,R,5.)
 SPEDAT(SET,MFM,VISCON,R,0.)
 SPEDAT(SET,MATERIAL,199,L,T)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=5893. ;EX(V1)=6.939
 EX(W1)=24.77 ;EX(KE)=23.219999
 EX(EP)=1.066E+04 ;EX(H1)=1.765E+06
 EX(PRPS)=0.9385 ;EX(PRDM)=0.3405
 EX(OXDM)=0.5908 ;EX(FRFM)=7.153E-03
 EX(TMPM)=1015. ;EX(SMOS)=3.601E-06
 EX(SMOM)=7.359E-06 ;EX(F1)=0.4534
 EX(F2)=0.1159 ;EX(F3)=0.04251
 EX(F4)=0.03242 ;EX(F5)=0.03352
 EX(F6)=0.03745 ;EX(F7)=0.03901
 EX(F8)=0.05495 ;EX(F9)=0.05943
 EX(F10)=0.02534 ;EX(F11)=0.04501
 EX(MNSQ)=0.2738 ;EX(AVEF)=0.0262
 EX(MIXL)=0. ;EX(RATE)=1176.
 EX(DEGF)=0. ;EX(RRAT)=0.
 EX(EPKE)=1176. ;EX(HPOR)=0.
 EX(NPOR)=0. ;EX(EPOR)=0.
 EX(VPOR)=0. ;EX(ENUT)=5.209E-03
 EX(RHO1)=2.618 ;EX(TMP1)=983.799988
 EX(PROD)=0.3687 ;EX(OXID)=0.5638
 EX(FUEL)=5.911E-03 ;EX(MIXF)=0.02613
 ************************************************************
  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(V1,Y,N,Y,Y,Y,Y)
 OUTPUT(W1,Y,N,Y,Y,Y,Y)
 OUTPUT(KE,Y,N,Y,Y,Y,Y)
 OUTPUT(EP,Y,N,Y,Y,Y,Y)
 OUTPUT(H1,Y,N,Y,Y,Y,Y)
 OUTPUT(PRPS,Y,N,Y,N,N,N)
 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(MNSQ,Y,N,Y,N,N,N)
 OUTPUT(AVEF,Y,N,Y,N,N,N)
 OUTPUT(MIXL,Y,N,Y,N,N,N)
 OUTPUT(RATE,Y,N,Y,N,N,N)
 OUTPUT(DEGF,Y,N,Y,N,N,N)
 OUTPUT(RRAT,Y,N,Y,N,N,N)
 OUTPUT(EPKE,Y,N,Y,N,N,N)
 OUTPUT(HPOR,N,N,N,N,N,N)
 OUTPUT(NPOR,N,N,N,N,N,N)
 OUTPUT(EPOR,N,N,N,N,N,N)
 OUTPUT(VPOR,N,N,N,N,N,N)
 OUTPUT(ENUT,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,Y,Y,Y)
 OUTPUT(MIXF,Y,N,Y,Y,Y,Y)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 3 ;IYMON = 5 ;IZMON = 5
 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 500 ;NUMCLS = 5
 NYPRIN = 1 ;IYPRF = 1 ;IYPRL = 10000
 NZPRIN = 2 ;IZPRF = 1 ;IZPRL = 10000
 XZPR = F ;YZPR = T
 IPLTF = 2 ;IPLTL = 100 ;NPLT = 5
 ISWPRF = 1 ;ISWPRL = 100000
 ITABL = 1 ;IPROF = 1
 ABSIZ =0.5 ;ORSIZ =0.4
 NTZPRF = 1 ;NCOLPF = 50
 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 40
 
 PATCH(IZ4 ,CONTUR, 1, 6, 1, 10, 4, 4, 1, 1)
 PLOT(IZ4 ,TMP1,1. ,10. )
 
 PATCH(XSECIN1 ,CONTUR, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(XSECIN1 ,W1 ,0. ,10. )
 PLOT(XSECIN1 ,TMP1,0. ,10. )
 PLOT(XSECIN1 ,FUEL,0. ,10. )
 PLOT(XSECIN1 ,MIXF,0. ,10. )
 
 PATCH(XSECIN2 ,CONTUR, 4, 4, 1, 10, 1, 13, 1, 1)
 PLOT(XSECIN2 ,W1 ,0. ,10. )
 PLOT(XSECIN2 ,TMP1,0. ,10. )
 PLOT(XSECIN2 ,FUEL,0. ,10. )
 PLOT(XSECIN2 ,MIXF,0. ,10. )
 
 PATCH(PROF1 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(PROF1 ,F1 ,0. ,0. )
 
 PATCH(PROF2 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(PROF2 ,F2 ,0. ,0. )
 
 PATCH(PROF3 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(PROF3 ,F3 ,0. ,0. )
 
 PATCH(PROF4 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(PROF4 ,F4 ,0. ,0. )
 
 PATCH(PROF5 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(PROF5 ,F5 ,0. ,0. )
 
 PATCH(PROF6 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(PROF6 ,F6 ,0. ,0. )
 
 PATCH(PROF7 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(PROF7 ,F7 ,0. ,0. )
 
 PATCH(PROF8 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(PROF8 ,F8 ,0. ,0. )
 
 PATCH(PROF9 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(PROF9 ,F9 ,0. ,0. )
 
 PATCH(PROF10 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(PROF10 ,F10 ,0. ,0. )
 
 PATCH(PROF11 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(PROF11 ,F11 ,0. ,0. )
 
 PATCH(MIDDLE ,PROFIL, 1, 1, 5, 5, 1, 13, 1, 1)
 PLOT(MIDDLE ,MNSQ,0. ,0. )
 PLOT(MIDDLE ,AVEF,0. ,0. )
 PLOT(MIDDLE ,ENUT,0. ,0. )
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP