TALK=T;RUN(1,1)
  PHOTON USE
  P;phi; 10 1;;

  msg(  Single fluid central-jet concentration contours
  con h1 x 1 fil;.001
  pause
  con cl; red
  msg(  Averaged 17-fluid contours
  con cav x 1 fil;.001
  pause
  con cl; red
  msg( Concentration fluctuations by transport equation
  con gg  x 1 fil;.001
  pause
  con cl; red
  msg( Averaged 17-fluid concentration fluctuation
  con gav x 1 fil;.001
  msg
  msg Hit Enter for FPD hystogram
  msg
  pause
   p

   20 1

   con fpd x 1 fil;10.
   msg Hit Enter to continue
  ENDUSE
  DISPLAY
     In this  case,  PLANT is used to introduce a turbulence
     model for concentration fluctuations which  employs  no
     conservation  equations  for  statistical properties of
     the  fluctuations.  It  is  a  variant  of  Multi-Fluid
     concept of Brian Spalding as employed by Sergei Zhubrin.

     The 17-fluid  model  is considered here to simulate the
     turbulent mixing resulting from the  admission  of  two
     separate,   isothermal   coaxial   jets   of  different
     composition  into  a  concentric   duct   as   depicted
     diagrammatically below.
  ENDDIS
    PLANTBEGIN
  **Source term for g
PATCH(SORG,VOLUME,1,NX,1,NY,1,NZ,1,1)
   CO=2.0*:RHO1:*EPKE
   VAL=GENG/(2.0*:RHO1:*EPKE+TINY)
COVAL(SORG,G  , GRND  ,GRND )
    Above statements  contain  the  formulae  for   combined
    source/sink  term  for the production and dissipation of
    the concentration fluctuations.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
PATCH(WG,VOLUME,1,NX,NY,NY,1,NZ,1,1)
   VAL=GENG/(2.0*:RHO1:*EPKE+TINY)
COVAL(WG,G  , FIXVAL  ,GRND )
    Production of G is made  equal  to  its  dissipation  at
    North-Wall  boundary  by  FIXVALing  its  value  to  the
    production rate divided by twice product of density  and
    EPKE.  The  latter  is  built-in  variable  standing for
    EP/KE.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
  ** Provide the re-calculation of reference residuals for G
   RES=SUM(VOL*(GENG-2.*:RHO1:*EPKE*G)/(NY*NZ))
   RESREF(G)=RES
    By above two statements the reference residuals for G is
    calculated at  the  and  of  z-slab  as  a  sum  of  its
    generation rate per cell.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
   **Calculation of GENG
     * Auxiliary variables
STORE(DFZ,DFY,DFZH,DFYN)
FIINIT(DFZ) =0.0;FIINIT(DFY) =0.0
FIINIT(DFZH)=0.0;FIINIT(DFYN)=0.0
   DFZ=((H1[,,+1]-H1)/DZGNZ)**2
   REGION(1,NX,1,NY,1,NZ-1,1,1)
   DFY=((H1[,+1,]-H1)/DYG2D)**2
   REGION(1,1,1,NY-1,1,NZ,1,1)
   DFZH=((H1-H1[,,-1])/DZGNZ[,,-1])**2
   REGION(1,NX,1,NY,NZ,NZ,1,1)
   DFYN=((H1-H1[,-1,])/DYG2D[,-1,])**2
   REGION(1,NX,NY,NY,1,NZ,1,1)
   The above statements calculate the  square  concentration
   derivatives  separately  for  internal  and  near  domain
   bounadary regions.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
   GENG=2.8*:RHO1:*ENUT*(DFZ+DFY+DFZH+DFYN)
   The sum of the radial  and  longitudinal  derivatives  is
   multiplied  by  density and turbulent viscosity times 2.8
   to get the generation term.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
  **Output calculations
    GG  - concentration fluctuation;
    GGF - concentration fluctuation normalised by local
          concentration of central jet fluid.
STORE(GG,GGF)
FIINIT(GGF)=0.0
   GG=SQRT(G)
  IF(ISWEEP.EQ.LSWEEP)
   GGF=GG/(H1+TINY)
  IF(ISWEEP.EQ.LSWEEP)
    At the   end   of   z-slab   for   the  last  sweep  the
    concentration fluctuation is calculated  and  normalized
    by the local average concentration.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<

     ********************     MFTM section    ***************
     ** Number of fluids in population
INTEGER(NFLUIDS)
NFLUIDS=17
     ** Micro-mixing constant
REAL(MMC)
MMC = 5.0 ; RG(1) = MMC
     ** Solve for fluid mass fractions F1, F2, ..., F17
DO II=1,NFLUIDS
 SOLVE(F:II:)
 TERMS(F:II:,N,Y,y,y,y,y)
 PRT(F:II:)= 0.86;PRNDTL(F:II:)= 0.71
 RELAX(F:II:,linrlx,0.15)
 VARMIN(F:II:)=0.0;VARMAX(F:II:)=1.0
 PATCH(PROF:II:,PROFIL,1,1,1,1,1,20,1,1)
 PLOT(PROF:II:,F:II:,0.000E+00, 0.000E+00)
ENDDO
ABSIZ=0.5; ORSIZ=0.2
     ** Fluid population boundary conditions
INLET(IN1,LOW,1,1,1,3,1,1,1,1)
INLET(IN2,LOW,1,1,4,NY,1,1,1,1)
DO II=1,NFLUIDS
 VALUE(IN1,F:II:,0.0)
 VALUE(IN2,F:II:,0.0)
ENDDO
VALUE(IN1,F1 , 1.0); fiinit(f1)=0
VALUE(IN2,F:NFLUIDS:, 1.0)
     ** Coupling/splitting rates
PATCH(MIX,PHASEM,1,NX,1,NY,1,NZ,1,1)
      *  Fluid 1
   CO=RG(1)*EPKE*(F3+F5+F7+F9+F11+F13+F15+F17)
COVAL(MIX,F1   , GRND  ,0.0 )
      *  Fluid 2
   VAL=2.*RG(1)*EPKE*(F1*F3)-           RG(1)*EPKE*(F4+F6+F$
8+F10+F12+F14+F16)*F2
COVAL(MIX,F2, FIXFLU,GRND)
       *  Fluid 3
   VAL=2.*RG(1)*EPKE*(F2*F4+F1*F5)-           RG(1)*EPKE*(F$
1+F17+F5+F7+F9+F11+F13+F15)*F3
COVAL(MIX,F3, FIXFLU,GRND)
       *  Fluid 4
   VAL=2.*RG(1)*EPKE*(F3*F5+F2*F6+F1*F7)-           RG(1)*E$
PKE*(F2+F6+F8+F10+F12+F14+F16)*F4
COVAL(MIX,F4, FIXFLU,GRND)
       *  Fluid 5
   VAL=2.*RG(1)*EPKE*(F4*F6+F3*F7+F2*F8+F1*F9)-           R$
G(1)*EPKE*(F1+F3+F17+F7+F9+F11+F13+F15)*F5
COVAL(MIX,F5, FIXFLU,GRND)
       *  Fluid 6
   VAL=2.*RG(1)*EPKE*(F5*F7+F4*F8+F3*F9+$
              F2*F10+F1*F11)-           RG(1)*EPKE*(F2+F4+F8+F10+F1$
2+F14+F16)*F6
           RG(1)*EPKE*(F2+F4+F8+F10+F12+F14+F16)*F6
COVAL(MIX,F6, FIXFLU,GRND)
        * Fluid 7
   VAL=2.*RG(1)*EPKE*(F6*F8+F5*F9+F4*F10+$
         F3*F11+F2*F12+F1*F13)-           RG(1)*EPKE*(F1+F3+F5+F17+$
F9+F11+F13+F15)*F7
           RG(1)*EPKE*(F1+F3+F5+F17+F9+F11+F13+F15)*F7
COVAL(MIX,F7, FIXFLU,GRND)
         * Fluid 8
   VAL=2.*RG(1)*EPKE*(F7*F9+F6*F10+F5*F11+$
    F4*F12+F3*F13+F2*F14+F1*F15)-           RG(1)*EPKE*(F2+F4+F6+F1$
0+F12+F14+F16)*F8
           RG(1)*EPKE*(F2+F4+F6+F10+F12+F14+F16)*F8
COVAL(MIX,F8, FIXFLU,GRND)
         * Fluid 9
   VAL=2.*RG(1)*EPKE*(F8*F10+F7*F11+F6*F12+               F$
5*F13+F4*F14+F3*F15+F2*F16+F1*F17)-           RG(1)*EPKE*(F1+F3+F5+$
F7+F17+F11+F13+F15)*F9
           RG(1)*EPKE*(F1+F3+F5+F7+F17+F11+F13+F15)*F9
COVAL(MIX,F9, FIXFLU,GRND)
         * Fluid 10
   VAL=2.*RG(1)*EPKE*(F9*F11+F8*F12+F7*F13+$
      F6*F14+F5*F15+F4*F16+F3*F17)-           RG(1)*EPKE*(F2+F4+F6+$
F8+F12+F14+F16)*F10
           RG(1)*EPKE*(F2+F4+F6+F8+F12+F14+F16)*F10
COVAL(MIX,F10, FIXFLU,GRND)
         * Fluid 11
   VAL=2.*RG(1)*EPKE*(F10*F12+F9*F13+F8*F14+$
               F7*F15+F6*F16+F5*F17)-           RG(1)*EPKE*(F1+F3+F$
5+F7+F9+F17+F13+F15)*F11
           RG(1)*EPKE*(F1+F3+F5+F7+F9+F17+F13+F15)*F11
COVAL(MIX,F11, FIXFLU,GRND)
         * Fluid 12
   VAL=2.*RG(1)*EPKE*(F11*F13+F10*F14+F9*F15+$
                        F8*F16+F7*F17)-           RG(1)*EPKE*(F2+F4$
+F6+F8+F10+F14+F16)*F12
           RG(1)*EPKE*(F2+F4+F6+F8+F10+F14+F16)*F12
COVAL(MIX,F12, FIXFLU,GRND)
         * Fluid 13
   VAL=2.*RG(1)*EPKE*(F12*F14+F11*F15+$
         F10*F16+F9*F17)-           RG(1)*EPKE*(F1+F3+F5+F7+F9+F11+$
F15+F17)*F13
           RG(1)*EPKE*(F1+F3+F5+F7+F9+F11+F15+F17)*F13
COVAL(MIX,F13, FIXFLU,GRND)
         * Fluid 14
   VAL=2.*RG(1)*EPKE*(F13*F15+F12*F16+F11*F17)-           R$
G(1)*EPKE*(F2+F4+F6+F8+F10+F12+F16)*F14
COVAL(MIX,F14, FIXFLU,GRND)
         * Fluid 15
   VAL=2.*RG(1)*EPKE*(F14*F16+F13*F17)-           RG(1)*EPK$
E*(F1+F3+F5+F7+F9+F11+F13+F17)*F15
COVAL(MIX,F15, FIXFLU,GRND)
         * Fluid 16
   VAL=2.*RG(1)*EPKE*(F15*F17)-           RG(1)*EPKE*(F2+F4$
+F6+F8+F10+F12+F14)*F16
COVAL(MIX,F16, FIXFLU,GRND)
         * Fluid 17
   CO=RG(1)*EPKE*(F1+F3+F5+F7+F9+F11+F13+F15)
COVAL(MIX,F17  , GRND  ,0.0 )
    The above source/sink terms in  the  fluid-mass-fraction
    equations  are  shared according to a coupling/splitting
    scheme derived from Spalding concept.

    The scheme hypotheses is that the coupling may only take
    place  between  those  parent fluids which would produce
    the appropriate offsprings inheriting the ATTRIBUTES  of
    either parent in EQUAL proportion.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
  ** Output calculations
STORE(CAV,MAS,GAV,GF)
FIINIT(GF)=0.0
   CAV=16./16.*F1 + 15./16.*F2 + 14./16.*F3 +$
 13./16.*F4 + 12./16.*F5 + 11./16.*F6 +               10./16.*F7 +$
 9./16.*F8 +  8./16.*F9 +                7./16.*F10+  6./16.*F11+$
5./16.*F12+$
               10./16.*F7 +  9./16.*F8 +  8./16.*F9 +$
  7./16.*F10+  6./16.*F11+  5./16.*F12+                4./16.*F13+$
 3./16.*F14+  2./16.*F15+                1./16.*F16+  0./16.*F17
                4./16.*F13+  3./16.*F14+  2./16.*F15+$
  1./16.*F16+  0./16.*F17
  IF(ISWEEP.EQ.LSWEEP)
   At the end of  the  iz-slab  for  the  last  sweep,  CAV,
    averaged   concentration   of   central  jet  fluid,  is
    calculated from the individual  fluid  mass-fractionsand
    their arrributes;
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
   MAS=F1+F2+F3+F4+F5+F6+F7+F8+F9+F10+F11+F12+$
  F13+F14+F15+F16+F17
  IF(ISWEEP.EQ.LSWEEP)
   At the end of the iz-slab for the last sweep, MAS, sum of
   fluid mass fractions, is calculated to check its equality
   to unity;
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
   GAV=ABS(CAV-16./16)*F1 + ABS(CAV-15./16.)*F2 +$
     ABS(CAV-14./16)*F3 + ABS(CAV-13./16.)*F4 +               ABS(C$
AV-12./16)*F5 + ABS(CAV-11./16.)*F6 +               ABS(CAV-10./16)$
*F7 + ABS(CAV- 9./16.)*F8 +$
               ABS(CAV-12./16)*F5 + ABS(CAV-11./16.)*F6 +$
     ABS(CAV-10./16)*F7 + ABS(CAV- 9./16.)*F8 +               ABS(C$
AV- 8./16)*F9 + ABS(CAV- 7./16.)*F10+               ABS(CAV- 6./16)$
*F11+ ABS(CAV- 5./16.)*F12+$
               ABS(CAV- 8./16)*F9 + ABS(CAV- 7./16.)*F10+$
     ABS(CAV- 6./16)*F11+ ABS(CAV- 5./16.)*F12+               ABS(C$
AV- 4./16)*F13+ ABS(CAV- 3./16.)*F14+               ABS(CAV- 2./16)$
*F15+ ABS(CAV- 1./16.)*F16+$
               ABS(CAV- 4./16)*F13+ ABS(CAV- 3./16.)*F14+$
     ABS(CAV- 2./16)*F15+ ABS(CAV- 1./16.)*F16+               ABS(C$
AV- 0./16)*F17
               ABS(CAV- 0./16)*F17
  IF(ISWEEP.EQ.LSWEEP)
   At the end of  the  iz-slab  for  the  last  sweep,  GAV,
   averaged concentration fluctuation,  is calculated as the
   sum of local deviations of averaged  concentrations  from
   the individual concentration attributes;
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
   GF =GAV/(CAV+TINY)
  IF(ISWEEP.EQ.LSWEEP)
   At the end  of  the  iz-slab  for  the  last  sweep,  GF,
   averaged  concentration  fluctuation  normalised by local
   averaged  concentration  of   central   jet   fluid,   is
   calculated.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
    ** Output data processing for plotting PDF
     * Specify the cell in question: IY=IG(1), IZ=IG(2)
IG(1)=4; IG(2)=4
STORE(FPD);FIINIT(FPD)=0.0
   FPD=F1[1,IG(1),IG(2)]*AMAX1(ABS(F1[1,IG(1),IG(2)]$
              -YV2D)/(F1[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
   REGION(1,1,1,NY,1,1) /ISWEEP.EQ.LSWEEP
   FPD=F2[1,IG(1),IG(2)]*AMAX1(ABS(F2[1,IG(1),IG(2)]$
              -YV2D)/(F2[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
   REGION(1,1,1,NY,2,2) /ISWEEP.EQ.LSWEEP
   FPD=F3[1,IG(1),IG(2)]*AMAX1(ABS(F3[1,IG(1),IG(2)]$
              -YV2D)/(F3[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
   REGION(1,1,1,NY,3,3) /ISWEEP.EQ.LSWEEP
   FPD=F4[1,IG(1),IG(2)]*AMAX1(ABS(F4[1,IG(1),IG(2)]$
              -YV2D)/(F4[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
   REGION(1,1,1,NY,4,4) /ISWEEP.EQ.LSWEEP
   FPD=F5[1,IG(1),IG(2)]*AMAX1(ABS(F5[1,IG(1),IG(2)]$
              -YV2D)/(F5[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
   REGION(1,1,1,NY,5,5) /ISWEEP.EQ.LSWEEP
   FPD=F6[1,IG(1),IG(2)]*AMAX1(ABS(F6[1,IG(1),IG(2)]$
              -YV2D)/(F6[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
   REGION(1,1,1,NY,6,6) /ISWEEP.EQ.LSWEEP
   FPD=F7[1,IG(1),IG(2)]*AMAX1(ABS(F7[1,IG(1),IG(2)]$
              -YV2D)/(F7[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
   REGION(1,1,1,NY,7,7) /ISWEEP.EQ.LSWEEP
   FPD=F8[1,IG(1),IG(2)]*AMAX1(ABS(F8[1,IG(1),IG(2)]$
              -YV2D)/(F8[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
   REGION(1,1,1,NY,8,8) /ISWEEP.EQ.LSWEEP
   FPD=F9[1,IG(1),IG(2)]*AMAX1(ABS(F9[1,IG(1),IG(2)]$
              -YV2D)/(F9[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
   REGION(1,1,1,NY,9,9) /ISWEEP.EQ.LSWEEP
   FPD=F10[1,IG(1),IG(2)]*AMAX1(ABS(F10[1,IG(1),IG(2)]$
                 -YV2D)/(F10[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
   REGION(1,1,1,NY,10,10) /ISWEEP.EQ.LSWEEP
   FPD=F11[1,IG(1),IG(2)]*AMAX1(ABS(F11[1,IG(1),IG(2)]$
                 -YV2D)/(F11[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
   REGION(1,1,1,NY,11,11) /ISWEEP.EQ.LSWEEP
   FPD=F12[1,IG(1),IG(2)]*AMAX1(ABS(F12[1,IG(1),IG(2)]$
                 -YV2D)/(F12[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
   REGION(1,1,1,NY,12,12) /ISWEEP.EQ.LSWEEP
   FPD=F13[1,IG(1),IG(2)]*AMAX1(ABS(F13[1,IG(1),IG(2)]$
                 -YV2D)/(F13[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
   REGION(1,1,1,NY,13,13) /ISWEEP.EQ.LSWEEP
   FPD=F14[1,IG(1),IG(2)]*AMAX1(ABS(F14[1,IG(1),IG(2)]$
                 -YV2D)/(F14[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
   REGION(1,1,1,NY,14,14) /ISWEEP.EQ.LSWEEP
   FPD=F15[1,IG(1),IG(2)]*AMAX1(ABS(F15[1,IG(1),IG(2)]$
                 -YV2D)/(F15[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
   REGION(1,1,1,NY,15,15) /ISWEEP.EQ.LSWEEP
   FPD=F16[1,IG(1),IG(2)]*AMAX1(ABS(F16[1,IG(1),IG(2)]$
                 -YV2D)/(F16[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
   REGION(1,1,1,NY,16,16) /ISWEEP.EQ.LSWEEP
   FPD=F17[1,IG(1),IG(2)]*AMAX1(ABS(F17[1,IG(1),IG(2)]$
                 -YV2D)/(F17[1,IG(1),IG(2)]+0.-YV2D) ,0.0)
   REGION(1,1,1,NY,17,17) /ISWEEP.EQ.LSWEEP
   The above  operations  are made at the end of the iz-slab
   for the last sweep to fill each IY-column of  the  domain
   by Fi value.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
    PLANTEND
  PHOTON USE
  P;phi; 10 1;;

  msg(  Single fluid central-jet concentration contours
  con h1 x 1 fil;.001
  pause
  con cl; red
  msg(  Averaged 17-fluid contours
  con cav x 1 fil;.001
  pause
  con cl; red
  msg( Concentration fluctuations by transport equation
  con gg  x 1 fil;.001
  pause
  con cl; red
  msg( Averaged 17-fluid concentration fluctuation
  con gav x 1 fil;.001
  msg
  msg Hit Enter for FPD hystogram
  msg
  pause
   p

   20 1

   con fpd x 1 fil;10.
   msg Hit Enter to continue
  ENDUSE
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(17-fluid turbulence model               )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 781
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = F
 NX = 1
 XULAST =0.1
 XFRAC(1)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 15
 YVLAST =0.15
 YFRAC(1)=0.066667 ;YFRAC(2)=0.133333
 YFRAC(3)=0.2 ;YFRAC(4)=0.266667
 YFRAC(5)=0.333333 ;YFRAC(6)=0.4
 YFRAC(7)=0.466667 ;YFRAC(8)=0.533333
 YFRAC(9)=0.6 ;YFRAC(10)=0.666667
 YFRAC(11)=0.733333 ;YFRAC(12)=0.8
 YFRAC(13)=0.866667 ;YFRAC(14)=0.933333
 YFRAC(15)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 20
 ZWLAST =3.
 ZFRAC(1)=0.05 ;ZFRAC(2)=0.1
 ZFRAC(3)=0.15 ;ZFRAC(4)=0.2
 ZFRAC(5)=0.25 ;ZFRAC(6)=0.3
 ZFRAC(7)=0.35 ;ZFRAC(8)=0.4
 ZFRAC(9)=0.45 ;ZFRAC(10)=0.5
 ZFRAC(11)=0.55 ;ZFRAC(12)=0.6
 ZFRAC(13)=0.65 ;ZFRAC(14)=0.7
 ZFRAC(15)=0.75 ;ZFRAC(16)=0.8
 ZFRAC(17)=0.85 ;ZFRAC(18)=0.9
 ZFRAC(19)=0.95 ;ZFRAC(20)=1.
 ************************************************************
  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(117)=FPD ;NAME(118)=GF
 NAME(119)=GAV ;NAME(120)=MAS
 NAME(121)=CAV ;NAME(122)=F17
 NAME(123)=F16 ;NAME(124)=F15
 NAME(125)=F14 ;NAME(126)=F13
 NAME(127)=F12 ;NAME(128)=F11
 NAME(129)=F10 ;NAME(130)=F9
 NAME(131)=F8 ;NAME(132)=F7
 NAME(133)=F6 ;NAME(134)=F5
 NAME(135)=F4 ;NAME(136)=F3
 NAME(137)=F2 ;NAME(138)=F1
 NAME(139)=GGF ;NAME(140)=GG
 NAME(141)=DFYN ;NAME(142)=DFZH
 NAME(143)=DFY ;NAME(144)=DFZ
 NAME(145)=GENG ;NAME(146)=EPKE
 NAME(147)=GEN1 ;NAME(148)=LEN1
 NAME(149)=ENUT ;NAME(150)=G
    * 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(FPD,Y,N,N,N,N,Y)
 SOLUTN(GF,Y,N,N,N,N,Y)
 SOLUTN(GAV,Y,N,N,N,N,Y)
 SOLUTN(MAS,Y,N,N,N,N,Y)
 SOLUTN(CAV,Y,N,N,N,N,Y)
 SOLUTN(F17,Y,Y,N,N,N,Y)
 SOLUTN(F16,Y,Y,N,N,N,Y)
 SOLUTN(F15,Y,Y,N,N,N,Y)
 SOLUTN(F14,Y,Y,N,N,N,Y)
 SOLUTN(F13,Y,Y,N,N,N,Y)
 SOLUTN(F12,Y,Y,N,N,N,Y)
 SOLUTN(F11,Y,Y,N,N,N,Y)
 SOLUTN(F10,Y,Y,N,N,N,Y)
 SOLUTN(F9,Y,Y,N,N,N,Y)
 SOLUTN(F8,Y,Y,N,N,N,Y)
 SOLUTN(F7,Y,Y,N,N,N,Y)
 SOLUTN(F6,Y,Y,N,N,N,Y)
 SOLUTN(F5,Y,Y,N,N,N,Y)
 SOLUTN(F4,Y,Y,N,N,N,Y)
 SOLUTN(F3,Y,Y,N,N,N,Y)
 SOLUTN(F2,Y,Y,N,N,N,Y)
 SOLUTN(F1,Y,Y,N,N,N,Y)
 SOLUTN(GGF,Y,N,N,N,N,Y)
 SOLUTN(GG,Y,N,N,N,N,Y)
 SOLUTN(DFYN,Y,N,N,N,N,Y)
 SOLUTN(DFZH,Y,N,N,N,N,Y)
 SOLUTN(DFY,Y,N,N,N,N,Y)
 SOLUTN(DFZ,Y,N,N,N,N,Y)
 SOLUTN(GENG,Y,N,N,N,N,Y)
 SOLUTN(EPKE,Y,N,N,N,N,Y)
 SOLUTN(GEN1,Y,N,N,N,N,Y)
 SOLUTN(LEN1,Y,N,N,N,N,Y)
 SOLUTN(ENUT,Y,N,N,N,N,Y)
 SOLUTN(G,Y,Y,N,N,N,Y)
 VIST = 149
 LEN1 = 148
 ************************************************************
  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,Y,Y,Y)
 TERMS(F17,N,Y,Y,Y,Y,Y)
 TERMS(F16,N,Y,Y,Y,Y,Y)
 TERMS(F15,N,Y,Y,Y,Y,Y)
 TERMS(F14,N,Y,Y,Y,Y,Y)
 TERMS(F13,N,Y,Y,Y,Y,Y)
 TERMS(F12,N,Y,Y,Y,Y,Y)
 TERMS(F11,N,Y,Y,Y,Y,Y)
 TERMS(F10,N,Y,Y,Y,Y,Y)
 TERMS(F9,N,Y,Y,Y,Y,Y)
 TERMS(F8,N,Y,Y,Y,Y,Y)
 TERMS(F7,N,Y,Y,Y,Y,Y)
 TERMS(F6,N,Y,Y,Y,Y,Y)
 TERMS(F5,N,Y,Y,Y,Y,Y)
 TERMS(F4,N,Y,Y,Y,Y,Y)
 TERMS(F3,N,Y,Y,Y,Y,Y)
 TERMS(F2,N,Y,Y,Y,Y,Y)
 TERMS(F1,N,Y,Y,Y,Y,Y)
 TERMS(G,N,Y,Y,Y,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWENT = T
 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 = GRND4
 TSURR =0. ;TEMP0 =0.
 PRESS0 =0.
 DVO1DT =0. ;DRH1DP =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 EL1A =0. ;EL1B =0.
 EL1C =0.
 ENUL =3.0E-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)=0.71 ;PRNDTL(F17)=0.71
 PRNDTL(F16)=0.71 ;PRNDTL(F15)=0.71
 PRNDTL(F14)=0.71 ;PRNDTL(F13)=0.71
 PRNDTL(F12)=0.71 ;PRNDTL(F11)=0.71
 PRNDTL(F10)=0.71 ;PRNDTL(F9)=0.71
 PRNDTL(F8)=0.71 ;PRNDTL(F7)=0.71
 PRNDTL(F6)=0.71 ;PRNDTL(F5)=0.71
 PRNDTL(F4)=0.71 ;PRNDTL(F3)=0.71
 PRNDTL(F2)=0.71 ;PRNDTL(F1)=0.71
 PRNDTL(G)=0.7
 PRT(V1)=1. ;PRT(W1)=1.
 PRT(KE)=1. ;PRT(EP)=1.314
 PRT(H1)=0.86 ;PRT(F17)=0.86
 PRT(F16)=0.86 ;PRT(F15)=0.86
 PRT(F14)=0.86 ;PRT(F13)=0.86
 PRT(F12)=0.86 ;PRT(F11)=0.86
 PRT(F10)=0.86 ;PRT(F9)=0.86
 PRT(F8)=0.86 ;PRT(F7)=0.86
 PRT(F6)=0.86 ;PRT(F5)=0.86
 PRT(F4)=0.86 ;PRT(F3)=0.86
 PRT(F2)=0.86 ;PRT(F1)=0.86
 PRT(G)=0.7
 CP1 =1. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1.3E-04 ;FIINIT(V1)=1.0E-10
 FIINIT(W1)=6. ;FIINIT(KE)=0.234
 FIINIT(EP)=15.149659 ;FIINIT(H1)=1.
 FIINIT(FPD)=0. ;FIINIT(GF)=0.
 FIINIT(GAV)=1.0E-10 ;FIINIT(MAS)=1.0E-10
 FIINIT(CAV)=1.0E-10 ;FIINIT(F17)=1.0E-10
 FIINIT(F16)=1.0E-10 ;FIINIT(F15)=1.0E-10
 FIINIT(F14)=1.0E-10 ;FIINIT(F13)=1.0E-10
 FIINIT(F12)=1.0E-10 ;FIINIT(F11)=1.0E-10
 FIINIT(F10)=1.0E-10 ;FIINIT(F9)=1.0E-10
 FIINIT(F8)=1.0E-10 ;FIINIT(F7)=1.0E-10
 FIINIT(F6)=1.0E-10 ;FIINIT(F5)=1.0E-10
 FIINIT(F4)=1.0E-10 ;FIINIT(F3)=1.0E-10
 FIINIT(F2)=1.0E-10 ;FIINIT(F1)=0.
 FIINIT(GGF)=0. ;FIINIT(GG)=1.0E-10
 FIINIT(DFYN)=0. ;FIINIT(DFZH)=0.
 FIINIT(DFY)=0. ;FIINIT(DFZ)=0.
 FIINIT(GENG)=1.0E-10 ;FIINIT(EPKE)=1.0E-10
 FIINIT(GEN1)=1.0E-10 ;FIINIT(LEN1)=0.015
 FIINIT(ENUT)=0.015 ;FIINIT(G)=1.0E-10
   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(KESOURCE,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1)
 COVAL(KESOURCE,KE , GRND4 , GRND4 )
 COVAL(KESOURCE,EP , GRND4 , GRND4 )
 
 PATCH(IN1 ,LOW , 1, 1, 1, 3, 1, 1, 1, 1)
 COVAL(IN1 ,P1 , FIXFLU ,10. )
 COVAL(IN1 ,V1 ,0. ,0. )
 COVAL(IN1 ,W1 ,0. ,10. )
 COVAL(IN1 ,KE ,0. ,0.45 )
 COVAL(IN1 ,EP ,0. ,30.058847 )
 COVAL(IN1 ,H1 ,0. ,1. )
 COVAL(IN1 ,F17 ,0. ,0. )
 COVAL(IN1 ,F16 ,0. ,0. )
 COVAL(IN1 ,F15 ,0. ,0. )
 COVAL(IN1 ,F14 ,0. ,0. )
 COVAL(IN1 ,F13 ,0. ,0. )
 COVAL(IN1 ,F12 ,0. ,0. )
 COVAL(IN1 ,F11 ,0. ,0. )
 COVAL(IN1 ,F10 ,0. ,0. )
 COVAL(IN1 ,F9 ,0. ,0. )
 COVAL(IN1 ,F8 ,0. ,0. )
 COVAL(IN1 ,F7 ,0. ,0. )
 COVAL(IN1 ,F6 ,0. ,0. )
 COVAL(IN1 ,F5 ,0. ,0. )
 COVAL(IN1 ,F4 ,0. ,0. )
 COVAL(IN1 ,F3 ,0. ,0. )
 COVAL(IN1 ,F2 ,0. ,0. )
 COVAL(IN1 ,F1 ,0. ,1. )
 COVAL(IN1 ,G ,0. ,0. )
 
 PATCH(IN2 ,LOW , 1, 1, 4, 15, 1, 1, 1, 1)
 COVAL(IN2 ,P1 , FIXFLU ,2. )
 COVAL(IN2 ,V1 ,0. ,0. )
 COVAL(IN2 ,W1 ,0. ,2. )
 COVAL(IN2 ,KE ,0. ,0.018 )
 COVAL(IN2 ,EP ,0. ,0.240471 )
 COVAL(IN2 ,H1 ,0. ,0. )
 COVAL(IN2 ,F17 ,0. ,1. )
 COVAL(IN2 ,F16 ,0. ,0. )
 COVAL(IN2 ,F15 ,0. ,0. )
 COVAL(IN2 ,F14 ,0. ,0. )
 COVAL(IN2 ,F13 ,0. ,0. )
 COVAL(IN2 ,F12 ,0. ,0. )
 COVAL(IN2 ,F11 ,0. ,0. )
 COVAL(IN2 ,F10 ,0. ,0. )
 COVAL(IN2 ,F9 ,0. ,0. )
 COVAL(IN2 ,F8 ,0. ,0. )
 COVAL(IN2 ,F7 ,0. ,0. )
 COVAL(IN2 ,F6 ,0. ,0. )
 COVAL(IN2 ,F5 ,0. ,0. )
 COVAL(IN2 ,F4 ,0. ,0. )
 COVAL(IN2 ,F3 ,0. ,0. )
 COVAL(IN2 ,F2 ,0. ,0. )
 COVAL(IN2 ,F1 ,0. ,0. )
 COVAL(IN2 ,G ,0. ,0. )
 
 PATCH(OUTLET ,HIGH , 1, 1, 1, 15, 20, 20, 1, 1)
 COVAL(OUTLET ,P1 ,1.0E+05 ,0. )
 COVAL(OUTLET ,V1 ,0. ,0. )
 COVAL(OUTLET ,W1 ,0. ,0. )
 COVAL(OUTLET ,KE ,0. ,0. )
 COVAL(OUTLET ,EP ,0. ,0. )
 
 PATCH(WFNN ,NWALL , 1, 1, 15, 15, 1, 20, 1, 1)
 COVAL(WFNN ,W1 , GRND2 ,0. )
 COVAL(WFNN ,KE , GRND2 , GRND2 )
 COVAL(WFNN ,EP , GRND2 , GRND2 )
 
 PATCH(ISORG ,VOLUME, 1, 1, 1, 15, 1, 20, 1, 1)
 COVAL(ISORG ,G ,In-Form:source - see Grp 19)
 
 PATCH(IWG ,VOLUME, 1, 1, 15, 15, 1, 20, 1, 1)
 COVAL(IWG ,G ,In-Form:source - see Grp 19)
 
 PATCH(PAT1 ,CELL , 1, 1, 1, 15, 1, 19, 1, 1)
 
 PATCH(PAT2 ,CELL , 1, 1, 1, 14, 1, 20, 1, 1)
 
 PATCH(PAT3 ,CELL , 1, 1, 1, 15, 20, 20, 1, 1)
 
 PATCH(PAT4 ,CELL , 1, 1, 15, 15, 1, 20, 1, 1)
 
 PATCH(IMIX ,PHASEM, 1, 1, 1, 15, 1, 20, 1, 1)
 COVAL(IMIX ,F17 ,In-Form:source - see Grp 19)
 COVAL(IMIX ,F16 ,In-Form:source - see Grp 19)
 COVAL(IMIX ,F15 ,In-Form:source - see Grp 19)
 COVAL(IMIX ,F14 ,In-Form:source - see Grp 19)
 COVAL(IMIX ,F13 ,In-Form:source - see Grp 19)
 COVAL(IMIX ,F12 ,In-Form:source - see Grp 19)
 COVAL(IMIX ,F11 ,In-Form:source - see Grp 19)
 COVAL(IMIX ,F10 ,In-Form:source - see Grp 19)
 COVAL(IMIX ,F9 ,In-Form:source - see Grp 19)
 COVAL(IMIX ,F8 ,In-Form:source - see Grp 19)
 COVAL(IMIX ,F7 ,In-Form:source - see Grp 19)
 COVAL(IMIX ,F6 ,In-Form:source - see Grp 19)
 COVAL(IMIX ,F5 ,In-Form:source - see Grp 19)
 COVAL(IMIX ,F4 ,In-Form:source - see Grp 19)
 COVAL(IMIX ,F3 ,In-Form:source - see Grp 19)
 COVAL(IMIX ,F2 ,In-Form:source - see Grp 19)
 COVAL(IMIX ,F1 ,In-Form:source - see Grp 19)
 
 PATCH(PATCH1 ,CELL , 1, 1, 1, 15, 1, 1, 1, 1)
 
 PATCH(PATCH2 ,CELL , 1, 1, 1, 15, 2, 2, 1, 1)
 
 PATCH(PATCH3 ,CELL , 1, 1, 1, 15, 3, 3, 1, 1)
 
 PATCH(PATCH4 ,CELL , 1, 1, 1, 15, 4, 4, 1, 1)
 
 PATCH(PATCH5 ,CELL , 1, 1, 1, 15, 5, 5, 1, 1)
 
 PATCH(PATCH6 ,CELL , 1, 1, 1, 15, 6, 6, 1, 1)
 
 PATCH(PATCH7 ,CELL , 1, 1, 1, 15, 7, 7, 1, 1)
 
 PATCH(PATCH8 ,CELL , 1, 1, 1, 15, 8, 8, 1, 1)
 
 PATCH(PATCH9 ,CELL , 1, 1, 1, 15, 9, 9, 1, 1)
 
 PATCH(PATCH10 ,CELL , 1, 1, 1, 15, 10, 10, 1, 1)
 
 PATCH(PATCH11 ,CELL , 1, 1, 1, 15, 11, 11, 1, 1)
 
 PATCH(PATCH12 ,CELL , 1, 1, 1, 15, 12, 12, 1, 1)
 
 PATCH(PATCH13 ,CELL , 1, 1, 1, 15, 13, 13, 1, 1)
 
 PATCH(PATCH14 ,CELL , 1, 1, 1, 15, 14, 14, 1, 1)
 
 PATCH(PATCH15 ,CELL , 1, 1, 1, 15, 15, 15, 1, 1)
 
 PATCH(PATCH16 ,CELL , 1, 1, 1, 15, 16, 16, 1, 1)
 
 PATCH(PATCH17 ,CELL , 1, 1, 1, 15, 17, 17, 1, 1)
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 250 ;ISWC1 = 1
 LITHYD = 10 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-02
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(V1)=10
 LITER(W1)=10 ;LITER(KE)=20
 LITER(EP)=20 ;LITER(H1)=20
 LITER(F17)=20 ;LITER(F16)=20
 LITER(F15)=20 ;LITER(F14)=20
 LITER(F13)=20 ;LITER(F12)=20
 LITER(F11)=20 ;LITER(F10)=20
 LITER(F9)=20 ;LITER(F8)=20
 LITER(F7)=20 ;LITER(F6)=20
 LITER(F5)=20 ;LITER(F4)=20
 LITER(F3)=20 ;LITER(F2)=20
 LITER(F1)=20 ;LITER(G)=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(F17)=1.0E-03 ;ENDIT(F16)=1.0E-03
 ENDIT(F15)=1.0E-03 ;ENDIT(F14)=1.0E-03
 ENDIT(F13)=1.0E-03 ;ENDIT(F12)=1.0E-03
 ENDIT(F11)=1.0E-03 ;ENDIT(F10)=1.0E-03
 ENDIT(F9)=1.0E-03 ;ENDIT(F8)=1.0E-03
 ENDIT(F7)=1.0E-03 ;ENDIT(F6)=1.0E-03
 ENDIT(F5)=1.0E-03 ;ENDIT(F4)=1.0E-03
 ENDIT(F3)=1.0E-03 ;ENDIT(F2)=1.0E-03
 ENDIT(F1)=1.0E-03 ;ENDIT(G)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,0.25)
 RELAX(V1,FALSDT,0.025)
 RELAX(W1,FALSDT,0.025)
 RELAX(KE,FALSDT,0.025)
 RELAX(EP,FALSDT,0.025)
 RELAX(H1,FALSDT,1.0E+09)
 RELAX(FPD,LINRLX,1.)
 RELAX(GF,LINRLX,1.)
 RELAX(GAV,LINRLX,1.)
 RELAX(MAS,LINRLX,1.)
 RELAX(CAV,LINRLX,1.)
 RELAX(F17,LINRLX,0.15)
 RELAX(F16,LINRLX,0.15)
 RELAX(F15,LINRLX,0.15)
 RELAX(F14,LINRLX,0.15)
 RELAX(F13,LINRLX,0.15)
 RELAX(F12,LINRLX,0.15)
 RELAX(F11,LINRLX,0.15)
 RELAX(F10,LINRLX,0.15)
 RELAX(F9,LINRLX,0.15)
 RELAX(F8,LINRLX,0.15)
 RELAX(F7,LINRLX,0.15)
 RELAX(F6,LINRLX,0.15)
 RELAX(F5,LINRLX,0.15)
 RELAX(F4,LINRLX,0.15)
 RELAX(F3,LINRLX,0.15)
 RELAX(F2,LINRLX,0.15)
 RELAX(F1,LINRLX,0.15)
 RELAX(GGF,LINRLX,1.)
 RELAX(GG,LINRLX,1.)
 RELAX(DFYN,LINRLX,1.)
 RELAX(DFZH,LINRLX,1.)
 RELAX(DFY,LINRLX,1.)
 RELAX(DFZ,LINRLX,1.)
 RELAX(GENG,LINRLX,1.)
 RELAX(EPKE,LINRLX,1.)
 RELAX(GEN1,LINRLX,1.)
 RELAX(LEN1,LINRLX,1.)
 RELAX(ENUT,LINRLX,1.)
 RELAX(G,FALSDT,0.025)
 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(FPD)=1.0E+10 ;VARMIN(FPD)=-1.0E+10
 VARMAX(GF)=1.0E+10 ;VARMIN(GF)=-1.0E+10
 VARMAX(GAV)=1.0E+10 ;VARMIN(GAV)=-1.0E+10
 VARMAX(MAS)=1.0E+10 ;VARMIN(MAS)=-1.0E+10
 VARMAX(CAV)=1.0E+10 ;VARMIN(CAV)=-1.0E+10
 VARMAX(F17)=1. ;VARMIN(F17)=0.
 VARMAX(F16)=1. ;VARMIN(F16)=0.
 VARMAX(F15)=1. ;VARMIN(F15)=0.
 VARMAX(F14)=1. ;VARMIN(F14)=0.
 VARMAX(F13)=1. ;VARMIN(F13)=0.
 VARMAX(F12)=1. ;VARMIN(F12)=0.
 VARMAX(F11)=1. ;VARMIN(F11)=0.
 VARMAX(F10)=1. ;VARMIN(F10)=0.
 VARMAX(F9)=1. ;VARMIN(F9)=0.
 VARMAX(F8)=1. ;VARMIN(F8)=0.
 VARMAX(F7)=1. ;VARMIN(F7)=0.
 VARMAX(F6)=1. ;VARMIN(F6)=0.
 VARMAX(F5)=1. ;VARMIN(F5)=0.
 VARMAX(F4)=1. ;VARMIN(F4)=0.
 VARMAX(F3)=1. ;VARMIN(F3)=0.
 VARMAX(F2)=1. ;VARMIN(F2)=0.
 VARMAX(F1)=1. ;VARMIN(F1)=0.
 VARMAX(GGF)=1.0E+10 ;VARMIN(GGF)=-1.0E+10
 VARMAX(GG)=1.0E+10 ;VARMIN(GG)=-1.0E+10
 VARMAX(DFYN)=1.0E+10 ;VARMIN(DFYN)=-1.0E+10
 VARMAX(DFZH)=1.0E+10 ;VARMIN(DFZH)=-1.0E+10
 VARMAX(DFY)=1.0E+10 ;VARMIN(DFY)=-1.0E+10
 VARMAX(DFZ)=1.0E+10 ;VARMIN(DFZ)=-1.0E+10
 VARMAX(GENG)=1.0E+10 ;VARMIN(GENG)=-1.0E+10
 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
 VARMAX(GEN1)=1.0E+10 ;VARMIN(GEN1)=-1.0E+10
 VARMAX(LEN1)=1.0E+10 ;VARMIN(LEN1)=-1.0E+10
 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
 VARMAX(G)=1.0E+10 ;VARMIN(G)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 GENK = T
 PARSOL = F
 ISG62 = 1
 SPEDAT(SET,SOURCE,G!ISORG,C,=2.0*1.*EPKE*(GENG/(2.0*1.*EPKE+TINY)$)
 SPEDAT(SET,SOURCE,G!ISORG,C,-G))
 SPEDAT(SET,SOURCE,G!IWG,C,=GENG/(2.0*1.*EPKE+TINY)!FIXVAL)
 SPEDAT(SET,MAKE,RES,C,=1.E-10)
 SPEDAT(SET,STORED,RES,C,=SUM(VOL*(GENG-2.*1.*EPKE*G)/(15*20))!ZSL$)
 SPEDAT(SET,STORED,RES,C,FIN)
 SPEDAT(SET,STORED,RESREF(G),C,=RES!ZSLFIN)
 SPEDAT(SET,STORED,DFZ!PAT1,C,=((H1[&&+1]-H1)/DZG)^2!ZSLFIN)
 SPEDAT(SET,STORED,DFY!PAT2,C,=((H1[&+1&]-H1)/DYG)^2!ZSLFIN)
 SPEDAT(SET,STORED,DFZH!PAT3,C,=((H1-H1[&&-1])/DZG[&&-1])^2!ZSLFIN)
 SPEDAT(SET,STORED,DFYN!PAT4,C,=((H1-H1[&-1&])/DYG[&-1&])^2!ZSLFIN)
 SPEDAT(SET,STORED,GENG,C,=2.8*1.*ENUT*(DFZ+DFY+DFZH+DFYN)!ZSLFIN)
 SPEDAT(SET,STORED,GG,C,=SQRT(G)!IF(ISWEEP.EQ.250)!ZSLFIN)
 SPEDAT(SET,STORED,GGF,C,=GG/(H1+1.0E-20)!IF(ISWEEP.EQ.250)!ZSLFIN)
 SPEDAT(SET,SOURCE,F1!IMIX,C,=5.*EPKE*(F3+F5+F7+F9+F11+F13+F15+F17$)
 SPEDAT(SET,SOURCE,F1!IMIX,C,)*(0.-F1)!LINE)
 SPEDAT(SET,SOURCE,F2!IMIX,C,=2.*5.*EPKE*(F1*F3)-5.*EPKE*(F4+F6+F8$)
 SPEDAT(SET,SOURCE,F2!IMIX,C,+F10+F12+F14+F16)*F2)
 SPEDAT(SET,SOURCE,F3!IMIX,C,=2.*5.*EPKE*(F2*F4+F1*F5)-5.*EPKE*(F1$)
 SPEDAT(SET,SOURCE,F3!IMIX,C,+F17+F5+F7+F9+F11+F13+F15)*F3)
 SPEDAT(SET,SOURCE,F4!IMIX,C,=2.*5.*EPKE*(F3*F5+F2*F6+F1*F7)-5.*EP$)
 SPEDAT(SET,SOURCE,F4!IMIX,C,KE*(F2+F6+F8+F10+F12+F14+F16)*F4)
 SPEDAT(SET,SOURCE,F5!IMIX,C,=2.*5.*EPKE*(F4*F6+F3*F7+F2*F8+F1*F9)$)
 SPEDAT(SET,SOURCE,F5!IMIX,C,-5.*EPKE*(F1+F3+F17+F7+F9+F11+F13+F15$)
 SPEDAT(SET,SOURCE,F5!IMIX,C,)*F5)
 SPEDAT(SET,SOURCE,F6!IMIX,C,=2.*5.*EPKE*(F5*F7+F4*F8+F3*F9+F2*F10$)
 SPEDAT(SET,SOURCE,F6!IMIX,C,+F1*F11)-5.*EPKE*(F2+F4+F8+F10+F12+F1$)
 SPEDAT(SET,SOURCE,F6!IMIX,C,4+F16)*F6)
 SPEDAT(SET,SOURCE,F7!IMIX,C,=2.*5.*EPKE*(F6*F8+F5*F9+F4*F10+F3*F1$)
 SPEDAT(SET,SOURCE,F7!IMIX,C,1+F2*F12+F1*F13)-5.*EPKE*(F1+F3+F5+F1$)
 SPEDAT(SET,SOURCE,F7!IMIX,C,7+F9+F11+F13+F15)*F7)
 SPEDAT(SET,SOURCE,F8!IMIX,C,=2.*5.*EPKE*(F7*F9+F6*F10+F5*F11+F4*F$)
 SPEDAT(SET,SOURCE,F8!IMIX,C,12+F3*F13+F2*F14+F1*F15)-5.*EPKE*(F2+$)
 SPEDAT(SET,SOURCE,F8!IMIX,C,F4+F6+F10+F12+F14+F16)*F8)
 SPEDAT(SET,SOURCE,F9!IMIX,C,=2.*5.*EPKE*(F8*F10+F7*F11+F6*F12+F5*$)
 SPEDAT(SET,SOURCE,F9!IMIX,C,F13+F4*F14+F3*F15+F2*F16+F1*F17)-5.*E$)
 SPEDAT(SET,SOURCE,F9!IMIX,C,PKE*(F1+F3+F5+F7+F17+F11+F13+F15)*F9)
 SPEDAT(SET,SOURCE,F10!IMIX,C,=2.*5.*EPKE*(F9*F11+F8*F12+F7*F13+F6$)
 SPEDAT(SET,SOURCE,F10!IMIX,C,*F14+F5*F15+F4*F16+F3*F17)-5.*EPKE*($)
 SPEDAT(SET,SOURCE,F10!IMIX,C,F2+F4+F6+F8+F12+F14+F16)*F10)
 SPEDAT(SET,SOURCE,F11!IMIX,C,=2.*5.*EPKE*(F10*F12+F9*F13+F8*F14+F$)
 SPEDAT(SET,SOURCE,F11!IMIX,C,7*F15+F6*F16+F5*F17)-5.*EPKE*(F1+F3+$)
 SPEDAT(SET,SOURCE,F11!IMIX,C,F5+F7+F9+F17+F13+F15)*F11)
 SPEDAT(SET,SOURCE,F12!IMIX,C,=2.*5.*EPKE*(F11*F13+F10*F14+F9*F15+$)
 SPEDAT(SET,SOURCE,F12!IMIX,C,F8*F16+F7*F17)-5.*EPKE*(F2+F4+F6+F8+$)
 SPEDAT(SET,SOURCE,F12!IMIX,C,F10+F14+F16)*F12)
 SPEDAT(SET,SOURCE,F13!IMIX,C,=2.*5.*EPKE*(F12*F14+F11*F15+F10*F16$)
 SPEDAT(SET,SOURCE,F13!IMIX,C,+F9*F17)-5.*EPKE*(F1+F3+F5+F7+F9+F11$)
 SPEDAT(SET,SOURCE,F13!IMIX,C,+F15+F17)*F13)
 SPEDAT(SET,SOURCE,F14!IMIX,C,=2.*5.*EPKE*(F13*F15+F12*F16+F11*F17$)
 SPEDAT(SET,SOURCE,F14!IMIX,C,)-5.*EPKE*(F2+F4+F6+F8+F10+F12+F16)*$)
 SPEDAT(SET,SOURCE,F14!IMIX,C,F14)
 SPEDAT(SET,SOURCE,F15!IMIX,C,=2.*5.*EPKE*(F14*F16+F13*F17)-5.*EPK$)
 SPEDAT(SET,SOURCE,F15!IMIX,C,E*(F1+F3+F5+F7+F9+F11+F13+F17)*F15)
 SPEDAT(SET,SOURCE,F16!IMIX,C,=2.*5.*EPKE*(F15*F17)-5.*EPKE*(F2+F4$)
 SPEDAT(SET,SOURCE,F16!IMIX,C,+F6+F8+F10+F12+F14)*F16)
 SPEDAT(SET,SOURCE,F17!IMIX,C,=5.*EPKE*(F1+F3+F5+F7+F9+F11+F13+F15$)
 SPEDAT(SET,SOURCE,F17!IMIX,C,)*(0.-F17)!LINE)
 SPEDAT(SET,STORED,CAV,C,=16./16.*F1+15./16.*F2+14./16.*F3+13./16.$)
 SPEDAT(SET,STORED,CAV,C,*F4+12./16.*F5+11./16.*F6+10./16.*F7+9./1$)
 SPEDAT(SET,STORED,CAV,C,6.*F8+8./16.*F9!IF(ISWEEP.EQ.250)!ZSLFIN)
 SPEDAT(SET,STORED,CAV,C,=CAV+7./16.*F10+6./16.*F11+5./16.*F12+4./$)
 SPEDAT(SET,STORED,CAV,C,16.*F13+3./16.*F14+2./16.*F15+1./16.*F16+$)
 SPEDAT(SET,STORED,CAV,C,0./16.*F17!IF(ISWEEP.EQ.250)!ZSLFIN)
 SPEDAT(SET,STORED,MAS,C,=F1+F2+F3+F4+F5+F6+F7+F8+F9+F10+F11+F12+F$)
 SPEDAT(SET,STORED,MAS,C,13+F14+F15+F16+F17!IF(ISWEEP.EQ.250)!ZSLF$)
 SPEDAT(SET,STORED,MAS,C,IN)
 SPEDAT(SET,STORED,GAV,C,=ABS(CAV-16./16)*F1+ABS(CAV-15./16.)*F2+A$)
 SPEDAT(SET,STORED,GAV,C,BS(CAV-14./16)*F3+ABS(CAV-13./16.)*F4+ABS$)
 SPEDAT(SET,STORED,GAV,C,(CAV-12./16)*F5+ABS(CAV-11./16.)*F6!ZSLFI$)
 SPEDAT(SET,STORED,GAV,C,N)
 SPEDAT(SET,STORED,GAV,C,=GAV+ABS(CAV-10./16)*F7+ABS(CAV-9./16.)*F$)
 SPEDAT(SET,STORED,GAV,C,8+ABS(CAV-8./16)*F9+ABS(CAV-7./16.)*F10+A$)
 SPEDAT(SET,STORED,GAV,C,BS(CAV-6./16)*F11+ABS(CAV-5./16.)*F12!ZSL$)
 SPEDAT(SET,STORED,GAV,C,FIN)
 SPEDAT(SET,STORED,GAV,C,=GAV+ABS(CAV-4./16)*F13+ABS(CAV-3./16.)*F$)
 SPEDAT(SET,STORED,GAV,C,14+ABS(CAV-2./16)*F15+ABS(CAV-1./16.)*F16$)
 SPEDAT(SET,STORED,GAV,C,+ABS(CAV-0./16)*F17!ZSLFIN)
 SPEDAT(SET,STORED,GF,C,=GAV/(CAV+1.0E-20)!IF(ISWEEP.EQ.250)!ZSLFIN)
 SPEDAT(SET,STORED,FPD!PATCH1,C,=F1[1&4&4]*MAX(ABS(F1[1&4&4]-YV)/($)
 SPEDAT(SET,STORED,FPD!PATCH1,C,F1[1&4&4]-YV)&0.)!IF(ISWEEP.EQ.250$)
 SPEDAT(SET,STORED,FPD!PATCH1,C,)!ZSLFIN)
 SPEDAT(SET,STORED,FPD!PATCH2,C,=F2[1&4&4]*MAX(ABS(F2[1&4&4]-YV)/($)
 SPEDAT(SET,STORED,FPD!PATCH2,C,F2[1&4&4]-YV)&0.0)!IF(ISWEEP.EQ.25$)
 SPEDAT(SET,STORED,FPD!PATCH2,C,0)!ZSLFIN)
 SPEDAT(SET,STORED,FPD!PATCH3,C,=F3[1&4&4]*MAX(ABS(F3[1&4&4]-YV)/($)
 SPEDAT(SET,STORED,FPD!PATCH3,C,F3[1&4&4]-YV)&0.0)!IF(ISWEEP.EQ.25$)
 SPEDAT(SET,STORED,FPD!PATCH3,C,0)!ZSLFIN)
 SPEDAT(SET,STORED,FPD!PATCH4,C,=F4[1&4&4]*MAX(ABS(F4[1&4&4]-YV)/($)
 SPEDAT(SET,STORED,FPD!PATCH4,C,F4[1&4&4]-YV)&0.0)!IF(ISWEEP.EQ.25$)
 SPEDAT(SET,STORED,FPD!PATCH4,C,0)!ZSLFIN)
 SPEDAT(SET,STORED,FPD!PATCH5,C,=F5[1&4&4]*MAX(ABS(F5[1&4&4]-YV)/($)
 SPEDAT(SET,STORED,FPD!PATCH5,C,F5[1&4&4]-YV)&0.0)!IF(ISWEEP.EQ.25$)
 SPEDAT(SET,STORED,FPD!PATCH5,C,0)!ZSLFIN)
 SPEDAT(SET,STORED,FPD!PATCH6,C,=F6[1&4&4]*MAX(ABS(F6[1&4&4]-YV)/($)
 SPEDAT(SET,STORED,FPD!PATCH6,C,F6[1&4&4]-YV)&0.0)!IF(ISWEEP.EQ.25$)
 SPEDAT(SET,STORED,FPD!PATCH6,C,0)!ZSLFIN)
 SPEDAT(SET,STORED,FPD!PATCH7,C,=F7[1&4&4]*MAX(ABS(F7[1&4&4]-YV)/($)
 SPEDAT(SET,STORED,FPD!PATCH7,C,F7[1&4&4]-YV)&0.0)!IF(ISWEEP.EQ.25$)
 SPEDAT(SET,STORED,FPD!PATCH7,C,0)!ZSLFIN)
 SPEDAT(SET,STORED,FPD!PATCH8,C,=F8[1&4&4]*MAX(ABS(F8[1&4&4]-YV)/($)
 SPEDAT(SET,STORED,FPD!PATCH8,C,F8[1&4&4]-YV)&0.0)!IF(ISWEEP.EQ.25$)
 SPEDAT(SET,STORED,FPD!PATCH8,C,0)!ZSLFIN)
 SPEDAT(SET,STORED,FPD!PATCH9,C,=F9[1&4&4]*MAX(ABS(F9[1&4&4]-YV)/($)
 SPEDAT(SET,STORED,FPD!PATCH9,C,F9[1&4&4]-YV)&0.0)!IF(ISWEEP.EQ.25$)
 SPEDAT(SET,STORED,FPD!PATCH9,C,0)!ZSLFIN)
 SPEDAT(SET,STORED,FPD!PATCH10,C,=F10[1&4&4]*MAX(ABS(F10[1&4&4]-YV$)
 SPEDAT(SET,STORED,FPD!PATCH10,C,)/(F10[1&4&4]-YV)&0.0)!IF(ISWEEP.$)
 SPEDAT(SET,STORED,FPD!PATCH10,C,EQ.250)!ZSLFIN)
 SPEDAT(SET,STORED,FPD!PATCH11,C,=F11[1&4&4]*MAX(ABS(F11[1&4&4]-YV$)
 SPEDAT(SET,STORED,FPD!PATCH11,C,)/(F11[1&4&4]-YV)&0.0)!IF(ISWEEP.$)
 SPEDAT(SET,STORED,FPD!PATCH11,C,EQ.250)!ZSLFIN)
 SPEDAT(SET,STORED,FPD!PATCH12,C,=F12[1&4&4]*MAX(ABS(F12[1&4&4]-YV$)
 SPEDAT(SET,STORED,FPD!PATCH12,C,)/(F12[1&4&4]-YV)&0.0)!IF(ISWEEP.$)
 SPEDAT(SET,STORED,FPD!PATCH12,C,EQ.250)!ZSLFIN)
 SPEDAT(SET,STORED,FPD!PATCH13,C,=F13[1&4&4]*MAX(ABS(F13[1&4&4]-YV$)
 SPEDAT(SET,STORED,FPD!PATCH13,C,)/(F13[1&4&4]-YV)&0.0)!IF(ISWEEP.$)
 SPEDAT(SET,STORED,FPD!PATCH13,C,EQ.250)!ZSLFIN)
 SPEDAT(SET,STORED,FPD!PATCH14,C,=F14[1&4&4]*MAX(ABS(F14[1&4&4]-YV$)
 SPEDAT(SET,STORED,FPD!PATCH14,C,)/(F14[1&4&4]-YV)&0.0)!IF(ISWEEP.$)
 SPEDAT(SET,STORED,FPD!PATCH14,C,EQ.250)!ZSLFIN)
 SPEDAT(SET,STORED,FPD!PATCH15,C,=F15[1&4&4]*MAX(ABS(F15[1&4&4]-YV$)
 SPEDAT(SET,STORED,FPD!PATCH15,C,)/(F15[1&4&4]-YV)&0.0)!IF(ISWEEP.$)
 SPEDAT(SET,STORED,FPD!PATCH15,C,EQ.250)!ZSLFIN)
 SPEDAT(SET,STORED,FPD!PATCH16,C,=F16[1&4&4]*MAX(ABS(F16[1&4&4]-YV$)
 SPEDAT(SET,STORED,FPD!PATCH16,C,)/(F16[1&4&4]-YV)&0.0)!IF(ISWEEP.$)
 SPEDAT(SET,STORED,FPD!PATCH16,C,EQ.250)!ZSLFIN)
 SPEDAT(SET,STORED,FPD!PATCH17,C,=F17[1&4&4]*MAX(ABS(F17[1&4&4]-YV$)
 SPEDAT(SET,STORED,FPD!PATCH17,C,)/(F17[1&4&4]-YV)&0.0)!IF(ISWEEP.$)
 SPEDAT(SET,STORED,FPD!PATCH17,C,EQ.250)!ZSLFIN)
 IG( 1) = 4
 IG( 2) = 4
 RG( 1) =5.
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=0.5613 ;EX(V1)=0.02532
 EX(W1)=2.84 ;EX(KE)=0.4442
 EX(EP)=6.611 ;EX(H1)=0.1939
 EX(FPD)=0.02724 ;EX(GF)=0.4936
 EX(GAV)=0.05976 ;EX(MAS)=1.
 EX(CAV)=0.1939 ;EX(F17)=0.1402
 EX(F16)=0.07228 ;EX(F15)=0.2095
 EX(F14)=0.3191 ;EX(F13)=0.1075
 EX(F12)=0.04737 ;EX(F11)=0.0269
 EX(F10)=0.01746 ;EX(F9)=0.01188
 EX(F8)=8.836E-03 ;EX(F7)=6.076E-03
 EX(F6)=5.269E-03 ;EX(F5)=3.765E-03
 EX(F4)=3.483E-03 ;EX(F3)=2.355E-03
 EX(F2)=3.138E-03 ;EX(F1)=0.01499
 EX(GGF)=2.83E+13 ;EX(GG)=0.05768
 EX(DFYN)=0.3957 ;EX(DFZH)=3.238E-06
 EX(DFY)=17.65 ;EX(DFZ)=0.04226
 EX(GENG)=0.2216 ;EX(EPKE)=8.359
 EX(GEN1)=1724. ;EX(LEN1)=0.01567
 EX(ENUT)=4.963E-03 ;EX(G)=5.591E-03
 ************************************************************
  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,Y,Y,Y,Y,Y)
 OUTPUT(EP,Y,N,Y,Y,Y,Y)
 OUTPUT(H1,Y,Y,Y,Y,Y,Y)
 OUTPUT(FPD,Y,N,Y,N,N,N)
 OUTPUT(GF,Y,N,Y,N,N,N)
 OUTPUT(GAV,Y,N,Y,N,N,N)
 OUTPUT(MAS,Y,N,Y,N,N,N)
 OUTPUT(CAV,Y,N,Y,N,N,N)
 OUTPUT(F17,Y,N,Y,Y,Y,Y)
 OUTPUT(F16,Y,N,Y,Y,Y,Y)
 OUTPUT(F15,Y,N,Y,Y,Y,Y)
 OUTPUT(F14,Y,N,Y,Y,Y,Y)
 OUTPUT(F13,Y,N,Y,Y,Y,Y)
 OUTPUT(F12,Y,N,Y,Y,Y,Y)
 OUTPUT(F11,Y,N,Y,Y,Y,Y)
 OUTPUT(F10,Y,N,Y,Y,Y,Y)
 OUTPUT(F9,Y,N,Y,Y,Y,Y)
 OUTPUT(F8,Y,N,Y,Y,Y,Y)
 OUTPUT(F7,Y,N,Y,Y,Y,Y)
 OUTPUT(F6,Y,N,Y,Y,Y,Y)
 OUTPUT(F5,Y,N,Y,Y,Y,Y)
 OUTPUT(F4,Y,N,Y,Y,Y,Y)
 OUTPUT(F3,Y,N,Y,Y,Y,Y)
 OUTPUT(F2,Y,N,Y,Y,Y,Y)
 OUTPUT(F1,Y,N,Y,Y,Y,Y)
 OUTPUT(GGF,Y,N,Y,N,N,N)
 OUTPUT(GG,Y,N,Y,N,N,N)
 OUTPUT(DFYN,Y,N,Y,N,N,N)
 OUTPUT(DFZH,Y,N,Y,N,N,N)
 OUTPUT(DFY,Y,N,Y,N,N,N)
 OUTPUT(DFZ,Y,N,Y,N,N,N)
 OUTPUT(GENG,Y,N,Y,N,N,N)
 OUTPUT(EPKE,Y,N,Y,N,N,N)
 OUTPUT(GEN1,Y,N,Y,N,N,N)
 OUTPUT(LEN1,Y,N,Y,N,N,N)
 OUTPUT(ENUT,Y,N,Y,N,N,N)
 OUTPUT(G,Y,N,Y,Y,Y,Y)
 WALPRN = T
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 1 ;IYMON = 14 ;IZMON = 19
 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NYPRIN = 1 ;IYPRF = 1 ;IYPRL = 30
 NZPRIN = 1 ;IZPRF = 1 ;IZPRL = 10000
 XZPR = F ;YZPR = F
 IPLTF = 1 ;IPLTL = -1 ;NPLT = 1
 ISWPRF = 1 ;ISWPRL = 100000
 ITABL = 3 ;IPROF = 1
 ABSIZ =0.5 ;ORSIZ =0.2
 NTZPRF = 1 ;NCOLPF = 50
 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
 
 PATCH(PROF1 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
 PLOT(PROF1 ,F1 ,0. ,0. )
 
 PATCH(PROF2 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
 PLOT(PROF2 ,F2 ,0. ,0. )
 
 PATCH(PROF3 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
 PLOT(PROF3 ,F3 ,0. ,0. )
 
 PATCH(PROF4 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
 PLOT(PROF4 ,F4 ,0. ,0. )
 
 PATCH(PROF5 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
 PLOT(PROF5 ,F5 ,0. ,0. )
 
 PATCH(PROF6 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
 PLOT(PROF6 ,F6 ,0. ,0. )
 
 PATCH(PROF7 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
 PLOT(PROF7 ,F7 ,0. ,0. )
 
 PATCH(PROF8 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
 PLOT(PROF8 ,F8 ,0. ,0. )
 
 PATCH(PROF9 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
 PLOT(PROF9 ,F9 ,0. ,0. )
 
 PATCH(PROF10 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
 PLOT(PROF10 ,F10 ,0. ,0. )
 
 PATCH(PROF11 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
 PLOT(PROF11 ,F11 ,0. ,0. )
 
 PATCH(PROF12 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
 PLOT(PROF12 ,F12 ,0. ,0. )
 
 PATCH(PROF13 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
 PLOT(PROF13 ,F13 ,0. ,0. )
 
 PATCH(PROF14 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
 PLOT(PROF14 ,F14 ,0. ,0. )
 
 PATCH(PROF15 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
 PLOT(PROF15 ,F15 ,0. ,0. )
 
 PATCH(PROF16 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
 PLOT(PROF16 ,F16 ,0. ,0. )
 
 PATCH(PROF17 ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 1)
 PLOT(PROF17 ,F17 ,0. ,0. )
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP