**** PIL macro, introducing IMMERSOL library csaes ****
TEXT(Gateway to IMMERSOL library cases
integer(choice)
char(runcase,blank)
blank=
#cls
label zero
TITLE
  DISPLAY
  The IMMERSOL-using cases in the PHOENICS input library include:
  
  1. One-dimensional cases having exact solutions
  2. Two-dimensional cases
  3. Three-dimensional cases
  4. A multi-physics case illustrating also the simulation of 
     stresses in solids
  5. Cases which include chemical reaction  
     
  Please press key 1, 2, 3, 4 or 5, according to your choice   
  
  
  ENDDIS
  -----------------------------------------------------------------
readvdu(choice,int,0)
if(choice.eq.0) then
 goto end
if(choice.eq.1) then
 goto one
else  
if(choice.eq.2) then
 goto two
else
if(choice.eq.3) then
 goto three
else
if(choice.eq.4) then
 goto four
else
if(choice.eq.5) then
 goto five
endif
endif  
endif  
endif  
endif  
  -----------------------------------------------------------------
label one
#cls
mesg(  The one-dimensional cases all concern radiative heat transfer 
mesg(  between wide uniform parallel plates, separated by a medium 
mesg(  which may participate by way of: 
mesg(     (1) emission and absorption af radiation, and
mesg(     (2) conduction.
mesg(     
mesg(  The cases in which conduction plays no part have exact 
mesg(  solutions, which are computed by the input files for 
mesg(  comparison with the PHOENICS predictions.
mesg(     
mesg(  Each case may be operated in the x, y or z direction.
mesg(  T3, the 'radiation temperature' is computed in all cases
mesg(  as are the distance between the walls, WGAP, and the LTLS
mesg(  variable from which WGAP is deduced.
#pause 
mesg(  The case numbers in question are:
mesg(  r201, in which the medium plays no part whatsoever
mesg(  r202, in which the medium, absorbs & emits; 
mesg(       It may conduct; but the conductivity is zero by default
mesg(  r203, which is like 202 but solves for enthalpy, not TEM1
mesg(  r206, in which the radiating walls are represented by patches
mesg(       rather than solids. The temperature T1 is solved for.
mesg(  r207, which is like 206, but solves for enthalpy
mesg(  r208, in which the heat flux is prescribed and the wall
mesg(       temperatures are to be computed
mesg(
goto choose
  -----------------------------------------------------------------
label two
#cls
mesg(  Two-dimensional cases are as follows:
mesg(  r204, steady channel flow, with plates in radiating-wall duct
mesg(  r209, wherein convection-cooled plates are heated by radiation
mesg(  r210, similar to 209, different geometry; TEM1 solved
mesg(  r211, similar to 210, with H1 (enthalpy) solved
mesg(  r400, a heat-treatment furnace
mesg(  r401, array of cylindrical heating elements
mesg(  r402, cooling of radiating rod by free convection
mesg(  r405, Convection affected by radiation            
mesg(  r406, Lid-driven flow in a radiating cavity       
mesg(
mesg(

goto choose
  -----------------------------------------------------------------
label three
#cls  
mesg( Three-dimensional cases are:
mesg(  r205, wherein a solid block in a hot gas stream loses heat by
mesg(       radiation to cooled walls
mesg(  r407, a heat-treatment oven
mesg(  h406, an electronics-cooling example
mesg( 
mesg( 
mesg( 
mesg( 
mesg( 
mesg( 
mesg( 
mesg( 
goto choose
  -----------------------------------------------------------------
label four
#cls
mesg(  The case in question is s400
goto choose       
  -----------------------------------------------------------------
label five
#cls
mesg(  The cases involving chemical reaction are:-
mesg(  r212, a model burner employing SCRS
mesg(  r403, Heat treatment in a direct-fired furnace
mesg(
mesg(
mesg(
mesg(
mesg(
mesg(
mesg(
mesg(
mesg(
mesg(
mesg(
goto choose
  -----------------------------------------------------------------
label choose
mesg(  Please enter the number of your choice, or s to start again
readvdu(runcase,char,s)
runcase
if(:runcase:.eq.s) then
 goto zero
else
 load($:runcase:)
endif
label end
q1quit=t 
stop