PHOTON USE
  p;;;;
 
  msg temperature contours
  con tem1 y 1 fi;0.01
  con tem1 y 5 fi;0.01
  con tem1 y 10 fi;0.01
  con tem1 y 15 fi;0.01
  con tem1 y 20 fi;0.01
  con tem1 y 25 fi;0.01
  con tem1 y 30 fi;0.01
  pause
  msg velocity vectors at the top
  vec y m
  enduse
 ************************************************************
  Group 1. Run Title
 TEXT(Free Convectn And Radiatn Heat Sink 
 TITLE
 LIBREF=608
   DISPLAY
 
   A vertical metal plate is fitted with vertical fins, which cool
   it by free convection and radiation. The flow is steady. The
   grid size is 20 * 30 * 30; so the computer time is not small.
 
   PHOTON USE commands are included in the Q1 file.
 
   ENDDIS
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
 REAL(TOL,LX,LY,LZ,X0,Y0,Z0)
 INTEGER(NXREG,NYREG,NZREG)
 
    * Overall number of cells
 NX=30; NY=30; NZ=20; TOL=1.E-5
 RSET(M,NX,NY,NZ,TOL)
 
    * Overall domain extent
 XULAST=0.05;YVLAST=0.2;ZWLAST=0.1
 RSET(D,CHAM,XULAST,YVLAST,ZWLAST)
 
    * Set objects: name  x0       y0        z0
    *                    lx       ly        lz
  object 1
X0=0.0; Y0=0.05; Z0=0.0455; LX=0.035; LY=0.1; LZ=0.003
RSET(B,BASE,X0,Y0,Z0,LX,LY,LZ)
 
  object 2
X0=0.0; Y0=0.08; Z0=0.0485; LX=0.009; LY=0.04; LZ=0.002
RSET(B,HEATB,X0,Y0,Z0,LX,LY,LZ)
 
  object 3
X0=0.0; Y0=0.08; Z0=0.0505; LX=0.0046; LY=0.04; LZ=0.008
RSET(B,HEATT,X0,Y0,Z0,LX,LY,LZ)
 
  object 4
X0=0.0145; Y0=0.05; Z0=0.0485; LX=0.001475; LY=0.1; LZ=0.024
RSET(B,FIN1,X0,Y0,Z0,LX,LY,LZ)
 
  object 5
X0=0.0237; Y0=0.05; Z0=0.0485; LX=0.001475; LY=0.1; LZ=0.024
RSET(B,FIN2,X0,Y0,Z0,LX,LY,LZ)
 
  object 6
X0=0.0329; Y0=0.05; Z0=0.0485; LX=0.001475; LY=0.1; LZ=0.024
RSET(B,FIN3,X0,Y0,Z0,LX,LY,LZ)
 
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
#solvel
#soltem
 STORE(DEN1,BLOK)
 ************************************************************
  Group 8. Terms & Devices
 CSG3=CNGR
 ************************************************************
  Group 9. Properties
#use_props
 PRESS0  = 1.0E+05; TEMP0   = 2.730E+02
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(V1  ) =  1.0E-01 ; FIINIT(BLOK) =  1.0
 :fluid: =airideal ; FIINIT(DEN1) =  1.189E+00
 FIINIT(TEM1) =  2.0E+01
 
 CONPOR(BASE    , -1.00,CELL  ,#1,#9,#2,#4,#2,#2)
 INIT(BASE    ,BLOK, 0.0, 2.0); INIT(BASE    ,PRPS, 0.0, 1.0E+02)
 
 CONPOR(HEATT   , -1.00,CELL  ,#1,#1,#3,#3,#4,#4)
 INIT(HEATT   ,BLOK, 0.0, 3.0); INIT(HEATT   ,PRPS, 0.0, 1.0E+02)
 
 CONPOR(HEATB   , -1.00,CELL  ,#1,#2,#3,#3,#3,#3)
 INIT(HEATB   ,BLOK, 0.0, 4.0); INIT(HEATB   ,PRPS, 0.0, 1.0E+02)
 
 CONPOR(FIN1    , -1.00,CELL  ,#4,#4,#2,#4,#3,#5)
 INIT(FIN1    ,BLOK, 0.0, 5.0); INIT(FIN1    ,PRPS, 0.0, 1.0E+02)
 
 CONPOR(FIN2    , -1.00,CELL  ,#6,#6,#2,#4,#3,#5)
 INIT(FIN2    ,BLOK, 0.0, 6.0); INIT(FIN2    ,PRPS, 0.0, 1.0E+02)
 
 CONPOR(FIN3    , -1.00,CELL  ,#8,#8,#2,#4,#3,#5)
 INIT(FIN3    ,BLOK, 0.0, 7.0); INIT(FIN3    ,PRPS, 0.0, 1.0E+02)
 INIADD  =    F
 ************************************************************
  Group 13. Boundary & Special Sources
 
 PATCH (X2      ,EAST  ,nx,nx,1,ny,1,nz,1,1)
 COVAL (X2      ,P1  , 1.0, 0.0)
 COVAL (X2      ,U1  , 0.0, 0.0     )
 COVAL (X2      ,V1  , 0.0, 0.0     )
 COVAL (X2      ,W1  , 0.0, 0.0     )
 COVAL (X2      ,TEM1, 0.0, 2.0E+01)
 
 PATCH (Y1      ,SOUTH ,1,nx,1,1,1,nz,1,1)
 COVAL (Y1      ,P1  , 1.0, 0.0)
 COVAL (Y1      ,U1  , 0.0, 0.0     )
 COVAL (Y1      ,V1  , 0.0, 0.0     )
 COVAL (Y1      ,W1  , 0.0, 0.0     )
 COVAL (Y1      ,TEM1, 0.0, 2.0E+01)
 
 PATCH (Y2      ,NORTH ,1,nx,ny,ny,1,nz,1,1)
 COVAL (Y2      ,P1  , 1.0, 0.0)
 COVAL (Y2      ,U1  , 0.0, 0.0     )
 COVAL (Y2      ,V1  , 0.0, 0.0     )
 COVAL (Y2      ,W1  , 0.0, 0.0     )
 COVAL (Y2      ,TEM1, 0.0, 2.0E+01)
 
 PATCH (Z1      ,LOW   ,1,nx,1,ny,1,1,1,1)
 COVAL (Z1      ,P1  , 1.0, 0.0)
 COVAL (Z1      ,U1  , 0.0, 0.0     )
 COVAL (Z1      ,V1  , 0.0, 0.0     )
 COVAL (Z1      ,W1  , 0.0, 0.0     )
 COVAL (Z1      ,TEM1, 0.0, 2.0E+01)
 
 PATCH (Z2      ,HIGH  ,1,nx,1,ny,nz,nz,1,1)
 COVAL (Z2      ,P1  , 1.0, 0.0)
 COVAL (Z2      ,U1  , 0.0, 0.0     )
 COVAL (Z2      ,V1  , 0.0, 0.0     )
 COVAL (Z2      ,W1  , 0.0, 0.0     )
 COVAL (Z2      ,TEM1, 0.0, 2.0E+01)
 
 PATCH (HOTHEATT,VOLUME,#1,#1,#3,#3,#4,#4,1,1)
 COVAL (HOTHEATT,TEM1, FIXFLU   , 5.989E+06)
 
 PATCH (HOTHEATB,VOLUME,#1,#2,#3,#3,#3,#3,1,1)
 COVAL (HOTHEATB,TEM1, FIXFLU   , 5.989E+06)
 
#gravity
gravdir=2; rhoref=1.189
#densdiff
 
 PATCH (*RADBASE,LOW   ,#1,#9,#2,#4,#3,#3,1,1)
 COVAL (*RADBASE,TEM1, 5.670E-09, 2.0E+01)
 
 PATCH (*RADFIN ,EAST  ,#8,#8,#2,#4,#3,#5,1,1)
 COVAL (*RADFIN ,TEM1, 5.670E-09, 2.0E+01)
 
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =      100; SELREF  =    T; RESFAC  = 1.0E-02
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 RELAX(U1  ,FALSDT, 2.0E-02); RELAX(V1  ,FALSDT, 2.0E-02)
 RELAX(W1  ,FALSDT, 2.0E-02); RELAX(TEM1,FALSDT, 1.0E+02)
 ************************************************************
  Group 18. Limits on variables or increments to them
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  Group 21. Print-out of Variables
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =      1 ; IYMON  =      NY-1 ; IZMON  =      11
 TSTSWP  =   -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NXPRIN  =  6; NYPRIN  =  6; NZPRIN  =  4; NPLT    =  1
 ************************************************************
  Group 24. Dumps For Restarts