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,GXMONI,PLOTALL,L,T)
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,Y,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