```
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
************************************************************
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.
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 )
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
************************************************************
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
```