TALK=T;RUN(1,1)
  DISPLAY

  This case is based on a standard fire model test case used for
  the validation of CFD field models. Its primary purpose is to
  test the fire models predictive capability in predicting
  temperature and flow distributions in a small compartment
  subject to a steady non-spreading fire. This case was studied
  experimentally by Steckler et al (NBSIR 82-2520, NBS, 1982).

  A steady non-spreading fire of 100kW is located centrally in a
  compartment measuring 2.8m square by 2.18m high. The walls and
  ceiling were 0.1m thick and covered with a ceramic fibre
  insulation board to establish near steady-state conditions within
  30 minutes. The door measures 0.74m wide and 1.83m high and is
  located centrally in one of the walls. The CFD model assumes all
  the walls are adiabatic and perfect radiative reflectors. The
  fire is simulated by means of a simple volumetric heat source and
  radiative transfer is modelled by means of the P-1 radiation
  model.

  ix=11 locates mid-plane of doorway

  The total sweep number is 800 and the case is run for 600 sweeps,
  and therafter GENTRA is activated so as to represent the spray of
  a steady sprinkler system located on the ceiling above the fire.

  ENDDIS
   VRV USE
  * Start of frame
  * Setting display switches
  AXIS ON
  POSITION TITLE   2.975420E-01  8.998145E-01
  POSITION PROBE   7.671410E-01  9.276438E-02
  * Settings for current slice
  PROBE  5.040420E-01  1.389167E+00  9.259090E-01; PROBE ON
  SLICE Y
  SLICE OUTLINE ON
  * View and up directions
  VIEW -3.537522E-01 -9.002048E-01  2.539546E-01
  UP  6.479594E-02  2.472723E-01  9.667771E-01
  * View perspective
  VIEW DEPTH  3.000000E+00;VIEW TILT 0.8
  * Setting object visibility and painting status
  OBJECT SHOW TYPE BLOCKAGE
  OBJECT HIDE LIST
  LIST FIRE
  OBJECT SHOW TYPE OPENING
  OBJECT HIDE LIST
  LIST OPN2
  OBJECT PAINT TYPE OPENING OFF
  OBJECT SHOW TYPE SPRAY-HEAD
  OBJECT PAINT TYPE SPRAY-HEAD OFF
  OBJECT HIDE TYPE PLATE
  OBJECT PAINT TYPE PLATE OFF
  VARIABLE Pressure;VECTOR ON
  msg Velocity vectors in Y-plane
  PAUSE
  * Start of frame
  VARIABLE Temperature; VARIABLE RANGE  3.027381E+02  2.564551E+03
  VECTOR OFF
  CONTOUR ON
  CONTOUR OPAQUENESS 100
  msg Temperature contours in Y-plane
  PAUSE
  * Start of frame
  VARIABLE  Pressure; VECTOR OFF;CONTOUR ON
  msg Pressure contours in Y-plane
  PAUSE
  * Start of frame
  VARIABLE Temperature; VARIABLE RANGE  3.027381E+02  2.564551E+03
  VECTOR ON
  CONTOUR ON
  msg Temprature contours and velocity vectors
  PAUSE
  * Start of frame
  VARIABLE Temperature; VARIABLE RANGE  2.938006E+02  2.564551E+03
  VECTOR OFF
  CONTOUR OFF
  * Generating streamlines
  STREAM CLEAR
  STREAM ANIMATE OFF
  HISTORY READ
  GHIS
  msg Spray tracks
  PAUSE
  * Start of frame
  PROBE  1.314041E+00  2.800000E+00  0.000000E+00; PROBE ON
  SLICE X
  SLICE OUTLINE ON
  * View and up directions
  VIEW -2.541715E-01 -9.079202E-01  3.332854E-01
  UP  3.135502E-02  3.366860E-01  9.410949E-01
  * View centre
  VIEW CENTRE  1.905373E+00  1.588076E+00  1.068990E+00
  * View size
  VIEW SIZE  2.946434E+00
  * View perspective
  VIEW DEPTH  3.000000E+00;VIEW TILT 0.8
  VARIABLE Temperature; VARIABLE RANGE  2.938006E+02  2.564551E+03
  VECTOR OFF
  CONTOUR ON
  HISTORY READ
  GHIS
  msg Temperature contours and spray tracks in X-plane
  PAUSE
  * Start of frame
  POSITION CONTOURKEY   2.587322E-02  9.276438E-02
  POSITION TITLE   2.975420E-01  8.998145E-01
  POSITION PROBE   7.671410E-01  9.276438E-02
  * Settings for current slice
  PROBE  1.140410E-01  1.320001E+00  0.000000E+00; PROBE ON
  SLICE Y
  SLICE OUTLINE ON
  * View and up directions
  VIEW -8.648933E-01 -2.411327E-01  4.402422E-01
  UP  3.938385E-01  2.177698E-01  8.930104E-01
  * View centre
  VIEW CENTRE  1.905373E+00  1.588076E+00  1.068990E+00
  * View size
  VIEW SIZE  2.946434E+00
  * View perspective
  VIEW DEPTH  3.000000E+00;VIEW TILT 0.8
  VARIABLE Temperature; VARIABLE RANGE  3.051951E+02  2.650010E+03
  VECTOR OFF
  CONTOUR ON
  * Generating streamlines
  HISTORY READ
  GHIS
  msg Temperature contours and spray tracks in Y-plane
   PAUSE
  ENDUSE
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(Fire-Spray in a compartment             )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 1
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 31
 XULAST =5.
 XFRAC(1)=0.029692 ;XFRAC(2)=0.058729
 XFRAC(3)=0.087049 ;XFRAC(4)=0.114568
 XFRAC(5)=0.141181 ;XFRAC(6)=0.166745
 XFRAC(7)=0.19105 ;XFRAC(8)=0.213761
 XFRAC(9)=0.234226 ;XFRAC(10)=0.25
 XFRAC(11)=0.268 ;XFRAC(12)=0.281333
 XFRAC(13)=0.294667 ;XFRAC(14)=0.308
 XFRAC(15)=0.31 ;XFRAC(16)=0.325774
 XFRAC(17)=0.346239 ;XFRAC(18)=0.36895
 XFRAC(19)=0.393255 ;XFRAC(20)=0.418819
 XFRAC(21)=0.445432 ;XFRAC(22)=0.472951
 XFRAC(23)=0.50127 ;XFRAC(24)=0.530308
 XFRAC(25)=0.56 ;XFRAC(26)=0.58
 XFRAC(27)=0.617566 ;XFRAC(28)=0.686253
 XFRAC(29)=0.775198 ;XFRAC(30)=0.880527
 XFRAC(31)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 20
 YVLAST =2.8
 YFRAC(1)=0.06131 ;YFRAC(2)=0.122619
 YFRAC(3)=0.183929 ;YFRAC(4)=0.245238
 YFRAC(5)=0.306548 ;YFRAC(6)=0.367857
 YFRAC(7)=0.407143 ;YFRAC(8)=0.446429
 YFRAC(9)=0.478571 ;YFRAC(10)=0.514286
 YFRAC(11)=0.55 ;YFRAC(12)=0.553571
 YFRAC(13)=0.589286 ;YFRAC(14)=0.625
 YFRAC(15)=0.6875 ;YFRAC(16)=0.75
 YFRAC(17)=0.8125 ;YFRAC(18)=0.875
 YFRAC(19)=0.9375 ;YFRAC(20)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 17
 ZWLAST =2.2
 ZFRAC(1)=0.045455 ;ZFRAC(2)=0.090909
 ZFRAC(3)=0.136364 ;ZFRAC(4)=0.199587
 ZFRAC(5)=0.26281 ;ZFRAC(6)=0.326033
 ZFRAC(7)=0.389256 ;ZFRAC(8)=0.452479
 ZFRAC(9)=0.515702 ;ZFRAC(10)=0.578926
 ZFRAC(11)=0.642149 ;ZFRAC(12)=0.705372
 ZFRAC(13)=0.768595 ;ZFRAC(14)=0.831818
 ZFRAC(15)=0.893182 ;ZFRAC(16)=0.954545
 ZFRAC(17)=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(127)=MH2O ;NAME(128)=PRPS
 NAME(129)=QRZ ;NAME(130)=ENUL
 NAME(131)=QRY ;NAME(132)=TDEC
 NAME(133)=SPH1 ;NAME(134)=QRX
 NAME(135)=LTLS ;NAME(136)=WDIS
 NAME(137)=SCAT ;NAME(138)=EMIS
 NAME(139)=WGAP ;NAME(140)=T3
 NAME(141)=MASS ;NAME(142)=HEAT
 NAME(143)=MOMZ ;NAME(144)=MOMY
 NAME(145)=MOMX ;NAME(146)=TEM1
 NAME(147)=EPKE ;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(MH2O,Y,Y,Y,N,N,Y)
 SOLUTN(PRPS,Y,N,N,N,N,N)
 SOLUTN(QRZ,Y,N,N,N,N,Y)
 SOLUTN(ENUL,Y,N,N,N,N,Y)
 SOLUTN(QRY,Y,N,N,N,N,Y)
 SOLUTN(TDEC,Y,N,N,N,N,Y)
 SOLUTN(SPH1,Y,N,N,N,N,Y)
 SOLUTN(QRX,Y,N,N,N,N,Y)
 SOLUTN(LTLS,Y,Y,Y,N,N,Y)
 SOLUTN(WDIS,Y,N,N,N,N,Y)
 SOLUTN(SCAT,Y,N,N,N,N,Y)
 SOLUTN(EMIS,Y,N,N,N,N,Y)
 SOLUTN(WGAP,Y,N,N,N,N,Y)
 SOLUTN(T3,Y,Y,N,N,N,Y)
 SOLUTN(MASS,Y,N,N,N,N,Y)
 SOLUTN(HEAT,Y,N,N,N,N,Y)
 SOLUTN(MOMZ,Y,N,N,N,N,Y)
 SOLUTN(MOMY,Y,N,N,N,N,Y)
 SOLUTN(MOMX,Y,N,N,N,N,Y)
 SOLUTN(TEM1,Y,Y,Y,N,N,Y)
 SOLUTN(EPKE,Y,N,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 = 130
 VIST = 150
 LEN1 = 149
 PRPS = 128
 ************************************************************
  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(MH2O,N,Y,Y,Y,N,Y)
 TERMS(LTLS,N,N,Y,N,Y,Y)
 TERMS(T3,N,N,Y,N,Y,N)
 TERMS(TEM1,N,Y,Y,Y,Y,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 =0. ;PRESS0 =1.0E+05
 DVO1DT =3.33E-03 ;DRH1DP = GRND5
 RHO1A =287.591705 ;RHO1B =461.744507 ;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 = 0
 PRNDTL(U1)=1. ;PRNDTL(V1)=1.
 PRNDTL(W1)=1. ;PRNDTL(KE)=1.
 PRNDTL(EP)=1. ;PRNDTL(MH2O)=1.
 PRNDTL(LTLS)=1. ;PRNDTL(T3)=-0.0263
 PRNDTL(TEM1)=-0.0263
 PRT(U1)=1. ;PRT(V1)=1.
 PRT(W1)=1. ;PRT(KE)=1.
 PRT(EP)=1.314 ;PRT(MH2O)=1.
 PRT(LTLS)=1. ;PRT(T3)=1.
 PRT(TEM1)=1.
 CP1 =1004. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1. ;FIINIT(U1)=1.0E-02
 FIINIT(V1)=1.0E-02 ;FIINIT(W1)=1.0E-02
 FIINIT(KE)=1.0E-10 ;FIINIT(EP)=1.0E-10
 FIINIT(MH2O)=0. ;FIINIT(PRPS)=-1.
 FIINIT(QRZ)=1.0E-10 ;FIINIT(ENUL)=1.0E-10
 FIINIT(QRY)=1.0E-10 ;FIINIT(TDEC)=1.0E-10
 FIINIT(SPH1)=1.0E-10 ;FIINIT(QRX)=1.0E-10
 FIINIT(LTLS)=1.0E-10 ;FIINIT(WDIS)=0.22
 FIINIT(SCAT)=0. ;FIINIT(EMIS)=0.
 FIINIT(WGAP)=1.0E-10 ;FIINIT(T3)=298.
 FIINIT(MASS)=1.0E-10 ;FIINIT(HEAT)=1.0E-10
 FIINIT(MOMZ)=1.0E-10 ;FIINIT(MOMY)=1.0E-10
 FIINIT(MOMX)=1.0E-10 ;FIINIT(TEM1)=293.
 FIINIT(EPKE)=1.0E-10 ;FIINIT(DEN1)=1.0E-10
 FIINIT(EL1)=1.0E-10 ;FIINIT(ENUT)=1.0E-10
 
 PATCH(OB1 ,INIVAL, 26, 26, 1, 6, 1, 17, 1, 1)
 INIT(OB1 ,PRPS,0. ,199. )
 
 PATCH(OB2 ,INIVAL, 26, 26, 15, 20, 1, 17, 1, 1)
 INIT(OB2 ,PRPS,0. ,199. )
 
 PATCH(OB3 ,INIVAL, 26, 26, 7, 14, 15, 17, 1, 1)
 INIT(OB3 ,PRPS,0. ,199. )
 
 PATCH(OBC ,INIVAL, 11, 15, 9, 12, 1, 3, 1, 1)
 INIT(OBC ,P1 ,0. ,0. )
 INIT(OBC ,MH2O,0. ,0. )
 INIT(OBC ,PRPS,0. ,-1. )
 INIADD = F
 FSWEEP = 1
 NAMFI =stdy
 ************************************************************
  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(BUOYANCY,PHASEM, 1, 31, 1, 20, 1, 17, 1, 1)
 COVAL(BUOYANCY,U1 , FIXFLU , GRND2 )
 COVAL(BUOYANCY,V1 , FIXFLU , GRND2 )
 COVAL(BUOYANCY,W1 , FIXFLU , GRND2 )
 
 PATCH(GENPAT ,CELL , 1, 31, 1, 20, 1, 17, 1, 1)
 COVAL(GENPAT ,U1 , FIXFLU , GRND )
 COVAL(GENPAT ,V1 , FIXFLU , GRND )
 COVAL(GENPAT ,W1 , FIXFLU , GRND )
 COVAL(GENPAT ,MH2O, FIXFLU , GRND )
 COVAL(GENPAT ,TEM1, FIXFLU , GRND )
 
 PATCH(GENMAS ,CELL , 1, 31, 1, 20, 1, 17, 1, 1)
 COVAL(GENMAS ,P1 , FIXFLU , GRND )
 
 PATCH(KESOURCE,PHASEM, 1, 31, 1, 20, 1, 17, 1, 1)
 COVAL(KESOURCE,KE , GRND4 , GRND4 )
 COVAL(KESOURCE,EP , GRND4 , GRND4 )
 
 PATCH(GX4 ,NORTH , 27, 31, 20, 20, 1, 17, 1, 1)
 COVAL(GX4 ,P1 ,1. ,0. )
 COVAL(GX4 ,U1 ,0. ,0. )
 COVAL(GX4 ,V1 ,0. ,0. )
 COVAL(GX4 ,W1 ,0. ,0. )
 COVAL(GX4 ,KE ,0. ,1.0E-10 )
 COVAL(GX4 ,EP ,0. ,1.0E-10 )
 COVAL(GX4 ,MH2O,0. ,0. )
 COVAL(GX4 ,TEM1,0. ,303. )
 
 PATCH(*RAD4 ,NORTH , 27, 31, 20, 20, 1, 17, 1, 1)
 COVAL(*RAD4 ,T3 ,5.6697E-08 ,300. )
 
 PATCH(GX5 ,SOUTH , 27, 31, 1, 1, 1, 17, 1, 1)
 COVAL(GX5 ,P1 ,1. ,0. )
 COVAL(GX5 ,U1 ,0. ,0. )
 COVAL(GX5 ,V1 ,0. ,0. )
 COVAL(GX5 ,W1 ,0. ,0. )
 COVAL(GX5 ,KE ,0. ,1.0E-10 )
 COVAL(GX5 ,EP ,0. ,1.0E-10 )
 COVAL(GX5 ,MH2O,0. ,0. )
 COVAL(GX5 ,TEM1,0. ,303. )
 
 PATCH(*RAD5 ,SOUTH , 27, 31, 1, 1, 1, 17, 1, 1)
 COVAL(*RAD5 ,T3 ,5.6697E-08 ,300. )
 
 PATCH(GX6 ,EAST , 31, 31, 1, 20, 1, 17, 1, 1)
 COVAL(GX6 ,P1 ,1. ,0. )
 COVAL(GX6 ,U1 ,0. ,0. )
 COVAL(GX6 ,V1 ,0. ,0. )
 COVAL(GX6 ,W1 ,0. ,0. )
 COVAL(GX6 ,KE ,0. ,1.0E-10 )
 COVAL(GX6 ,EP ,0. ,1.0E-10 )
 COVAL(GX6 ,MH2O,0. ,0. )
 COVAL(GX6 ,TEM1,0. ,303. )
 
 PATCH(*RAD6 ,EAST , 31, 31, 1, 20, 1, 17, 1, 1)
 COVAL(*RAD6 ,T3 ,5.6697E-08 ,300. )
 
 PATCH(GX7 ,HIGH , 27, 31, 1, 20, 17, 17, 1, 1)
 COVAL(GX7 ,P1 ,1. ,0. )
 COVAL(GX7 ,U1 ,0. ,0. )
 COVAL(GX7 ,V1 ,0. ,0. )
 COVAL(GX7 ,W1 ,0. ,0. )
 COVAL(GX7 ,KE ,0. ,1.0E-10 )
 COVAL(GX7 ,EP ,0. ,1.0E-10 )
 COVAL(GX7 ,MH2O,0. ,0. )
 COVAL(GX7 ,TEM1,0. ,303. )
 
 PATCH(*RAD7 ,HIGH , 27, 31, 1, 20, 17, 17, 1, 1)
 COVAL(*RAD7 ,T3 ,5.6697E-08 ,300. )
 
 PATCH(IM8 ,NWALL , 1, 25, 20, 20, 1, 17, 1, 1)
 COVAL(IM8 ,U1 , GRND2 ,0. )
 COVAL(IM8 ,W1 , GRND2 ,0. )
 COVAL(IM8 ,KE , GRND2 , GRND2 )
 COVAL(IM8 ,EP , GRND2 , GRND2 )
 COVAL(IM8 ,T3 , GRND4 ,300. )
 COVAL(IM8 ,TEM1, GRND2 ,300. )
 
 PATCH(IM9 ,SWALL , 1, 25, 1, 1, 1, 17, 1, 1)
 COVAL(IM9 ,U1 , GRND2 ,0. )
 COVAL(IM9 ,W1 , GRND2 ,0. )
 COVAL(IM9 ,KE , GRND2 , GRND2 )
 COVAL(IM9 ,EP , GRND2 , GRND2 )
 COVAL(IM9 ,T3 , GRND4 ,300. )
 COVAL(IM9 ,TEM1, GRND2 ,300. )
 
 PATCH(IMA ,WWALL , 1, 1, 1, 20, 1, 17, 1, 1)
 COVAL(IMA ,V1 , GRND2 ,0. )
 COVAL(IMA ,W1 , GRND2 ,0. )
 COVAL(IMA ,KE , GRND2 , GRND2 )
 COVAL(IMA ,EP , GRND2 , GRND2 )
 COVAL(IMA ,T3 , GRND4 ,300. )
 COVAL(IMA ,TEM1, GRND2 ,300. )
 
 PATCH(IMB ,LWALL , 1, 31, 1, 20, 1, 1, 1, 1)
 COVAL(IMB ,U1 , GRND2 ,0. )
 COVAL(IMB ,V1 , GRND2 ,0. )
 COVAL(IMB ,KE , GRND2 , GRND2 )
 COVAL(IMB ,EP , GRND2 , GRND2 )
 COVAL(IMB ,T3 , GRND4 ,300. )
 COVAL(IMB ,TEM1, GRND2 ,300. )
 
 PATCH(OCC ,VOLUME, 11, 15, 9, 12, 1, 3, 1, 1)
 COVAL(OCC ,TEM1,In-Form:source - see Grp 19)
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 BUOYA =0. ; BUOYB =0. ; BUOYC =-9.81
 BUOYD =1.2763
 BUOYE =27.
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 800 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-03
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=50 ;LITER(U1)=10
 LITER(V1)=10 ;LITER(W1)=10
 LITER(KE)=10 ;LITER(EP)=10
 LITER(MH2O)=20 ;LITER(LTLS)=20
 LITER(T3)=20 ;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(MH2O)=1.0E-03 ;ENDIT(LTLS)=1.0E-03
 ENDIT(T3)=1.0E-03 ;ENDIT(TEM1)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(U1,FALSDT,1.)
 RELAX(V1,FALSDT,1.)
 RELAX(W1,FALSDT,1.)
 RELAX(KE,LINRLX,0.5)
 RELAX(EP,LINRLX,0.5)
 RELAX(MH2O,LINRLX,0.3)
 RELAX(PRPS,LINRLX,1.)
 RELAX(QRZ,LINRLX,1.)
 RELAX(ENUL,LINRLX,1.)
 RELAX(QRY,LINRLX,1.)
 RELAX(TDEC,LINRLX,0.3)
 RELAX(SPH1,LINRLX,1.)
 RELAX(QRX,LINRLX,1.)
 RELAX(LTLS,FALSDT,1.0E+09)
 RELAX(WDIS,LINRLX,1.)
 RELAX(SCAT,LINRLX,1.)
 RELAX(EMIS,LINRLX,1.)
 RELAX(WGAP,LINRLX,1.)
 RELAX(T3,FALSDT,1.0E+09)
 RELAX(MASS,LINRLX,0.3)
 RELAX(HEAT,LINRLX,0.3)
 RELAX(MOMZ,LINRLX,0.3)
 RELAX(MOMY,LINRLX,0.3)
 RELAX(MOMX,LINRLX,0.3)
 RELAX(TEM1,FALSDT,1.0E+09)
 RELAX(EPKE,LINRLX,1.)
 RELAX(DEN1,LINRLX,0.3)
 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)=-7.5E+04
 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(MH2O)=1. ;VARMIN(MH2O)=0.
 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
 VARMAX(QRZ)=1.0E+10 ;VARMIN(QRZ)=-1.0E+10
 VARMAX(ENUL)=1.0E+10 ;VARMIN(ENUL)=-1.0E+10
 VARMAX(QRY)=1.0E+10 ;VARMIN(QRY)=-1.0E+10
 VARMAX(TDEC)=1.0E+10 ;VARMIN(TDEC)=-1.0E+10
 VARMAX(SPH1)=1.0E+10 ;VARMIN(SPH1)=-1.0E+10
 VARMAX(QRX)=1.0E+10 ;VARMIN(QRX)=-1.0E+10
 VARMAX(LTLS)=1.0E+10 ;VARMIN(LTLS)=-1.0E+10
 VARMAX(WDIS)=1.0E+10 ;VARMIN(WDIS)=-1.0E+10
 VARMAX(SCAT)=1.0E+10 ;VARMIN(SCAT)=-1.0E+10
 VARMAX(EMIS)=1.0E+10 ;VARMIN(EMIS)=-1.0E+10
 VARMAX(WGAP)=1.0E+10 ;VARMIN(WGAP)=-1.0E+10
 VARMAX(T3)=1.0E+10 ;VARMIN(T3)=-30.
 VARMAX(MASS)=1.0E+10 ;VARMIN(MASS)=-1.0E+10
 VARMAX(HEAT)=1.0E+10 ;VARMIN(HEAT)=-1.0E+10
 VARMAX(MOMZ)=1.0E+10 ;VARMIN(MOMZ)=-1.0E+10
 VARMAX(MOMY)=1.0E+10 ;VARMIN(MOMY)=-1.0E+10
 VARMAX(MOMX)=1.0E+10 ;VARMIN(MOMX)=-1.0E+10
 VARMAX(TEM1)=3000. ;VARMIN(TEM1)=-204.75
 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
 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
 CONWIZ = T
 ISG62 = 1
 SPEDAT(SET,DOMAIN,PHASE_1_MAT,I,2)
 SPEDAT(SET,FACETDAT,NUMOBJ,I,13)
 SPEDAT(SET,GENTRA,RESTART,C,NONE)
 SPEDAT(SET,MAXINC,T3,R,100.)
 SPEDAT(SET,MAXINC,TEM1,R,100.)
 SPEDAT(SET,STORED,TDEC,C,=TEM1-273.!ZSLFIN)
 SPEDAT(SET,PROPERTY,RHO1,C,=(P1+1.0E+05)/((175.2*MH2O+286.7)*TEM1$)
 SPEDAT(SET,PROPERTY,RHO1,C,)!IMAT<100)
 SPEDAT(SET,SOURCE,TEM1!FIRE,C,=2.E6*EXP(-10.*MH2O*TIM)!VOLU)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 SPEDAT(SET,OBJNAM,!OB1,C,BLK1)
 SPEDAT(SET,OBJTYP,!OB1,C,BLOCKAGE)
 SPEDAT(SET,BLK1,MATERIAL,R,199.)
 SPEDAT(SET,OBJNAM,!OB2,C,BLK2)
 SPEDAT(SET,OBJTYP,!OB2,C,BLOCKAGE)
 SPEDAT(SET,BLK2,MATERIAL,R,199.)
 SPEDAT(SET,OBJNAM,!OB3,C,BLK3)
 SPEDAT(SET,OBJTYP,!OB3,C,BLOCKAGE)
 SPEDAT(SET,BLK3,MATERIAL,R,199.)
 SPEDAT(SET,OBJNAM,!GX4,C,OPN1)
 SPEDAT(SET,OBJTYP,!GX4,C,OPENING)
 SPEDAT(SET,OBJNAM,!*RAD4,C,OPN1)
 SPEDAT(SET,OBJTYP,!*RAD4,C,OPENING)
 SPEDAT(SET,ARATIO,!GX4,R,1.)
 SPEDAT(SET,OBJNAM,!GX5,C,OPN2)
 SPEDAT(SET,OBJTYP,!GX5,C,OPENING)
 SPEDAT(SET,OBJNAM,!*RAD5,C,OPN2)
 SPEDAT(SET,OBJTYP,!*RAD5,C,OPENING)
 SPEDAT(SET,ARATIO,!GX5,R,1.)
 SPEDAT(SET,OBJNAM,!GX6,C,OPN3)
 SPEDAT(SET,OBJTYP,!GX6,C,OPENING)
 SPEDAT(SET,OBJNAM,!*RAD6,C,OPN3)
 SPEDAT(SET,OBJTYP,!*RAD6,C,OPENING)
 SPEDAT(SET,ARATIO,!GX6,R,1.)
 SPEDAT(SET,OBJNAM,!GX7,C,OPN4)
 SPEDAT(SET,OBJTYP,!GX7,C,OPENING)
 SPEDAT(SET,OBJNAM,!*RAD7,C,OPN4)
 SPEDAT(SET,OBJTYP,!*RAD7,C,OPENING)
 SPEDAT(SET,ARATIO,!GX7,R,1.)
 SPEDAT(SET,OBJNAM,!OB8,C,WALL1)
 SPEDAT(SET,OBJTYP,!OB8,C,PLATE)
 SPEDAT(SET,EMISSIVITY,OF!IM8,R,1.)
 SPEDAT(SET,OBJNAM,!IM8,C,WALL1)
 SPEDAT(SET,OBJTYP,!IM8,C,PLATE)
 SPEDAT(SET,OBJNAM,!OB9,C,WALL2)
 SPEDAT(SET,OBJTYP,!OB9,C,PLATE)
 SPEDAT(SET,EMISSIVITY,OF!IM9,R,1.)
 SPEDAT(SET,OBJNAM,!IM9,C,WALL2)
 SPEDAT(SET,OBJTYP,!IM9,C,PLATE)
 SPEDAT(SET,OBJNAM,!OBA,C,WALL3)
 SPEDAT(SET,OBJTYP,!OBA,C,PLATE)
 SPEDAT(SET,EMISSIVITY,OF!IMA,R,1.)
 SPEDAT(SET,OBJNAM,!IMA,C,WALL3)
 SPEDAT(SET,OBJTYP,!IMA,C,PLATE)
 SPEDAT(SET,OBJNAM,!OBB,C,FLOOR)
 SPEDAT(SET,OBJTYP,!OBB,C,PLATE)
 SPEDAT(SET,EMISSIVITY,OF!IMB,R,1.)
 SPEDAT(SET,OBJNAM,!IMB,C,FLOOR)
 SPEDAT(SET,OBJTYP,!IMB,C,PLATE)
 SPEDAT(SET,OBJNAM,!OBC,C,FIRE)
 SPEDAT(SET,OBJTYP,!OBC,C,BLOCKAGE)
 SPEDAT(SET,FIRE,MATERIAL,R,-1.)
 SPEDAT(SET,OBJNAM,!OCC,C,FIRE)
 SPEDAT(SET,OBJTYP,!OCC,C,BLOCKAGE)
 SPEDAT(SET,OBJNAM,!OBD,C,SPRAY)
 SPEDAT(SET,OBJTYP,!OBD,C,USER_DEFINED)
 SPEDAT(SET,MATERIAL,199,L,T)
 LG( 20) = T
 IG( 9) = 1
 IG( 10) = 1
 IG( 14) = 1
 IG( 15) = 100
 IG( 16) = 5
 IG( 17) = 1
 IG( 18) = 600
 IG( 19) = 1
 IG( 20) = 60
 RG( 71) = GRND1
 RG( 72) =4131.799805
 RG( 73) =8.617739E-05
 RG( 74) = GRND1
 RG( 75) = GRND1
 RG( 76) = GRND1
 RG( 77) =1.
 RG( 80) =1.
 RG( 82) =4131.799805
 RG( 83) =100.
 RG( 84) =1.
 RG( 86) =18.
 RG( 87) =28.9
 RG( 88) = GRND1
 RG( 89) =1007.
 RG( 90) = GRND1
 RG( 91) =273.149994
 RG( 92) =273.149994
 RG( 93) =1.
 RG( 94) =3.335E+05
 RG( 95) = GRND1
 RG( 96) =0.0263
 RG( 97) = GRND1
 RG( 98) =-9.81
 CG(  6) ='GHIS'
 CG(  8) ='spin'
 CG(  9) ='NONE'
 CG( 10) ='NONE'
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  Group 21. Print-out of Variables
 INIFLD = F ;SUBWGR = F
    * Y in OUTPUT argument list denotes:
    * 1-field 2-correction-eq. monitor 3-selective dumping      
    * 4-whole-field residual 5-spot-value table 6-residual table
 OUTPUT(P1,Y,N,Y,Y,Y,Y)
 OUTPUT(U1,Y,N,Y,Y,Y,Y)
 OUTPUT(V1,Y,N,Y,Y,Y,Y)
 OUTPUT(W1,Y,N,Y,Y,Y,Y)
 OUTPUT(KE,Y,N,Y,Y,Y,Y)
 OUTPUT(EP,Y,N,Y,Y,Y,Y)
 OUTPUT(MH2O,Y,N,Y,Y,Y,Y)
 OUTPUT(PRPS,Y,N,Y,N,N,N)
 OUTPUT(QRZ,Y,N,Y,N,N,N)
 OUTPUT(ENUL,Y,N,N,N,N,N)
 OUTPUT(QRY,Y,N,Y,N,N,N)
 OUTPUT(TDEC,Y,N,Y,N,N,N)
 OUTPUT(SPH1,Y,N,Y,N,N,N)
 OUTPUT(QRX,Y,N,Y,N,N,N)
 OUTPUT(LTLS,Y,N,Y,Y,Y,Y)
 OUTPUT(WDIS,Y,N,N,N,N,N)
 OUTPUT(SCAT,Y,N,Y,N,N,N)
 OUTPUT(EMIS,Y,N,Y,N,N,N)
 OUTPUT(WGAP,Y,N,Y,N,N,N)
 OUTPUT(T3,Y,N,Y,Y,Y,Y)
 OUTPUT(MASS,Y,N,Y,N,N,N)
 OUTPUT(HEAT,Y,N,Y,N,N,N)
 OUTPUT(MOMZ,Y,N,Y,N,N,N)
 OUTPUT(MOMY,Y,N,Y,N,N,N)
 OUTPUT(MOMX,Y,N,Y,N,N,N)
 OUTPUT(TEM1,Y,N,Y,Y,Y,Y)
 OUTPUT(EPKE,Y,N,Y,N,N,N)
 OUTPUT(DEN1,Y,N,Y,N,N,N)
 OUTPUT(EL1,Y,N,Y,N,N,N)
 OUTPUT(ENUT,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 4 ;IYMON = 3 ;IZMON = 8
 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
STOP