TALK=T;RUN(1,1)
DISPLAY
T$
RANSPORT OF A TOXIC SPILL
An application of multicomponent chemical reaction system model
is presented aimed at the at the simulation of the relevant
physicochemical phenomena taking place in a real-life
environmental case study.
ENDDIS
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(Transport Of Toxic Spill )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 77
XULAST =77.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 30
YVLAST =30.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 1
ZWLAST =1.
************************************************************
Group 6. Body-Fitted Coordinates
BFC = T ;NONORT = T
NCRT = -1
RSTGEO = F ;SAVGEO = F
UUP = F ;VUP = F ;WUP = F
NGEOM =CHAM
NAMXYZ =CHAM
ANGMIN =20.
DOMAIN( 1, 78, 1, 31, 1, 2)
* Set fixed sub-domain
FIXDOM( 1, 0, 0, 0, 0, 0, 0)
FIXDOM( 2, 0, 0, 0, 0, 0, 0)
FIXDOM( 3, 0, 0, 0, 0, 0, 0)
FIXDOM( 4, 0, 0, 0, 0, 0, 0)
FIXDOM( 5, 0, 0, 0, 0, 0, 0)
FIXDOM( 6, 0, 0, 0, 0, 0, 0)
FIXDOM( 7, 0, 0, 0, 0, 0, 0)
FIXDOM( 8, 0, 0, 0, 0, 0, 0)
FIXDOM( 9, 0, 0, 0, 0, 0, 0)
FIXDOM( 10, 0, 0, 0, 0, 0, 0)
LIJ = F ;LJK = T ;LIK = F
GCV = T
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(3)=U1
NAME(5)=V1 ;NAME(16)=C1
NAME(18)=UC1 ;NAME(19)=VC1
NAME(139)=DM ;NAME(140)=EM
NAME(141)=CM ;NAME(144)=PRPS
NAME(145)=BM ;NAME(147)=AM
NAME(148)=WCRT ;NAME(149)=VCRT
NAME(150)=UCRT
* Y in SOLUTN argument list denotes:
* 1-stored 2-solved 3-whole-field
* 4-point-by-point 5-explicit 6-harmonic averaging
SOLUTN(P1,Y,Y,Y,N,N,N)
SOLUTN(U1,Y,Y,N,Y,N,N)
SOLUTN(V1,Y,Y,N,Y,N,N)
SOLUTN(C1,Y,Y,N,N,N,Y)
SOLUTN(UC1,Y,Y,N,N,N,N)
SOLUTN(VC1,Y,Y,N,N,N,N)
SOLUTN(DM,Y,N,N,N,N,Y)
SOLUTN(EM,Y,Y,N,N,N,Y)
SOLUTN(CM,Y,Y,N,N,N,Y)
SOLUTN(PRPS,Y,N,N,N,N,N)
SOLUTN(BM,Y,Y,N,N,N,Y)
SOLUTN(AM,Y,Y,N,N,N,Y)
SOLUTN(WCRT,Y,N,N,N,N,Y)
SOLUTN(VCRT,Y,N,N,N,N,Y)
SOLUTN(UCRT,Y,N,N,N,N,Y)
PRPS = 144
************************************************************
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(C1,N,Y,Y,Y,Y,Y)
TERMS(UC1,N,Y,Y,Y,Y,Y)
TERMS(VC1,N,Y,Y,Y,Y,Y)
TERMS(EM,N,Y,Y,Y,Y,Y)
TERMS(CM,N,Y,Y,Y,N,Y)
TERMS(BM,N,Y,Y,Y,N,Y)
TERMS(AM,N,Y,Y,Y,N,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
UCONV = T
UDIFF = T
USOLVE = 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 =998.22998 ;TMP1 =0. ;EL1 =0.
TSURR =0. ;TEMP0 =0. ;PRESS0 =1.0E+05
DVO1DT =1.18E-04 ;DRH1DP =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
ENUL =1.006E-06 ;ENUT =1.0E-02
PRNDTL(U1)=1. ;PRNDTL(V1)=1.
PRNDTL(C1)=1. ;PRNDTL(UC1)=1.
PRNDTL(VC1)=1. ;PRNDTL(EM)=1.
PRNDTL(CM)=1. ;PRNDTL(BM)=1.
PRNDTL(AM)=1.
PRT(U1)=1. ;PRT(V1)=1.
PRT(C1)=1. ;PRT(UC1)=1.
PRT(VC1)=1. ;PRT(EM)=1.
PRT(CM)=1. ;PRT(BM)=1.
PRT(AM)=1.
CP1 =4181.799805 ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=1.0E-10 ;FIINIT(U1)=1.
FIINIT(V1)=1.0E-10 ;FIINIT(C1)=1.
FIINIT(UC1)=1. ;FIINIT(VC1)=1.0E-10
FIINIT(DM)=1.0E-10 ;FIINIT(EM)=0.
FIINIT(CM)=0. ;FIINIT(PRPS)=-1.
FIINIT(BM)=0. ;FIINIT(AM)=0.
FIINIT(WCRT)=1.0E-10 ;FIINIT(VCRT)=1.0E-10
FIINIT(UCRT)=1.0E-10
PATCH(OB1 ,INIVAL, 6, 7, 3, 3, 1, 1, 1, 1)
INIT(OB1 ,PRPS,0. ,199. )
PATCH(OB2 ,INIVAL, 9, 25, 10, 10, 1, 1, 1, 1)
INIT(OB2 ,PRPS,0. ,199. )
PATCH(OB3 ,INIVAL, 22, 32, 7, 7, 1, 1, 1, 1)
INIT(OB3 ,PRPS,0. ,199. )
PATCH(OB4 ,INIVAL, 34, 57, 6, 6, 1, 1, 1, 1)
INIT(OB4 ,PRPS,0. ,199. )
PATCH(OB5 ,INIVAL, 52, 65, 11, 11, 1, 1, 1, 1)
INIT(OB5 ,PRPS,0. ,199. )
PATCH(OB6 ,INIVAL, 63, 75, 7, 7, 1, 1, 1, 1)
INIT(OB6 ,PRPS,0. ,199. )
PATCH(OB26 ,INIVAL, 2, 2, 9, 9, 1, 1, 1, 1)
INIT(OB26 ,P1 ,0. ,0. )
INIT(OB26 ,C1 ,0. ,0. )
INIT(OB26 ,PRPS,0. ,-1. )
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(CHEM ,VOLUME, 1, 77, 1, 30, 1, 1, 1, 1)
COVAL(CHEM ,EM ,In-Form:source - see Grp 19)
COVAL(CHEM ,CM ,In-Form:source - see Grp 19)
COVAL(CHEM ,BM ,In-Form:source - see Grp 19)
COVAL(CHEM ,AM ,In-Form:source - see Grp 19)
PATCH(BFCOB7 ,WEST , 1, 1, 1, 1, 1, 1, 1, 1)
COVAL(BFCOB7 ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB7 ,U1 ,0. ,0. )
COVAL(BFCOB7 ,V1 ,0. ,0. )
COVAL(BFCOB7 ,C1 ,0. ,1. )
COVAL(BFCOB7 ,UC1 ,0. ,0.04 )
COVAL(BFCOB7 ,VC1 ,0. ,0. )
COVAL(BFCOB7 ,WCRT,0. ,0. )
COVAL(BFCOB7 ,VCRT,0. ,0. )
COVAL(BFCOB7 ,UCRT,0. ,0.04 )
PATCH(BFCOB8 ,WEST , 1, 1, 2, 2, 1, 1, 1, 1)
COVAL(BFCOB8 ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB8 ,U1 ,0. ,0. )
COVAL(BFCOB8 ,V1 ,0. ,0. )
COVAL(BFCOB8 ,C1 ,0. ,1. )
COVAL(BFCOB8 ,UC1 ,0. ,0.16 )
COVAL(BFCOB8 ,VC1 ,0. ,0. )
COVAL(BFCOB8 ,WCRT,0. ,0. )
COVAL(BFCOB8 ,VCRT,0. ,0. )
COVAL(BFCOB8 ,UCRT,0. ,0.16 )
PATCH(BFCOB9 ,WEST , 1, 1, 3, 3, 1, 1, 1, 1)
COVAL(BFCOB9 ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB9 ,U1 ,0. ,0. )
COVAL(BFCOB9 ,V1 ,0. ,0. )
COVAL(BFCOB9 ,C1 ,0. ,1. )
COVAL(BFCOB9 ,UC1 ,0. ,0.24 )
COVAL(BFCOB9 ,VC1 ,0. ,0. )
COVAL(BFCOB9 ,WCRT,0. ,0. )
COVAL(BFCOB9 ,VCRT,0. ,0. )
COVAL(BFCOB9 ,UCRT,0. ,0.24 )
PATCH(BFCOBA ,WEST , 1, 1, 4, 4, 1, 1, 1, 1)
COVAL(BFCOBA ,P1 , FIXFLU , GRND3 )
COVAL(BFCOBA ,U1 ,0. ,0. )
COVAL(BFCOBA ,V1 ,0. ,0. )
COVAL(BFCOBA ,C1 ,0. ,1. )
COVAL(BFCOBA ,UC1 ,0. ,0.34 )
COVAL(BFCOBA ,VC1 ,0. ,0. )
COVAL(BFCOBA ,WCRT,0. ,0. )
COVAL(BFCOBA ,VCRT,0. ,0. )
COVAL(BFCOBA ,UCRT,0. ,0.34 )
PATCH(BFCOBB ,WEST , 1, 1, 5, 5, 1, 1, 1, 1)
COVAL(BFCOBB ,P1 , FIXFLU , GRND3 )
COVAL(BFCOBB ,U1 ,0. ,0. )
COVAL(BFCOBB ,V1 ,0. ,0. )
COVAL(BFCOBB ,C1 ,0. ,1. )
COVAL(BFCOBB ,UC1 ,0. ,0.42 )
COVAL(BFCOBB ,VC1 ,0. ,0. )
COVAL(BFCOBB ,WCRT,0. ,0. )
COVAL(BFCOBB ,VCRT,0. ,0. )
COVAL(BFCOBB ,UCRT,0. ,0.42 )
PATCH(BFCOBC ,WEST , 1, 1, 6, 6, 1, 1, 1, 1)
COVAL(BFCOBC ,P1 , FIXFLU , GRND3 )
COVAL(BFCOBC ,U1 ,0. ,0. )
COVAL(BFCOBC ,V1 ,0. ,0. )
COVAL(BFCOBC ,C1 ,0. ,1. )
COVAL(BFCOBC ,UC1 ,0. ,0.5 )
COVAL(BFCOBC ,VC1 ,0. ,0. )
COVAL(BFCOBC ,WCRT,0. ,0. )
COVAL(BFCOBC ,VCRT,0. ,0. )
COVAL(BFCOBC ,UCRT,0. ,0.5 )
PATCH(BFCOBD ,WEST , 1, 1, 7, 7, 1, 1, 1, 1)
COVAL(BFCOBD ,P1 , FIXFLU , GRND3 )
COVAL(BFCOBD ,U1 ,0. ,0. )
COVAL(BFCOBD ,V1 ,0. ,0. )
COVAL(BFCOBD ,C1 ,0. ,1. )
COVAL(BFCOBD ,UC1 ,0. ,0.58 )
COVAL(BFCOBD ,VC1 ,0. ,0. )
COVAL(BFCOBD ,WCRT,0. ,0. )
COVAL(BFCOBD ,VCRT,0. ,0. )
COVAL(BFCOBD ,UCRT,0. ,0.58 )
PATCH(BFCOBE ,WEST , 1, 1, 8, 8, 1, 1, 1, 1)
COVAL(BFCOBE ,P1 , FIXFLU , GRND3 )
COVAL(BFCOBE ,U1 ,0. ,0. )
COVAL(BFCOBE ,V1 ,0. ,0. )
COVAL(BFCOBE ,C1 ,0. ,1. )
COVAL(BFCOBE ,UC1 ,0. ,0.65 )
COVAL(BFCOBE ,VC1 ,0. ,0. )
COVAL(BFCOBE ,WCRT,0. ,0. )
COVAL(BFCOBE ,VCRT,0. ,0. )
COVAL(BFCOBE ,UCRT,0. ,0.65 )
PATCH(BFCOBF ,WEST , 1, 1, 9, 9, 1, 1, 1, 1)
COVAL(BFCOBF ,P1 , FIXFLU , GRND3 )
COVAL(BFCOBF ,U1 ,0. ,0. )
COVAL(BFCOBF ,V1 ,0. ,0. )
COVAL(BFCOBF ,C1 ,0. ,1. )
COVAL(BFCOBF ,UC1 ,0. ,0.71 )
COVAL(BFCOBF ,VC1 ,0. ,0. )
COVAL(BFCOBF ,WCRT,0. ,0. )
COVAL(BFCOBF ,VCRT,0. ,0. )
COVAL(BFCOBF ,UCRT,0. ,0.71 )
PATCH(BFCOB10 ,WEST , 1, 1, 10, 10, 1, 1, 1, 1)
COVAL(BFCOB10 ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB10 ,U1 ,0. ,0. )
COVAL(BFCOB10 ,V1 ,0. ,0. )
COVAL(BFCOB10 ,C1 ,0. ,1. )
COVAL(BFCOB10 ,UC1 ,0. ,0.78 )
COVAL(BFCOB10 ,VC1 ,0. ,0. )
COVAL(BFCOB10 ,WCRT,0. ,0. )
COVAL(BFCOB10 ,VCRT,0. ,0. )
COVAL(BFCOB10 ,UCRT,0. ,0.78 )
PATCH(BFCOB11 ,WEST , 1, 1, 11, 11, 1, 1, 1, 1)
COVAL(BFCOB11 ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB11 ,U1 ,0. ,0. )
COVAL(BFCOB11 ,V1 ,0. ,0. )
COVAL(BFCOB11 ,C1 ,0. ,1. )
COVAL(BFCOB11 ,UC1 ,0. ,0.82 )
COVAL(BFCOB11 ,VC1 ,0. ,0. )
COVAL(BFCOB11 ,WCRT,0. ,0. )
COVAL(BFCOB11 ,VCRT,0. ,0. )
COVAL(BFCOB11 ,UCRT,0. ,0.82 )
PATCH(BFCOB12 ,WEST , 1, 1, 12, 12, 1, 1, 1, 1)
COVAL(BFCOB12 ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB12 ,U1 ,0. ,0. )
COVAL(BFCOB12 ,V1 ,0. ,0. )
COVAL(BFCOB12 ,C1 ,0. ,1. )
COVAL(BFCOB12 ,UC1 ,0. ,0.84 )
COVAL(BFCOB12 ,VC1 ,0. ,0. )
COVAL(BFCOB12 ,WCRT,0. ,0. )
COVAL(BFCOB12 ,VCRT,0. ,0. )
COVAL(BFCOB12 ,UCRT,0. ,0.84 )
PATCH(BFCOB13 ,WEST , 1, 1, 13, 13, 1, 1, 1, 1)
COVAL(BFCOB13 ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB13 ,U1 ,0. ,0. )
COVAL(BFCOB13 ,V1 ,0. ,0. )
COVAL(BFCOB13 ,C1 ,0. ,1. )
COVAL(BFCOB13 ,UC1 ,0. ,0.85 )
COVAL(BFCOB13 ,VC1 ,0. ,0. )
COVAL(BFCOB13 ,WCRT,0. ,0. )
COVAL(BFCOB13 ,VCRT,0. ,0. )
COVAL(BFCOB13 ,UCRT,0. ,0.85 )
PATCH(BFCOB14 ,WEST , 1, 1, 14, 14, 1, 1, 1, 1)
COVAL(BFCOB14 ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB14 ,U1 ,0. ,0. )
COVAL(BFCOB14 ,V1 ,0. ,0. )
COVAL(BFCOB14 ,C1 ,0. ,1. )
COVAL(BFCOB14 ,UC1 ,0. ,0.855 )
COVAL(BFCOB14 ,VC1 ,0. ,0. )
COVAL(BFCOB14 ,WCRT,0. ,0. )
COVAL(BFCOB14 ,VCRT,0. ,0. )
COVAL(BFCOB14 ,UCRT,0. ,0.855 )
PATCH(BFCOB15 ,WEST , 1, 1, 15, 15, 1, 1, 1, 1)
COVAL(BFCOB15 ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB15 ,U1 ,0. ,0. )
COVAL(BFCOB15 ,V1 ,0. ,0. )
COVAL(BFCOB15 ,C1 ,0. ,1. )
COVAL(BFCOB15 ,UC1 ,0. ,0.85 )
COVAL(BFCOB15 ,VC1 ,0. ,0. )
COVAL(BFCOB15 ,WCRT,0. ,0. )
COVAL(BFCOB15 ,VCRT,0. ,0. )
COVAL(BFCOB15 ,UCRT,0. ,0.85 )
PATCH(BFCOB16 ,WEST , 1, 1, 16, 16, 1, 1, 1, 1)
COVAL(BFCOB16 ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB16 ,U1 ,0. ,0. )
COVAL(BFCOB16 ,V1 ,0. ,0. )
COVAL(BFCOB16 ,C1 ,0. ,1. )
COVAL(BFCOB16 ,UC1 ,0. ,0.83 )
COVAL(BFCOB16 ,VC1 ,0. ,0. )
COVAL(BFCOB16 ,WCRT,0. ,0. )
COVAL(BFCOB16 ,VCRT,0. ,0. )
COVAL(BFCOB16 ,UCRT,0. ,0.83 )
PATCH(BFCOB17 ,WEST , 1, 1, 17, 17, 1, 1, 1, 1)
COVAL(BFCOB17 ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB17 ,U1 ,0. ,0. )
COVAL(BFCOB17 ,V1 ,0. ,0. )
COVAL(BFCOB17 ,C1 ,0. ,1. )
COVAL(BFCOB17 ,UC1 ,0. ,0.81 )
COVAL(BFCOB17 ,VC1 ,0. ,0. )
COVAL(BFCOB17 ,WCRT,0. ,0. )
COVAL(BFCOB17 ,VCRT,0. ,0. )
COVAL(BFCOB17 ,UCRT,0. ,0.81 )
PATCH(BFCOB18 ,WEST , 1, 1, 18, 18, 1, 1, 1, 1)
COVAL(BFCOB18 ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB18 ,U1 ,0. ,0. )
COVAL(BFCOB18 ,V1 ,0. ,0. )
COVAL(BFCOB18 ,C1 ,0. ,1. )
COVAL(BFCOB18 ,UC1 ,0. ,0.79 )
COVAL(BFCOB18 ,VC1 ,0. ,0. )
COVAL(BFCOB18 ,WCRT,0. ,0. )
COVAL(BFCOB18 ,VCRT,0. ,0. )
COVAL(BFCOB18 ,UCRT,0. ,0.79 )
PATCH(BFCOB19 ,WEST , 1, 1, 19, 19, 1, 1, 1, 1)
COVAL(BFCOB19 ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB19 ,U1 ,0. ,0. )
COVAL(BFCOB19 ,V1 ,0. ,0. )
COVAL(BFCOB19 ,C1 ,0. ,1. )
COVAL(BFCOB19 ,UC1 ,0. ,0.76 )
COVAL(BFCOB19 ,VC1 ,0. ,0. )
COVAL(BFCOB19 ,WCRT,0. ,0. )
COVAL(BFCOB19 ,VCRT,0. ,0. )
COVAL(BFCOB19 ,UCRT,0. ,0.76 )
PATCH(BFCOB1A ,WEST , 1, 1, 20, 20, 1, 1, 1, 1)
COVAL(BFCOB1A ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB1A ,U1 ,0. ,0. )
COVAL(BFCOB1A ,V1 ,0. ,0. )
COVAL(BFCOB1A ,C1 ,0. ,1. )
COVAL(BFCOB1A ,UC1 ,0. ,0.73 )
COVAL(BFCOB1A ,VC1 ,0. ,0. )
COVAL(BFCOB1A ,WCRT,0. ,0. )
COVAL(BFCOB1A ,VCRT,0. ,0. )
COVAL(BFCOB1A ,UCRT,0. ,0.73 )
PATCH(BFCOB1B ,WEST , 1, 1, 21, 21, 1, 1, 1, 1)
COVAL(BFCOB1B ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB1B ,U1 ,0. ,0. )
COVAL(BFCOB1B ,V1 ,0. ,0. )
COVAL(BFCOB1B ,C1 ,0. ,1. )
COVAL(BFCOB1B ,UC1 ,0. ,0.7 )
COVAL(BFCOB1B ,VC1 ,0. ,0. )
COVAL(BFCOB1B ,WCRT,0. ,0. )
COVAL(BFCOB1B ,VCRT,0. ,0. )
COVAL(BFCOB1B ,UCRT,0. ,0.7 )
PATCH(BFCOB1C ,WEST , 1, 1, 22, 22, 1, 1, 1, 1)
COVAL(BFCOB1C ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB1C ,U1 ,0. ,0. )
COVAL(BFCOB1C ,V1 ,0. ,0. )
COVAL(BFCOB1C ,C1 ,0. ,1. )
COVAL(BFCOB1C ,UC1 ,0. ,0.67 )
COVAL(BFCOB1C ,VC1 ,0. ,0. )
COVAL(BFCOB1C ,WCRT,0. ,0. )
COVAL(BFCOB1C ,VCRT,0. ,0. )
COVAL(BFCOB1C ,UCRT,0. ,0.67 )
PATCH(BFCOB1D ,WEST , 1, 1, 23, 23, 1, 1, 1, 1)
COVAL(BFCOB1D ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB1D ,U1 ,0. ,0. )
COVAL(BFCOB1D ,V1 ,0. ,0. )
COVAL(BFCOB1D ,C1 ,0. ,1. )
COVAL(BFCOB1D ,UC1 ,0. ,0.64 )
COVAL(BFCOB1D ,VC1 ,0. ,0. )
COVAL(BFCOB1D ,WCRT,0. ,0. )
COVAL(BFCOB1D ,VCRT,0. ,0. )
COVAL(BFCOB1D ,UCRT,0. ,0.64 )
PATCH(BFCOB1E ,WEST , 1, 1, 24, 24, 1, 1, 1, 1)
COVAL(BFCOB1E ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB1E ,U1 ,0. ,0. )
COVAL(BFCOB1E ,V1 ,0. ,0. )
COVAL(BFCOB1E ,C1 ,0. ,1. )
COVAL(BFCOB1E ,UC1 ,0. ,0.62 )
COVAL(BFCOB1E ,VC1 ,0. ,0. )
COVAL(BFCOB1E ,WCRT,0. ,0. )
COVAL(BFCOB1E ,VCRT,0. ,0. )
COVAL(BFCOB1E ,UCRT,0. ,0.62 )
PATCH(BFCOB1F ,WEST , 1, 1, 25, 25, 1, 1, 1, 1)
COVAL(BFCOB1F ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB1F ,U1 ,0. ,0. )
COVAL(BFCOB1F ,V1 ,0. ,0. )
COVAL(BFCOB1F ,C1 ,0. ,1. )
COVAL(BFCOB1F ,UC1 ,0. ,0.57 )
COVAL(BFCOB1F ,VC1 ,0. ,0. )
COVAL(BFCOB1F ,WCRT,0. ,0. )
COVAL(BFCOB1F ,VCRT,0. ,0. )
COVAL(BFCOB1F ,UCRT,0. ,0.57 )
PATCH(BFCOB20 ,WEST , 1, 1, 26, 26, 1, 1, 1, 1)
COVAL(BFCOB20 ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB20 ,U1 ,0. ,0. )
COVAL(BFCOB20 ,V1 ,0. ,0. )
COVAL(BFCOB20 ,C1 ,0. ,1. )
COVAL(BFCOB20 ,UC1 ,0. ,0.54 )
COVAL(BFCOB20 ,VC1 ,0. ,0. )
COVAL(BFCOB20 ,WCRT,0. ,0. )
COVAL(BFCOB20 ,VCRT,0. ,0. )
COVAL(BFCOB20 ,UCRT,0. ,0.54 )
PATCH(BFCOB21 ,WEST , 1, 1, 27, 27, 1, 1, 1, 1)
COVAL(BFCOB21 ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB21 ,U1 ,0. ,0. )
COVAL(BFCOB21 ,V1 ,0. ,0. )
COVAL(BFCOB21 ,C1 ,0. ,1. )
COVAL(BFCOB21 ,UC1 ,0. ,0.48 )
COVAL(BFCOB21 ,VC1 ,0. ,0. )
COVAL(BFCOB21 ,WCRT,0. ,0. )
COVAL(BFCOB21 ,VCRT,0. ,0. )
COVAL(BFCOB21 ,UCRT,0. ,0.48 )
PATCH(BFCOB22 ,WEST , 1, 1, 28, 28, 1, 1, 1, 1)
COVAL(BFCOB22 ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB22 ,U1 ,0. ,0. )
COVAL(BFCOB22 ,V1 ,0. ,0. )
COVAL(BFCOB22 ,C1 ,0. ,1. )
COVAL(BFCOB22 ,UC1 ,0. ,0.4 )
COVAL(BFCOB22 ,VC1 ,0. ,0. )
COVAL(BFCOB22 ,WCRT,0. ,0. )
COVAL(BFCOB22 ,VCRT,0. ,0. )
COVAL(BFCOB22 ,UCRT,0. ,0.4 )
PATCH(BFCOB23 ,WEST , 1, 1, 29, 29, 1, 1, 1, 1)
COVAL(BFCOB23 ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB23 ,U1 ,0. ,0. )
COVAL(BFCOB23 ,V1 ,0. ,0. )
COVAL(BFCOB23 ,C1 ,0. ,1. )
COVAL(BFCOB23 ,UC1 ,0. ,0.3 )
COVAL(BFCOB23 ,VC1 ,0. ,0. )
COVAL(BFCOB23 ,WCRT,0. ,0. )
COVAL(BFCOB23 ,VCRT,0. ,0. )
COVAL(BFCOB23 ,UCRT,0. ,0.3 )
PATCH(BFCOB24 ,WEST , 1, 1, 30, 30, 1, 1, 1, 1)
COVAL(BFCOB24 ,P1 , FIXFLU , GRND3 )
COVAL(BFCOB24 ,U1 ,0. ,0. )
COVAL(BFCOB24 ,V1 ,0. ,0. )
COVAL(BFCOB24 ,C1 ,0. ,1. )
COVAL(BFCOB24 ,UC1 ,0. ,0.08 )
COVAL(BFCOB24 ,VC1 ,0. ,0. )
COVAL(BFCOB24 ,WCRT,0. ,0. )
COVAL(BFCOB24 ,VCRT,0. ,0. )
COVAL(BFCOB24 ,UCRT,0. ,0.08 )
PATCH(OB25 ,EAST , 77, 77, 1, 30, 1, 1, 1, 1)
COVAL(OB25 ,P1 ,1000. ,0. )
COVAL(OB25 ,U1 ,0. ,0. )
COVAL(OB25 ,V1 ,0. ,0. )
COVAL(OB25 ,UC1 ,0. ,0. )
COVAL(OB25 ,VC1 ,0. ,0. )
PATCH(OC26 ,VOLUME, 2, 2, 9, 9, 1, 1, 1, 1)
COVAL(OC26 ,C1 , FIXVAL ,0.99 )
COVAL(OC26 ,EM , FIXVAL ,0. )
COVAL(OC26 ,CM , FIXVAL ,0. )
COVAL(OC26 ,BM , FIXVAL ,7.5E-03 )
COVAL(OC26 ,AM , FIXVAL ,2.5E-03 )
PATCH(OB27 ,LWALL , 1, 77, 1, 30, 1, 1, 1, 1)
COVAL(OB27 ,UC1 , GRND2 ,0. )
COVAL(OB27 ,VC1 , GRND2 ,0. )
XCYCLE = F
EGWF = T
WALLCO = GRND2
BFCA =1000.
************************************************************
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-03
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(U1)=10
LITER(V1)=10 ;LITER(C1)=20
LITER(UC1)=10 ;LITER(VC1)=10
LITER(EM)=20 ;LITER(CM)=20
LITER(BM)=20 ;LITER(AM)=20
ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
ENDIT(V1)=1.0E-03 ;ENDIT(C1)=1.0E-03
ENDIT(UC1)=1.0E-03 ;ENDIT(VC1)=1.0E-03
ENDIT(EM)=1.0E-03 ;ENDIT(CM)=1.0E-03
ENDIT(BM)=1.0E-03 ;ENDIT(AM)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(U1,FALSDT,1000.)
RELAX(V1,FALSDT,1000.)
RELAX(C1,FALSDT,1.0E+09)
RELAX(UC1,FALSDT,1000.)
RELAX(VC1,FALSDT,1000.)
RELAX(DM,LINRLX,1.)
RELAX(EM,FALSDT,1000.)
RELAX(CM,FALSDT,1000.)
RELAX(PRPS,LINRLX,1.)
RELAX(BM,FALSDT,1000.)
RELAX(AM,FALSDT,1000.)
RELAX(WCRT,LINRLX,1.)
RELAX(VCRT,LINRLX,1.)
RELAX(UCRT,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(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(C1)=1.0E+10 ;VARMIN(C1)=-1.0E+10
VARMAX(UC1)=1.0E+06 ;VARMIN(UC1)=-1.0E+06
VARMAX(VC1)=1.0E+06 ;VARMIN(VC1)=-1.0E+06
VARMAX(DM)=1. ;VARMIN(DM)=0.
VARMAX(EM)=1. ;VARMIN(EM)=0.
VARMAX(CM)=1. ;VARMIN(CM)=0.
VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
VARMAX(BM)=1. ;VARMIN(BM)=0.
VARMAX(AM)=1. ;VARMIN(AM)=0.
VARMAX(WCRT)=1.0E+10 ;VARMIN(WCRT)=-1.0E+10
VARMAX(VCRT)=1.0E+10 ;VARMIN(VCRT)=-1.0E+10
VARMAX(UCRT)=1.0E+10 ;VARMIN(UCRT)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
PARSOL = F
ISG62 = 1
SPEDAT(SET,DOMAIN,PHASE_1_MAT,I,67)
SPEDAT(SET,SOURCE,AM!CHEM,C,=COVAL(2.5*BM&0.0))
SPEDAT(SET,SOURCE,BM!CHEM,C,=COVAL(5.55*AM&0.0))
SPEDAT(SET,SOURCE,CM!CHEM,C,=COVAL(0.215&24.9*AM*BM))
SPEDAT(SET,SOURCE,EM!CHEM,C,=COVAL(2.0E-10&0.18*CM))
SPEDAT(SET,STORED,DM,C,=MAX(0.0&(1.-AM-BM-CM-EM-C1))!ZSLFIN)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,OBJNAM,!OB1,C,CMP0)
SPEDAT(SET,OBJTYP,!OB1,C,BLOCKAGE)
SPEDAT(SET,CMP0,MATERIAL,R,199.)
SPEDAT(SET,OBJNAM,!OB2,C,CMP1)
SPEDAT(SET,OBJTYP,!OB2,C,BLOCKAGE)
SPEDAT(SET,CMP1,MATERIAL,R,199.)
SPEDAT(SET,OBJNAM,!OB3,C,CMP2)
SPEDAT(SET,OBJTYP,!OB3,C,BLOCKAGE)
SPEDAT(SET,CMP2,MATERIAL,R,199.)
SPEDAT(SET,OBJNAM,!OB4,C,CMP3)
SPEDAT(SET,OBJTYP,!OB4,C,BLOCKAGE)
SPEDAT(SET,CMP3,MATERIAL,R,199.)
SPEDAT(SET,OBJNAM,!OB5,C,CMP4)
SPEDAT(SET,OBJTYP,!OB5,C,BLOCKAGE)
SPEDAT(SET,CMP4,MATERIAL,R,199.)
SPEDAT(SET,OBJNAM,!OB6,C,CMP5)
SPEDAT(SET,OBJTYP,!OB6,C,BLOCKAGE)
SPEDAT(SET,CMP5,MATERIAL,R,199.)
SPEDAT(SET,OBJNAM,!BFCOB7,C,INLE1)
SPEDAT(SET,OBJTYP,!BFCOB7,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB8,C,INLE2)
SPEDAT(SET,OBJTYP,!BFCOB8,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB9,C,INLE3)
SPEDAT(SET,OBJTYP,!BFCOB9,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOBA,C,INLE4)
SPEDAT(SET,OBJTYP,!BFCOBA,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOBB,C,INLE5)
SPEDAT(SET,OBJTYP,!BFCOBB,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOBC,C,INLE6)
SPEDAT(SET,OBJTYP,!BFCOBC,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOBD,C,INLE7)
SPEDAT(SET,OBJTYP,!BFCOBD,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOBE,C,INLE8)
SPEDAT(SET,OBJTYP,!BFCOBE,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOBF,C,INLE9)
SPEDAT(SET,OBJTYP,!BFCOBF,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB10,C,INLE10)
SPEDAT(SET,OBJTYP,!BFCOB10,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB11,C,INLE11)
SPEDAT(SET,OBJTYP,!BFCOB11,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB12,C,INLE12)
SPEDAT(SET,OBJTYP,!BFCOB12,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB13,C,INLE13)
SPEDAT(SET,OBJTYP,!BFCOB13,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB14,C,INLE14)
SPEDAT(SET,OBJTYP,!BFCOB14,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB15,C,INLE15)
SPEDAT(SET,OBJTYP,!BFCOB15,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB16,C,INLE16)
SPEDAT(SET,OBJTYP,!BFCOB16,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB17,C,INLE17)
SPEDAT(SET,OBJTYP,!BFCOB17,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB18,C,INLE18)
SPEDAT(SET,OBJTYP,!BFCOB18,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB19,C,INLE19)
SPEDAT(SET,OBJTYP,!BFCOB19,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB1A,C,INLE20)
SPEDAT(SET,OBJTYP,!BFCOB1A,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB1B,C,INLE21)
SPEDAT(SET,OBJTYP,!BFCOB1B,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB1C,C,INLE22)
SPEDAT(SET,OBJTYP,!BFCOB1C,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB1D,C,INLE23)
SPEDAT(SET,OBJTYP,!BFCOB1D,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB1E,C,INLE24)
SPEDAT(SET,OBJTYP,!BFCOB1E,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB1F,C,INLE25)
SPEDAT(SET,OBJTYP,!BFCOB1F,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB20,C,INLE26)
SPEDAT(SET,OBJTYP,!BFCOB20,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB21,C,INLE27)
SPEDAT(SET,OBJTYP,!BFCOB21,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB22,C,INLE28)
SPEDAT(SET,OBJTYP,!BFCOB22,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB23,C,INLE29)
SPEDAT(SET,OBJTYP,!BFCOB23,C,INLET)
SPEDAT(SET,OBJNAM,!BFCOB24,C,INLE30)
SPEDAT(SET,OBJTYP,!BFCOB24,C,INLET)
SPEDAT(SET,OBJNAM,!OB25,C,OUTLET)
SPEDAT(SET,OBJTYP,!OB25,C,OUTLET)
SPEDAT(SET,ARATIO,!OB25,R,1.)
SPEDAT(SET,OBJNAM,!OB26,C,B39)
SPEDAT(SET,OBJTYP,!OB26,C,BLOCKAGE)
SPEDAT(SET,B39,MATERIAL,R,-1.)
SPEDAT(SET,OBJNAM,!OC26,C,B39)
SPEDAT(SET,OBJTYP,!OC26,C,BLOCKAGE)
SPEDAT(SET,OBJNAM,!OB27,C,B40)
SPEDAT(SET,OBJTYP,!OB27,C,PLATE)
SPEDAT(SET,FACETDAT,NUMOBJ,I,39)
SPEDAT(SET,MATERIAL,199,L,T)
************************************************************
Group 20. Preliminary Printout
************************************************************
Group 21. Print-out of Variables
INIFLD = T ;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,N,N,N)
OUTPUT(V1,Y,N,Y,N,N,N)
OUTPUT(C1,Y,N,Y,Y,Y,Y)
OUTPUT(UC1,Y,N,Y,Y,Y,Y)
OUTPUT(VC1,Y,N,Y,Y,Y,Y)
OUTPUT(DM,Y,N,Y,N,N,N)
OUTPUT(EM,Y,N,Y,Y,Y,Y)
OUTPUT(CM,Y,N,Y,Y,Y,Y)
OUTPUT(PRPS,Y,N,Y,N,N,N)
OUTPUT(BM,Y,N,Y,Y,Y,Y)
OUTPUT(AM,Y,N,Y,Y,Y,Y)
OUTPUT(WCRT,Y,N,Y,N,N,N)
OUTPUT(VCRT,Y,N,Y,N,N,N)
OUTPUT(UCRT,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 76 ;IYMON = 3 ;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 = 0 ;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