TALK=T;RUN(1,1)
  DISPLAY
  Homogeneous turbulence decays behind a grid in a stream of uniform
  velocity. The transport equations of the k-epsilon model reduce to
  ordinary differential equations with analytical solutions of which
  that for k is:

    k=( k0**a - a*e0*(x-x0)/{u0*k0**(1-a)} ) ** (1/a)

  where k0 and e0 are k and e for x=0, and u0 is the flow velocity;
  a=1.0-C2E. The value of C2E that fits the experimental data is 1.9

  The analytical solution is computed in the Q1 file, by way of a DO
  loop, and stored in the variable ANKE. The small differences
  between the computed KE and ANKE result from the upwind
  differencing: ANKE is the cell-centre value, whereas KE is the
  convected-out value.
  ENDDIS
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(Decay Of Turbulence. K-Epsilon Model    )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 14
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 200
 XULAST =1.
 XFRAC(1)=5.0E-03 ;XFRAC(41)=0.205
 XFRAC(81)=0.405 ;XFRAC(121)=0.605
 XFRAC(161)=0.805
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 1
 YVLAST =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(3)=U1 ;NAME(12)=KE
 NAME(13)=EP ;NAME(149)=ANKE
 NAME(150)=VIST
    * Y in SOLUTN argument list denotes:
    * 1-stored 2-solved 3-whole-field
    * 4-point-by-point 5-explicit 6-harmonic averaging 
 SOLUTN(U1,Y,N,N,N,N,Y)
 SOLUTN(KE,Y,Y,N,N,N,N)
 SOLUTN(EP,Y,Y,N,N,N,N)
 SOLUTN(ANKE,Y,N,N,N,N,Y)
 SOLUTN(VIST,Y,N,N,N,N,Y)
 VIST = 150
 ************************************************************
  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(KE,N,Y,Y,Y,Y,N)
 TERMS(EP,N,Y,Y,Y,Y,N)
 DIFCUT =0. ;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 =5.0E-06 ;ENUT = GRND3
 ENUTA =0. ;ENUTB =0. ;ENUTC =0.
 IENUTA = 0
 PRNDTL(KE)=1. ;PRNDTL(EP)=1.
 PRT(KE)=1. ;PRT(EP)=1.314
 CP1 =1. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(U1)=10. ;FIINIT(KE)=1.
 FIINIT(EP)=7.5 ;FIINIT(ANKE)=1.0E-10
 FIINIT(VIST)=1.0E-10
 
 PATCH(IN1 ,INIVAL, 1, 1, 1, 1, 1, 1, 1, 1)
 INIT(IN1 ,ANKE,0. ,0.998128 )
 
 PATCH(IN2 ,INIVAL, 2, 2, 1, 1, 1, 1, 1, 1)
 INIT(IN2 ,ANKE,0. ,0.994405 )
 
 PATCH(IN3 ,INIVAL, 3, 3, 1, 1, 1, 1, 1, 1)
 INIT(IN3 ,ANKE,0. ,0.990709 )
 
 PATCH(IN4 ,INIVAL, 4, 4, 1, 1, 1, 1, 1, 1)
 INIT(IN4 ,ANKE,0. ,0.987038 )
 
 PATCH(IN5 ,INIVAL, 5, 5, 1, 1, 1, 1, 1, 1)
 INIT(IN5 ,ANKE,0. ,0.983394 )
 
 PATCH(IN6 ,INIVAL, 6, 6, 1, 1, 1, 1, 1, 1)
 INIT(IN6 ,ANKE,0. ,0.979776 )
 
 PATCH(IN7 ,INIVAL, 7, 7, 1, 1, 1, 1, 1, 1)
 INIT(IN7 ,ANKE,0. ,0.976182 )
 
 PATCH(IN8 ,INIVAL, 8, 8, 1, 1, 1, 1, 1, 1)
 INIT(IN8 ,ANKE,0. ,0.972615 )
 
 PATCH(IN9 ,INIVAL, 9, 9, 1, 1, 1, 1, 1, 1)
 INIT(IN9 ,ANKE,0. ,0.969072 )
 
 PATCH(IN10 ,INIVAL, 10, 10, 1, 1, 1, 1, 1, 1)
 INIT(IN10 ,ANKE,0. ,0.965554 )
 
 PATCH(IN11 ,INIVAL, 11, 11, 1, 1, 1, 1, 1, 1)
 INIT(IN11 ,ANKE,0. ,0.96206 )
 
 PATCH(IN12 ,INIVAL, 12, 12, 1, 1, 1, 1, 1, 1)
 INIT(IN12 ,ANKE,0. ,0.95859 )
 
 PATCH(IN13 ,INIVAL, 13, 13, 1, 1, 1, 1, 1, 1)
 INIT(IN13 ,ANKE,0. ,0.955145 )
 
 PATCH(IN14 ,INIVAL, 14, 14, 1, 1, 1, 1, 1, 1)
 INIT(IN14 ,ANKE,0. ,0.951723 )
 
 PATCH(IN15 ,INIVAL, 15, 15, 1, 1, 1, 1, 1, 1)
 INIT(IN15 ,ANKE,0. ,0.948324 )
 
 PATCH(IN16 ,INIVAL, 16, 16, 1, 1, 1, 1, 1, 1)
 INIT(IN16 ,ANKE,0. ,0.944949 )
 
 PATCH(IN17 ,INIVAL, 17, 17, 1, 1, 1, 1, 1, 1)
 INIT(IN17 ,ANKE,0. ,0.941597 )
 
 PATCH(IN18 ,INIVAL, 18, 18, 1, 1, 1, 1, 1, 1)
 INIT(IN18 ,ANKE,0. ,0.938267 )
 
 PATCH(IN19 ,INIVAL, 19, 19, 1, 1, 1, 1, 1, 1)
 INIT(IN19 ,ANKE,0. ,0.93496 )
 
 PATCH(IN20 ,INIVAL, 20, 20, 1, 1, 1, 1, 1, 1)
 INIT(IN20 ,ANKE,0. ,0.931676 )
 
 PATCH(IN21 ,INIVAL, 21, 21, 1, 1, 1, 1, 1, 1)
 INIT(IN21 ,ANKE,0. ,0.928413 )
 
 PATCH(IN22 ,INIVAL, 22, 22, 1, 1, 1, 1, 1, 1)
 INIT(IN22 ,ANKE,0. ,0.925173 )
 
 PATCH(IN23 ,INIVAL, 23, 23, 1, 1, 1, 1, 1, 1)
 INIT(IN23 ,ANKE,0. ,0.921954 )
 
 PATCH(IN24 ,INIVAL, 24, 24, 1, 1, 1, 1, 1, 1)
 INIT(IN24 ,ANKE,0. ,0.918756 )
 
 PATCH(IN25 ,INIVAL, 25, 25, 1, 1, 1, 1, 1, 1)
 INIT(IN25 ,ANKE,0. ,0.91558 )
 
 PATCH(IN26 ,INIVAL, 26, 26, 1, 1, 1, 1, 1, 1)
 INIT(IN26 ,ANKE,0. ,0.912424 )
 
 PATCH(IN27 ,INIVAL, 27, 27, 1, 1, 1, 1, 1, 1)
 INIT(IN27 ,ANKE,0. ,0.90929 )
 
 PATCH(IN28 ,INIVAL, 28, 28, 1, 1, 1, 1, 1, 1)
 INIT(IN28 ,ANKE,0. ,0.906176 )
 
 PATCH(IN29 ,INIVAL, 29, 29, 1, 1, 1, 1, 1, 1)
 INIT(IN29 ,ANKE,0. ,0.903082 )
 
 PATCH(IN30 ,INIVAL, 30, 30, 1, 1, 1, 1, 1, 1)
 INIT(IN30 ,ANKE,0. ,0.900009 )
 
 PATCH(IN31 ,INIVAL, 31, 31, 1, 1, 1, 1, 1, 1)
 INIT(IN31 ,ANKE,0. ,0.896956 )
 
 PATCH(IN32 ,INIVAL, 32, 32, 1, 1, 1, 1, 1, 1)
 INIT(IN32 ,ANKE,0. ,0.893922 )
 
 PATCH(IN33 ,INIVAL, 33, 33, 1, 1, 1, 1, 1, 1)
 INIT(IN33 ,ANKE,0. ,0.890908 )
 
 PATCH(IN34 ,INIVAL, 34, 34, 1, 1, 1, 1, 1, 1)
 INIT(IN34 ,ANKE,0. ,0.887914 )
 
 PATCH(IN35 ,INIVAL, 35, 35, 1, 1, 1, 1, 1, 1)
 INIT(IN35 ,ANKE,0. ,0.884939 )
 
 PATCH(IN36 ,INIVAL, 36, 36, 1, 1, 1, 1, 1, 1)
 INIT(IN36 ,ANKE,0. ,0.881983 )
 
 PATCH(IN37 ,INIVAL, 37, 37, 1, 1, 1, 1, 1, 1)
 INIT(IN37 ,ANKE,0. ,0.879046 )
 
 PATCH(IN38 ,INIVAL, 38, 38, 1, 1, 1, 1, 1, 1)
 INIT(IN38 ,ANKE,0. ,0.876127 )
 
 PATCH(IN39 ,INIVAL, 39, 39, 1, 1, 1, 1, 1, 1)
 INIT(IN39 ,ANKE,0. ,0.873227 )
 
 PATCH(IN40 ,INIVAL, 40, 40, 1, 1, 1, 1, 1, 1)
 INIT(IN40 ,ANKE,0. ,0.870346 )
 
 PATCH(IN41 ,INIVAL, 41, 41, 1, 1, 1, 1, 1, 1)
 INIT(IN41 ,ANKE,0. ,0.867483 )
 
 PATCH(IN42 ,INIVAL, 42, 42, 1, 1, 1, 1, 1, 1)
 INIT(IN42 ,ANKE,0. ,0.864637 )
 
 PATCH(IN43 ,INIVAL, 43, 43, 1, 1, 1, 1, 1, 1)
 INIT(IN43 ,ANKE,0. ,0.86181 )
 
 PATCH(IN44 ,INIVAL, 44, 44, 1, 1, 1, 1, 1, 1)
 INIT(IN44 ,ANKE,0. ,0.859 )
 
 PATCH(IN45 ,INIVAL, 45, 45, 1, 1, 1, 1, 1, 1)
 INIT(IN45 ,ANKE,0. ,0.856208 )
 
 PATCH(IN46 ,INIVAL, 46, 46, 1, 1, 1, 1, 1, 1)
 INIT(IN46 ,ANKE,0. ,0.853433 )
 
 PATCH(IN47 ,INIVAL, 47, 47, 1, 1, 1, 1, 1, 1)
 INIT(IN47 ,ANKE,0. ,0.850676 )
 
 PATCH(IN48 ,INIVAL, 48, 48, 1, 1, 1, 1, 1, 1)
 INIT(IN48 ,ANKE,0. ,0.847935 )
 
 PATCH(IN49 ,INIVAL, 49, 49, 1, 1, 1, 1, 1, 1)
 INIT(IN49 ,ANKE,0. ,0.845212 )
 
 PATCH(IN50 ,INIVAL, 50, 50, 1, 1, 1, 1, 1, 1)
 INIT(IN50 ,ANKE,0. ,0.842505 )
 
 PATCH(IN51 ,INIVAL, 51, 51, 1, 1, 1, 1, 1, 1)
 INIT(IN51 ,ANKE,0. ,0.839814 )
 
 PATCH(IN52 ,INIVAL, 52, 52, 1, 1, 1, 1, 1, 1)
 INIT(IN52 ,ANKE,0. ,0.837141 )
 
 PATCH(IN53 ,INIVAL, 53, 53, 1, 1, 1, 1, 1, 1)
 INIT(IN53 ,ANKE,0. ,0.834483 )
 
 PATCH(IN54 ,INIVAL, 54, 54, 1, 1, 1, 1, 1, 1)
 INIT(IN54 ,ANKE,0. ,0.831842 )
 
 PATCH(IN55 ,INIVAL, 55, 55, 1, 1, 1, 1, 1, 1)
 INIT(IN55 ,ANKE,0. ,0.829216 )
 
 PATCH(IN56 ,INIVAL, 56, 56, 1, 1, 1, 1, 1, 1)
 INIT(IN56 ,ANKE,0. ,0.826607 )
 
 PATCH(IN57 ,INIVAL, 57, 57, 1, 1, 1, 1, 1, 1)
 INIT(IN57 ,ANKE,0. ,0.824013 )
 
 PATCH(IN58 ,INIVAL, 58, 58, 1, 1, 1, 1, 1, 1)
 INIT(IN58 ,ANKE,0. ,0.821435 )
 
 PATCH(IN59 ,INIVAL, 59, 59, 1, 1, 1, 1, 1, 1)
 INIT(IN59 ,ANKE,0. ,0.818872 )
 
 PATCH(IN60 ,INIVAL, 60, 60, 1, 1, 1, 1, 1, 1)
 INIT(IN60 ,ANKE,0. ,0.816325 )
 
 PATCH(IN61 ,INIVAL, 61, 61, 1, 1, 1, 1, 1, 1)
 INIT(IN61 ,ANKE,0. ,0.813792 )
 
 PATCH(IN62 ,INIVAL, 62, 62, 1, 1, 1, 1, 1, 1)
 INIT(IN62 ,ANKE,0. ,0.811275 )
 
 PATCH(IN63 ,INIVAL, 63, 63, 1, 1, 1, 1, 1, 1)
 INIT(IN63 ,ANKE,0. ,0.808773 )
 
 PATCH(IN64 ,INIVAL, 64, 64, 1, 1, 1, 1, 1, 1)
 INIT(IN64 ,ANKE,0. ,0.806285 )
 
 PATCH(IN65 ,INIVAL, 65, 65, 1, 1, 1, 1, 1, 1)
 INIT(IN65 ,ANKE,0. ,0.803812 )
 
 PATCH(IN66 ,INIVAL, 66, 66, 1, 1, 1, 1, 1, 1)
 INIT(IN66 ,ANKE,0. ,0.801354 )
 
 PATCH(IN67 ,INIVAL, 67, 67, 1, 1, 1, 1, 1, 1)
 INIT(IN67 ,ANKE,0. ,0.79891 )
 
 PATCH(IN68 ,INIVAL, 68, 68, 1, 1, 1, 1, 1, 1)
 INIT(IN68 ,ANKE,0. ,0.79648 )
 
 PATCH(IN69 ,INIVAL, 69, 69, 1, 1, 1, 1, 1, 1)
 INIT(IN69 ,ANKE,0. ,0.794065 )
 
 PATCH(IN70 ,INIVAL, 70, 70, 1, 1, 1, 1, 1, 1)
 INIT(IN70 ,ANKE,0. ,0.791663 )
 
 PATCH(IN71 ,INIVAL, 71, 71, 1, 1, 1, 1, 1, 1)
 INIT(IN71 ,ANKE,0. ,0.789275 )
 
 PATCH(IN72 ,INIVAL, 72, 72, 1, 1, 1, 1, 1, 1)
 INIT(IN72 ,ANKE,0. ,0.786901 )
 
 PATCH(IN73 ,INIVAL, 73, 73, 1, 1, 1, 1, 1, 1)
 INIT(IN73 ,ANKE,0. ,0.784541 )
 
 PATCH(IN74 ,INIVAL, 74, 74, 1, 1, 1, 1, 1, 1)
 INIT(IN74 ,ANKE,0. ,0.782195 )
 
 PATCH(IN75 ,INIVAL, 75, 75, 1, 1, 1, 1, 1, 1)
 INIT(IN75 ,ANKE,0. ,0.779862 )
 
 PATCH(IN76 ,INIVAL, 76, 76, 1, 1, 1, 1, 1, 1)
 INIT(IN76 ,ANKE,0. ,0.777542 )
 
 PATCH(IN77 ,INIVAL, 77, 77, 1, 1, 1, 1, 1, 1)
 INIT(IN77 ,ANKE,0. ,0.775235 )
 
 PATCH(IN78 ,INIVAL, 78, 78, 1, 1, 1, 1, 1, 1)
 INIT(IN78 ,ANKE,0. ,0.772941 )
 
 PATCH(IN79 ,INIVAL, 79, 79, 1, 1, 1, 1, 1, 1)
 INIT(IN79 ,ANKE,0. ,0.770661 )
 
 PATCH(IN80 ,INIVAL, 80, 80, 1, 1, 1, 1, 1, 1)
 INIT(IN80 ,ANKE,0. ,0.768393 )
 
 PATCH(IN81 ,INIVAL, 81, 81, 1, 1, 1, 1, 1, 1)
 INIT(IN81 ,ANKE,0. ,0.766138 )
 
 PATCH(IN82 ,INIVAL, 82, 82, 1, 1, 1, 1, 1, 1)
 INIT(IN82 ,ANKE,0. ,0.763896 )
 
 PATCH(IN83 ,INIVAL, 83, 83, 1, 1, 1, 1, 1, 1)
 INIT(IN83 ,ANKE,0. ,0.761666 )
 
 PATCH(IN84 ,INIVAL, 84, 84, 1, 1, 1, 1, 1, 1)
 INIT(IN84 ,ANKE,0. ,0.759449 )
 
 PATCH(IN85 ,INIVAL, 85, 85, 1, 1, 1, 1, 1, 1)
 INIT(IN85 ,ANKE,0. ,0.757244 )
 
 PATCH(IN86 ,INIVAL, 86, 86, 1, 1, 1, 1, 1, 1)
 INIT(IN86 ,ANKE,0. ,0.755052 )
 
 PATCH(IN87 ,INIVAL, 87, 87, 1, 1, 1, 1, 1, 1)
 INIT(IN87 ,ANKE,0. ,0.752871 )
 
 PATCH(IN88 ,INIVAL, 88, 88, 1, 1, 1, 1, 1, 1)
 INIT(IN88 ,ANKE,0. ,0.750703 )
 
 PATCH(IN89 ,INIVAL, 89, 89, 1, 1, 1, 1, 1, 1)
 INIT(IN89 ,ANKE,0. ,0.748547 )
 
 PATCH(IN90 ,INIVAL, 90, 90, 1, 1, 1, 1, 1, 1)
 INIT(IN90 ,ANKE,0. ,0.746402 )
 
 PATCH(IN91 ,INIVAL, 91, 91, 1, 1, 1, 1, 1, 1)
 INIT(IN91 ,ANKE,0. ,0.744269 )
 
 PATCH(IN92 ,INIVAL, 92, 92, 1, 1, 1, 1, 1, 1)
 INIT(IN92 ,ANKE,0. ,0.742148 )
 
 PATCH(IN93 ,INIVAL, 93, 93, 1, 1, 1, 1, 1, 1)
 INIT(IN93 ,ANKE,0. ,0.740039 )
 
 PATCH(IN94 ,INIVAL, 94, 94, 1, 1, 1, 1, 1, 1)
 INIT(IN94 ,ANKE,0. ,0.737941 )
 
 PATCH(IN95 ,INIVAL, 95, 95, 1, 1, 1, 1, 1, 1)
 INIT(IN95 ,ANKE,0. ,0.735854 )
 
 PATCH(IN96 ,INIVAL, 96, 96, 1, 1, 1, 1, 1, 1)
 INIT(IN96 ,ANKE,0. ,0.733779 )
 
 PATCH(IN97 ,INIVAL, 97, 97, 1, 1, 1, 1, 1, 1)
 INIT(IN97 ,ANKE,0. ,0.731714 )
 
 PATCH(IN98 ,INIVAL, 98, 98, 1, 1, 1, 1, 1, 1)
 INIT(IN98 ,ANKE,0. ,0.729661 )
 
 PATCH(IN99 ,INIVAL, 99, 99, 1, 1, 1, 1, 1, 1)
 INIT(IN99 ,ANKE,0. ,0.727619 )
 
 PATCH(IN100 ,INIVAL, 100, 100, 1, 1, 1, 1, 1, 1)
 INIT(IN100 ,ANKE,0. ,0.725588 )
 
 PATCH(IN101 ,INIVAL, 101, 101, 1, 1, 1, 1, 1, 1)
 INIT(IN101 ,ANKE,0. ,0.723568 )
 
 PATCH(IN102 ,INIVAL, 102, 102, 1, 1, 1, 1, 1, 1)
 INIT(IN102 ,ANKE,0. ,0.721559 )
 
 PATCH(IN103 ,INIVAL, 103, 103, 1, 1, 1, 1, 1, 1)
 INIT(IN103 ,ANKE,0. ,0.71956 )
 
 PATCH(IN104 ,INIVAL, 104, 104, 1, 1, 1, 1, 1, 1)
 INIT(IN104 ,ANKE,0. ,0.717572 )
 
 PATCH(IN105 ,INIVAL, 105, 105, 1, 1, 1, 1, 1, 1)
 INIT(IN105 ,ANKE,0. ,0.715594 )
 
 PATCH(IN106 ,INIVAL, 106, 106, 1, 1, 1, 1, 1, 1)
 INIT(IN106 ,ANKE,0. ,0.713627 )
 
 PATCH(IN107 ,INIVAL, 107, 107, 1, 1, 1, 1, 1, 1)
 INIT(IN107 ,ANKE,0. ,0.71167 )
 
 PATCH(IN108 ,INIVAL, 108, 108, 1, 1, 1, 1, 1, 1)
 INIT(IN108 ,ANKE,0. ,0.709724 )
 
 PATCH(IN109 ,INIVAL, 109, 109, 1, 1, 1, 1, 1, 1)
 INIT(IN109 ,ANKE,0. ,0.707787 )
 
 PATCH(IN110 ,INIVAL, 110, 110, 1, 1, 1, 1, 1, 1)
 INIT(IN110 ,ANKE,0. ,0.705861 )
 
 PATCH(IN111 ,INIVAL, 111, 111, 1, 1, 1, 1, 1, 1)
 INIT(IN111 ,ANKE,0. ,0.703945 )
 
 PATCH(IN112 ,INIVAL, 112, 112, 1, 1, 1, 1, 1, 1)
 INIT(IN112 ,ANKE,0. ,0.702039 )
 
 PATCH(IN113 ,INIVAL, 113, 113, 1, 1, 1, 1, 1, 1)
 INIT(IN113 ,ANKE,0. ,0.700142 )
 
 PATCH(IN114 ,INIVAL, 114, 114, 1, 1, 1, 1, 1, 1)
 INIT(IN114 ,ANKE,0. ,0.698256 )
 
 PATCH(IN115 ,INIVAL, 115, 115, 1, 1, 1, 1, 1, 1)
 INIT(IN115 ,ANKE,0. ,0.696379 )
 
 PATCH(IN116 ,INIVAL, 116, 116, 1, 1, 1, 1, 1, 1)
 INIT(IN116 ,ANKE,0. ,0.694512 )
 
 PATCH(IN117 ,INIVAL, 117, 117, 1, 1, 1, 1, 1, 1)
 INIT(IN117 ,ANKE,0. ,0.692654 )
 
 PATCH(IN118 ,INIVAL, 118, 118, 1, 1, 1, 1, 1, 1)
 INIT(IN118 ,ANKE,0. ,0.690806 )
 
 PATCH(IN119 ,INIVAL, 119, 119, 1, 1, 1, 1, 1, 1)
 INIT(IN119 ,ANKE,0. ,0.688968 )
 
 PATCH(IN120 ,INIVAL, 120, 120, 1, 1, 1, 1, 1, 1)
 INIT(IN120 ,ANKE,0. ,0.687138 )
 
 PATCH(IN121 ,INIVAL, 121, 121, 1, 1, 1, 1, 1, 1)
 INIT(IN121 ,ANKE,0. ,0.685319 )
 
 PATCH(IN122 ,INIVAL, 122, 122, 1, 1, 1, 1, 1, 1)
 INIT(IN122 ,ANKE,0. ,0.683508 )
 
 PATCH(IN123 ,INIVAL, 123, 123, 1, 1, 1, 1, 1, 1)
 INIT(IN123 ,ANKE,0. ,0.681706 )
 
 PATCH(IN124 ,INIVAL, 124, 124, 1, 1, 1, 1, 1, 1)
 INIT(IN124 ,ANKE,0. ,0.679914 )
 
 PATCH(IN125 ,INIVAL, 125, 125, 1, 1, 1, 1, 1, 1)
 INIT(IN125 ,ANKE,0. ,0.678131 )
 
 PATCH(IN126 ,INIVAL, 126, 126, 1, 1, 1, 1, 1, 1)
 INIT(IN126 ,ANKE,0. ,0.676356 )
 
 PATCH(IN127 ,INIVAL, 127, 127, 1, 1, 1, 1, 1, 1)
 INIT(IN127 ,ANKE,0. ,0.674591 )
 
 PATCH(IN128 ,INIVAL, 128, 128, 1, 1, 1, 1, 1, 1)
 INIT(IN128 ,ANKE,0. ,0.672834 )
 
 PATCH(IN129 ,INIVAL, 129, 129, 1, 1, 1, 1, 1, 1)
 INIT(IN129 ,ANKE,0. ,0.671086 )
 
 PATCH(IN130 ,INIVAL, 130, 130, 1, 1, 1, 1, 1, 1)
 INIT(IN130 ,ANKE,0. ,0.669347 )
 
 PATCH(IN131 ,INIVAL, 131, 131, 1, 1, 1, 1, 1, 1)
 INIT(IN131 ,ANKE,0. ,0.667616 )
 
 PATCH(IN132 ,INIVAL, 132, 132, 1, 1, 1, 1, 1, 1)
 INIT(IN132 ,ANKE,0. ,0.665894 )
 
 PATCH(IN133 ,INIVAL, 133, 133, 1, 1, 1, 1, 1, 1)
 INIT(IN133 ,ANKE,0. ,0.66418 )
 
 PATCH(IN134 ,INIVAL, 134, 134, 1, 1, 1, 1, 1, 1)
 INIT(IN134 ,ANKE,0. ,0.662475 )
 
 PATCH(IN135 ,INIVAL, 135, 135, 1, 1, 1, 1, 1, 1)
 INIT(IN135 ,ANKE,0. ,0.660779 )
 
 PATCH(IN136 ,INIVAL, 136, 136, 1, 1, 1, 1, 1, 1)
 INIT(IN136 ,ANKE,0. ,0.65909 )
 
 PATCH(IN137 ,INIVAL, 137, 137, 1, 1, 1, 1, 1, 1)
 INIT(IN137 ,ANKE,0. ,0.65741 )
 
 PATCH(IN138 ,INIVAL, 138, 138, 1, 1, 1, 1, 1, 1)
 INIT(IN138 ,ANKE,0. ,0.655738 )
 
 PATCH(IN139 ,INIVAL, 139, 139, 1, 1, 1, 1, 1, 1)
 INIT(IN139 ,ANKE,0. ,0.654074 )
 
 PATCH(IN140 ,INIVAL, 140, 140, 1, 1, 1, 1, 1, 1)
 INIT(IN140 ,ANKE,0. ,0.652419 )
 
 PATCH(IN141 ,INIVAL, 141, 141, 1, 1, 1, 1, 1, 1)
 INIT(IN141 ,ANKE,0. ,0.650771 )
 
 PATCH(IN142 ,INIVAL, 142, 142, 1, 1, 1, 1, 1, 1)
 INIT(IN142 ,ANKE,0. ,0.649131 )
 
 PATCH(IN143 ,INIVAL, 143, 143, 1, 1, 1, 1, 1, 1)
 INIT(IN143 ,ANKE,0. ,0.6475 )
 
 PATCH(IN144 ,INIVAL, 144, 144, 1, 1, 1, 1, 1, 1)
 INIT(IN144 ,ANKE,0. ,0.645876 )
 
 PATCH(IN145 ,INIVAL, 145, 145, 1, 1, 1, 1, 1, 1)
 INIT(IN145 ,ANKE,0. ,0.64426 )
 
 PATCH(IN146 ,INIVAL, 146, 146, 1, 1, 1, 1, 1, 1)
 INIT(IN146 ,ANKE,0. ,0.642651 )
 
 PATCH(IN147 ,INIVAL, 147, 147, 1, 1, 1, 1, 1, 1)
 INIT(IN147 ,ANKE,0. ,0.641051 )
 
 PATCH(IN148 ,INIVAL, 148, 148, 1, 1, 1, 1, 1, 1)
 INIT(IN148 ,ANKE,0. ,0.639458 )
 
 PATCH(IN149 ,INIVAL, 149, 149, 1, 1, 1, 1, 1, 1)
 INIT(IN149 ,ANKE,0. ,0.637872 )
 
 PATCH(IN150 ,INIVAL, 150, 150, 1, 1, 1, 1, 1, 1)
 INIT(IN150 ,ANKE,0. ,0.636294 )
 
 PATCH(IN151 ,INIVAL, 151, 151, 1, 1, 1, 1, 1, 1)
 INIT(IN151 ,ANKE,0. ,0.634724 )
 
 PATCH(IN152 ,INIVAL, 152, 152, 1, 1, 1, 1, 1, 1)
 INIT(IN152 ,ANKE,0. ,0.633161 )
 
 PATCH(IN153 ,INIVAL, 153, 153, 1, 1, 1, 1, 1, 1)
 INIT(IN153 ,ANKE,0. ,0.631605 )
 
 PATCH(IN154 ,INIVAL, 154, 154, 1, 1, 1, 1, 1, 1)
 INIT(IN154 ,ANKE,0. ,0.630057 )
 
 PATCH(IN155 ,INIVAL, 155, 155, 1, 1, 1, 1, 1, 1)
 INIT(IN155 ,ANKE,0. ,0.628516 )
 
 PATCH(IN156 ,INIVAL, 156, 156, 1, 1, 1, 1, 1, 1)
 INIT(IN156 ,ANKE,0. ,0.626982 )
 
 PATCH(IN157 ,INIVAL, 157, 157, 1, 1, 1, 1, 1, 1)
 INIT(IN157 ,ANKE,0. ,0.625455 )
 
 PATCH(IN158 ,INIVAL, 158, 158, 1, 1, 1, 1, 1, 1)
 INIT(IN158 ,ANKE,0. ,0.623936 )
 
 PATCH(IN159 ,INIVAL, 159, 159, 1, 1, 1, 1, 1, 1)
 INIT(IN159 ,ANKE,0. ,0.622423 )
 
 PATCH(IN160 ,INIVAL, 160, 160, 1, 1, 1, 1, 1, 1)
 INIT(IN160 ,ANKE,0. ,0.620918 )
 
 PATCH(IN161 ,INIVAL, 161, 161, 1, 1, 1, 1, 1, 1)
 INIT(IN161 ,ANKE,0. ,0.619419 )
 
 PATCH(IN162 ,INIVAL, 162, 162, 1, 1, 1, 1, 1, 1)
 INIT(IN162 ,ANKE,0. ,0.617928 )
 
 PATCH(IN163 ,INIVAL, 163, 163, 1, 1, 1, 1, 1, 1)
 INIT(IN163 ,ANKE,0. ,0.616443 )
 
 PATCH(IN164 ,INIVAL, 164, 164, 1, 1, 1, 1, 1, 1)
 INIT(IN164 ,ANKE,0. ,0.614965 )
 
 PATCH(IN165 ,INIVAL, 165, 165, 1, 1, 1, 1, 1, 1)
 INIT(IN165 ,ANKE,0. ,0.613494 )
 
 PATCH(IN166 ,INIVAL, 166, 166, 1, 1, 1, 1, 1, 1)
 INIT(IN166 ,ANKE,0. ,0.61203 )
 
 PATCH(IN167 ,INIVAL, 167, 167, 1, 1, 1, 1, 1, 1)
 INIT(IN167 ,ANKE,0. ,0.610572 )
 
 PATCH(IN168 ,INIVAL, 168, 168, 1, 1, 1, 1, 1, 1)
 INIT(IN168 ,ANKE,0. ,0.609121 )
 
 PATCH(IN169 ,INIVAL, 169, 169, 1, 1, 1, 1, 1, 1)
 INIT(IN169 ,ANKE,0. ,0.607677 )
 
 PATCH(IN170 ,INIVAL, 170, 170, 1, 1, 1, 1, 1, 1)
 INIT(IN170 ,ANKE,0. ,0.606239 )
 
 PATCH(IN171 ,INIVAL, 171, 171, 1, 1, 1, 1, 1, 1)
 INIT(IN171 ,ANKE,0. ,0.604808 )
 
 PATCH(IN172 ,INIVAL, 172, 172, 1, 1, 1, 1, 1, 1)
 INIT(IN172 ,ANKE,0. ,0.603383 )
 
 PATCH(IN173 ,INIVAL, 173, 173, 1, 1, 1, 1, 1, 1)
 INIT(IN173 ,ANKE,0. ,0.601965 )
 
 PATCH(IN174 ,INIVAL, 174, 174, 1, 1, 1, 1, 1, 1)
 INIT(IN174 ,ANKE,0. ,0.600553 )
 
 PATCH(IN175 ,INIVAL, 175, 175, 1, 1, 1, 1, 1, 1)
 INIT(IN175 ,ANKE,0. ,0.599147 )
 
 PATCH(IN176 ,INIVAL, 176, 176, 1, 1, 1, 1, 1, 1)
 INIT(IN176 ,ANKE,0. ,0.597748 )
 
 PATCH(IN177 ,INIVAL, 177, 177, 1, 1, 1, 1, 1, 1)
 INIT(IN177 ,ANKE,0. ,0.596355 )
 
 PATCH(IN178 ,INIVAL, 178, 178, 1, 1, 1, 1, 1, 1)
 INIT(IN178 ,ANKE,0. ,0.594968 )
 
 PATCH(IN179 ,INIVAL, 179, 179, 1, 1, 1, 1, 1, 1)
 INIT(IN179 ,ANKE,0. ,0.593587 )
 
 PATCH(IN180 ,INIVAL, 180, 180, 1, 1, 1, 1, 1, 1)
 INIT(IN180 ,ANKE,0. ,0.592213 )
 
 PATCH(IN181 ,INIVAL, 181, 181, 1, 1, 1, 1, 1, 1)
 INIT(IN181 ,ANKE,0. ,0.590844 )
 
 PATCH(IN182 ,INIVAL, 182, 182, 1, 1, 1, 1, 1, 1)
 INIT(IN182 ,ANKE,0. ,0.589482 )
 
 PATCH(IN183 ,INIVAL, 183, 183, 1, 1, 1, 1, 1, 1)
 INIT(IN183 ,ANKE,0. ,0.588125 )
 
 PATCH(IN184 ,INIVAL, 184, 184, 1, 1, 1, 1, 1, 1)
 INIT(IN184 ,ANKE,0. ,0.586775 )
 
 PATCH(IN185 ,INIVAL, 185, 185, 1, 1, 1, 1, 1, 1)
 INIT(IN185 ,ANKE,0. ,0.585431 )
 
 PATCH(IN186 ,INIVAL, 186, 186, 1, 1, 1, 1, 1, 1)
 INIT(IN186 ,ANKE,0. ,0.584092 )
 
 PATCH(IN187 ,INIVAL, 187, 187, 1, 1, 1, 1, 1, 1)
 INIT(IN187 ,ANKE,0. ,0.582759 )
 
 PATCH(IN188 ,INIVAL, 188, 188, 1, 1, 1, 1, 1, 1)
 INIT(IN188 ,ANKE,0. ,0.581433 )
 
 PATCH(IN189 ,INIVAL, 189, 189, 1, 1, 1, 1, 1, 1)
 INIT(IN189 ,ANKE,0. ,0.580111 )
 
 PATCH(IN190 ,INIVAL, 190, 190, 1, 1, 1, 1, 1, 1)
 INIT(IN190 ,ANKE,0. ,0.578796 )
 
 PATCH(IN191 ,INIVAL, 191, 191, 1, 1, 1, 1, 1, 1)
 INIT(IN191 ,ANKE,0. ,0.577487 )
 
 PATCH(IN192 ,INIVAL, 192, 192, 1, 1, 1, 1, 1, 1)
 INIT(IN192 ,ANKE,0. ,0.576183 )
 
 PATCH(IN193 ,INIVAL, 193, 193, 1, 1, 1, 1, 1, 1)
 INIT(IN193 ,ANKE,0. ,0.574884 )
 
 PATCH(IN194 ,INIVAL, 194, 194, 1, 1, 1, 1, 1, 1)
 INIT(IN194 ,ANKE,0. ,0.573592 )
 
 PATCH(IN195 ,INIVAL, 195, 195, 1, 1, 1, 1, 1, 1)
 INIT(IN195 ,ANKE,0. ,0.572305 )
 
 PATCH(IN196 ,INIVAL, 196, 196, 1, 1, 1, 1, 1, 1)
 INIT(IN196 ,ANKE,0. ,0.571023 )
 
 PATCH(IN197 ,INIVAL, 197, 197, 1, 1, 1, 1, 1, 1)
 INIT(IN197 ,ANKE,0. ,0.569747 )
 
 PATCH(IN198 ,INIVAL, 198, 198, 1, 1, 1, 1, 1, 1)
 INIT(IN198 ,ANKE,0. ,0.568476 )
 
 PATCH(IN199 ,INIVAL, 199, 199, 1, 1, 1, 1, 1, 1)
 INIT(IN199 ,ANKE,0. ,0.567211 )
 
 PATCH(IN200 ,INIVAL, 200, 200, 1, 1, 1, 1, 1, 1)
 INIT(IN200 ,ANKE,0. ,0.565951 )
 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(INLET ,WEST , 1, 1, 1, 1, 1, 1, 1, 1)
 COVAL(INLET ,KE ,10. ,1. )
 COVAL(INLET ,EP ,10. ,7.5 )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 1000 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-04
 ************************************************************
  Group 16. Terminate Iterations
 LITER(KE)=20 ;LITER(EP)=20
 ENDIT(KE)=1.0E-03 ;ENDIT(EP)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(U1,LINRLX,1.)
 RELAX(KE,LINRLX,0.5)
 RELAX(EP,LINRLX,0.5)
 RELAX(ANKE,LINRLX,1.)
 RELAX(VIST,LINRLX,1.)
 KELIN = 0
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
 VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
 VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
 VARMAX(ANKE)=1.0E+10 ;VARMIN(ANKE)=-1.0E+10
 VARMAX(VIST)=1.0E+10 ;VARMIN(VIST)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 GENK = T
 PARSOL = F
 ISG62 = 1
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  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(U1,N,N,N,N,N,N)
 OUTPUT(KE,Y,N,Y,Y,Y,Y)
 OUTPUT(EP,Y,N,Y,Y,Y,Y)
 OUTPUT(ANKE,Y,N,Y,N,N,N)
 OUTPUT(VIST,N,N,N,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 200 ;IYMON = 1 ;IZMON = 1
 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = 1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NXPRIN = 40 ;IXPRF = 1 ;IXPRL = 10000
 IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
 ISWPRF = 1 ;ISWPRL = 100000
 ITABL = 3 ;IPROF = 1
 ABSIZ =0.5 ;ORSIZ =0.4
 NTZPRF = 1 ;NCOLPF = 50
 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
 
 PATCH(PROFILE ,PROFIL, 1, 200, 1, 1, 1, 1, 1, 1)
 PLOT(PROFILE ,KE ,0. ,1.5 )
 PLOT(PROFILE ,EP ,0. ,15. )
 PLOT(PROFILE ,VIST,0. ,0. )
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP