talk=f;run(1,1)
  DISPLAY

  Library case H108:

  ELECTRONIC SYSTEM RACK

  This case represents a real-life sytem rack.  The model is
  set up of 17 plates, 6 solids, 2 apertures and  two  fans;
  fans are part of the equipment within the rack. Turbulence
  model used is k-e model. System  units  generate heat and
  the model is used to check maximum temperature.

  Current grid set-up is nx=20, ny=48, nz=18. 
  However, when  you  run  this case you have the
  opportunity to reduce cell number by a factor of two.

  (Use PHOTON to display  results.  Enter  USE Q1  in PHOTON
  command mode,  or create a file called U,  which  contains
  the one line: USE Q1   before typing  RUNPHO.

  Note that USE file contains commands corresponding to  the
  finer grid.)

  ENDDIS
    *************************************************
  PHOTON USE
  p;;;;

  use patgeo
  msg Geometry of the problem
  msg
  use hotspot
  msg
  msg Press return to plot temperature field
  msg
  pause
  cont tem1 z 8 fil;.001
  msg Press return to plot velocity field
  pause
  vec x 10 sh
  msg Press E to end
  ENDUSE
 ************************************************************
  Group 1. Run Title
 TEXT(Electronic System Rack
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
mesg( Current cell set-up: nx=20, ny=48, nz=18
mesg( Do you want this reduced by a factor of two?
mesg( (Note that if coarser grid is used the accuracy of
mesg(  results may not be very good.)  (y/n)
ans=y
readvdu(ans,char,n)
if(:ans:.eq.y) then
 RSET(M,10,24,9)
else
 RSET(M,20,48,18)
endif
    * Set overall domain extent:
    *        xulast  yvlast  zwlast    name
 XSI= 6.000E-01;YSI= 1.847E+00;ZSI= 8.000E-01;RSET(D,CABT    )
    * Set objects: x0     y0     z0
    *              dx     dy     dz    name
 XPO= 7.500E-02;YPO= 5.000E-02;ZPO= 8.000E-01
 XSI= 4.500E-01;YSI= 2.250E-01;ZSI= 0.000E+00;RSET(B,AIR     )
 XPO= 7.500E-02;YPO= 2.750E-01;ZPO= 8.000E-01
 XSI= 4.500E-01;YSI= 4.500E-02;ZSI= 0.000E+00;RSET(B,BP1     )
 XPO= 7.500E-02;YPO= 3.200E-01;ZPO= 7.000E-01
 XSI= 4.500E-01;YSI= 9.000E-02;ZSI= 1.000E-01;RSET(B,SS      )
 XPO= 7.500E-02;YPO= 4.100E-01;ZPO= 8.000E-01
 XSI= 4.500E-01;YSI= 4.500E-02;ZSI= 0.000E+00;RSET(B,BP2     )
 XPO= 7.500E-02;YPO= 4.550E-01;ZPO= 4.800E-01
 XSI= 0.000E+00;YSI= 2.250E-01;ZSI= 3.200E-01;RSET(B,FDC1    )
 XPO= 7.500E-02;YPO= 4.550E-01;ZPO= 8.000E-01
 XSI= 4.500E-01;YSI= 2.250E-01;ZSI= 0.000E+00;RSET(B,FDC5    )
 XPO= 5.250E-01;YPO= 4.550E-01;ZPO= 4.800E-01
 XSI= 0.000E+00;YSI= 2.250E-01;ZSI= 3.200E-01;RSET(B,FDC2    )
 XPO= 7.500E-02;YPO= 4.550E-01;ZPO= 2.350E-01
 XSI= 4.500E-01;YSI= 2.250E-01;ZSI= 0.000E+00;RSET(B,FDC4    )
 XPO= 7.500E-02;YPO= 6.800E-01;ZPO= 2.350E-01
 XSI= 4.500E-01;YSI= 0.000E+00;ZSI= 5.650E-01;RSET(B,FDC3    )
 XPO= 7.500E-02;YPO= 4.550E-01;ZPO= 2.350E-01
 XSI= 4.500E-01;YSI= 0.000E+00;ZSI= 5.650E-01;RSET(B,FDC6    )
 XPO= 7.500E-02;YPO= 4.550E-01;ZPO= 2.350E-01
 XSI= 0.000E+00;YSI= 2.250E-01;ZSI= 2.450E-01;RSET(B,FDIN    )
 XPO= 5.250E-01;YPO= 4.550E-01;ZPO= 2.350E-01
 XSI= 0.000E+00;YSI= 2.250E-01;ZSI= 2.450E-01;RSET(B,FDOT    )
 XPO= 7.500E-02;YPO= 6.800E-01;ZPO= 5.590E-01
 XSI= 4.500E-01;YSI= 4.500E-02;ZSI= 2.410E-01;RSET(B,TOIN    )
 XPO= 7.500E-02;YPO= 7.250E-01;ZPO= 5.800E-01
 XSI= 4.500E-01;YSI= 1.800E-01;ZSI= 2.200E-01;RSET(B,PCU2    )
 XPO= 7.500E-02;YPO= 9.050E-01;ZPO= 8.000E-01
 XSI= 4.500E-01;YSI= 4.500E-02;ZSI= 0.000E+00;RSET(B,BP3     )
 XPO= 7.500E-02;YPO= 9.500E-01;ZPO= 8.000E-01
 XSI= 4.500E-01;YSI= 1.350E-01;ZSI= 0.000E+00;RSET(B,MSU5    )
 XPO= 7.500E-02;YPO= 9.500E-01;ZPO= 2.500E-01
 XSI= 0.000E+00;YSI= 1.350E-01;ZSI= 5.500E-01;RSET(B,MSU1    )
 XPO= 5.250E-01;YPO= 9.500E-01;ZPO= 2.500E-01
 XSI= 0.000E+00;YSI= 1.350E-01;ZSI= 5.500E-01;RSET(B,MSU2    )
 XPO= 7.500E-02;YPO= 9.500E-01;ZPO= 2.500E-01
 XSI= 4.500E-01;YSI= 1.350E-01;ZSI= 0.000E+00;RSET(B,MSU4    )
 XPO= 7.500E-02;YPO= 9.500E-01;ZPO= 2.500E-01
 XSI= 4.500E-01;YSI= 0.000E+00;ZSI= 5.500E-01;RSET(B,MSU3    )
 XPO= 7.500E-02;YPO= 1.085E+00;ZPO= 2.500E-01
 XSI= 4.500E-01;YSI= 0.000E+00;ZSI= 5.500E-01;RSET(B,MSU6    )
 XPO= 7.500E-02;YPO= 1.085E+00;ZPO= 8.000E-01
 XSI= 4.500E-01;YSI= 9.000E-02;ZSI= 0.000E+00;RSET(B,BP4     )
 XPO= 7.500E-02;YPO= 1.400E+00;ZPO= 8.000E-01
 XSI= 4.500E-01;YSI= 9.000E-02;ZSI= 0.000E+00;RSET(B,BP5     )
 XPO= 7.500E-02;YPO= 1.175E+00;ZPO= 8.000E-01
 XSI= 4.500E-01;YSI= 2.250E-01;ZSI= 0.000E+00;RSET(B,RDC5    )
 XPO= 7.500E-02;YPO= 1.175E+00;ZPO= 4.800E-01
 XSI= 0.000E+00;YSI= 2.250E-01;ZSI= 3.200E-01;RSET(B,RDC1    )
 XPO= 7.500E-02;YPO= 1.175E+00;ZPO= 2.350E-01
 XSI= 0.000E+00;YSI= 2.250E-01;ZSI= 2.450E-01;RSET(B,RDIN    )
 XPO= 7.500E-02;YPO= 1.400E+00;ZPO= 2.350E-01
 XSI= 4.500E-01;YSI= 0.000E+00;ZSI= 5.650E-01;RSET(B,RDC4    )
 XPO= 7.500E-02;YPO= 1.175E+00;ZPO= 2.350E-01
 XSI= 4.500E-01;YSI= 2.250E-01;ZSI= 0.000E+00;RSET(B,RDC3    )
 XPO= 5.250E-01;YPO= 1.175E+00;ZPO= 2.350E-01
 XSI= 0.000E+00;YSI= 2.250E-01;ZSI= 2.450E-01;RSET(B,RDOT    )
 XPO= 5.250E-01;YPO= 1.175E+00;ZPO= 4.800E-01
 XSI= 0.000E+00;YSI= 2.250E-01;ZSI= 3.200E-01;RSET(B,RDC2    )
 XPO= 7.500E-02;YPO= 1.175E+00;ZPO= 2.350E-01
 XSI= 4.500E-01;YSI= 0.000E+00;ZSI= 5.650E-01;RSET(B,RDC6    )
 XPO= 0.000E+00;YPO= 0.000E+00;ZPO= 0.000E+00
 XSI= 6.000E-01;YSI= 1.847E+00;ZSI= 8.000E-01;RSET(B,B32     )
 XPO= 0.000E+00;YPO= 0.000E+00;ZPO= 0.000E+00
 XSI= 6.000E-01;YSI= 1.847E+00;ZSI= 8.000E-01;RSET(B,B33     )
 XPO= 7.500E-02;YPO= 1.490E+00;ZPO= 5.590E-01
 XSI= 4.500E-01;YSI= 4.500E-02;ZSI= 2.410E-01;RSET(B,TIO2    )
 XPO= 7.500E-02;YPO= 1.535E+00;ZPO= 6.000E-01
 XSI= 4.500E-01;YSI= 9.000E-02;ZSI= 2.000E-01;RSET(B,UNIT    )
 XPO= 7.500E-02;YPO= 1.625E+00;ZPO= 8.000E-01
 XSI= 4.500E-01;YSI= 4.500E-02;ZSI= 0.000E+00;RSET(B,BP6     )
 XPO= 7.500E-02;YPO= 1.670E+00;ZPO= 5.500E-01
 XSI= 4.500E-01;YSI= 1.350E-01;ZSI= 2.500E-01;RSET(B,ACD     )
 XPO= 0.000E+00;YPO= 1.834E+00;ZPO= 8.000E-01
 XSI= 6.000E-01;YSI= 1.200E-02;ZSI= 0.000E+00;RSET(B,EX1     )
 XPO= 0.000E+00;YPO= 1.834E+00;ZPO= 0.000E+00
 XSI= 6.000E-01;YSI= 1.200E-02;ZSI= 0.000E+00;RSET(B,EX2     )
 XPO= 0.000E+00;YPO= 1.834E+00;ZPO= 0.000E+00
 XSI= 0.000E+00;YSI= 1.200E-02;ZSI= 8.000E-01;RSET(B,EX3     )
 XPO= 6.000E-01;YPO= 1.834E+00;ZPO= 0.000E+00
 XSI= 0.000E+00;YSI= 1.200E-02;ZSI= 8.000E-01;RSET(B,EX4     )
    * Modify default grid
 RSET(X,1,4,-1.500E+00)
 RSET(X,3,4,1.500E+00)
 RSET(Y,18,2,1.000E+00)
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 NAME( 43) =EPKE ; NAME( 44) =NPOR
 NAME( 45) =HPOR ; NAME( 46) =EPOR
 NAME( 47) =BLOK ; NAME( 48) =PRPS
 NAME( 49) =DEN1 ; NAME( 50) =TEM1
    * Solved variables list
 SOLVE(P1  ,U1  ,V1  ,W1  ,TEM1)
    * Stored variables list
 STORE(DEN1,PRPS,BLOK,EPOR,HPOR,NPOR,EPKE)
    * Additional solver options
 SOLUTN(P1  ,Y,Y,Y,N,N,N)
 SOLUTN(TEM1,Y,Y,Y,N,N,Y)
 IVARBK  =      -1 ;ISOLBK =       1
 TURMOD(KEMODL)
 ************************************************************
  Group 8. Terms & Devices
 NEWRH1  =    T; NEWENL  =    F; NEWENT  =    T
 ISOLX   =       0 ;ISOLY  =       0 ;ISOLZ  =       0
 ************************************************************
  Group 9. Properties
 RHO1    = FILE
 PRESS0  = 1.000E+05; TEMP0   = 2.730E+02
 CP1     = FILE; ENUL    = FILE; PRNDTL(TEM1) = CONDFILE
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(EPKE) =  1.000E+00 ;FIINIT(NPOR) =  1.000E+00
 FIINIT(HPOR) =  1.000E+00 ;FIINIT(EPOR) =  1.000E+00
 FIINIT(BLOK) =  1.000E+00 ;FIINIT(PRPS) =  2.000E+00
 FIINIT(TEM1) =  2.000E+01
 CONPOR(FDC1    ,  0.00,EAST  ,-#1,-#1,-#6,-#6,-#4,-#9)
 CONPOR(RDC2    ,  0.00,EAST  ,-#2,-#2,-#12,-#12,-#4,-#9)
 CONPOR(FDC2    ,  0.00,EAST  ,-#2,-#2,-#6,-#6,-#4,-#9)
 CONPOR(FDIN    ,  1.00,EAST  ,#1,#1,#6,#6,#2,#3)
 CONPOR(RDOT    ,  1.00,EAST  ,#2,#2,#12,#12,#2,#3)
 CONPOR(FDOT    ,  1.00,EAST  ,#2,#2,#6,#6,#2,#3)
 CONPOR(FDC4    ,  0.00,HIGH  ,-#2,-#2,-#6,-#6,-#1,-#1)
 CONPOR(RDC1    ,  0.00,EAST  ,-#1,-#1,-#12,-#12,-#4,-#9)
 CONPOR(RDC4    ,  0.00,NORTH ,-#2,-#2,-#12,-#12,-#2,-#9)
 CONPOR(RDC6    ,  0.00,NORTH ,-#2,-#2,-#11,-#11,-#2,-#9)
 CONPOR(FDC6    ,  0.00,NORTH ,-#2,-#2,-#5,-#5,-#2,-#9)
 CONPOR(FDC3    ,  0.00,NORTH ,-#2,-#2,-#6,-#6,-#2,-#9)
 CONPOR(MSU1    ,  0.00,EAST  ,-#1,-#1,-#10,-#10,-#3,-#9)
 CONPOR(MSU2    ,  0.00,EAST  ,-#2,-#2,-#10,-#10,-#3,-#9)
 CONPOR(MSU6    ,  0.25,NORTH ,#2,#2,#10,#10,#3,#9)
 CONPOR(MSU3    ,  0.25,NORTH ,#2,#2,#9,#9,#3,#9)
 CONPOR(MSU4    ,  0.00,HIGH  ,-#2,-#2,-#10,-#10,-#2,-#2)
 CONPOR(SS      , -1.00,CELL  ,-#2,-#2,-#4,-#4,-#9,-#9)
 INIT(SS      ,BLOK, 0.000E+00, 2.000E+00)
 INIT(SS      ,PRPS, 0.000E+00, 1.000E+02)
 CONPOR(TOIN    , -1.00,CELL  ,-#2,-#2,-#7,-#7,-#6,-#9)
 INIT(TOIN    ,BLOK, 0.000E+00, 3.000E+00)
 INIT(TOIN    ,PRPS, 0.000E+00, 1.000E+02)
 CONPOR(PCU2    , -1.00,CELL  ,-#2,-#2,-#8,-#8,-#7,-#9)
 INIT(PCU2    ,BLOK, 0.000E+00, 4.000E+00)
 INIT(PCU2    ,PRPS, 0.000E+00, 1.000E+02)
 CONPOR(TIO2    , -1.00,CELL  ,-#2,-#2,-#14,-#14,-#6,-#9)
 INIT(TIO2    ,BLOK, 0.000E+00, 5.000E+00)
 INIT(TIO2    ,PRPS, 0.000E+00, 1.000E+02)
 CONPOR(UNIT    , -1.00,CELL  ,-#2,-#2,-#15,-#15,-#8,-#9)
 INIT(UNIT    ,BLOK, 0.000E+00, 6.000E+00)
 INIT(UNIT    ,PRPS, 0.000E+00, 1.000E+02)
 CONPOR(ACD     , -1.00,CELL  ,-#2,-#2,-#17,-#17,-#5,-#9)
 INIT(ACD     ,BLOK, 0.000E+00, 7.000E+00)
 INIT(ACD     ,PRPS, 0.000E+00, 1.000E+02)
 INIADD  =    F
 ************************************************************
  Group 13. Boundary & Special Sources
 PATCH (BUOYANCY,PHASEM,#1,#NREGX,#1,#NREGY,#1,#NREGZ,#1,#NREGT)
 COVAL (BUOYANCY,V1  , FIXFLU   , DENSDIFF    )

 PATCH (HOTFDC1 ,WEST  ,#2,#2,#6,#6,#2,#9,1,1)
 COVAL (HOTFDC1 ,TEM1, FIXFLU   , 5.565E+02)

 PATCH (HOTRDC2 ,EAST  ,#2,#2,#12,#12,#2,#9,1,1)
 COVAL (HOTRDC2 ,TEM1, FIXFLU   , 5.565E+02)

 PATCH (HOTFDC2 ,EAST  ,#2,#2,#6,#6,#2,#9,1,1)
 COVAL (HOTFDC2 ,TEM1, FIXFLU   , 5.565E+02)

 PATCH (FDC5    ,HWALL ,#2,#2,#6,#6,#9,#9,1,1)
 COVAL (FDC5    ,U1  , LOGLAW    , 0.000E+00)
 COVAL (FDC5    ,V1  , LOGLAW    , 0.000E+00)
 COVAL (FDC5    ,KE  , LOGLAW    , LOGLAW    )
 COVAL (FDC5    ,EP  , LOGLAW    , LOGLAW    )

 PATCH (HOTFDC5 ,HIGH  ,#2,#2,#6,#6,#9,#9,1,1)
 COVAL (HOTFDC5 ,TEM1, FIXFLU   , 5.553E+02)

 PATCH (FFDIN   ,EAST  ,#2,#2,#6,#6,#2,#3,1,1)
 COVAL (FFDIN   ,U1  , FIXVAL   , 1.130E+00)

 PATCH (FRDIN   ,EAST  ,#1,#1,#12,#12,#2,#3,1,1)
 COVAL (FRDIN   ,U1  , FIXVAL   , 1.130E+00)

 PATCH (HOTFDC4 ,LOW   ,#2,#2,#6,#6,#2,#2,1,1)
 COVAL (HOTFDC4 ,TEM1, FIXFLU   , 5.553E+02)

 PATCH (HOTRDC1 ,WEST  ,#2,#2,#12,#12,#2,#9,1,1)
 COVAL (HOTRDC1 ,TEM1, FIXFLU   , 5.565E+02)

 PATCH (RDC5    ,HWALL ,#2,#2,#12,#12,#9,#9,1,1)
 COVAL (RDC5    ,U1  , LOGLAW    , 0.000E+00)
 COVAL (RDC5    ,V1  , LOGLAW    , 0.000E+00)
 COVAL (RDC5    ,KE  , LOGLAW    , LOGLAW    )
 COVAL (RDC5    ,EP  , LOGLAW    , LOGLAW    )

 PATCH (HOTRDC5 ,HIGH  ,#2,#2,#12,#12,#9,#9,1,1)
 COVAL (HOTRDC5 ,TEM1, FIXFLU   , 5.553E+02)

 PATCH (HOTRDC4 ,NORTH ,#2,#2,#12,#12,#2,#9,1,1)
 COVAL (HOTRDC4 ,TEM1, FIXFLU   , 5.553E+02)

 PATCH (HOTRDC6 ,SOUTH ,#2,#2,#12,#12,#2,#9,1,1)
 COVAL (HOTRDC6 ,TEM1, FIXFLU   , 5.521E+02)

 PATCH (HOTFDC6 ,SOUTH ,#2,#2,#6,#6,#2,#9,1,1)
 COVAL (HOTFDC6 ,TEM1, FIXFLU   , 5.521E+02)

 PATCH (HOTFDC3 ,NORTH ,#2,#2,#6,#6,#2,#9,1,1)
 COVAL (HOTFDC3 ,TEM1, FIXFLU   , 5.521E+02)

 PATCH (HOTMSU1 ,WEST  ,#2,#2,#10,#10,#3,#9,1,1)
 COVAL (HOTMSU1 ,TEM1, FIXFLU   , 2.859E+01)

 PATCH (HOTMSU2 ,EAST  ,#2,#2,#10,#10,#3,#9,1,1)
 COVAL (HOTMSU2 ,TEM1, FIXFLU   , 2.859E+01)

 PATCH (MSU5    ,HWALL ,#2,#2,#10,#10,#9,#9,1,1)
 COVAL (MSU5    ,U1  , LOGLAW    , 0.000E+00)
 COVAL (MSU5    ,V1  , LOGLAW    , 0.000E+00)
 COVAL (MSU5    ,KE  , LOGLAW    , LOGLAW    )
 COVAL (MSU5    ,EP  , LOGLAW    , LOGLAW    )

 PATCH (HOTMSU5 ,HIGH  ,#2,#2,#10,#10,#9,#9,1,1)
 COVAL (HOTMSU5 ,TEM1, FIXFLU   , 2.600E+01)

 PATCH (HOTMSU6 ,NORTH ,#2,#2,#10,#10,#3,#9,1,1)
 COVAL (HOTMSU6 ,TEM1, FIXFLU   , 9.540E+01)

 PATCH (HOTMSU3 ,SOUTH ,#2,#2,#10,#10,#3,#9,1,1)
 COVAL (HOTMSU3 ,TEM1, FIXFLU   , 9.540E+01)

 PATCH (HOTMSU4 ,LOW   ,#2,#2,#10,#10,#3,#3,1,1)
 COVAL (HOTMSU4 ,TEM1, FIXFLU   , 2.600E+01)

 PATCH (HOTTOIN ,VOLUME,#2,#2,#7,#7,#6,#9,#1,#1)
 COVAL (HOTTOIN ,TEM1, FIXFLU   , 3.074E+03)

 PATCH (HOTTIO2 ,VOLUME,#2,#2,#14,#14,#6,#9,#1,#1)
 COVAL (HOTTIO2 ,TEM1, FIXFLU   , 3.074E+03)

 PATCH (AIR     ,HIGH  ,#2,#2,#2,#2,#9,#9,1,1)
 COVAL (AIR     ,P1  , 1.000E+00, 0.000E+00)
 COVAL (AIR     ,KE  , 0.000E+00, SAME     )
 COVAL (AIR     ,EP  , 0.000E+00, SAME     )
 COVAL (AIR     ,TEM1, 0.000E+00, 2.000E+01)

 PATCH (EX1     ,HIGH  ,#1,#3,#19,#19,#9,#9,1,1)
 COVAL (EX1     ,P1  , 1.000E+00, 0.000E+00)
 COVAL (EX1     ,KE  , 0.000E+00, SAME     )
 COVAL (EX1     ,EP  , 0.000E+00, SAME     )
 COVAL (EX1     ,TEM1, 0.000E+00, 2.000E+01)

 PATCH (EX2     ,LOW   ,#1,#3,#19,#19,#1,#1,1,1)
 COVAL (EX2     ,P1  , 1.000E+00, 0.000E+00)
 COVAL (EX2     ,KE  , 0.000E+00, SAME     )
 COVAL (EX2     ,EP  , 0.000E+00, SAME     )
 COVAL (EX2     ,TEM1, 0.000E+00, 2.000E+01)

 PATCH (EX3     ,WEST  ,#1,#1,#19,#19,#1,#9,1,1)
 COVAL (EX3     ,P1  , 1.000E+00, 0.000E+00)
 COVAL (EX3     ,KE  , 0.000E+00, SAME     )
 COVAL (EX3     ,EP  , 0.000E+00, SAME     )
 COVAL (EX3     ,TEM1, 0.000E+00, 2.000E+01)

 PATCH (EX4     ,EAST  ,#3,#3,#19,#19,#1,#9,1,1)
 COVAL (EX4     ,P1  , 1.000E+00, 0.000E+00)
 COVAL (EX4     ,KE  , 0.000E+00, SAME     )
 COVAL (EX4     ,EP  , 0.000E+00, SAME     )
 COVAL (EX4     ,TEM1, 0.000E+00, 2.000E+01)

 PATCH (HIW     ,HWALL ,#1,#3,#3,#18,#9,#9,1,1)
 COVAL (HIW     ,U1  , LOGLAW    , 0.000E+00)
 COVAL (HIW     ,V1  , LOGLAW    , 0.000E+00)
 COVAL (HIW     ,KE  , LOGLAW    , LOGLAW    )
 COVAL (HIW     ,EP  , LOGLAW    , LOGLAW    )

 PATCH (LW      ,LWALL ,#1,#3,#1,#18,#1,#1,1,1)
 COVAL (LW      ,U1  , LOGLAW    , 0.000E+00)
 COVAL (LW      ,V1  , LOGLAW    , 0.000E+00)
 COVAL (LW      ,KE  , LOGLAW    , LOGLAW    )
 COVAL (LW      ,EP  , LOGLAW    , LOGLAW    )

 PATCH (WW      ,WWALL ,#1,#1,#1,#18,#1,#9,1,1)
 COVAL (WW      ,V1  , LOGLAW    , 0.000E+00)
 COVAL (WW      ,W1  , LOGLAW    , 0.000E+00)
 COVAL (WW      ,KE  , LOGLAW    , LOGLAW    )
 COVAL (WW      ,EP  , LOGLAW    , LOGLAW    )

 PATCH (EW      ,EWALL ,#3,#3,#1,#18,#1,#9,1,1)
 COVAL (EW      ,V1  , LOGLAW    , 0.000E+00)
 COVAL (EW      ,W1  , LOGLAW    , 0.000E+00)
 COVAL (EW      ,KE  , LOGLAW    , LOGLAW    )
 COVAL (EW      ,EP  , LOGLAW    , LOGLAW    )

 BUOYA   = 0.000E+00 ; BUOYB =-9.810E+00 ; BUOYC = 0.000E+00
 BUOYD   = 1.189E+00
 EGWF    =    T
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =      50; SELREF  =    T; RESFAC  = 1.000E-02
 ************************************************************
  Group 16. Terminate Iterations
 LITER (TEM1) =   15
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 1.000E+00); RELAX(U1  ,FALSDT, 1.635E-03)
 RELAX(V1  ,FALSDT, 1.635E-03); RELAX(W1  ,FALSDT, 1.635E-03)
 RELAX(KE  ,LINRLX, 4.000E-01); RELAX(EP  ,LINRLX, 4.000E-01)
 RELAX(TEM1,FALSDT, 1.635E+01)
 KELIN   =       1
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 NAMGRD  =HTBX
 ************************************************************
  Group 21. Print-out of Variables
 OUTPUT(BLOK,N,N,N,N,N,N); OUTPUT(PRPS,N,N,N,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =       5 ;IYMON  =      16 ;IZMON  =       7
 TSTSWP  =      -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NXPRIN  =       4; NYPRIN  =       4; NZPRIN  =       3
 NPLT    =       1; ITABL   =       3
STOP