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



     set prop off
     gr ou z 1 x 12 25 y 19 21
     gr ou z 1 x 12 25 y 4 18
     gr ou z 1 x 14 23 y 7 18
     gr ou z 1 x 12 25 y 22 31
     gr ou z 1 x 14 23 y 22 29
     *gr ou z 1 x 11 11 y 25 28
     *gr ou z 1 x 26 26 y 25 28
     gr z 1 x 11 11 y 20 21
     gr z 1 x 26 26 y 20 21
     *gr ou z 1 x 26 26 y 7 10
     *gr ou z 1 x 11 11 y 15 18
     gr ou z 1 x 14 16 y 1 3
     gr ou z 1 x 21 23 y 1 3
     gr ou z 1 x 11 11 y 16 17
     gr ou x 11 11 y 15 18
     gr ou z 1 x 26 26 y 8 9
     gr ou x 27 y 7 10
     gr ou z 1 x 11 11 y 26 27
     gr ou x 11 11 y 25 28
     gr ou z 1 x 26 26 y 26 27
     gr ou x 27 y 25 28
     gr ou x 15 15 y 19 21
     gr ou x 16 16 y 19 21
     gr ou x 17 17 y 19 21
     gr ou x 18 18 y 19 21
     gr ou x 19 19 y 19 21
     gr ou x 20 20 y 19 21
     gr ou x 21 21 y 19 21
     gr ou x 22 22 y 19 21
     gr ou x 23 23 y 19 21
     gr ou z 1 x 15 22 y 8 18
     gr ou z 1 x 16 21 y 9 18
     gr ou z 1 x 17 20 y 10 18
     gr ou z 1 x 18 19 y 11 18
     gr ou y 11 11 x 17 23
     gr ou y 15 15 x 14 20
     gr ou x 19 19 y 12 29
     set vec comp uu1 vv1 -
     vec z 1 sh
     set vec comp u2 v2 -
     vec z 1
     msg( Shell and tube fluid velocity vectors
     pause
     vec cl
     red
     set vec comp u1 v1 -
     vec z 1 x 11 m sh
     msg( Thermal displacement vectors in shell and tube bundle
  ENDUSE
  DISPLAY
    A 2D shell-and-tube heat-exchanger is used to  exemplify
    essential ideas of HEXAGON model,  which is probably the
    first  to  show  how  the   thermo-hydraulics   of   the
    shell-side  and tube-side fluids could be simultaneously
    computed with the displacements and thermal stresses  in
    tubes   and   shell   to   be   included   in   a   SFT,
    Solid-Fluid-Thermal, heat-exchanger analysis.

    The case illustrates the use a single  computer  program
    to  calculate  from  the  partial-differential equations
    governing relevant fluid processes the distributions of:
    *  shell-side fluid velocity components;
    *  the corresponding temperatures and pressures;
    *  the tube-side fluid velocity components;
    *  the corresponding temperatures and pressures;
    *  the tube metal temperatures; and
    *  the displacements and stresses in the tubes and the shell.

    The heat  exchanger  considered  is  an  imaginary  one,
    having  two  baffles  within the shell,  with the U-bend
    tubes arranged in  array  and  header  distributing  the
    in-fluid between the tubes and collecting out-fluid.

    No attempt  has  been  made to pick-up and implement the
    actual resistance formulae  which  are  widely  used  in
    thermal  engineering.  But  because  PLANT  is  used  to
    represent them,  the artificial formulae can  be  easily
    replaced by required ones.

    The heat exchanger is a rectangular box,  2.0m high,  1m
    wide and 1m long.  It consists of the header,  the hight
    of which is 0.8m and shell closed at the bottom and open
    at the top.  The header is divided into two halfes by  a
    vertical plate.

    The shell  is uniformly filled with the tubes.  The tube
    fluid (water) enters the header through the inlet at its
    west  wall,  flows  downwards in west half of the shell,
    turns through the U bend at the bottom and rises  upward
    in  the  other  shell half to enter the east half of the
    header going out through the outlet at header east side.

    The shell fluid ( air ) entering the shell  through  the
    inlet  at  the  east  wall  is  made to pass between two
    baffles in a zig-zag manner,  until it goes out  through
    the outlet at the top of the west wall of shell.

    Only X-Y  plane  of  the  exchanger  is  included in the
    calculation domain, because of 2D-nature of analysis.

    A uniform 26*32*1 grid is used, to cover computational
    space.

    PLANT is used to:
     * set and/or compute the fluid properties in sub-domains;
     * introduce the non-linear flow resistances ;
     * calculate the distribution of overall heat transfer
       coefficient;
     * calculate the distribution of the tube metal
       temperatures;
     * link the sub-domains for data transfer and
       manipulations and
     * output data processing.
   ENDDIS
    PLANTBEGIN
  **  Set the fluid densities:
               Header    : water
               Tube side : water
               Shell side: air
    DEN1=1000.
   REGION() 1
    DEN1=1000.
   REGION() 2
    DEN1=1.2
   REGION() 3
     The above three statements,  followed  by  the  pointer
     RHO1=GRND  and parameterized REGION commands,  instruct
     PLANT  to  make  the  density  distributions   as   the
     distribition of in-cell marker values dictates.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
   ** Set fluid viscosities
        Header    : effective viscosity proprtional to
                    local velocity magnitude and distance to
                    nearest wall.
        Tube side : Constant=0.01
        Shell side: as for header
    VISL=1.*SQRT(U1**2+V1**2)*WDIS
   REGION() 1
    VISL=0.01
   REGION() 2
    VISL=1.*SQRT(U1**2+V1**2)*WDIS
   REGION() 3
     The above  three statements do the same for viscosities
     as previous three has done  for  densities.  Note  that
     viscosities  in  the  domains  marked  1 and 3 are made
     proprtional  to  the   products   of   local   velocity
     magnitudes and distances to the nearest wall.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
   **  Non-linear resistance to tube-fluid flow exerted
       by tubes, throughout the U-tube array.
PATCH(SS002U,PHASEM,1,10,1,NY,1,NZ,1,1)
   CO=.2*(U1**2+V1**2)**0.15
COVAL(SS002U,U1,GRND,0.0)
   CO=.2*(U1**2+V1**2)**0.15
COVAL(SS002U,V1,GRND,0.0)
    Momentum sinks are introduced by above formulae over all
    cells  having  marker  value  appearing in the number of
    PATCH name, 002.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
   **  Non-linear resistance to shell-fluid flow exerted
       by tubes, throughout the shell-side.
PATCH(SS003H,PHASEM,1,10,1,NY,1,NZ,1,1)
   CO=2.2*(U1**2+V1**2)**0.25
COVAL(SS003H,U1,GRND,0.0)
   CO=2.2*(U1**2+V1**2)**0.25
COVAL(SS003H,V1,GRND,0.0)
    Momentum sinks are introduced by above formulae over all
    cells  having  marker  value  appearing in the number of
    PATCH name, 003.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
    **  Tube fluid heat transfer coefficient
STORE(ALF2);FIINIT(ALF2)=0.0
    ALF2=1.+1.*SQRT(U1**2+V1**2+TINY)
   REGION() 2
    **  Shell fluid heat transfer coefficient
STORE(ALF3);FIINIT(ALF3)=0.0
    ALF3=1.+3.*SQRT(U1**2+V1**2+TINY)
   REGION() 3
    Tube and  shell  fluid heat transfer coefficients,  ALF2
    and  ALF3,  are  made  dependent   on   local   velocity
    magnitudes   over   the   REGIONs   marked   2   and   3
    correspondingly.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
    **  Overall heat transfer coefficient
STORE(HTC);FIINIT(HTC)=0.0
    HTC=1./(1/ALF2+1/ALF3[,-IG(1),])
   REGION() 2
    HTC=1./(1/ALF3 +1/ALF2[,+IG(1),])
   REGION() 3
    Overall heat transfer coefficient, HTC, distribution are
    calculated  by  reference  to  appropriate  local   heat
    transfer coefficients in REGIONs 2 and 3.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
    **  Heat-exchange with shell-fluid, throughout the shell.
PATCH(SS002T,PHASEM,1,NX,1,NY,1,NZ,1,1)
   CO =HTC
   VAL=TEMP[,-IG(1),]
COVAL(SS002T,TEMP,GRND,GRND)
    **  Heat-exchange with tube-fluid, throughout the shell.
PATCH(SS003S,PHASEM,1,NX,1,NY,1,NZ,1,1)
   CO =HTC
   VAL=TEMP[,+IG(1),]
COVAL(SS003S,TEMP,GRND,GRND)
    The PATCH names indicate the sub-domain cell markers , 2
    and 3, over which the heat-exchange sources are applied.
    The  indicial  operations  for  TEMP  are  arranged   in
    appropriate manner.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<

    ===========================================================
       Data preparation for conjugate and stress analysis
    ===========================================================
    ** Transfer shell fluid temperatures
PATCH(SS005T,CELL,1,NX,1,NY,1,NZ,1,1)
   CO=1.e10
   VAL=TEMP[-13,-6,]
COVAL(SS005T,TEM1,GRND,GRND)
    The temperatures   of   the   shell   fluid,  TEMP,  are
    transfered into the stress analysis sub-domain,  MARK=5,
    to be used as TEM1.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
    **  Tube-wall temperature
PATCH(SS100T,CELL,1,NX,1,NY,1,NZ,1,1)
   CO=1.e10
   VAL=(ALF2[-13,+6,]*TEMP[-13,+6,]+                ALF3[-1$
3,-6,]*TEMP[-13,-6,])          /(ALF2[-13,+6,]+ALF3[-13,-6,]+TINY)
          /(ALF2[-13,+6,]+ALF3[-13,-6,]+TINY)
COVAL(SS100T,TEM1,GRND,GRND)
    Here the tube wall temperatures, TEM1, are calculated in
    the sub-domain  indicated  by  MARK=100  as  PATCH  name
    number specifies.  TEM1s are computed via shell and tube
    fluid  temperatures  and  heat   transfer   coefficients
    transfered  from  cooresponding  sub-domains as indicial
    numbers show.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
    ** Transfer the header temperatures
PATCH(SS004T,CELL,1,NX,1,NY,1,NZ,1,1)
   CO=1.e10
   VAL=TEMP[-13,+3,]
COVAL(SS004T,TEM1,GRND,GRND)
    The temperatures of the header  tube  fluid,  TEMP,  are
    transfered into the stress analysis sub-domain,  MARK=4,
    to be used as TEM1.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<


    ===========================================================
       Output data processing
    ===========================================================
     ** Tube fluid velocities transfer
STORE(UU1,VV1)
    UU1=U1[-13,+3,]
   REGION(14,23,22,29,1,NZ)
    VV1=V1[-13,+3,]
   REGION(14,23,22,29,1,NZ)
    UU1=U1[-13,+6,]
   REGION(14,23,7,18,1,NZ)
    VV1=V1[-13,+6,]
   REGION(14,23,7,18,1,NZ)
    Tube fluid velocities are  transfered  from  where  they
    have been calculated for easy visualisation.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
     ** Shell fluid velocities transfer
    U2=U1[-13,-6,]
   REGION(14,23,7,18,1,NZ)
    V2=V1[-13,-6,]
   REGION(14,23,7,18,1,NZ)
    Shell fluid velocities are transfered  from  where  they
    have been calculated for easy visualisation.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
    PLANTEND
  PHOTON USE
     p



     set prop off
     gr ou z 1 x 12 25 y 19 21
     gr ou z 1 x 12 25 y 4 18
     gr ou z 1 x 14 23 y 7 18
     gr ou z 1 x 12 25 y 22 31
     gr ou z 1 x 14 23 y 22 29
     *gr ou z 1 x 11 11 y 25 28
     *gr ou z 1 x 26 26 y 25 28
     gr z 1 x 11 11 y 20 21
     gr z 1 x 26 26 y 20 21
     *gr ou z 1 x 26 26 y 7 10
     *gr ou z 1 x 11 11 y 15 18
     gr ou z 1 x 14 16 y 1 3
     gr ou z 1 x 21 23 y 1 3
     gr ou z 1 x 11 11 y 16 17
     gr ou x 11 11 y 15 18
     gr ou z 1 x 26 26 y 8 9
     gr ou x 27 y 7 10
     gr ou z 1 x 11 11 y 26 27
     gr ou x 11 11 y 25 28
     gr ou z 1 x 26 26 y 26 27
     gr ou x 27 y 25 28
     gr ou x 15 15 y 19 21
     gr ou x 16 16 y 19 21
     gr ou x 17 17 y 19 21
     gr ou x 18 18 y 19 21
     gr ou x 19 19 y 19 21
     gr ou x 20 20 y 19 21
     gr ou x 21 21 y 19 21
     gr ou x 22 22 y 19 21
     gr ou x 23 23 y 19 21
     gr ou z 1 x 15 22 y 8 18
     gr ou z 1 x 16 21 y 9 18
     gr ou z 1 x 17 20 y 10 18
     gr ou z 1 x 18 19 y 11 18
     gr ou y 11 11 x 17 23
     gr ou y 15 15 x 14 20
     gr ou x 19 19 y 12 29
     set vec comp uu1 vv1 -
     vec z 1 sh
     set vec comp u2 v2 -
     vec z 1
     msg( Shell and tube fluid velocity vectors
     pause
     vec cl
     red
     set vec comp u1 v1 -
     vec z 1 x 11 m sh
     msg( Thermal displacement vectors in shell and tube bundle
  ENDUSE
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(In-Form equivalent of PLANT case Z604   )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 780
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 26
 XULAST =2.6
 XFRAC(1)=0.038462 ;XFRAC(2)=0.076923
 XFRAC(3)=0.115385 ;XFRAC(4)=0.153846
 XFRAC(5)=0.192308 ;XFRAC(6)=0.230769
 XFRAC(7)=0.269231 ;XFRAC(8)=0.307692
 XFRAC(9)=0.346154 ;XFRAC(10)=0.384615
 XFRAC(11)=0.423077 ;XFRAC(12)=0.461538
 XFRAC(13)=0.5 ;XFRAC(14)=0.538462
 XFRAC(15)=0.576923 ;XFRAC(16)=0.615385
 XFRAC(17)=0.653846 ;XFRAC(18)=0.692308
 XFRAC(19)=0.730769 ;XFRAC(20)=0.769231
 XFRAC(21)=0.807692 ;XFRAC(22)=0.846154
 XFRAC(23)=0.884615 ;XFRAC(24)=0.923077
 XFRAC(25)=0.961538 ;XFRAC(26)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 32
 YVLAST =3.2
 YFRAC(1)=0.03125 ;YFRAC(2)=0.0625
 YFRAC(3)=0.09375 ;YFRAC(4)=0.125
 YFRAC(5)=0.15625 ;YFRAC(6)=0.1875
 YFRAC(7)=0.21875 ;YFRAC(8)=0.25
 YFRAC(9)=0.28125 ;YFRAC(10)=0.3125
 YFRAC(11)=0.34375 ;YFRAC(12)=0.375
 YFRAC(13)=0.40625 ;YFRAC(14)=0.4375
 YFRAC(15)=0.46875 ;YFRAC(16)=0.5
 YFRAC(17)=0.53125 ;YFRAC(18)=0.5625
 YFRAC(19)=0.59375 ;YFRAC(20)=0.625
 YFRAC(21)=0.65625 ;YFRAC(22)=0.6875
 YFRAC(23)=0.71875 ;YFRAC(24)=0.75
 YFRAC(25)=0.78125 ;YFRAC(26)=0.8125
 YFRAC(27)=0.84375 ;YFRAC(28)=0.875
 YFRAC(29)=0.90625 ;YFRAC(30)=0.9375
 YFRAC(31)=0.96875 ;YFRAC(32)=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(1)=P1 ;NAME(3)=U1
 NAME(4)=U2 ;NAME(5)=V1
 NAME(6)=V2 ;NAME(14)=TEMP
 NAME(131)=EPSZ ;NAME(132)=VV1
 NAME(133)=UU1 ;NAME(134)=HTC
 NAME(135)=ALF3 ;NAME(136)=ALF2
 NAME(137)=EPST ;NAME(138)=STRX
 NAME(139)=EPSX ;NAME(140)=STRY
 NAME(141)=EPSY ;NAME(142)=LTLS
 NAME(143)=WDIS ;NAME(144)=TEM1
 NAME(145)=MARK ;NAME(146)=VISL
 NAME(147)=DEN1 ;NAME(148)=PRPS
 NAME(149)=NPOR ;NAME(150)=EPOR
    * 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,N,N,N,Y)
 SOLUTN(U1,Y,Y,N,N,N,Y)
 SOLUTN(U2,Y,N,N,N,N,Y)
 SOLUTN(V1,Y,Y,N,N,N,Y)
 SOLUTN(V2,Y,N,N,N,N,Y)
 SOLUTN(TEMP,Y,Y,N,N,N,Y)
 SOLUTN(EPSZ,Y,N,N,N,N,Y)
 SOLUTN(VV1,Y,N,N,N,N,Y)
 SOLUTN(UU1,Y,N,N,N,N,Y)
 SOLUTN(HTC,Y,N,N,N,N,Y)
 SOLUTN(ALF3,Y,N,N,N,N,Y)
 SOLUTN(ALF2,Y,N,N,N,N,Y)
 SOLUTN(EPST,Y,N,N,N,N,Y)
 SOLUTN(STRX,Y,N,N,N,N,Y)
 SOLUTN(EPSX,Y,N,N,N,N,Y)
 SOLUTN(STRY,Y,N,N,N,N,Y)
 SOLUTN(EPSY,Y,N,N,N,N,Y)
 SOLUTN(LTLS,Y,Y,Y,N,N,Y)
 SOLUTN(WDIS,Y,N,N,N,N,N)
 SOLUTN(TEM1,Y,Y,Y,N,N,Y)
 SOLUTN(MARK,Y,N,N,N,N,Y)
 SOLUTN(VISL,Y,N,N,N,N,Y)
 SOLUTN(DEN1,Y,N,N,N,N,Y)
 SOLUTN(PRPS,Y,N,N,N,N,Y)
 SOLUTN(NPOR,Y,N,N,N,N,Y)
 SOLUTN(EPOR,Y,N,N,N,N,Y)
 DEN1 = 147
 VISL = 146
 EPOR = 150 ;HPOR = 0 ;NPOR = 149 ;VPOR = 0
 PRPS = 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(U1,Y,Y,Y,Y,Y,Y)
 TERMS(V1,Y,Y,Y,Y,Y,Y)
 TERMS(TEMP,N,Y,Y,N,Y,N)
 TERMS(LTLS,N,N,Y,N,Y,Y)
 TERMS(TEM1,N,Y,Y,Y,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 ISOLX = 0 ;ISOLY = 0 ;ISOLZ = 0
 ************************************************************
  Group 9. Properties used if PRPS is not
  stored, and where PRPS = -1.0 if it is!
 RHO1 =1. ;TMP1 =0.
 EL1 =0.
 TSURR =0. ;TEMP0 =0.
 PRESS0 =0.
 DVO1DT =1.0E-03 ;DRH1DP =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 ENUL =1. ;ENUT =0.
 PRNDTL(U1)=1. ;PRNDTL(V1)=1.
 PRNDTL(TEMP)=0.702 ;PRNDTL(LTLS)=1.
 PRNDTL(TEM1)=1.
 PRT(U1)=1. ;PRT(V1)=1.
 PRT(TEMP)=1. ;PRT(LTLS)=1.
 PRT(TEM1)=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.
 FIINIT(U2)=1.0E-10 ;FIINIT(V1)=1.0E-10
 FIINIT(V2)=1.0E-10 ;FIINIT(TEMP)=1.0E-10
 FIINIT(EPSZ)=1.0E-10 ;FIINIT(VV1)=1.0E-10
 FIINIT(UU1)=1.0E-10 ;FIINIT(HTC)=0.
 FIINIT(ALF3)=0. ;FIINIT(ALF2)=0.
 FIINIT(EPST)=0. ;FIINIT(STRX)=0.
 FIINIT(EPSX)=0. ;FIINIT(STRY)=0.
 FIINIT(EPSY)=0. ;FIINIT(LTLS)=1.0E-10
 FIINIT(WDIS)=0.1 ;FIINIT(TEM1)=0.
 FIINIT(MARK)=1.0E-10 ;FIINIT(VISL)=1.0E-10
 FIINIT(DEN1)=1.0E-10 ;FIINIT(PRPS)=111.
 FIINIT(NPOR)=1. ;FIINIT(EPOR)=1.
 
 PATCH(BACGRND ,INIVAL, 12, 25, 1, 31, 1, 1, 1, 1)
 INIT(BACGRND ,PRPS,0. ,111. )
 
 PATCH(HEADSIDE,INIVAL, 1, 10, 25, 32, 1, 1, 1, 1)
 INIT(HEADSIDE,MARK,0. ,1. )
 INIT(HEADSIDE,PRPS,0. ,0. )
 
 PATCH(TUBESIDE,INIVAL, 1, 10, 13, 24, 1, 1, 1, 1)
 INIT(TUBESIDE,MARK,0. ,2. )
 INIT(TUBESIDE,PRPS,0. ,0. )
 
 PATCH(SHELSIDE,INIVAL, 1, 10, 1, 12, 1, 1, 1, 1)
 INIT(SHELSIDE,MARK,0. ,3. )
 INIT(SHELSIDE,PRPS,0. ,0. )
 
 PATCH(TUBSHEET,INIVAL, 11, 26, 19, 21, 1, 1, 1, 1)
 INIT(TUBSHEET,PRPS,0. ,103. )
 
 PATCH(ADDHEAD ,INIVAL, 14, 23, 22, 29, 1, 1, 1, 1)
 INIT(ADDHEAD ,MARK,0. ,4. )
 INIT(ADDHEAD ,PRPS,0. ,0. )
 
 PATCH(ADDSHELL,INIVAL, 14, 23, 7, 18, 1, 1, 1, 1)
 INIT(ADDSHELL,MARK,0. ,5. )
 INIT(ADDSHELL,PRPS,0. ,0. )
 
 PATCH(TUBEWALL,INIVAL, 16, 21, 9, 18, 1, 1, 1, 1)
 INIT(TUBEWALL,MARK,0. ,100. )
 INIT(TUBEWALL,PRPS,0. ,100. )
 
 PATCH(SUPOR1 ,INIVAL, 11, 13, 1, 3, 1, 1, 1, 1)
 INIT(SUPOR1 ,PRPS,0. ,0. )
 
 PATCH(SUPOR2 ,INIVAL, 17, 20, 1, 3, 1, 1, 1, 1)
 INIT(SUPOR2 ,PRPS,0. ,0. )
 
 PATCH(SUPOR3 ,INIVAL, 24, 26, 1, 3, 1, 1, 1, 1)
 INIT(SUPOR3 ,PRPS,0. ,0. )
 
 PATCH(WESLAYER,INIVAL, 11, 11, 1, 32, 1, 1, 1, 1)
 INIT(WESLAYER,PRPS,0. ,0. )
 
 PATCH(EASLAYER,INIVAL, 26, 26, 1, 32, 1, 1, 1, 1)
 INIT(EASLAYER,PRPS,0. ,0. )
 
 PATCH(NORLAYER,INIVAL, 11, 26, 32, 32, 1, 1, 1, 1)
 INIT(NORLAYER,PRPS,0. ,0. )
 
 PATCH(CMP14 ,INIVAL, 1, 10, 12, 12, 1, 1, 1, 1)
 INIT(CMP14 ,V1 ,0. ,0. )
 INIT(CMP14 ,V2 ,0. ,0. )
 INIT(CMP14 ,NPOR,0. ,0. )
 
 PATCH(CMP15 ,INIVAL, 10, 10, 1, 32, 1, 1, 1, 1)
 INIT(CMP15 ,U1 ,0. ,0. )
 INIT(CMP15 ,U2 ,0. ,0. )
 INIT(CMP15 ,EPOR,0. ,0. )
 
 PATCH(CMP16 ,INIVAL, 5, 5, 18, 32, 1, 1, 1, 1)
 INIT(CMP16 ,U1 ,0. ,0. )
 INIT(CMP16 ,U2 ,0. ,0. )
 INIT(CMP16 ,EPOR,0. ,0. )
 
 PATCH(CMP17 ,INIVAL, 1, 7, 8, 8, 1, 1, 1, 1)
 INIT(CMP17 ,V1 ,0. ,0. )
 INIT(CMP17 ,V2 ,0. ,0. )
 INIT(CMP17 ,NPOR,0. ,0. )
 
 PATCH(CMP18 ,INIVAL, 4, 10, 4, 4, 1, 1, 1, 1)
 INIT(CMP18 ,V1 ,0. ,0. )
 INIT(CMP18 ,V2 ,0. ,0. )
 INIT(CMP18 ,NPOR,0. ,0. )
 
 PATCH(CMP19 ,INIVAL, 2, 9, 13, 24, 1, 1, 1, 1)
 INIT(CMP19 ,NPOR,0. ,0.5 )
 
 PATCH(CMP20 ,INIVAL, 2, 9, 14, 24, 1, 1, 1, 1)
 INIT(CMP20 ,EPOR,0. ,0.5 )
 
 PATCH(CMP21 ,INIVAL, 5, 6, 16, 16, 1, 1, 1, 1)
 INIT(CMP21 ,V1 ,0. ,0. )
 INIT(CMP21 ,V2 ,0. ,0. )
 INIT(CMP21 ,NPOR,0. ,0. )
 
 PATCH(CMP22 ,INIVAL, 4, 7, 15, 15, 1, 1, 1, 1)
 INIT(CMP22 ,V1 ,0. ,0. )
 INIT(CMP22 ,V2 ,0. ,0. )
 INIT(CMP22 ,NPOR,0. ,0. )
 
 PATCH(CMP23 ,INIVAL, 3, 8, 14, 14, 1, 1, 1, 1)
 INIT(CMP23 ,V1 ,0. ,0. )
 INIT(CMP23 ,V2 ,0. ,0. )
 INIT(CMP23 ,NPOR,0. ,0. )
 
 PATCH(CMP24 ,INIVAL, 2, 9, 13, 13, 1, 1, 1, 1)
 INIT(CMP24 ,V1 ,0. ,0. )
 INIT(CMP24 ,V2 ,0. ,0. )
 INIT(CMP24 ,NPOR,0. ,0. )
 
 PATCH(CMP25 ,INIVAL, 1, 1, 14, 24, 1, 1, 1, 1)
 INIT(CMP25 ,U1 ,0. ,0. )
 INIT(CMP25 ,U2 ,0. ,0. )
 INIT(CMP25 ,EPOR,0. ,0. )
 
 PATCH(CMP26 ,INIVAL, 2, 2, 15, 24, 1, 1, 1, 1)
 INIT(CMP26 ,U1 ,0. ,0. )
 INIT(CMP26 ,U2 ,0. ,0. )
 INIT(CMP26 ,EPOR,0. ,0. )
 
 PATCH(CMP27 ,INIVAL, 3, 3, 16, 24, 1, 1, 1, 1)
 INIT(CMP27 ,U1 ,0. ,0. )
 INIT(CMP27 ,U2 ,0. ,0. )
 INIT(CMP27 ,EPOR,0. ,0. )
 
 PATCH(CMP28 ,INIVAL, 4, 4, 17, 24, 1, 1, 1, 1)
 INIT(CMP28 ,U1 ,0. ,0. )
 INIT(CMP28 ,U2 ,0. ,0. )
 INIT(CMP28 ,EPOR,0. ,0. )
 
 PATCH(CMP29 ,INIVAL, 5, 5, 18, 24, 1, 1, 1, 1)
 INIT(CMP29 ,U1 ,0. ,0. )
 INIT(CMP29 ,U2 ,0. ,0. )
 INIT(CMP29 ,EPOR,0. ,0. )
 
 PATCH(CMP30 ,INIVAL, 6, 6, 17, 24, 1, 1, 1, 1)
 INIT(CMP30 ,U1 ,0. ,0. )
 INIT(CMP30 ,U2 ,0. ,0. )
 INIT(CMP30 ,EPOR,0. ,0. )
 
 PATCH(CMP31 ,INIVAL, 7, 7, 16, 24, 1, 1, 1, 1)
 INIT(CMP31 ,U1 ,0. ,0. )
 INIT(CMP31 ,U2 ,0. ,0. )
 INIT(CMP31 ,EPOR,0. ,0. )
 
 PATCH(CMP32 ,INIVAL, 8, 8, 15, 24, 1, 1, 1, 1)
 INIT(CMP32 ,U1 ,0. ,0. )
 INIT(CMP32 ,U2 ,0. ,0. )
 INIT(CMP32 ,EPOR,0. ,0. )
 
 PATCH(CMP33 ,INIVAL, 9, 9, 14, 24, 1, 1, 1, 1)
 INIT(CMP33 ,U1 ,0. ,0. )
 INIT(CMP33 ,U2 ,0. ,0. )
 INIT(CMP33 ,EPOR,0. ,0. )
 
 PATCH(NL1N ,INIVAL, 1, 7, 8, 8, 1, 1, 1, 1)
 INIT(NL1N ,V1 ,0. ,0. )
 INIT(NL1N ,V2 ,0. ,0. )
 INIT(NL1N ,LTLS,1. ,0. )
 INIT(NL1N ,NPOR,0. ,0. )
 
 PATCH(NL1S ,INIVAL, 1, 7, 8, 8, 1, 1, 1, 1)
 INIT(NL1S ,V1 ,0. ,0. )
 INIT(NL1S ,V2 ,0. ,0. )
 INIT(NL1S ,LTLS,1. ,0. )
 INIT(NL1S ,NPOR,0. ,0. )
 
 PATCH(NL2N ,INIVAL, 4, 10, 4, 4, 1, 1, 1, 1)
 INIT(NL2N ,V1 ,0. ,0. )
 INIT(NL2N ,V2 ,0. ,0. )
 INIT(NL2N ,LTLS,1. ,0. )
 INIT(NL2N ,NPOR,0. ,0. )
 
 PATCH(NL2S ,INIVAL, 4, 10, 4, 4, 1, 1, 1, 1)
 INIT(NL2S ,V1 ,0. ,0. )
 INIT(NL2S ,V2 ,0. ,0. )
 INIT(NL2S ,LTLS,1. ,0. )
 INIT(NL2S ,NPOR,0. ,0. )
 
 PATCH(NL3N ,INIVAL, 1, 10, 12, 12, 1, 1, 1, 1)
 INIT(NL3N ,V1 ,0. ,0. )
 INIT(NL3N ,V2 ,0. ,0. )
 INIT(NL3N ,LTLS,1. ,0. )
 INIT(NL3N ,NPOR,0. ,0. )
 
 PATCH(NL4E1 ,INIVAL, 10, 10, 4, 12, 1, 1, 1, 1)
 INIT(NL4E1 ,U1 ,0. ,0. )
 INIT(NL4E1 ,U2 ,0. ,0. )
 INIT(NL4E1 ,LTLS,1. ,0. )
 INIT(NL4E1 ,EPOR,0. ,0. )
 
 PATCH(NL5W3 ,INIVAL, 5, 5, 25, 32, 1, 1, 1, 1)
 INIT(NL5W3 ,U1 ,0. ,0. )
 INIT(NL5W3 ,U2 ,0. ,0. )
 INIT(NL5W3 ,LTLS,1. ,0. )
 INIT(NL5W3 ,EPOR,0. ,0. )
 
 PATCH(NL5E1 ,INIVAL, 10, 10, 25, 28, 1, 1, 1, 1)
 INIT(NL5E1 ,U1 ,0. ,0. )
 INIT(NL5E1 ,U2 ,0. ,0. )
 INIT(NL5E1 ,LTLS,1. ,0. )
 INIT(NL5E1 ,EPOR,0. ,0. )
 
 PATCH(NL5E2 ,INIVAL, 10, 10, 31, 32, 1, 1, 1, 1)
 INIT(NL5E2 ,U1 ,0. ,0. )
 INIT(NL5E2 ,U2 ,0. ,0. )
 INIT(NL5E2 ,LTLS,1. ,0. )
 INIT(NL5E2 ,EPOR,0. ,0. )
 
 PATCH(NL5E3 ,INIVAL, 5, 5, 25, 32, 1, 1, 1, 1)
 INIT(NL5E3 ,U1 ,0. ,0. )
 INIT(NL5E3 ,U2 ,0. ,0. )
 INIT(NL5E3 ,LTLS,1. ,0. )
 INIT(NL5E3 ,EPOR,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(NL1N-NW ,NWALL , 1, 7, 8, 8, 1, 1, 1, 1)
 COVAL(NL1N-NW ,U1 ,1. ,0. )
 
 PATCH(NL1N-SW ,SWALL , 1, 7, 9, 9, 1, 1, 1, 1)
 COVAL(NL1N-SW ,U1 ,1. ,0. )
 
 PATCH(NL1S-NW ,NWALL , 1, 7, 8, 8, 1, 1, 1, 1)
 COVAL(NL1S-NW ,U1 ,1. ,0. )
 
 PATCH(NL1S-SW ,SWALL , 1, 7, 9, 9, 1, 1, 1, 1)
 COVAL(NL1S-SW ,U1 ,1. ,0. )
 
 PATCH(NL2N-NW ,NWALL , 4, 10, 4, 4, 1, 1, 1, 1)
 COVAL(NL2N-NW ,U1 ,1. ,0. )
 
 PATCH(NL2N-SW ,SWALL , 4, 10, 5, 5, 1, 1, 1, 1)
 COVAL(NL2N-SW ,U1 ,1. ,0. )
 
 PATCH(NL2S-NW ,NWALL , 4, 10, 4, 4, 1, 1, 1, 1)
 COVAL(NL2S-NW ,U1 ,1. ,0. )
 
 PATCH(NL2S-SW ,SWALL , 4, 10, 5, 5, 1, 1, 1, 1)
 COVAL(NL2S-SW ,U1 ,1. ,0. )
 
 PATCH(NL3N-NW ,NWALL , 1, 10, 12, 12, 1, 1, 1, 1)
 COVAL(NL3N-NW ,U1 ,1. ,0. )
 
 PATCH(NL3N-SW ,SWALL , 1, 10, 13, 13, 1, 1, 1, 1)
 COVAL(NL3N-SW ,U1 ,1. ,0. )
 
 PATCH(NL3S ,SWALL , 1, 10, 1, 1, 1, 1, 1, 1)
 COVAL(NL3S ,U1 ,1. ,0. )
 COVAL(NL3S ,LTLS,1. ,0. )
 
 PATCH(NL4W1 ,WWALL , 1, 1, 1, 9, 1, 1, 1, 1)
 COVAL(NL4W1 ,V1 ,1. ,0. )
 COVAL(NL4W1 ,LTLS,1. ,0. )
 
 PATCH(NL4W2 ,WWALL , 1, 1, 12, 12, 1, 1, 1, 1)
 COVAL(NL4W2 ,V1 ,1. ,0. )
 COVAL(NL4W2 ,LTLS,1. ,0. )
 
 PATCH(NL4E1-EW,EWALL , 10, 10, 4, 12, 1, 1, 1, 1)
 COVAL(NL4E1-EW,V1 ,1. ,0. )
 
 PATCH(NL4E1-WW,WWALL , 11, 11, 4, 12, 1, 1, 1, 1)
 COVAL(NL4E1-WW,V1 ,1. ,0. )
 
 PATCH(NL5N1 ,NWALL , 1, 10, 32, 32, 1, 1, 1, 1)
 COVAL(NL5N1 ,U1 ,1. ,0. )
 COVAL(NL5N1 ,LTLS,1. ,0. )
 
 PATCH(NL5W1 ,WWALL , 1, 1, 25, 28, 1, 1, 1, 1)
 COVAL(NL5W1 ,V1 ,1. ,0. )
 COVAL(NL5W1 ,LTLS,1. ,0. )
 
 PATCH(NL5W2 ,WWALL , 1, 1, 31, 32, 1, 1, 1, 1)
 COVAL(NL5W2 ,V1 ,1. ,0. )
 COVAL(NL5W2 ,LTLS,1. ,0. )
 
 PATCH(NL5W3-EW,EWALL , 5, 5, 25, 32, 1, 1, 1, 1)
 COVAL(NL5W3-EW,V1 ,1. ,0. )
 
 PATCH(NL5W3-WW,WWALL , 6, 6, 25, 32, 1, 1, 1, 1)
 COVAL(NL5W3-WW,V1 ,1. ,0. )
 
 PATCH(NL5E1-EW,EWALL , 10, 10, 25, 28, 1, 1, 1, 1)
 COVAL(NL5E1-EW,V1 ,1. ,0. )
 
 PATCH(NL5E1-WW,WWALL , 11, 11, 25, 28, 1, 1, 1, 1)
 COVAL(NL5E1-WW,V1 ,1. ,0. )
 
 PATCH(NL5E2-EW,EWALL , 10, 10, 31, 32, 1, 1, 1, 1)
 COVAL(NL5E2-EW,V1 ,1. ,0. )
 
 PATCH(NL5E2-WW,WWALL , 11, 11, 31, 32, 1, 1, 1, 1)
 COVAL(NL5E2-WW,V1 ,1. ,0. )
 
 PATCH(NL5E3-EW,EWALL , 5, 5, 25, 32, 1, 1, 1, 1)
 COVAL(NL5E3-EW,V1 ,1. ,0. )
 
 PATCH(NL5E3-WW,WWALL , 6, 6, 25, 32, 1, 1, 1, 1)
 COVAL(NL5E3-WW,V1 ,1. ,0. )
 
 PATCH(FIXL ,CELL , 1, 10, 13, 24, 1, 1, 1, 1)
 COVAL(FIXL ,LTLS, FIXVAL ,0. )
 
 PATCH(INTUBE ,WEST , 1, 1, 29, 30, 1, 1, 1, 1)
 COVAL(INTUBE ,P1 , FIXFLU ,1000. )
 COVAL(INTUBE ,U1 ,0. ,1. )
 COVAL(INTUBE ,TEMP,0. ,0. )
 
 PATCH(OUTUBE ,EAST , 10, 10, 29, 30, 1, 1, 1, 1)
 COVAL(OUTUBE ,P1 ,1000. ,0. )
 
 PATCH(INSHEL ,EAST , 10, 10, 2, 3, 1, 1, 1, 1)
 COVAL(INSHEL ,P1 , FIXFLU ,1.2 )
 COVAL(INSHEL ,U1 ,0. ,-1. )
 COVAL(INSHEL ,TEMP,0. ,1. )
 
 PATCH(OUSHEL ,WEST , 1, 1, 10, 11, 1, 1, 1, 1)
 COVAL(OUSHEL ,P1 ,1000. ,0. )
 
 PATCH(BASEW ,NORTH , 12, 13, 3, 3, 1, 1, 1, 1)
 COVAL(BASEW ,V1 , FIXFLU ,0. )
 
 PATCH(BASEM ,NORTH , 17, 20, 3, 3, 1, 1, 1, 1)
 COVAL(BASEM ,V1 , FIXFLU ,0. )
 
 PATCH(BASEE ,NORTH , 24, 25, 3, 3, 1, 1, 1, 1)
 COVAL(BASEE ,V1 , FIXFLU ,0. )
 
 PATCH(SUP1W ,EAST , 13, 13, 1, 3, 1, 1, 1, 1)
 COVAL(SUP1W ,U1 , FIXFLU ,0. )
 
 PATCH(SUP1E ,EAST , 16, 16, 1, 3, 1, 1, 1, 1)
 COVAL(SUP1E ,U1 , FIXFLU ,0. )
 
 PATCH(SUP2W ,EAST , 20, 20, 1, 3, 1, 1, 1, 1)
 COVAL(SUP2W ,U1 , FIXFLU ,0. )
 
 PATCH(SUP2E ,EAST , 23, 23, 1, 3, 1, 1, 1, 1)
 COVAL(SUP2E ,U1 , FIXFLU ,0. )
 
 PATCH(OUWW ,EAST , 11, 11, 4, 31, 1, 1, 1, 1)
 COVAL(OUWW ,U1 , FIXFLU ,0. )
 
 PATCH(OUWE ,EAST , 25, 25, 4, 31, 1, 1, 1, 1)
 COVAL(OUWE ,U1 , FIXFLU ,0. )
 
 PATCH(OUWN ,NORTH , 12, 25, 31, 31, 1, 1, 1, 1)
 COVAL(OUWN ,V1 , FIXFLU ,0. )
 
 PATCH(HEDWN ,NORTH , 14, 23, 29, 29, 1, 1, 1, 1)
 COVAL(HEDWN ,V1 , FIXFLU ,0. )
 
 PATCH(HEDWS ,NORTH , 14, 23, 21, 21, 1, 1, 1, 1)
 COVAL(HEDWS ,V1 , FIXFLU ,0. )
 
 PATCH(HEDWW ,EAST , 13, 13, 22, 29, 1, 1, 1, 1)
 COVAL(HEDWW ,U1 , FIXFLU ,0. )
 
 PATCH(HEDWE ,EAST , 23, 23, 22, 29, 1, 1, 1, 1)
 COVAL(HEDWE ,U1 , FIXFLU ,0. )
 
 PATCH(SHWNW ,NORTH , 14, 15, 18, 18, 1, 1, 1, 1)
 COVAL(SHWNW ,V1 , FIXFLU ,0. )
 
 PATCH(SHWNE ,NORTH , 22, 23, 18, 18, 1, 1, 1, 1)
 COVAL(SHWNE ,V1 , FIXFLU ,0. )
 
 PATCH(SHWS ,NORTH , 14, 23, 6, 6, 1, 1, 1, 1)
 COVAL(SHWS ,V1 , FIXFLU ,0. )
 
 PATCH(SHWW ,EAST , 13, 13, 7, 18, 1, 1, 1, 1)
 COVAL(SHWW ,U1 , FIXFLU ,0. )
 
 PATCH(SHWE ,EAST , 23, 23, 7, 18, 1, 1, 1, 1)
 COVAL(SHWE ,U1 , FIXFLU ,0. )
 
 PATCH(TUBS ,NORTH , 16, 21, 8, 8, 1, 1, 1, 1)
 COVAL(TUBS ,V1 , FIXFLU ,0. )
 
 PATCH(TUBW ,EAST , 15, 15, 9, 18, 1, 1, 1, 1)
 COVAL(TUBW ,U1 , FIXFLU ,0. )
 
 PATCH(TUBE ,EAST , 21, 21, 9, 18, 1, 1, 1, 1)
 COVAL(TUBE ,U1 , FIXFLU ,0. )
 
 PATCH(FIXW ,EAST , 11, 11, 19, 21, 1, 1, 1, 1)
 COVAL(FIXW ,U1 , FIXVAL ,0. )
 
 PATCH(FIXE ,EAST , 25, 25, 19, 21, 1, 1, 1, 1)
 COVAL(FIXE ,U1 , FIXVAL ,0. )
 
 PATCH(HEADFLOW,PHASEM, 1, 10, 25, 32, 1, 1, 1, 1)
 
 PATCH(TBSDFLOW,PHASEM, 1, 10, 13, 24, 1, 1, 1, 1)
 COVAL(TBSDFLOW,U1 ,In-Form:source - see Grp 19)
 COVAL(TBSDFLOW,V1 ,In-Form:source - see Grp 19)
 COVAL(TBSDFLOW,TEMP,In-Form:source - see Grp 19)
 
 PATCH(SHSDFLOW,PHASEM, 1, 10, 1, 12, 1, 1, 1, 1)
 COVAL(SHSDFLOW,U1 ,In-Form:source - see Grp 19)
 COVAL(SHSDFLOW,V1 ,In-Form:source - see Grp 19)
 COVAL(SHSDFLOW,TEMP,In-Form:source - see Grp 19)
 
 PATCH(HEADTHRM,CELL , 14, 23, 22, 29, 1, 1, 1, 1)
 COVAL(HEADTHRM,TEM1,In-Form:source - see Grp 19)
 
 PATCH(SHSDTHRM,CELL , 14, 23, 7, 18, 1, 1, 1, 1)
 COVAL(SHSDTHRM,TEM1,In-Form:source - see Grp 19)
 
 PATCH(TBWLTHRM,CELL , 16, 21, 9, 18, 1, 1, 1, 1)
 COVAL(TBWLTHRM,TEM1,In-Form:source - see Grp 19)
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 400 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-08
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(U1)=10
 LITER(V1)=10 ;LITER(TEMP)=20
 LITER(LTLS)=20 ;LITER(TEM1)=20
 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
 ENDIT(V1)=1.0E-03 ;ENDIT(TEMP)=1.0E-03
 ENDIT(LTLS)=1.0E-03 ;ENDIT(TEM1)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,0.25)
 RELAX(U1,FALSDT,0.075)
 RELAX(U2,LINRLX,1.)
 RELAX(V1,FALSDT,0.075)
 RELAX(V2,LINRLX,1.)
 RELAX(TEMP,FALSDT,3.)
 RELAX(EPSZ,LINRLX,1.)
 RELAX(VV1,LINRLX,1.)
 RELAX(UU1,LINRLX,1.)
 RELAX(HTC,LINRLX,1.)
 RELAX(ALF3,LINRLX,1.)
 RELAX(ALF2,LINRLX,1.)
 RELAX(EPST,LINRLX,1.)
 RELAX(STRX,LINRLX,1.)
 RELAX(EPSX,LINRLX,1.)
 RELAX(STRY,LINRLX,1.)
 RELAX(EPSY,LINRLX,1.)
 RELAX(LTLS,FALSDT,1.0E+09)
 RELAX(WDIS,LINRLX,1.)
 RELAX(TEM1,FALSDT,1.0E+04)
 RELAX(MARK,LINRLX,1.)
 RELAX(VISL,LINRLX,1.)
 RELAX(DEN1,LINRLX,1.)
 RELAX(PRPS,LINRLX,1.)
 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(U2)=1.0E+06 ;VARMIN(U2)=-1.0E+06
 VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
 VARMAX(V2)=1.0E+06 ;VARMIN(V2)=-1.0E+06
 VARMAX(TEMP)=1.0E+10 ;VARMIN(TEMP)=-1.0E+10
 VARMAX(EPSZ)=1.0E+10 ;VARMIN(EPSZ)=-1.0E+10
 VARMAX(VV1)=1.0E+10 ;VARMIN(VV1)=-1.0E+10
 VARMAX(UU1)=1.0E+10 ;VARMIN(UU1)=-1.0E+10
 VARMAX(HTC)=1.0E+10 ;VARMIN(HTC)=-1.0E+10
 VARMAX(ALF3)=1.0E+10 ;VARMIN(ALF3)=-1.0E+10
 VARMAX(ALF2)=1.0E+10 ;VARMIN(ALF2)=-1.0E+10
 VARMAX(EPST)=1.0E+10 ;VARMIN(EPST)=-1.0E+10
 VARMAX(STRX)=1.0E+10 ;VARMIN(STRX)=-1.0E+10
 VARMAX(EPSX)=1.0E+10 ;VARMIN(EPSX)=-1.0E+10
 VARMAX(STRY)=1.0E+10 ;VARMIN(STRY)=-1.0E+10
 VARMAX(EPSY)=1.0E+10 ;VARMIN(EPSY)=-1.0E+10
 VARMAX(LTLS)=1.0E+10 ;VARMIN(LTLS)=-1.0E+10
 VARMAX(WDIS)=1.0E+10 ;VARMIN(WDIS)=-1.0E+10
 VARMAX(TEM1)=1.0E+10 ;VARMIN(TEM1)=-1.0E+10
 VARMAX(MARK)=1.0E+10 ;VARMIN(MARK)=-1.0E+10
 VARMAX(VISL)=1.0E+10 ;VARMIN(VISL)=-1.0E+10
 VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=-1.0E+10
 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
 VARMAX(NPOR)=1.0E+10 ;VARMIN(NPOR)=-1.0E+10
 VARMAX(EPOR)=1.0E+10 ;VARMIN(EPOR)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 STRA = T
 PARSOL = F
 ISG62 = 1
 POISSN =0.3333
 SPEDAT(SET,STRAIN,CALSTR,L,T)
 SPEDAT(SET,STRAIN,POISSN,R,0.3333)
 SPEDAT(SET,STRAIN,EXCOLI,R,1.0E-03)
 SPEDAT(SET,STRAIN,EXCOC1,R,0.)
 SPEDAT(SET,STRAIN,EXCOC2,R,0.)
 SPEDAT(SET,STRAIN,STIFFN,R,2.0E+11)
 SPEDAT(SET,STRAIN,STIFC1,R,0.)
 SPEDAT(SET,STRAIN,STIFC2,R,0.)
 SPEDAT(SET,PROPERTY,RHO1!HEADFLOW,C,=1000.)
 SPEDAT(SET,PROPERTY,RHO1!TBSDFLOW,C,=1000.)
 SPEDAT(SET,PROPERTY,RHO1!SHSDFLOW,C,=1.2)
 SPEDAT(SET,PROPERTY,ENUL!HEADFLOW,C,=1.*SQRT(U1^2+V1^2)*WDIS)
 SPEDAT(SET,PROPERTY,ENUL!TBSDFLOW,C,=0.01)
 SPEDAT(SET,PROPERTY,ENUL!SHSDFLOW,C,=1.*SQRT(U1^2+V1^2)*WDIS)
 SPEDAT(SET,SOURCE,U1!TBSDFLOW,C,=-0.2*(U1^2+V1^2)^0.15*U1)
 SPEDAT(SET,SOURCE,V1!TBSDFLOW,C,=-0.2*(U1^2+V1^2)^0.15*V1)
 SPEDAT(SET,SOURCE,U1!SHSDFLOW,C,=-U1*2.2*(U1^2+V1^2)^0.25)
 SPEDAT(SET,SOURCE,V1!SHSDFLOW,C,=-V1*2.2*(U1^2+V1^2)^0.25)
 SPEDAT(SET,STORED,ALF2!TBSDFLOW,C,=1.+1.*SQRT(U1^2+V1^2+1.0E-20)!$)
 SPEDAT(SET,STORED,ALF2!TBSDFLOW,C,ZSLFIN)
 SPEDAT(SET,STORED,ALF3!SHSDFLOW,C,=1.+3.*SQRT(U1^2+V1^2+1.0E-20)!$)
 SPEDAT(SET,STORED,ALF3!SHSDFLOW,C,ZSLFIN)
 SPEDAT(SET,STORED,HTC!TBSDFLOW,C,=1./(1/ALF2+1/ALF3[&-12&])!ZSLFIN)
 SPEDAT(SET,STORED,HTC!SHSDFLOW,C,=1./(1/ALF3+1/ALF2[&+12&])!ZSLFIN)
 SPEDAT(SET,SOURCE,TEMP!TBSDFLOW,C,=HTC*(TEMP[&-12&]-TEMP)!LINE)
 SPEDAT(SET,SOURCE,TEMP!SHSDFLOW,C,=HTC*(TEMP[&+12&]-TEMP)!LINE)
 SPEDAT(SET,SOURCE,TEM1!SHSDTHRM,C,=TEMP[-13&-6]!FIXVAL)
 SPEDAT(SET,SOURCE,TEM1!TBWLTHRM,C,=(ALF2[-13&+6]*TEMP[-13&+6]+ALF$)
 SPEDAT(SET,SOURCE,TEM1!TBWLTHRM,C,3[-13&-6]*TEMP[-13&-6])/(ALF2[-$)
 SPEDAT(SET,SOURCE,TEM1!TBWLTHRM,C,13&+6]+ALF3[-13&-6]+1.0E-20)!FI$)
 SPEDAT(SET,SOURCE,TEM1!TBWLTHRM,C,XVAL!IMAT=100)
 SPEDAT(SET,SOURCE,TEM1!HEADTHRM,C,=TEMP[-13&+3]!FIXVAL)
 SPEDAT(SET,STORED,UU1!HEADTHRM,C,=U1[-13&+3&]!ZSLFIN)
 SPEDAT(SET,STORED,VV1!HEADTHRM,C,=V1[-13&+3&]!ZSLFIN)
 SPEDAT(SET,STORED,UU1!SHSDTHRM,C,=U1[-13&+6&]!ZSLFIN)
 SPEDAT(SET,STORED,UU1!TBWLTHRM,C,=U1[-13&+6&]!IMAT=100!ZSLFIN)
 SPEDAT(SET,STORED,VV1!SHSDTHRM,C,=V1[-13&+6&]!ZSLFIN)
 SPEDAT(SET,STORED,VV1!TBWLTHRM,C,=V1[-13&+6&]!IMAT=100!ZSLFIN)
 SPEDAT(SET,STORED,U2!SHSDTHRM,C,=U1[-13&-6&]!ZSLFIN)
 SPEDAT(SET,STORED,U2!TBWLTHRM,C,=U1[-13&-6&]!IMAT=100!ZSLFIN)
 SPEDAT(SET,STORED,V2!SHSDTHRM,C,=V1[-13&-6&]!ZSLFIN)
 SPEDAT(SET,STORED,V2!TBWLTHRM,C,=V1[-13&-6&]!IMAT=100!ZSLFIN)
 SPEDAT(SET,MATERIAL,111,L,T)
 SPEDAT(SET,MATERIAL,0,L,T)
 SPEDAT(SET,MATERIAL,103,L,T)
 SPEDAT(SET,MATERIAL,100,L,T)
 IG( 1) = 12
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=732.5 ;EX(U1)=0.09472
 EX(U2)=0.0573 ;EX(V1)=0.1366
 EX(V2)=0.02435 ;EX(TEMP)=0.1982
 EX(EPSZ)=0. ;EX(VV1)=0.1086
 EX(UU1)=0.03593 ;EX(HTC)=0.2839
 EX(ALF3)=0.3517 ;EX(ALF2)=0.244
 EX(EPST)=1.679E-06 ;EX(STRX)=1.086E+05
 EX(EPSX)=2.141E-06 ;EX(STRY)=2.916E+05
 EX(EPSY)=3.287E-06 ;EX(LTLS)=0.01575
 EX(WDIS)=0.05966 ;EX(TEM1)=0.2956
 EX(MARK)=8.774 ;EX(VISL)=0.01852
 EX(DEN1)=2463. ;EX(PRPS)=34.82
 EX(NPOR)=0.9014 ;EX(EPOR)=0.8425
 ************************************************************
  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(U1,Y,N,Y,Y,Y,Y)
 OUTPUT(U2,Y,N,Y,N,N,N)
 OUTPUT(V1,Y,N,Y,Y,Y,Y)
 OUTPUT(V2,Y,N,Y,N,N,N)
 OUTPUT(TEMP,Y,N,Y,Y,Y,Y)
 OUTPUT(EPSZ,Y,N,Y,N,N,N)
 OUTPUT(VV1,Y,N,Y,N,N,N)
 OUTPUT(UU1,Y,N,Y,N,N,N)
 OUTPUT(HTC,Y,N,Y,N,N,N)
 OUTPUT(ALF3,Y,N,Y,N,N,N)
 OUTPUT(ALF2,Y,N,Y,N,N,N)
 OUTPUT(EPST,Y,N,N,N,N,N)
 OUTPUT(STRX,Y,N,N,N,N,N)
 OUTPUT(EPSX,Y,N,N,N,N,N)
 OUTPUT(STRY,Y,N,N,N,N,N)
 OUTPUT(EPSY,Y,N,N,N,N,N)
 OUTPUT(LTLS,Y,N,Y,Y,Y,Y)
 OUTPUT(WDIS,Y,N,N,N,N,N)
 OUTPUT(TEM1,Y,N,Y,Y,Y,Y)
 OUTPUT(MARK,Y,N,Y,N,N,N)
 OUTPUT(VISL,Y,N,Y,N,N,N)
 OUTPUT(DEN1,Y,N,Y,N,N,N)
 OUTPUT(PRPS,Y,N,Y,N,N,N)
 OUTPUT(NPOR,Y,N,Y,N,N,N)
 OUTPUT(EPOR,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 18 ;IYMON = 5 ;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 = 1 ;IXPRF = 1 ;IXPRL = 10000
 NYPRIN = 1 ;IYPRF = 1 ;IYPRL = 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
   No PATCHes yet used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP