TALK=T;RUN(1,1)
  DISPLAY
  Library Case I407: Aerosol dispersion and deposition in a small
  room.
  The case considered is steady, isothermal airflow with aerosol
  transport and deposition of three particle sizes in a small room.
  The Eulerian drift flux model is used to simulate aerosol
  dispersion with surface deposition via the 3-layer deposition
  model of Chen & Lai (2004), which accounts for the deposition
  mechanisms of gravity, Brownian & turbulent diffusion.
  ENDDIS
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(I407: Aerosol deposition in a small room)
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 0
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 25
 XULAST =2.
 XFRAC(1)=0.03257 ;XFRAC(2)=0.068397
 XFRAC(3)=0.107807 ;XFRAC(4)=0.151157
 XFRAC(5)=0.198843 ;XFRAC(6)=0.242193
 XFRAC(7)=0.281603 ;XFRAC(8)=0.31743
 XFRAC(9)=0.35 ;XFRAC(10)=0.374194
 XFRAC(11)=0.400806 ;XFRAC(12)=0.425
 XFRAC(13)=0.45 ;XFRAC(14)=0.482258
 XFRAC(15)=0.517742 ;XFRAC(16)=0.55
 XFRAC(17)=0.575 ;XFRAC(18)=0.6125
 XFRAC(19)=0.65 ;XFRAC(20)=0.70287
 XFRAC(21)=0.761027 ;XFRAC(22)=0.825
 XFRAC(23)=0.888973 ;XFRAC(24)=0.94713
 XFRAC(25)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 39
 YVLAST =2.
 YFRAC(1)=0.017989 ;YFRAC(2)=0.037777
 YFRAC(3)=0.059543 ;YFRAC(4)=0.083487
 YFRAC(5)=0.109824 ;YFRAC(6)=0.138796
 YFRAC(7)=0.170664 ;YFRAC(8)=0.20572
 YFRAC(9)=0.24428 ;YFRAC(10)=0.279336
 YFRAC(11)=0.311204 ;YFRAC(12)=0.340176
 YFRAC(13)=0.366513 ;YFRAC(14)=0.390457
 YFRAC(15)=0.412223 ;YFRAC(16)=0.432011
 YFRAC(17)=0.45 ;YFRAC(18)=0.475
 YFRAC(19)=0.491129 ;YFRAC(20)=0.508871
 YFRAC(21)=0.525 ;YFRAC(22)=0.55
 YFRAC(23)=0.567989 ;YFRAC(24)=0.587777
 YFRAC(25)=0.609543 ;YFRAC(26)=0.633487
 YFRAC(27)=0.659824 ;YFRAC(28)=0.688796
 YFRAC(29)=0.720664 ;YFRAC(30)=0.75572
 YFRAC(31)=0.79428 ;YFRAC(32)=0.829336
 YFRAC(33)=0.861204 ;YFRAC(34)=0.890176
 YFRAC(35)=0.916513 ;YFRAC(36)=0.940457
 YFRAC(37)=0.962223 ;YFRAC(38)=0.982011
 YFRAC(39)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 58
 ZWLAST =2.8
 ZFRAC(1)=4.6845E-03 ;ZFRAC(2)=9.83745E-03
 ZFRAC(3)=0.015506 ;ZFRAC(4)=0.021741
 ZFRAC(5)=0.028599 ;ZFRAC(6)=0.036144
 ZFRAC(7)=0.044443 ;ZFRAC(8)=0.053571
 ZFRAC(9)=0.0627 ;ZFRAC(10)=0.070999
 ZFRAC(11)=0.078543 ;ZFRAC(12)=0.085402
 ZFRAC(13)=0.091637 ;ZFRAC(14)=0.097305
 ZFRAC(15)=0.102458 ;ZFRAC(16)=0.107143
 ZFRAC(17)=0.114118 ;ZFRAC(18)=0.121791
 ZFRAC(19)=0.130231 ;ZFRAC(20)=0.139515
 ZFRAC(21)=0.149727 ;ZFRAC(22)=0.160961
 ZFRAC(23)=0.173317 ;ZFRAC(24)=0.18691
 ZFRAC(25)=0.201862 ;ZFRAC(26)=0.218309
 ZFRAC(27)=0.236401 ;ZFRAC(28)=0.256302
 ZFRAC(29)=0.278193 ;ZFRAC(30)=0.302273
 ZFRAC(31)=0.328761 ;ZFRAC(32)=0.357898
 ZFRAC(33)=0.389949 ;ZFRAC(34)=0.425205
 ZFRAC(35)=0.463986 ;ZFRAC(36)=0.506646
 ZFRAC(37)=0.553571 ;ZFRAC(38)=0.600497
 ZFRAC(39)=0.643156 ;ZFRAC(40)=0.681938
 ZFRAC(41)=0.717194 ;ZFRAC(42)=0.749244
 ZFRAC(43)=0.778381 ;ZFRAC(44)=0.80487
 ZFRAC(45)=0.82895 ;ZFRAC(46)=0.850841
 ZFRAC(47)=0.870742 ;ZFRAC(48)=0.888834
 ZFRAC(49)=0.905281 ;ZFRAC(50)=0.920233
 ZFRAC(51)=0.933825 ;ZFRAC(52)=0.946182
 ZFRAC(53)=0.957416 ;ZFRAC(54)=0.967628
 ZFRAC(55)=0.976912 ;ZFRAC(56)=0.985352
 ZFRAC(57)=0.993025 ;ZFRAC(58)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(1)=P1 ;NAME(3)=U1
 NAME(5)=V1 ;NAME(7)=W1
 NAME(12)=KE ;NAME(13)=EP
 NAME(21)=C6 ;NAME(22)=C7
 NAME(23)=C8 ;NAME(126)=EPKE
 NAME(127)=V3DH ;NAME(128)=VD3
 NAME(129)=VD2 ;NAME(130)=VS3
 NAME(131)=VS1 ;NAME(132)=VDPH
 NAME(133)=DI1 ;NAME(134)=VP1
 NAME(135)=VSTR ;NAME(136)=TP1
 NAME(137)=V2DH ;NAME(138)=PRPS
 NAME(139)=VABS ;NAME(140)=YPLS
 NAME(141)=VD1 ;NAME(142)=STRS
 NAME(143)=V1DH ;NAME(144)=ENUL
 NAME(145)=SPH1 ;NAME(146)=VS2
 NAME(147)=TEM1 ;NAME(148)=DEN1
 NAME(149)=EL1 ;NAME(150)=ENUT
    * 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,Y)
 SOLUTN(U1,Y,Y,N,N,N,Y)
 SOLUTN(V1,Y,Y,N,N,N,Y)
 SOLUTN(W1,Y,Y,N,N,N,Y)
 SOLUTN(KE,Y,Y,N,N,N,N)
 SOLUTN(EP,Y,Y,N,N,N,N)
 SOLUTN(C6,Y,Y,Y,N,N,Y)
 SOLUTN(C7,Y,Y,Y,N,N,Y)
 SOLUTN(C8,Y,Y,Y,N,N,Y)
 SOLUTN(EPKE,Y,N,N,N,N,Y)
 SOLUTN(V3DH,Y,N,N,N,N,Y)
 SOLUTN(VD3,Y,N,N,N,N,Y)
 SOLUTN(VD2,Y,N,N,N,N,Y)
 SOLUTN(VS3,Y,N,N,N,N,Y)
 SOLUTN(VS1,Y,N,N,N,N,Y)
 SOLUTN(VDPH,Y,N,N,N,N,Y)
 SOLUTN(DI1,Y,N,N,N,N,Y)
 SOLUTN(VP1,Y,N,N,N,N,Y)
 SOLUTN(VSTR,Y,N,N,N,N,Y)
 SOLUTN(TP1,Y,N,N,N,N,Y)
 SOLUTN(V2DH,Y,N,N,N,N,Y)
 SOLUTN(PRPS,Y,N,N,N,N,Y)
 SOLUTN(VABS,Y,N,N,N,N,Y)
 SOLUTN(YPLS,Y,N,N,N,N,Y)
 SOLUTN(VD1,Y,N,N,N,N,Y)
 SOLUTN(STRS,Y,N,N,N,N,Y)
 SOLUTN(V1DH,Y,N,N,N,N,Y)
 SOLUTN(ENUL,Y,N,N,N,N,Y)
 SOLUTN(SPH1,Y,N,N,N,N,Y)
 SOLUTN(VS2,Y,N,N,N,N,Y)
 SOLUTN(TEM1,Y,Y,N,N,N,Y)
 SOLUTN(DEN1,Y,N,N,N,N,Y)
 SOLUTN(EL1,Y,N,N,N,N,Y)
 SOLUTN(ENUT,Y,N,N,N,N,Y)
 DEN1 = 148
 VISL = 144
 VIST = 150
 LEN1 = 149
 PRPS = 138
 ************************************************************
  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(W1,Y,Y,Y,Y,Y,Y)
 TERMS(KE,N,Y,Y,Y,Y,N)
 TERMS(EP,N,Y,Y,Y,Y,N)
 TERMS(C6,N,Y,Y,Y,N,Y)
 TERMS(C7,N,Y,Y,Y,Y,Y)
 TERMS(C8,N,Y,Y,Y,N,Y)
 TERMS(TEM1,N,Y,Y,Y,N,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWRH1 = T
 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 = GRND5 ;TMP1 =0. ;EL1 = GRND4
 TSURR =0. ;TEMP0 =273. ;PRESS0 =1.01325E+05
 DVO1DT =3.33E-03 ;DRH1DP = GRND5
 RHO1A =0. ;RHO1B =3.4843E-03 ;RHO1C =0.71429
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 EL1A =0. ;EL1B =0. ;EL1C =0.
 ENUL =1.0E-05 ;ENUT = GRND3
 ENUTA =0. ;ENUTB =0. ;ENUTC =0.
 IENUTA = 1
 PRNDTL(U1)=1. ;PRNDTL(V1)=1.
 PRNDTL(W1)=1. ;PRNDTL(KE)=1.
 PRNDTL(EP)=1. ;PRNDTL(C6)=0.7
 PRNDTL(C7)=1. ;PRNDTL(C8)=1.
 PRNDTL(TEM1)=-0.0263
 PRT(U1)=1. ;PRT(V1)=1.
 PRT(W1)=1. ;PRT(KE)=0.7194
 PRT(EP)=0.7194 ;PRT(C6)=1.
 PRT(C7)=1. ;PRT(C8)=1.
 PRT(TEM1)=1.
 CP1 =1004. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=0. ;FIINIT(U1)=1.0E-10
 FIINIT(V1)=1.0E-10 ;FIINIT(W1)=0.225
 FIINIT(KE)=1.0E-03 ;FIINIT(EP)=9.0E-03
 FIINIT(C6)=0. ;FIINIT(C7)=0.
 FIINIT(C8)=0. ;FIINIT(EPKE)=1.0E-10
 FIINIT(V3DH)=1.0E-10 ;FIINIT(VD3)=1.0E-10
 FIINIT(VD2)=1.0E-10 ;FIINIT(VS3)=1.0E-10
 FIINIT(VS1)=1.0E-10 ;FIINIT(VDPH)=1.0E-10
 FIINIT(DI1)=1.0E-10 ;FIINIT(VP1)=1.0E-10
 FIINIT(VSTR)=1.0E-10 ;FIINIT(TP1)=1.0E-10
 FIINIT(V2DH)=1.0E-10 ;FIINIT(PRPS)=-1.
 FIINIT(VABS)=1.0E-10 ;FIINIT(YPLS)=1.0E-10
 FIINIT(VD1)=1.0E-10 ;FIINIT(STRS)=1.0E-10
 FIINIT(V1DH)=1.0E-10 ;FIINIT(ENUL)=1.0E-10
 FIINIT(SPH1)=1.0E-10 ;FIINIT(VS2)=1.0E-10
 FIINIT(TEM1)=20. ;FIINIT(DEN1)=1.0E-10
 FIINIT(EL1)=1.0E-10 ;FIINIT(ENUT)=1.0E-10
   No PATCHes yet used for this Group
 INIADD = F
 FSWEEP = 1
 NAMFI =phida
 ************************************************************
  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(DFLUX ,CELL , 1, 25, 1, 39, 1, 58, 1, 1)
 COVAL(DFLUX ,C6 , GRND4 , GRND4 )
 COVAL(DFLUX ,C7 , GRND4 , GRND4 )
 COVAL(DFLUX ,C8 , GRND4 , GRND4 )
 
 PATCH(BUOYANCY,PHASEM, 1, 25, 1, 39, 1, 58, 1, 1)
 COVAL(BUOYANCY,U1 , FIXFLU , GRND2 )
 COVAL(BUOYANCY,V1 , FIXFLU , GRND2 )
 COVAL(BUOYANCY,W1 , FIXFLU , GRND2 )
 
 PATCH(KESOURCE,PHASEM, 1, 25, 1, 39, 1, 58, 1, 1)
 COVAL(KESOURCE,KE , GRND4 , GRND4 )
 COVAL(KESOURCE,EP , GRND4 , GRND4 )
 
 PATCH(RNGMKE ,PHASEM, 1, 25, 1, 39, 1, 58, 1, 1)
 COVAL(RNGMKE ,EP , GRND4 , GRND4 )
 
 PATCH(OB1 ,LOW , 13, 17, 19, 21, 1, 1, 1, 1)
 COVAL(OB1 ,P1 , FIXFLU ,0.117682 )
 COVAL(OB1 ,U1 ,0. ,0. )
 COVAL(OB1 ,V1 ,0. ,0. )
 COVAL(OB1 ,W1 ,0. ,0.1 )
 COVAL(OB1 ,KE ,0. ,1.0E-03 )
 COVAL(OB1 ,EP ,0. ,9.0E-03 )
 COVAL(OB1 ,C6 ,0. ,1. )
 COVAL(OB1 ,C7 ,0. ,1. )
 COVAL(OB1 ,C8 ,0. ,1. )
 COVAL(OB1 ,TEM1,0. ,20. )
 
 PATCH(OB2 ,HIGH , 14, 16, 18, 22, 58, 58, 1, 1)
 COVAL(OB2 ,P1 ,In-Form:source - see Grp 19)
 COVAL(OB2 ,KE ,0. ,5.625E-05 )
 COVAL(OB2 ,EP ,0. ,6.6119E-06 )
 COVAL(OB2 ,TEM1,0. ,20. )
 
 PATCH(OB3 ,NWALL , 1, 25, 39, 39, 1, 58, 1, 1)
 COVAL(OB3 ,U1 , GRND2 ,0. )
 COVAL(OB3 ,W1 , GRND2 ,0. )
 COVAL(OB3 ,KE , GRND2 , GRND2 )
 COVAL(OB3 ,EP , GRND2 , GRND2 )
 
 PATCH(PW3 ,NORTH , 1, 25, 39, 39, 1, 58, 1, 1)
 COVAL(PW3 ,C6 , GRND4 ,0. )
 COVAL(PW3 ,C7 , GRND4 ,0. )
 COVAL(PW3 ,C8 , GRND4 ,0. )
 
 PATCH(OB4 ,LWALL , 1, 25, 1, 39, 1, 1, 1, 1)
 COVAL(OB4 ,U1 , GRND2 ,0. )
 COVAL(OB4 ,V1 , GRND2 ,0. )
 COVAL(OB4 ,KE , GRND2 , GRND2 )
 COVAL(OB4 ,EP , GRND2 , GRND2 )
 
 PATCH(PW4 ,LOW , 1, 25, 1, 39, 1, 1, 1, 1)
 COVAL(PW4 ,C6 , GRND4 ,0. )
 COVAL(PW4 ,C7 , GRND4 ,0. )
 COVAL(PW4 ,C8 , GRND4 ,0. )
 
 PATCH(OB5 ,HWALL , 1, 25, 1, 39, 58, 58, 1, 1)
 COVAL(OB5 ,U1 , GRND2 ,0. )
 COVAL(OB5 ,V1 , GRND2 ,0. )
 COVAL(OB5 ,KE , GRND2 , GRND2 )
 COVAL(OB5 ,EP , GRND2 , GRND2 )
 
 PATCH(PW5 ,HIGH , 1, 25, 1, 39, 58, 58, 1, 1)
 COVAL(PW5 ,C6 , GRND4 ,0. )
 COVAL(PW5 ,C7 , GRND4 ,0. )
 COVAL(PW5 ,C8 , GRND4 ,0. )
 
 PATCH(OB6 ,EWALL , 25, 25, 1, 39, 1, 58, 1, 1)
 COVAL(OB6 ,V1 , GRND2 ,0. )
 COVAL(OB6 ,W1 , GRND2 ,0. )
 COVAL(OB6 ,KE , GRND2 , GRND2 )
 COVAL(OB6 ,EP , GRND2 , GRND2 )
 
 PATCH(PW6 ,EAST , 25, 25, 1, 39, 1, 58, 1, 1)
 COVAL(PW6 ,C6 , GRND4 ,0. )
 COVAL(PW6 ,C7 , GRND4 ,0. )
 COVAL(PW6 ,C8 , GRND4 ,0. )
 
 PATCH(OB7 ,WWALL , 1, 1, 1, 39, 1, 58, 1, 1)
 COVAL(OB7 ,V1 , GRND2 ,0. )
 COVAL(OB7 ,W1 , GRND2 ,0. )
 COVAL(OB7 ,KE , GRND2 , GRND2 )
 COVAL(OB7 ,EP , GRND2 , GRND2 )
 
 PATCH(PW7 ,WEST , 1, 1, 1, 39, 1, 58, 1, 1)
 COVAL(PW7 ,C6 , GRND4 ,0. )
 COVAL(PW7 ,C7 , GRND4 ,0. )
 COVAL(PW7 ,C8 , GRND4 ,0. )
 
 PATCH(OB8 ,NORTH , 10, 19, 39, 39, 1, 16, 1, 1)
 COVAL(OB8 ,P1 ,1. ,0. )
 COVAL(OB8 ,U1 ,0. ,0. )
 COVAL(OB8 ,V1 ,0. ,0. )
 COVAL(OB8 ,W1 ,0. ,0. )
 COVAL(OB8 ,KE ,0. , SAME )
 COVAL(OB8 ,EP ,0. , SAME )
 COVAL(OB8 ,TEM1,0. ,20. )
 
 PATCH(OB9 ,SWALL , 1, 25, 1, 1, 1, 58, 1, 1)
 COVAL(OB9 ,U1 , GRND2 ,0. )
 COVAL(OB9 ,W1 , GRND2 ,0. )
 COVAL(OB9 ,KE , GRND2 , GRND2 )
 COVAL(OB9 ,EP , GRND2 , GRND2 )
 
 PATCH(PW9 ,SOUTH , 1, 25, 1, 1, 1, 58, 1, 1)
 COVAL(PW9 ,C6 , GRND4 ,0. )
 COVAL(PW9 ,C7 , GRND4 ,0. )
 COVAL(PW9 ,C8 , GRND4 ,0. )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 BUOYA =0. ; BUOYB =0. ; BUOYC =-9.81
 BUOYD =1.204938
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 2000 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-04
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(U1)=10
 LITER(V1)=10 ;LITER(W1)=10
 LITER(KE)=10 ;LITER(EP)=10
 LITER(C6)=50 ;LITER(C7)=50
 LITER(C8)=50 ;LITER(TEM1)=50
 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
 ENDIT(V1)=1.0E-03 ;ENDIT(W1)=1.0E-03
 ENDIT(KE)=1.0E-03 ;ENDIT(EP)=1.0E-03
 ENDIT(C6)=1.0E-03 ;ENDIT(C7)=1.0E-03
 ENDIT(C8)=1.0E-03 ;ENDIT(TEM1)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(U1,FALSDT,10.)
 RELAX(V1,FALSDT,10.)
 RELAX(W1,FALSDT,10.)
 RELAX(KE,LINRLX,0.5)
 RELAX(EP,LINRLX,0.5)
 RELAX(C6,FALSDT,10.)
 RELAX(C7,FALSDT,10.)
 RELAX(C8,FALSDT,10.)
 RELAX(EPKE,LINRLX,1.)
 RELAX(V3DH,LINRLX,1.)
 RELAX(VD3,LINRLX,1.)
 RELAX(VD2,LINRLX,1.)
 RELAX(VS3,LINRLX,1.)
 RELAX(VS1,LINRLX,1.)
 RELAX(VDPH,LINRLX,1.)
 RELAX(DI1,LINRLX,1.)
 RELAX(VP1,LINRLX,1.)
 RELAX(VSTR,LINRLX,1.)
 RELAX(TP1,LINRLX,1.)
 RELAX(V2DH,LINRLX,1.)
 RELAX(PRPS,LINRLX,1.)
 RELAX(VABS,LINRLX,1.)
 RELAX(YPLS,LINRLX,1.)
 RELAX(VD1,LINRLX,1.)
 RELAX(STRS,LINRLX,1.)
 RELAX(V1DH,LINRLX,1.)
 RELAX(ENUL,LINRLX,1.)
 RELAX(SPH1,LINRLX,1.)
 RELAX(VS2,LINRLX,1.)
 RELAX(TEM1,LINRLX,0.5)
 RELAX(DEN1,LINRLX,1.)
 RELAX(EL1,LINRLX,1.)
 RELAX(ENUT,LINRLX,1.)
 KELIN = 3
 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(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
 VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
 VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
 VARMAX(C6)=1.1 ;VARMIN(C6)=0.
 VARMAX(C7)=1.1 ;VARMIN(C7)=0.
 VARMAX(C8)=1.1 ;VARMIN(C8)=0.
 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
 VARMAX(V3DH)=1.0E+10 ;VARMIN(V3DH)=-1.0E+10
 VARMAX(VD3)=1.0E+10 ;VARMIN(VD3)=-1.0E+10
 VARMAX(VD2)=1.0E+10 ;VARMIN(VD2)=-1.0E+10
 VARMAX(VS3)=1.0E+10 ;VARMIN(VS3)=-1.0E+10
 VARMAX(VS1)=1.0E+10 ;VARMIN(VS1)=-1.0E+10
 VARMAX(VDPH)=1.0E+10 ;VARMIN(VDPH)=-1.0E+10
 VARMAX(DI1)=1.0E+10 ;VARMIN(DI1)=-1.0E+10
 VARMAX(VP1)=1.0E+10 ;VARMIN(VP1)=-1.0E+10
 VARMAX(VSTR)=1.0E+10 ;VARMIN(VSTR)=-1.0E+10
 VARMAX(TP1)=1.0E+10 ;VARMIN(TP1)=-1.0E+10
 VARMAX(V2DH)=1.0E+10 ;VARMIN(V2DH)=-1.0E+10
 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
 VARMAX(VABS)=1.0E+10 ;VARMIN(VABS)=-1.0E+10
 VARMAX(YPLS)=1.0E+10 ;VARMIN(YPLS)=-1.0E+10
 VARMAX(VD1)=1.0E+10 ;VARMIN(VD1)=-1.0E+10
 VARMAX(STRS)=1.0E+10 ;VARMIN(STRS)=-1.0E+10
 VARMAX(V1DH)=1.0E+10 ;VARMIN(V1DH)=-1.0E+10
 VARMAX(ENUL)=1.0E+10 ;VARMIN(ENUL)=-1.0E+10
 VARMAX(SPH1)=1.0E+10 ;VARMIN(SPH1)=-1.0E+10
 VARMAX(VS2)=1.0E+10 ;VARMIN(VS2)=-1.0E+10
 VARMAX(TEM1)=3000. ;VARMIN(TEM1)=-30.
 VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=-1.0E+10
 VARMAX(EL1)=1.0E+10 ;VARMIN(EL1)=-1.0E+10
 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 NAMGRD =FLAR
 GENK = T
 PARSOL = F
 IENUTA = 1
 ISG62 = 1
 SPEDAT(SET,DOMAIN,PHASE_1_MAT,I,2)
 SPEDAT(SET,DFLUX,CONSLP,L,F)
 SPEDAT(SET,DFLUX,JWDEP,I,2)
 SPEDAT(SET,DFLUX,DFMODL,L,T)
 SPEDAT(SET,DFLUX,DEPOMOD,I,4)
 SPEDAT(SET,DFLUX,DENP1,R,1000.)
 SPEDAT(SET,DFLUX,DIAP1,R,1.0E-05)
 SPEDAT(SET,DFLUX,DENP2,R,1000.)
 SPEDAT(SET,DFLUX,DIAP2,R,6.0E-06)
 SPEDAT(SET,DFLUX,DENP3,R,1000.)
 SPEDAT(SET,DFLUX,DIAP3,R,2.0E-06)
 SPEDAT(SET,MAKE,SUMAF,C,=0.0)
 SPEDAT(SET,STORED,SUMAF!PW4,C,=SUM(AHIGH)!ZSLFIN)
 SPEDAT(SET,PRINT,FLOOR_AREA,C,=SUMAF)
 SPEDAT(SET,STORED,V1DH,C,=VD1*3600.0!ZSLFIN)
 SPEDAT(SET,MAKE,P1DFR,C,=0.0)
 SPEDAT(SET,MAKE,SUMVD1,C,=0.0)
 SPEDAT(SET,MAKE,SUMC6,C,=0.0)
 SPEDAT(SET,STORED,P1DFR,C,=ABS(NETS(C6&PW4))/NETS(C6&INL)!ZSLFIN)
 SPEDAT(SET,STORED,SUMVD1!PW4,C,=SUM(V1DH*AHIGH)!ZSLFIN)
 SPEDAT(SET,STORED,SUMC6!PW4,C,=SUM(C6*AHIGH)!ZSLFIN)
 SPEDAT(SET,PRINT,PAR1_DFRAC,C,=P1DFR)
 SPEDAT(SET,PRINT,FLOOR_AVG_VDEP1,C,=SUMVD1/SUMAF)
 SPEDAT(SET,PRINT,FLOOR_AVG_C6,C,=SUMC6/SUMAF)
 SPEDAT(SET,STORED,V2DH,C,=VD2*3600.0!ZSLFIN)
 SPEDAT(SET,MAKE,P2DFR,C,=0.0)
 SPEDAT(SET,MAKE,SUMVD2,C,=0.0)
 SPEDAT(SET,MAKE,SUMC7,C,=0.0)
 SPEDAT(SET,STORED,P2DFR,C,=ABS(NETS(C7&PW4))/NETS(C7&INL)!ZSLFIN)
 SPEDAT(SET,STORED,SUMVD2!PW4,C,=SUM(V2DH*AHIGH)!ZSLFIN)
 SPEDAT(SET,STORED,SUMC7!PW4,C,=SUM(C7*AHIGH)!ZSLFIN)
 SPEDAT(SET,PRINT,PAR2_DFRAC,C,=P2DFR)
 SPEDAT(SET,PRINT,FLOOR_AVG_VDEP2,C,=SUMVD2/SUMAF)
 SPEDAT(SET,PRINT,FLOOR_AVG_C7,C,=SUMC7/SUMAF)
 SPEDAT(SET,STORED,V3DH,C,=VD3*3600.0!ZSLFIN)
 SPEDAT(SET,MAKE,P3DFR,C,=0.0)
 SPEDAT(SET,MAKE,SUMVD3,C,=0.0)
 SPEDAT(SET,MAKE,SUMC8,C,=0.0)
 SPEDAT(SET,STORED,P3DFR,C,=ABS(NETS(C8&PW4))/NETS(C8&INL)!ZSLFIN)
 SPEDAT(SET,STORED,SUMVD3!PW4,C,=SUM(V3DH*AHIGH)!ZSLFIN)
 SPEDAT(SET,STORED,SUMC8!PW4,C,=SUM(C8*AHIGH)!ZSLFIN)
 SPEDAT(SET,PRINT,PAR3_DFRAC,C,=P3DFR)
 SPEDAT(SET,PRINT,FLOOR_AVG_VDEP3,C,=SUMVD3/SUMAF)
 SPEDAT(SET,PRINT,FLOOR_AVG_C8,C,=SUMC8/SUMAF)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 SPEDAT(SET,OBJNAM,!OB1,C,INL)
 SPEDAT(SET,OBJTYP,!OB1,C,INLET)
 SPEDAT(SET,OBJNAM,!OB2,C,EXHAUST)
 SPEDAT(SET,OBJTYP,!OB2,C,INLET)
 SPEDAT(SET,SOURCE,R1!OB2,C,=-0.15*DEN1!IMAT<100!VOLU)
 SPEDAT(SET,OBJNAM,!OB3,C,NW)
 SPEDAT(SET,OBJTYP,!OB3,C,PLATE)
 SPEDAT(SET,OBJNAM,!PW3,C,NW)
 SPEDAT(SET,OBJTYP,!PW3,C,USER_DEFINED)
 SPEDAT(SET,OBJNAM,!OB4,C,LW1)
 SPEDAT(SET,OBJTYP,!OB4,C,PLATE)
 SPEDAT(SET,OBJNAM,!PW4,C,LW1)
 SPEDAT(SET,OBJTYP,!PW4,C,USER_DEFINED)
 SPEDAT(SET,OBJNAM,!OB5,C,HW1)
 SPEDAT(SET,OBJTYP,!OB5,C,PLATE)
 SPEDAT(SET,OBJNAM,!PW5,C,HW1)
 SPEDAT(SET,OBJTYP,!PW5,C,USER_DEFINED)
 SPEDAT(SET,OBJNAM,!OB6,C,EW)
 SPEDAT(SET,OBJTYP,!OB6,C,PLATE)
 SPEDAT(SET,OBJNAM,!PW6,C,EW)
 SPEDAT(SET,OBJTYP,!PW6,C,USER_DEFINED)
 SPEDAT(SET,OBJNAM,!OB7,C,WW)
 SPEDAT(SET,OBJTYP,!OB7,C,PLATE)
 SPEDAT(SET,ROUGHNESS,OB7,R,0.)
 SPEDAT(SET,OBJNAM,!PW7,C,WW)
 SPEDAT(SET,OBJTYP,!PW7,C,USER_DEFINED)
 SPEDAT(SET,OBJNAM,!OB8,C,NOPEN)
 SPEDAT(SET,OBJTYP,!OB8,C,OPENING)
 SPEDAT(SET,ARATIO,!OB8,R,1.)
 SPEDAT(SET,OBJNAM,!OB9,C,SW)
 SPEDAT(SET,OBJTYP,!OB9,C,PLATE)
 SPEDAT(SET,OBJNAM,!PW9,C,SW)
 SPEDAT(SET,OBJTYP,!PW9,C,USER_DEFINED)
 SPEDAT(SET,FACETDAT,NUMOBJ,I,9)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=0.0208 ;EX(U1)=1.545E-03
 EX(V1)=2.758E-03 ;EX(W1)=4.951E-03
 EX(KE)=1.535E-05 ;EX(EP)=3.061E-06
 EX(C6)=0.09899 ;EX(C7)=0.1749
 EX(C8)=0.3205 ;EX(EPKE)=0.
 EX(V3DH)=0.01089 ;EX(VD3)=3.025E-06
 EX(VD2)=2.583E-05 ;EX(VS3)=1.957E-04
 EX(VS1)=4.593E-03 ;EX(VDPH)=1.0E-10
 EX(DI1)=2.433E+05 ;EX(VP1)=0.08505
 EX(VSTR)=1.885E-04 ;EX(TP1)=1.639E-05
 EX(V2DH)=0.09299 ;EX(PRPS)=1.
 EX(VABS)=7.045E-03 ;EX(YPLS)=0.4537
 EX(VD1)=7.098E-05 ;EX(STRS)=3.501E-07
 EX(V1DH)=0.2555 ;EX(ENUL)=1.0E-05
 EX(SPH1)=0. ;EX(VS2)=1.671E-03
 EX(TEM1)=20. ;EX(DEN1)=1.205
 EX(EL1)=0.02831 ;EX(ENUT)=4.428E-05
 ************************************************************
  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,N,N,Y,Y,Y,Y)
 OUTPUT(U1,N,N,Y,Y,Y,Y)
 OUTPUT(V1,N,N,Y,Y,Y,Y)
 OUTPUT(W1,N,N,Y,Y,Y,Y)
 OUTPUT(KE,N,N,Y,Y,Y,Y)
 OUTPUT(EP,N,N,Y,Y,Y,Y)
 OUTPUT(C6,N,N,Y,Y,Y,Y)
 OUTPUT(C7,N,N,Y,Y,Y,Y)
 OUTPUT(C8,N,N,Y,Y,Y,Y)
 OUTPUT(EPKE,N,N,Y,N,N,N)
 OUTPUT(V3DH,N,N,Y,N,N,N)
 OUTPUT(VD3,N,N,Y,N,N,N)
 OUTPUT(VD2,N,N,Y,N,N,N)
 OUTPUT(VS3,N,N,Y,N,N,N)
 OUTPUT(VS1,N,N,Y,N,N,N)
 OUTPUT(VDPH,N,N,Y,N,N,N)
 OUTPUT(DI1,N,N,Y,N,N,N)
 OUTPUT(VP1,N,N,Y,N,N,N)
 OUTPUT(VSTR,N,N,Y,N,N,N)
 OUTPUT(TP1,N,N,Y,N,N,N)
 OUTPUT(V2DH,N,N,Y,N,N,N)
 OUTPUT(PRPS,N,N,Y,N,N,N)
 OUTPUT(VABS,N,N,Y,N,N,N)
 OUTPUT(YPLS,N,N,Y,N,N,N)
 OUTPUT(VD1,N,N,Y,N,N,N)
 OUTPUT(STRS,N,N,Y,N,N,N)
 OUTPUT(V1DH,N,N,Y,N,N,N)
 OUTPUT(ENUL,N,N,Y,N,N,N)
 OUTPUT(SPH1,N,N,Y,N,N,N)
 OUTPUT(VS2,N,N,Y,N,N,N)
 OUTPUT(TEM1,N,N,Y,Y,Y,Y)
 OUTPUT(DEN1,N,N,Y,N,N,N)
 OUTPUT(EL1,N,N,Y,N,N,N)
 OUTPUT(ENUT,N,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 16 ;IYMON = 15 ;IZMON = 47
 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
 NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
 XZPR = F ;YZPR = F
 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
 IDISPA = 2000 ;IDISPB = 0 ;IDISPC = 0
 CSG1    ='SW'
STOP