** LOAD(090) from the PHOENICS Input Library
    GROUP 1. Run title and other preliminaries
TEXT(printing properties of materials
TITLE
  DISPLAY
  This case computes the properties of the materials listed in
  case 089, which also contains the formulae which express them
  in terms of temperature, measured in degrees Kelvin .
  
  Line-printer plots are also provided.
  ENDDIS
PRESS0=1.E5
REAL(TMIN,TMAX)
TMIN=0.0 + 273         ! temperatures must be in degrees Kelvin
TMAX=100.0 + 273       ! if formulae of case 089 are to be used
mesgm(Choose lowest temperature of range, now :TMIN:
readvdu(tmin,real,tmin)
TMIN
mesgm(Choose highest temperature of range, now :TMAX:
readvdu(tmax,real,tmax)
TMIN
TMAX
    GROUP 2. Transience; time-step specification
    GROUP 3. X-direction grid specification
NX=11
GRDPWR(X,NX,1.0,1.0)
    GROUP 7. Variables stored, solved & named
STORE(P1,RHO1,ENUL,CP1,COND,TEM1)
    GROUP 8. Terms (in differential equations) & devices
TERMS(TEM1,N,N,Y,P,P,P)
    GROUP 9. Properties of the medium (or media)
  Choose from list provided in case 089  
#$088  
TEXT(properties of :fluid_name:
    
    GROUP 11. Initial values
  inform11begin
REAL(XG1,FACTOR)
XG1=XULAST/(2*NX)
FACTOR=1/(XULAST-2.0*XG1)
(initial var TEM1 is TMIN + (TMAX-TMIN)*(xg-xg1)*factor)
  inform11end
    
    GROUP 13. Boundary conditions and special sources
   **Enthalpy is held to zero at the low-x end,
PATCH(XEQ1,CELL,1,1,1,1,1,1,1,1) 
COVAL(XEQ1,TEM1,FIXVAL,TMIN)
PATCH(XEQNX,CELL,NX,NX,1,1,1,1,1,1) 
COVAL(XEQNX,TEM1,FIXVAL,TMAX)
    GROUP 23. Field print-out and plot control
OUTPUT(TEM1,N,N,N,N,N,N) ! iprof=3 below provides sufficient   
OUTPUT(P1,N,N,N,N,N,N)   ! print-out from profile patch 
OUTPUT(RHO1,N,N,N,N,N,N)    
OUTPUT(CP1,N,N,N,N,N,N)    
OUTPUT(ENUL,N,N,N,N,N,N)    
OUTPUT(COND,N,N,N,N,N,N)    
   **Tabulate spot-values and residuals
IPROF=3
   **Plot profiles over the length of the domain
PATCH(XWISE,PROFIL,1,NX,1,1,1,1,1,LSTEP)
PLOT(XWISE,TEM1,0.0,0.0); 
PLOT(XWISE,RHO1,0.0,0.0)
PLOT(XWISE,ENUL,0.0,0.0)
PLOT(XWISE,CP1,0.0,0.0)
PLOT(XWISE,COND,0.0,0.0)
LIBREF=090
STOP