This set of PIL commands is useful for setting up radiation
  simulations based upon the six-flux model.
  It could be used as a separate file, eg with name: radidata .
  Then it can be "included" in the Q1 in group 7, by insertion of
  the line:
 
  INTRPT(R,RADIDATA)
 
  However, it has now been made into a library case, viz 493.htm
  which can be loaded into any other q1 by the lines: 
  #$493
  or
  L($493)
    -----------------------------------------------
  Declarations
REAL(RADW,ABSORB,SCAT,SIGMA,EMIW,ED2ME)
SIGMA=5.6697E-8;EMIW=0.3
RADW=SIGMA*TWALL**4;ABSORB=0.2;SCAT=0.01
RADIAT(FLUX,ABSORB,SCAT,H1)
FIINIT(RADX)=RADW;FIINIT(RADY)=RADW;FIINIT(RADZ)=RADW
ED2ME=EMIW/(2.0-EMIW)
RELAX(RADX,FALSDT,0.1);RELAX(RADY,FALSDT,1.0)
RELAX(RADZ,FALSDT,1.0)
  TEXT(Radiation Macro          
mesg(6-flux radiation model is in use
mesg(RADW    - intensity of the radiation from the wall = :radw:
mesg(EMIW    - emissivity of the wall                   = :emiw:
mesg(ABSORB  - absorption coefficient of the gases      = :absorb:
mesg(SCAT    - scattering     "       "   "    "        = :scat:
mesg(SIGMA   - Stefan-Boltzmann constant                = :sigma: