PHOTON USE
  p
 
 
 
 
  up -x
  set con fi dep 4
 
  vec z 1 sh
  t
  2
 
  Figure 1.  Velocity Vectors.
    .6E+02  .299E+04 CR
  pause
  vec cl;te cl;redr
  con tem1 z 1 fill
  te
  2
 
  Figure 2.  Temperature Contours (K).
    .6E+02  .299E+04 CR
  pause
  cl;text
  2
 
  LAMP HEATER
    .80874E+03  .29346E+04  C
  WAFER
    .88921E+03  .13181E+04  C
  inlet
    .94956E+03  .27209E+04  C
  outlet1
    .14485E+03  .43938E+03  C
  outlet2
    .16497E+04  .44341E+03  C
  Figure 3.  Reactor Geometry.
    .6E+02  .299E+04 CR
  use patgeo
  pause
  ENDUSE
    GROUP 7. Variables stored, solved & named
SOLVE(P1,V1,U1);SOLUTN(P1,Y,Y,Y,P,P,P)
STORE(DEN1,ENUL,KOND,SPH1)
STORE(PRPS)
SOLVE(TEM1);SOLUTN(TEM1,Y,Y,Y,P,P,P)
    GROUP 9. Properties of the medium (or media)
RHOIN = 1.189
PRESS0=1.E5
TEMP0=0.0 ! TEM1 in deg K
ENUT = 0.0
WVELIN = 0.01
    GROUP 10. Inter-phase-transfer processes and properties
    GROUP 11. Initialization of variable or porosity fields
FIINIT(PRPS) = 0.0
CONPOR(TOPL1,0.0,CELL,1,-NX1,1,-(NY/2-NY1),1,NZ)
CONPOR(TOPS1,0.0,CELL,NX1+1,-NX2A,-(NY/2-NY2A+1),-(NY/2-NY1),1,NZ)
CONPOR(TOPL2,0.0,CELL,1,-NX1,-(NY/2+NY1+1),NY,1,NZ)
CONPOR(TOPS2,0.0,CELL,NX1+1,-NX2A,-(NY/2+NY1+1),-(NY/2+NY2A),1,NZ)
CONPOR(BOTL,0.0,CELL,-(NX4A+1),NX,-(NY5+1),-(NY-NY5),1,NZ)
CONPOR(BOTS,0.0,CELL,-(NX3A+1),NX4A,-(NY5+NY4+1),-(NY/2+NY3A),1,NZ)
    GROUP 13. Boundary conditions and special sources
INLET(PLUG,WEST,1,1,NY/2-NY1+1,NY/2+NY1,1,NZ,1,LSTEP)
VALUE(PLUG,P1,WVELIN*RHOIN)
VALUE(PLUG,V1,0.0        )
VALUE(PLUG,U1,WVELIN     )
VALUE(PLUG,TEM1,GTIN     )
  ** outlets
PATCH(EXTP1,EAST ,NX,NX,1,NY5  ,1   ,NZ,1,LSTEP)
COVAL(EXTP1,P1,1000.0,0.0        )
PATCH(EXTP2,EAST ,NX,NX,NY/2+NY4A+1,NY  ,1   ,NZ,1,LSTEP)
COVAL(EXTP2,P1,1000.0,0.0        )
  ** wafer surface concentration b/c
PATCH(WAFCON,EAST,NX3A,NX3A,NY/2-NY3A+1,NY/2+NY3A,1,NZ,1,LSTEP)
PATCH(SIDE1,NWALL,NX1+1,NX,NY    ,NY  ,1,NZ,1,LSTEP)
COVAL(SIDE1,U1,1.0   ,0.0        )
PATCH(SIDE2,SWALL,NX1+1,NX,1     ,1   ,1,NZ,1,LSTEP)
COVAL(SIDE2,U1,1.0   ,0.0        )
  Radiation boundary conditions
  inlet patch radiating aperture temperature gtin(K)
patch(@RI001,west,1,1,ny/2-ny1+1,ny/2+ny1,1,1,nz,lstep)
coval(@RI001,tem1,0.0,gtin)
  reactor walls and susceptor plate - adiabatic
patch(@RI002,north,1,nx1,ny/2+ny1,ny/2+ny1,1,nz,1,lstep)
coval(@RI002,tem1,grnd1,0.0)
patch(@RI003,north,nx1+1,nx2a,ny/2+ny1,ny/2+ny1,1,nz,1,lstep)
coval(@RI003,tem1,grnd1,0.0)
patch(@RI004,west,nx2a+1,nx2a+1,ny/2+ny1+1,ny/2+ny2a,1,nz,1,lstep)
coval(@RI004,tem1,grnd1,0.0)
patch(@RI005,south,nx1+1,nx2a,ny/2+ny2a+1,ny/2+ny2a+1,1,nz,1,lstep)
coval(@RI005,tem1,grnd1,0.0)
patch(@RI006,west,nx1+1,nx1+1,ny/2+ny2a+1,ny,1,nz,1,lstep)
coval(@RI006,tem1,grnd1,0.0)
patch(@RI007,north,nx1+1,nx3a,ny,ny,1,nz,1,lstep)
coval(@RI007,tem1,grnd1,0.0)
patch(@RI008,north,nx3a+1,nx4a,ny,ny,1,nz,1,lstep)
coval(@RI008,tem1,grnd1,0.0)
patch(@RI009,north,nx4a+1,nx,ny,ny,1,nz,1,lstep)
coval(@RI009,tem1,grnd1,0.0)
patch(@RI010,south,nx4a+1,nx,ny/2+ny4a+1,ny/2+ny4a+1,1,nz,1,lstep)
coval(@RI010,tem1,grnd1,0.0)
patch(@RI011,east,nx4a,nx4a,ny/2+ny3a+1,ny/2+ny4a,1,nz,1,lstep)
coval(@RI011,tem1,grnd1,0.0)
patch(@RI012,south,nx3a+1,nx4a,ny/2+ny3a+1,ny/2+ny3a+1,1,nz,1,lstep)
coval(@RI012,tem1,grnd1,0.0)
patch(@RI013,south,1,nx1,ny/2-ny1+1,ny/2-ny1+1,1,nz,1,lstep)
coval(@RI013,tem1,grnd1,0.0)
patch(@RI014,south,nx1+1,nx2a,ny/2-ny1+1,ny/2-ny1+1,1,nz,1,lstep)
coval(@RI014,tem1,grnd1,0.0)
patch(@RI015,west,nx2a+1,nx2a+1,ny/2-ny2a+1,ny/2-ny1,1,nz,1,lstep)
coval(@RI015,tem1,grnd1,0.0)
patch(@RI016,north,nx1+1,nx2a,ny/2-ny2a,ny/2-ny2a,1,nz,1,lstep)
coval(@RI016,tem1,grnd1,0.0)
patch(@RI017,west,nx1+1,nx1+1,1,ny/2-ny2a,1,nz,1,lstep)
coval(@RI017,tem1,grnd1,0.0)
patch(@RI018,south,nx1+1,nx3a,1,1,1,nz,1,lstep)
coval(@RI018,tem1,grnd1,0.0)
patch(@RI019,south,nx3a+1,nx4a,1,1,1,nz,1,lstep)
coval(@RI019,tem1,grnd1,0.0)
patch(@RI020,south,nx4a+1,nx,1,1,1,nz,1,lstep)
coval(@RI020,tem1,grnd1,0.0)
patch(@RI021,north,nx4a+1,nx,ny5,ny5,1,nz,1,lstep)
coval(@RI021,tem1,grnd1,0.0)
patch(@RI022,east,nx4a,nx4a,ny5+1,ny5+ny4,1,nz,1,lstep)
coval(@RI022,tem1,grnd1,0.0)
patch(@RI023,north,nx3a+1,nx4a,ny5+ny4,ny5+ny4,1,nz,1,lstep)
coval(@RI023,tem1,grnd1,0.0)
    susceptor patches radiation flux gwaft(K)
patch(@RI024,east,nx3a,nx3a,ny/2+ny1+1,ny/2+ny3a,1,nz,1,lstep)
coval(@RI024,tem1,0.0,gwaft)
patch(@RI025,east,nx3a,nx3a,ny/2+1,ny/2+ny1,1,nz,1,lstep)
coval(@RI025,tem1,0.0,gwaft)
patch(@RI026,east,nx3a,nx3a,ny/2-ny3a+1,ny/2-ny1,1,nz,1,lstep)
coval(@RI026,tem1,0.0,gwaft)
patch(@RI027,east,nx3a,nx3a,ny/2-ny1+1,ny/2,1,nz,1,lstep)
coval(@RI027,tem1,0.0,gwaft)
  outlet patch
patch(@RI028,east,nx,nx,ny-ny5+1,ny,1,nz,1,lstep)
coval(@RI028,tem1,0.0,gtin)
patch(@RI029,east,nx,nx,1,ny5,1,nz,1,lstep)
coval(@RI029,tem1,0.0,gtin)
    GROUP 15. Termination of sweeps
LSWEEP = 200
    GROUP 16. Termination of iterations
    GROUP 17. Under-relaxation devices
DUM1  = XR1+XR2+XR3+XR4+XR5
DUM1  = DUM1/WVELIN/NZ
RELAX(P1,LINRLX,1.)
RELAX(V1,FALSDT,DUM1)
RELAX(U1,FALSDT,DUM1)
    GROUP 20. Preliminary print-out
    GROUP 21. Print-out of variables
TSTSWP = -1
  Monitor situated at stagnation region.
IYMON = NY;IXMON = NX3A;IZMON = 36
S2SR=t
SPEDAT(SET,CVD,RADCVD,L,T)
    GROUP 24. Dumps for restarts