#$r002
   GROUP 1. Run title and other preliminaries
TEXT(CVD REACTOR RADIATION EXAMPLE(BFC): R101
TITLE
  DISPLAY
    This case is identical to case R100, but it
    is set up for BFC with non-orthogonal grids.
  ENDDIS
REAL   (YR1,YR2,YR3,YR4,YR5,XR1,XR2,XR3,XR4,XR5)
INTEGER(NY1,NY2,NY3,NY4,NY5,NX1,NX2,NX3,NX4,NX5)
INTEGER(NY2A,NY3A,NY4A,NX2A,NX3A,NX4A,REFINX,REFINY)
REAL   (ZDTF,WVELIN,DUM1,RHOIN)
  ** Radiation settings
REAL(GWAFT,GTIN);GTIN=300.;GWAFT=500.
  For coarse or original grid set refine to 1.
  For fine grid set refine to 2/3/4...etc.
REFINX = 1
REFINY = 1
    GROUP 2. Transience; time-step specification
STEADY = T;PARAB = F;CARTES =T
    GROUP 3. X-direction grid specification
NZ     = 1;ZWLAST = .132
    GROUP 4. Y-direction grid specification
YR1    = 0.042;NY1 = 7*REFINY
YR2    = 0.003;NY2 = 1*REFINY
YR3    = 0.005;NY3 = 1*REFINY
YR4    = 0.015;NY4 = 3*REFINY
YR5    = 0.020;NY5 = 4*REFINY
NY     = 2*(NY1+NY2+NY3+NY4+NY5)
NY2A   = NY1 +NY2
NY3A   = NY2A+NY3
NY4A   = NY3A+NY4
NREGY  = 10
IREGY  = 1;GRDPWR(Y,-NY5,YR5,1.5)
IREGY  = 2;GRDPWR(Y,-NY4,YR4,1.5)
IREGY  = 3;GRDPWR(Y,-NY3,YR3,1.5)
IREGY  = 4;GRDPWR(Y, NY2,YR2,1.0)
IREGY  = 5;GRDPWR(Y,-NY1,YR1,1.2)
IREGY  = 6;GRDPWR(Y,-NY1,YR1,1.2)
IREGY  = 7;GRDPWR(Y, NY2,YR2,1.0)
IREGY  = 8;GRDPWR(Y,-NY3,YR3,1.5)
IREGY  = 9;GRDPWR(Y,-NY4,YR4,1.5)
IREGY  = 10;GRDPWR(Y,-NY5,YR5,1.5)
 
    GROUP 5. Z-direction grid specification
XR1    = 0.100;NX1 = 3*REFINX
XR2    = 0.015;NX2 = 1*REFINX
XR3    = 0.040;NX3 = 4*REFINX
XR4    = 0.005;NX4 = 1*REFINX
XR5    = 0.080;NX5 = 5*REFINX
NX     = NX1+NX2+NX3+NX4+NX5
NX2A   = NX1 +NX2
NX3A   = NX2A+NX3
NX4A   = NX3A+NX4
NREGX  = 5
IREGX  = 1;GRDPWR(X,-NX1,XR1,1.5)
IREGX  = 2;GRDPWR(X,-NX2,XR2,1.5)
IREGX  = 3;GRDPWR(X,-NX3,XR3,1.5)
IREGX  = 4;GRDPWR(X, NX4,XR4,1.0)
IREGX  = 5;GRDPWR(X, NX5,XR5,1.5)
 
   Use bfc for consistency test.
BFC    = T
FIXDOM(1,1,NX1+1,1,NY/2-NY1+1,1,NZ+1)
FIXDOM(2,NX1+1,NX2A+1,NY/2-NY2A+1,NY/2-NY1+1,1,NZ+1)
FIXDOM(3,1,NX1+1,NY/2+NY1+1,NY+1,1,NZ+1)
FIXDOM(4,NX1+1,NX2A+1,NY/2+NY1+1,NY/2+NY2A+1,1,NZ+1)
FIXDOM(5,NX4A+1,NX+1,NY5+1,NY-NY5+1,1,NZ+1)
FIXDOM(6,NX3A+1,NX4A+1,NY5+NY4+1,NY/2+NY3A+1,1,NZ+1)
DOMAIN(3,NX+1,1,NY+1,1,NZ+1)
MSWP=10;MAGIC(L)
LOAD($R333)