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

  up z
  msg velocity vectors
  * set vec ref;0.5
  vec x 1 sh;  gr ou x 1
  msg Press  to continue
  pause;  cl
  msg contour of pressure
  con p1 x 1 fi;0.002;  gr ou x 1
  msg Press  to continue
  pause;  cl
  msg contour of SLTM, deduced from MFM
  con SLTM x 1 fi;0.002;  gr ou x 1
  pause;  cl
  msg contour of SLTS, dedued from single-fluid model
  con SLTS x 1 fi;0.002;  gr ou x 1
  msg Press  to continue
  pause;  cl
  msg contour of c1
  con c1 x 1 fi;0.002;  gr ou x 1
  msg Press  to continue
  pause;  cl
  msg contour of enut
  con enut x 1 fi;0.002;  gr ou x 1
  msg Press  to continue
  pause;  cl
  msg contour of ep/ke
  con RATE x 1 fi;0.002;  gr ou x 1
  msg Press  to continue
  pause;  cl
  msg contour of ke
  con ke x 1 fi;0.002;  gr ou x 1
  pause;  cl
  msg contour of mnsq
  con mnsq x 1 fi;0.002;  gr ou x 1
  pause;  cl
  msg contour of avef
  con avef x 1 fi;0.002;  gr ou x 1
  pause;  cl
  msg contour of fluid 1 volume fraction
  con f1 x 1 fi;0.002;  gr ou x 1
  msg Press  to continue
  pause;  cl
  msg contour of fluid 2 volume fraction
  con f2 x 1 fi;0.002;  gr ou x 1
  msg Press  to continue
  pause;  cl
  msg contour of fluid 3 volume fraction
  con f3 x 1 fi;0.002;  gr ou x 1
  msg Press  to continue
  pause;  cl
  msg contour of fluid 4 volume fraction
  con f4 x 1 fi;0.002;  gr ou x 1
  msg Press  to continue
  pause;  cl
  msg contour of fluid 5 volume fraction
  con f5 x 1 fi;0.002;  gr ou x 1
  msg Press  to continue
  pause;  cl
  msg contour of fluid 6 volume fraction
  con f6 x 1 fi;0.002;  gr ou x 1
  msg Press  to continue
  pause;  cl
  msg contour of fluid 7 volume fraction
  con f7 x 1 fi;0.002;  gr ou x 1
  msg Press  to continue
  pause;  cl
  msg contour of fluid 8 volume fraction
  con f8 x 1 fi;0.002;  gr ou x 1
  msg Press  to continue
  pause;  cl
  msg contour of fluid 9 volume fraction
  con f9 x 1 fi;0.002;  gr ou x 1
  msg Press  to continue
  pause;  cl
  msg contour of fluid 10 volume fraction
  con f10 x 1 fi;0.002;  gr ou x 1
  pause;  cl
  msg contour of fluid 11 volume fraction
  con f11 x 1 fi;0.002;  gr ou x 1
  msg Press  to continue
  enduse
  DISPLAY
  The geometry of this case is similar to that of case W918; but,
  unlike that case, ONEPHS=T. There are therefore only 3 velocity
  components. On the other hand, 11 fluids are distinguished, formed
  by mixing an acidic fluid with an alkaline one in different
  proportions.

                                             __________|.|________
  The sketch illustrates the apparatus and  |          |||        |
  the initial state of the two liquids.     | upper    |.|        |
                                            | liquid   |||        |
  They are both at rest, and are in contact |          |.|  acid  |
  at a horizontal interface                 |..........|||........|
                                            | lower    |.| alkali |
  The paddle is supposed to be suddenly set | liquid   |||        |
  in motion.                                |       ---------     |
                                            |paddle ////\\\\\     |
  The computational task is to predict both |           .         |
  the macro-mixing, represented by the      ------------|---------
  subsequent distributions of velocity,                 .< axis of
  pressure and time-average concentration,                 symmetry
  but also the extent to which the two       The stirred mixing tank
  liquids are mixed together at any point.

  ENDDIS
  DISPLAY
  Notes:

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

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

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

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

  (4) the values are:
      mfm  = $L004
      mfm1 = $L003
      mfm2 = $L002
      mfm3 = $L001
  ENDDIS
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(MFM; Stirred Tank; Conmix=50.           )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 500
 ************************************************************
  Group 2. Time dependence
 STEADY = F
    * Set overall time and no. of steps
 TFIRST =0. ;TLAST =10.
 FSTEP = 1 ;LSTEP = 10
 TFRAC(1)=0.1 ;TFRAC(2)=0.2
 TFRAC(3)=0.3 ;TFRAC(4)=0.4
 TFRAC(5)=0.5 ;TFRAC(6)=0.6
 TFRAC(7)=0.7 ;TFRAC(8)=0.8
 TFRAC(9)=0.9 ;TFRAC(10)=1.
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = F
 NX = 1
 XULAST =0.1
 XFRAC(1)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 10
 YVLAST =0.5
 YFRAC(1)=0.1 ;YFRAC(2)=0.2
 YFRAC(3)=0.3 ;YFRAC(4)=0.4
 YFRAC(5)=0.5 ;YFRAC(6)=0.6
 YFRAC(7)=0.7 ;YFRAC(8)=0.8
 YFRAC(9)=0.9 ;YFRAC(10)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 10
 ZWLAST =1.
 ZFRAC(1)=0.1 ;ZFRAC(2)=0.2
 ZFRAC(3)=0.3 ;ZFRAC(4)=0.4
 ZFRAC(5)=0.5 ;ZFRAC(6)=0.6
 ZFRAC(7)=0.7 ;ZFRAC(8)=0.8
 ZFRAC(9)=0.9 ;ZFRAC(10)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(1)=P1 ;NAME(3)=U1
 NAME(5)=V1 ;NAME(7)=W1
 NAME(12)=KE ;NAME(13)=EP
 NAME(16)=C1 ;NAME(131)=F1
 NAME(132)=F2 ;NAME(133)=F3
 NAME(134)=F4 ;NAME(135)=F5
 NAME(136)=F6 ;NAME(137)=F7
 NAME(138)=F8 ;NAME(139)=F9
 NAME(140)=F10 ;NAME(141)=F11
 NAME(142)=MNSQ ;NAME(143)=AVEF
 NAME(144)=SLTS ;NAME(145)=SLTM
 NAME(146)=MIXL ;NAME(147)=EPKE
 NAME(148)=RATE ;NAME(149)=LEN1
 NAME(150)=ENUT
    * 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,Y)
 SOLUTN(U1,Y,Y,N,N,N,Y)
 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(C1,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(SLTS,Y,Y,N,N,N,Y)
 SOLUTN(SLTM,Y,Y,N,N,N,Y)
 SOLUTN(MIXL,Y,N,N,N,N,N)
 SOLUTN(EPKE,Y,N,N,N,N,N)
 SOLUTN(RATE,Y,N,N,N,N,Y)
 SOLUTN(LEN1,Y,N,N,N,N,Y)
 SOLUTN(ENUT,Y,N,N,N,N,Y)
 VIST = 150
 LEN1 = 149
 ************************************************************
  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(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(C1,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(SLTS,N,Y,Y,Y,Y,Y)
 TERMS(SLTM,N,Y,Y,Y,N,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 =1000. ;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 =1.0E-02 ;ENUT = GRND3
 ENUTA =0. ;ENUTB =0.
 ENUTC =0.
 IENUTA = 0
 PRNDTL(U1)=1. ;PRNDTL(V1)=1.
 PRNDTL(W1)=1. ;PRNDTL(KE)=1.
 PRNDTL(EP)=1. ;PRNDTL(C1)=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(SLTS)=1.
 PRNDTL(SLTM)=1.
 PRT(U1)=1. ;PRT(V1)=1.
 PRT(W1)=1. ;PRT(KE)=1.
 PRT(EP)=1.314 ;PRT(C1)=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(SLTS)=1.
 PRT(SLTM)=1.
 CP1 =1. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1.0E-10 ;FIINIT(U1)=1.0E-10
 FIINIT(V1)=1.0E-10 ;FIINIT(W1)=1.0E-10
 FIINIT(KE)=1. ;FIINIT(EP)=0.1
 FIINIT(C1)=0. ;FIINIT(F1)=0.
 FIINIT(F2)=0. ;FIINIT(F3)=0.
 FIINIT(F4)=0. ;FIINIT(F5)=0.
 FIINIT(F6)=0. ;FIINIT(F7)=0.
 FIINIT(F8)=0. ;FIINIT(F9)=0.
 FIINIT(F10)=0. ;FIINIT(F11)=1.
 FIINIT(MNSQ)=1.0E-10 ;FIINIT(AVEF)=1.0E-10
 FIINIT(SLTS)=1.0E-10 ;FIINIT(SLTM)=0.
 FIINIT(MIXL)=1.0E-10 ;FIINIT(EPKE)=1.0E-10
 FIINIT(RATE)=0.1 ;FIINIT(LEN1)=0.2
 FIINIT(ENUT)=0.1
 
 PATCH(START ,INIVAL, 1, 1, 1, 10, 1, 5, 1, 10)
 INIT(START ,C1 ,0. ,1. )
 INIT(START ,F1 ,0. ,1. )
 INIT(START ,F2 ,0. ,0. )
 INIT(START ,F3 ,0. ,0. )
 INIT(START ,F4 ,0. ,0. )
 INIT(START ,F5 ,0. ,0. )
 INIT(START ,F6 ,0. ,0. )
 INIT(START ,F7 ,0. ,0. )
 INIT(START ,F8 ,0. ,0. )
 INIT(START ,F9 ,0. ,0. )
 INIT(START ,F10 ,0. ,0. )
 INIT(START ,F11 ,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, 10)
 COVAL(KESOURCE,KE , GRND4 , GRND4 )
 COVAL(KESOURCE,EP , GRND4 , GRND4 )
 
 PATCH(CYLWALL ,NWALL , 1, 1, 10, 10, 1, 10, 1, 10)
 COVAL(CYLWALL ,U1 , GRND2 ,0. )
 COVAL(CYLWALL ,W1 , GRND2 ,0. )
 COVAL(CYLWALL ,KE , GRND2 , GRND2 )
 COVAL(CYLWALL ,EP , GRND2 , GRND2 )
 
 PATCH(TOP ,HWALL , 1, 1, 1, 10, 10, 10, 1, 10)
 COVAL(TOP ,U1 , GRND2 ,0. )
 COVAL(TOP ,V1 , GRND2 ,0. )
 COVAL(TOP ,KE , GRND2 , GRND2 )
 COVAL(TOP ,EP , GRND2 , GRND2 )
 
 PATCH(BOTTOM ,LWALL , 1, 1, 1, 10, 1, 1, 1, 10)
 COVAL(BOTTOM ,U1 , GRND2 ,0. )
 COVAL(BOTTOM ,V1 , GRND2 ,0. )
 COVAL(BOTTOM ,KE , GRND2 , GRND2 )
 COVAL(BOTTOM ,EP , GRND2 , GRND2 )
 
 PATCH(PRELIEF ,HIGH , 1, 1, 10, 10, 10, 10, 1, 100)
 COVAL(PRELIEF ,P1 ,1. ,0. )
 COVAL(PRELIEF ,C1 ,0. ,0. )
 COVAL(PRELIEF ,F1 ,0. , SAME )
 COVAL(PRELIEF ,F2 ,0. , SAME )
 COVAL(PRELIEF ,F3 ,0. , SAME )
 COVAL(PRELIEF ,F4 ,0. , SAME )
 COVAL(PRELIEF ,F5 ,0. , SAME )
 COVAL(PRELIEF ,F6 ,0. , SAME )
 COVAL(PRELIEF ,F7 ,0. , SAME )
 COVAL(PRELIEF ,F8 ,0. , SAME )
 COVAL(PRELIEF ,F9 ,0. , SAME )
 COVAL(PRELIEF ,F10 ,0. , SAME )
 COVAL(PRELIEF ,F11 ,0. , SAME )
 
 PATCH(PADDLE ,CELL , 1, 1, 5, 5, 2, 2, 1, 10)
 COVAL(PADDLE ,U1 , FIXVAL ,10. )
 
 PATCH(SALT1 ,PHASEM, 1, 1, 1, 10, 1, 10, 1, 10)
 COVAL(SALT1 ,SLTM, FIXFLU , GRND )
 
 PATCH(SALT2 ,PHASEM, 1, 1, 1, 10, 1, 10, 1, 10)
 COVAL(SALT2 ,SLTS, FIXFLU , GRND )
 IURVAL = -1
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 100 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SARAH =1.0E-03
 SELREF = T
 RESFAC =1.0E-03
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(U1)=10
 LITER(V1)=10 ;LITER(W1)=10
 LITER(KE)=20 ;LITER(EP)=20
 LITER(C1)=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(SLTS)=20 ;LITER(SLTM)=20
 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
 ENDIT(V1)=1.0E-03 ;ENDIT(W1)=1.0E-03
 ENDIT(KE)=1.0E-03 ;ENDIT(EP)=1.0E-03
 ENDIT(C1)=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(SLTS)=1.0E-03 ;ENDIT(SLTM)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(U1,LINRLX,0.1)
 RELAX(V1,FALSDT,9.999999E-04)
 RELAX(W1,FALSDT,9.999999E-04)
 RELAX(KE,LINRLX,0.5)
 RELAX(EP,LINRLX,0.5)
 RELAX(C1,FALSDT,1.0E+09)
 RELAX(F1,LINRLX,0.5)
 RELAX(F2,LINRLX,0.5)
 RELAX(F3,LINRLX,0.5)
 RELAX(F4,LINRLX,0.5)
 RELAX(F5,LINRLX,0.5)
 RELAX(F6,LINRLX,0.5)
 RELAX(F7,LINRLX,0.5)
 RELAX(F8,LINRLX,0.5)
 RELAX(F9,LINRLX,0.5)
 RELAX(F10,LINRLX,0.5)
 RELAX(F11,LINRLX,0.5)
 RELAX(MNSQ,LINRLX,1.)
 RELAX(AVEF,LINRLX,1.)
 RELAX(SLTS,FALSDT,1.0E+09)
 RELAX(SLTM,FALSDT,1.0E+09)
 RELAX(MIXL,LINRLX,1.)
 RELAX(EPKE,LINRLX,1.)
 RELAX(RATE,LINRLX,1.)
 RELAX(LEN1,LINRLX,1.)
 RELAX(ENUT,LINRLX,0.1)
 KELIN = 3
 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(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-04
 VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
 VARMAX(C1)=1.0E+10 ;VARMIN(C1)=-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(SLTS)=1.0E+10 ;VARMIN(SLTS)=-1.0E+10
 VARMAX(SLTM)=1.0E+10 ;VARMIN(SLTM)=-1.0E+10
 VARMAX(MIXL)=1.0E+10 ;VARMIN(MIXL)=-1.0E+10
 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
 VARMAX(RATE)=1.0E+04 ;VARMIN(RATE)=-1.0E+10
 VARMAX(LEN1)=1.0E+10 ;VARMIN(LEN1)=-1.0E+10
 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 GENK = T
 PARSOL = F
 ISG62 = 1
 SPEDAT(SET,STORED,MIXL,C,=LEN1!ZSLFIN)
 SPEDAT(SET,STORED,RATE,C,=EPKE!ZSLFIN)
 SPEDAT(SET,MFM,MFMMOD,C,MFM)
 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,50.)
 SPEDAT(SET,MFM,VISCON,R,0.)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=74.32 ;EX(U1)=0.
 EX(V1)=0.09456 ;EX(W1)=0.1466
 EX(KE)=0.05478 ;EX(EP)=0.2388
 EX(C1)=0.4929 ;EX(F1)=8.365E-13
 EX(F2)=2.605E-12 ;EX(F3)=8.119E-12
 EX(F4)=5.094E-08 ;EX(F5)=0.1629
 EX(F6)=0.4913 ;EX(F7)=0.2956
 EX(F8)=0.0499 ;EX(F9)=2.687E-05
 EX(F10)=1.927E-08 ;EX(F11)=2.752E-10
 EX(MNSQ)=0.04651 ;EX(AVEF)=0.5143
 EX(SLTS)=2.177E-06 ;EX(SLTM)=2.107E-06
 EX(MIXL)=0.01619 ;EX(EPKE)=2.002
 EX(RATE)=2.002 ;EX(LEN1)=0.01619
 EX(ENUT)=1.641E-03
 ************************************************************
  Group 21. Print-out of Variables
 IURPRN = -1
 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,Y,Y,Y,Y,Y)
 OUTPUT(U1,N,N,Y,Y,N,N)
 OUTPUT(V1,Y,Y,Y,Y,Y,Y)
 OUTPUT(W1,Y,Y,Y,Y,Y,Y)
 OUTPUT(KE,Y,N,Y,Y,Y,Y)
 OUTPUT(EP,Y,N,Y,Y,Y,Y)
 OUTPUT(C1,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(SLTS,Y,N,Y,Y,Y,Y)
 OUTPUT(SLTM,Y,N,Y,Y,Y,Y)
 OUTPUT(MIXL,Y,N,Y,N,N,N)
 OUTPUT(EPKE,Y,N,Y,N,N,N)
 OUTPUT(RATE,Y,N,Y,N,N,N)
 OUTPUT(LEN1,Y,N,Y,N,N,N)
 OUTPUT(ENUT,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 1 ;IYMON = 5 ;IZMON = 5
 NPRMON = 20 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NTPRIN = 100000 ;ISTPRF = 1 ;ISTPRL = 100000
 NYPRIN = 1 ;IYPRF = 1 ;IYPRL = 10000
 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(MAP ,CONTUR, 1, 1, 1, 10, 1, 10, 1, 10)
 PLOT(MAP ,U1 ,0. ,10. )
 PLOT(MAP ,V1 ,0. ,10. )
 PLOT(MAP ,W1 ,0. ,10. )
 PLOT(MAP ,C1 ,0. ,10. )
 
 PATCH(IZEQNZM1,PROFIL, 1, 1, 1, 10, 9, 9, 1, 10)
 PLOT(IZEQNZM1,W1 ,0. ,0. )
 
 PATCH(IZEQNZD2,PROFIL, 1, 1, 1, 10, 1, 5, 1, 10)
 PLOT(IZEQNZD2,W1 ,-1. ,-1. )
 
 PATCH(TIMEPLOT,PROFIL, 1, 1, 5, 5, 5, 5, 1, 10)
 PLOT(TIMEPLOT,U1 ,0. ,0. )
 PLOT(TIMEPLOT,V1 ,0. ,0. )
 PLOT(TIMEPLOT,W1 ,0. ,0. )
 PLOT(TIMEPLOT,C1 ,0. ,0. )
 
 PATCH(PROF1 ,PROFIL, 1, 1, 1, 10, 1, 10, 1, 1)
 PLOT(PROF1 ,F1 ,0. ,0. )
 
 PATCH(PROF2 ,PROFIL, 1, 1, 1, 10, 1, 10, 1, 1)
 PLOT(PROF2 ,F2 ,0. ,0. )
 
 PATCH(PROF3 ,PROFIL, 1, 1, 1, 10, 1, 10, 1, 1)
 PLOT(PROF3 ,F3 ,0. ,0. )
 
 PATCH(PROF4 ,PROFIL, 1, 1, 1, 10, 1, 10, 1, 1)
 PLOT(PROF4 ,F4 ,0. ,0. )
 
 PATCH(PROF5 ,PROFIL, 1, 1, 1, 10, 1, 10, 1, 1)
 PLOT(PROF5 ,F5 ,0. ,0. )
 
 PATCH(PROF6 ,PROFIL, 1, 1, 1, 10, 1, 10, 1, 1)
 PLOT(PROF6 ,F6 ,0. ,0. )
 
 PATCH(PROF7 ,PROFIL, 1, 1, 1, 10, 1, 10, 1, 1)
 PLOT(PROF7 ,F7 ,0. ,0. )
 
 PATCH(PROF8 ,PROFIL, 1, 1, 1, 10, 1, 10, 1, 1)
 PLOT(PROF8 ,F8 ,0. ,0. )
 
 PATCH(PROF9 ,PROFIL, 1, 1, 1, 10, 1, 10, 1, 1)
 PLOT(PROF9 ,F9 ,0. ,0. )
 
 PATCH(PROF10 ,PROFIL, 1, 1, 1, 10, 1, 10, 1, 1)
 PLOT(PROF10 ,F10 ,0. ,0. )
 
 PATCH(PROF11 ,PROFIL, 1, 1, 1, 10, 1, 10, 1, 1)
 PLOT(PROF11 ,F11 ,0. ,0. )
 
 PATCH(MIDDLE ,PROFIL, 1, 1, 5, 5, 1, 10, 1, 10)
 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