DISPLAY
   *********************************************************
   This case concerns the prediction of air flow and temperature
   distribution in and around a closed metal box, which contains
   two heat sources. These liberate 5W and 7W respectively. The
   box is cooled from the outside by air at 20 deg C being blown
   at 1.5 m/s.
 
   The case takes about 45 min to run on a 486/50. The conjugate
   gradient solver is used to improve convergence of the
   temperature equation.
 
   The results can be plotted in PHOTON by typing 'use q1'
  ENDDIS
        PHOTON USE
        p
 
 
 
 
 
        GR OU X 1;  GR OU X M;  GR OU Y 1;  GR OU Y M
        gr ou z 4 y 1 10 x 5 20; gr ou z 12 y 1 10 x 5 20
        gr ou x 5 y 1 10 z 4 11; gr ou x 21 y 1 10 z 4 11
        gr ou z 5 y 1 9 x 6 19; gr ou z 11 y 1 9 x 6 19
        gr ou x 6 y 1 9 z 5 10; gr ou x 20 y 1 9 z 5 10
        GR OU X    10 Y     2     3 Z     7     8 COL     13
        GR OU X     8 Y     2     3 Z     7     8 COL     13
        GR OU Y     4 X     8     9 Z     7     8 COL     13
        GR OU Y     2 X     8     9 Z     7     8 COL     13
        GR OU X    16 Y     2     3 Z     7     8 COL     13
        GR OU X    14 Y     2     3 Z     7     8 COL     13
        GR OU Y     4 X    14    15 Z     7     8 COL     13
        GR OU Y     2 X    14    15 Z     7     8 COL     13
        GR OU X     1 Y     1    16 Z     1    15 COL    5
        GR OU X    34 Y     1    16 Z     1    15 COL    12
        msg The geometry - Two heat sources in a closed metal box
        msg                which is force-cooled
        pause
        vec z 8 sh
        con tem1 z 8 sh;int 28
        msg Vectors and temperatures on the center plane
        pause
        vi z
        msg Viewed from the side
        pause
        con cl;vec cl;mag gr 3
        0.94204E+03 0.15046E+04 CR
        set vec ref .1
        vec z 8 x 6 19 y 1 9 sh
        msg Flow field inside the closed box
        pause
        vec cl; vi 1 1 1
        set con sca pl
        con tem1 x 25 fi;.001
        msg Temperature rise in wake of box
        ENDUSE
 IRUNN   =       1 ;LIBREF =       0
 ************************************************************
  Group 1. Run Title
 TEXT(Heat Sources; Air-Cooled Box; EGWF
 TITLE
   Originated by JCL Jan 1994
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,33,16,15)
    * Overall domain extent, RSET(D,name,XULAST,YVLAST,ZWLAST)
 RSET(D,AIR,1.500E+00,3.000E-01,5.000E-01)
    * Set objects: name  x0       y0        z0
    *                    dx       dy        dz
 RSET(B,BOX1    , 2.000E-01, 0.000E+00, 1.000E-01                , $
1.500E-02, 1.000E-01, 2.500E-01)
 RSET(B,BOX2    , 2.150E-01, 0.000E+00, 1.000E-01                , $
4.000E-01, 1.000E-01, 1.500E-02)
 RSET(B,BOX3    , 2.150E-01, 0.000E+00, 3.350E-01                , $
4.000E-01, 1.000E-01, 1.500E-02)
 RSET(B,BOX4    , 6.150E-01, 0.000E+00, 1.000E-01                , $
1.500E-02, 1.000E-01, 2.500E-01)
 RSET(B,TOP     , 2.000E-01, 1.000E-01, 1.000E-01                , $
4.300E-01, 5.000E-03, 2.500E-01)
 RSET(B,CHP1    , 2.774E-01, 1.000E-02, 1.889E-01                , $
5.570E-02, 2.000E-02, 6.334E-02)
 RSET(B,CHP2    , 4.641E-01, 1.000E-02, 1.889E-01                , $
5.570E-02, 2.000E-02, 6.334E-02)
 RSET(B,FAN1    , 0.000E+00, 0.000E+00, 0.000E+00                , $
0.000E+00, 3.000E-01, 5.000E-01)
 RSET(B,OUT1    , 1.500E+00, 0.000E+00, 0.000E+00                , $
0.000E+00, 3.000E-01, 5.000E-01)
    * Modify default grid
 RSET(X,1,4,-1.500E+00)
 RSET(X,3,2,1.000E+00)
 RSET(X,4,2,1.000E+00)
 RSET(X,5,4,1.000E+00)
 RSET(X,6,2,1.000E+00)
 RSET(X,7,4,1.000E+00)
 RSET(X,9,13,1.700E+00)
 RSET(Y,2,2,1.000E+00)
 RSET(Y,3,6,1.000E+00)
 RSET(Y,5,6,1.700E+00)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 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)
    * Additional solver options
 SOLUTN(P1  ,Y,Y,Y,N,N,N)
 SOLUTN(TEM1,Y,Y,Y,N,N,Y)
 IVARBK  =      -1 ; ISOLBK =       1
 ************************************************************
  Group 8. Terms & Devices
 NEWRH1  =    T
 NEWENL  =    T
 ISOLX   =       0 ; ISOLY  =       0 ; ISOLZ  =       0
 ************************************************************
  Group 9. Properties
 RHO1    = FILE
 PRESS0  = 1.000E+05
 TEMP0   = 2.730E+02
 ENUL    = FILE    ; ENUT   = 0.000E+00
 
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(BLOK) =  1.000E+00 ; FIINIT(PRPS) =  2.000E+00
 FIINIT(TEM1) =  2.000E+01
 
 CONPOR(BOX1    , -1.00,CELL  ,-#2,-#2,-#1,-#3,-#2,-#6)
 INIT(BOX1    ,BLOK, 0.000E+00, 2.000E+00)
 INIT(BOX1    ,PRPS, 0.000E+00, 1.110E+02)
 
 CONPOR(BOX2    , -1.00,CELL  ,-#3,-#7,-#1,-#3,-#2,-#2)
 INIT(BOX2    ,BLOK, 0.000E+00, 3.000E+00)
 INIT(BOX2    ,PRPS, 0.000E+00, 1.110E+02)
 
 CONPOR(BOX3    , -1.00,CELL  ,-#3,-#7,-#1,-#3,-#6,-#6)
 INIT(BOX3    ,BLOK, 0.000E+00, 4.000E+00)
 INIT(BOX3    ,PRPS, 0.000E+00, 1.110E+02)
 
 CONPOR(BOX4    , -1.00,CELL  ,-#8,-#8,-#1,-#3,-#2,-#6)
 INIT(BOX4    ,BLOK, 0.000E+00, 5.000E+00)
 INIT(BOX4    ,PRPS, 0.000E+00, 1.110E+02)
 
 CONPOR(TOP     , -1.00,CELL  ,-#2,-#8,-#4,-#4,-#2,-#6)
 INIT(TOP     ,BLOK, 0.000E+00, 6.000E+00)
 INIT(TOP     ,PRPS, 0.000E+00, 1.110E+02)
 
 CONPOR(CHP1    , -1.00,CELL  ,-#4,-#4,-#2,-#2,-#4,-#4)
 INIT(CHP1    ,BLOK, 0.000E+00, 7.000E+00)
 INIT(CHP1    ,PRPS, 0.000E+00, 1.060E+02)
 
 CONPOR(CHP2    , -1.00,CELL  ,-#6,-#6,-#2,-#2,-#4,-#4)
 INIT(CHP2    ,BLOK, 0.000E+00, 8.000E+00)
 INIT(CHP2    ,PRPS, 0.000E+00, 1.060E+02)
 
 INIADD  =    F
 ************************************************************
  Group 12. Convection and diffusion adjustments
 ************************************************************
  Group 13. Boundary & Special Sources
 PATCH (HOTCHP1 ,VOLUME,#4,#4,#2,#2,#4,#4,#1,#1)
 COVAL (HOTCHP1 ,TEM1, FIXFLU   , 7.087E+04)
 
 PATCH (HOTCHP2 ,VOLUME,#6,#6,#2,#2,#4,#4,#1,#1)
 COVAL (HOTCHP2 ,TEM1, FIXFLU   , 9.922E+04)
 
 INLET (FAN1    ,WEST  ,#1,#1,#1,#5,#1,#7,1,1);
 VALUE (FAN1    ,P1  , 1.783E+00)
 VALUE (FAN1    ,U1  , 1.500E+00)
 VALUE (FAN1    ,TEM1, 2.000E+01)
 
 PATCH (OUT1    ,EAST  ,#9,#9,#1,#5,#1,#7,1,1)
 COVAL (OUT1    ,P1  , 1.000E+00, 0.000E+00)
 COVAL (OUT1    ,TEM1, 0.000E+00, 2.000E+01)
 
#gravity
gravdir=2; rhoref=1.189
#densdiff
   PATCH (BUOYANCY,PHASEM,1,33,1,16,1,15,1,1)
   COVAL (BUOYANCY,V1  , FIXFLU   ,DENSDIFF)
 
   BUOYA   = 0.000E+00 ; BUOYB =-9.810E+00 ; BUOYC = 0.000E+00
   BUOYD   = 1.189E+00
 EGWF    =    T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =     200
 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.000E-01)
 RELAX(V1  ,FALSDT, 1.000E-01)
 RELAX(W1  ,FALSDT, 1.000E-01)
 RELAX(TEM1,FALSDT, 1.000E+04)
 ************************************************************
  Group 18. Limits
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
        NAMGRD  =HTBX
 CSG3=CRGR
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  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  =       5 ;IZMON  =       5
 TSTSWP  =      -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NXPRIN  =       6
 NYPRIN  =       3
 NZPRIN  =       3
 NPLT    =       1
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 mesg(Do you want to see the geometry? (N/y)
 readvdu(ans,char,n)
 if(:ans:.eq.y) then
+ rset(v,0); gview(p,1,1,1); gdom(1,nx+1,1,ny+1,1,nz+1,1,0)
+ gdom(5,21,1,11,4,12,6,0); gdom(6,20,1,10,5,11,6,0)
+ gpatch(chp1,15,0); gpatch(chp2,15,0)
+ gpatch(fan1,4,0); gpatch(out1,12,0); gdraw
 endif