TALK=T;RUN(1,1)
  PHOTON USE
  P;;;;

  gr ou y 1
  msg             Duct geometry
  msg  Press Enter to continue
  pause;  set vec ref .04;  red;  vec y 1 sh
  msg             Xylene velocity distribution, m/s
  msg Press ENTER to continue
  pause;  vec off;red;  con H1 y 1 fil;.000001
  msg             Enthalpy contours, J/kg
  msg Press ENTER to continue
  pause;  con off;red;  con T1 y 1 fil;.000001
  msg             Temperature contours, K
  msg Press ENTER to continue
  pause;  con off;red;  con den1 y 1 fil;.000001
  msg             Density contours, kg/m^3
  msg Press ENTER to continue
  pause;  con off;red;  con EMU y 1 fil;.000001
  msg             Dynamic viscosity contours, N.s/m^2
  msg Press ENTER to continue
  pause;  con off;red;  con visl y 1 fil;.000001
  msg             Kinematic viscosity contours, m^2/s
  msg Press ENTER to continue
  pause;  con off;red;  con COND y 1 fil;.000001
  msg             Thermal conductivity contours, W/(m.K)
  msg Press ENTER to continue
  pause;  con off;red;  con CP y 1 fil;.000001
  msg             Heat capacity contours, J/(kg.K)
  msg Press E to finish
  ENDUSE
  DISPLAY
   PLANT Input Library case Z101
   This case concerns flow and heat transfer of a fluid having
   properties which depend on temperature, incuding:
   * density,
   * thermal conductivity,
   * viscosity,
   * heat capacity and
   * Prandtl number.

   All are introduced by way of PLANT statements, which link them
   with the solved-for variable H1, the enthalpy.

  ENDDIS
  PLANTBEGIN
     ** Field properties variations
RHO1=GRND
REAL(DENC0,DENC1,DENC2,DENC3)
DENC0=1.1386e3;DENC1=1.0388;DENC2=5.81149E-4;DENC3=-6.7765E-7
   DEN1=:DENC0:+:DENC1:*T1+:DENC2:*T1**2+(:DENC3:)*T1**3
ENUL=GRND
    VISL=EMU/DEN1
PRNDTL(H1)=GRND
   LAMPR(H1)=EMU*CP/COND
   ** Property correlations
     * Temperature
   T1=H1/CP
     * Heat capacity
   CP=4186.8*(0.616-.0040428*T1+1.8333e-5*T1**2-2.38e-08*T1*$
*3)
     * Dynamic viscosity
   EMU=4.1868e-04*(69.647-0.5278*T1+1.358e-3*T1 **2-1.16e-06$
*T1**3)
     * Thermal conductivity
   COND=1.163*(-0.1909+.003155*T1-1.0341e-5*T1**2$
    +1.0666e-8*T1**3)
  PLANTEND
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(XYLEN FLOW IN CONV.-DIV, DUCT: Z101     )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 101
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 20
 XULAST =5.0E-03
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 1
 YVLAST =1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 30
 ZWLAST =0.026
 ************************************************************
  Group 6. Body-Fitted Coordinates
 BFC = T ;NONORT = T
 NCRT = 1
 RSTGEO = F ;SAVGEO = F
 UUP = F ;VUP = F ;WUP = F
 NGEOM =CHAM
 NAMXYZ =CHAM
 ANGMIN =20.
 DOMAIN(   1,  21,   1,   2,   3,  31)
    * Set fixed sub-domain
 FIXDOM(   1,   0,   0,   0,   0,   0,   0)
 FIXDOM(   2,   0,   0,   0,   0,   0,   0)
 FIXDOM(   3,   0,   0,   0,   0,   0,   0)
 FIXDOM(   4,   0,   0,   0,   0,   0,   0)
 FIXDOM(   5,   0,   0,   0,   0,   0,   0)
 FIXDOM(   6,   0,   0,   0,   0,   0,   0)
 FIXDOM(   7,   0,   0,   0,   0,   0,   0)
 FIXDOM(   8,   0,   0,   0,   0,   0,   0)
 FIXDOM(   9,   0,   0,   0,   0,   0,   0)
 FIXDOM(  10,   0,   0,   0,   0,   0,   0)
 LIJ = F ;LJK = T ;LIK = F
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(1)=P1 ;NAME(3)=U1
 NAME(7)=W1 ;NAME(14)=H1
 NAME(142)=CP ;NAME(143)=COND
 NAME(144)=EMU ;NAME(145)=T1
 NAME(146)=VISL ;NAME(147)=DEN1
 NAME(148)=WCRT ;NAME(149)=VCRT
 NAME(150)=UCRT
    * Y in SOLUTN argument list denotes:
    * 1-stored 2-solved 3-whole-field
    * 4-point-by-point 5-explicit 6-harmonic averaging 
 SOLUTN(P1,Y,Y,Y,N,N,N)
 SOLUTN(U1,Y,Y,N,N,N,Y)
 SOLUTN(W1,Y,Y,N,N,N,Y)
 SOLUTN(H1,Y,Y,N,N,N,Y)
 SOLUTN(CP,Y,N,N,N,N,Y)
 SOLUTN(COND,Y,N,N,N,N,Y)
 SOLUTN(EMU,Y,N,N,N,N,Y)
 SOLUTN(T1,Y,N,N,N,N,Y)
 SOLUTN(VISL,Y,N,N,N,N,Y)
 SOLUTN(DEN1,Y,N,N,N,N,Y)
 SOLUTN(WCRT,Y,N,N,N,N,N)
 SOLUTN(VCRT,Y,N,N,N,N,N)
 SOLUTN(UCRT,Y,N,N,N,N,N)
 DEN1 = 147
 VISL = 146
 ************************************************************
  Group 8. Terms & Devices
    * Y in TERMS argument list denotes:
    * 1-built-in source 2-convection 3-diffusion 4-transient
    * 5-first phase variable 6-interphase transport         
 TERMS(P1,Y,Y,Y,N,Y,Y)
 TERMS(U1,Y,Y,Y,Y,Y,Y)
 TERMS(W1,Y,Y,Y,Y,Y,Y)
 TERMS(H1,N,Y,Y,N,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWRH1 = T
 NEWENL = T
 ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
 ************************************************************
  Group 9. Properties used if PRPS is not
  stored, and where PRPS = -1.0 if it is!
 RHO1 = GRND ;TMP1 =0.
 EL1 =0.
 TSURR =0. ;TEMP0 =0.
 PRESS0 =0.
 DVO1DT =0. ;DRH1DP =0.
 RHO1A =0. ;RHO1B =0.
 RHO1C =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 ENUL = GRND ;ENUT =0.
 ENULA =0. ;ENULB =0.
 ENULC =0.
 PRNDTL(U1)=1. ;PRNDTL(W1)=1.
 PRNDTL(H1)=GRND
 PRLH1A =0. ;PRLH1B =0.
 PRLH1C =0.
 PRT(U1)=1. ;PRT(W1)=1.
 PRT(H1)=1.
 CP1 =1. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1.0E-10 ;FIINIT(U1)=1.0E-10
 FIINIT(W1)=0.1 ;FIINIT(H1)=5.5E+05
 FIINIT(CP)=1747.699951 ;FIINIT(COND)=0.1282
 FIINIT(EMU)=7.87E-04 ;FIINIT(T1)=1.0E-10
 FIINIT(VISL)=1.0E-10 ;FIINIT(DEN1)=1502.
 FIINIT(WCRT)=1.0E-10 ;FIINIT(VCRT)=1.0E-10
 FIINIT(UCRT)=1.0E-10
   No PATCHes yet used for this Group
 INIADD = F
 FSWEEP = 1
 NAMFI =CHAM
 ************************************************************
  Group 12. Patchwise adjustment of terms
  Patches for this group are printed with those
  for Group 13.
  Their names begin either with GP12 or &
 ************************************************************
  Group 13. Boundary & Special Sources
 
 PATCH(INLET ,LOW , 1, 20, 1, 1, 1, 1, 1, 1)
 COVAL(INLET ,P1 , FIXFLU ,30.039999 )
 COVAL(INLET ,W1 ,0. ,0.02 )
 COVAL(INLET ,H1 ,0. ,5.5E+05 )
 
 PATCH(OUTLET ,HIGH , 1, 20, 1, 1, 30, 30, 1, 1)
 COVAL(OUTLET ,P1 ,1000. ,0. )
 COVAL(OUTLET ,U1 ,0. ,0. )
 COVAL(OUTLET ,W1 ,0. ,0. )
 
 PATCH(WALL ,EWALL , 20, 20, 1, 1, 1, 30, 1, 1)
 COVAL(WALL ,W1 ,1. ,0. )
 COVAL(WALL ,H1 ,0.031438 ,3.34E+05 )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 120 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-04
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(U1)=10
 LITER(W1)=10 ;LITER(H1)=20
 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
 ENDIT(W1)=1.0E-03 ;ENDIT(H1)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,0.5)
 RELAX(U1,FALSDT,0.1)
 RELAX(W1,FALSDT,0.1)
 RELAX(H1,FALSDT,10.)
 RELAX(CP,LINRLX,1.)
 RELAX(COND,LINRLX,1.)
 RELAX(EMU,LINRLX,1.)
 RELAX(T1,LINRLX,1.)
 RELAX(VISL,LINRLX,1.)
 RELAX(DEN1,LINRLX,1.)
 RELAX(WCRT,LINRLX,1.)
 RELAX(VCRT,LINRLX,1.)
 RELAX(UCRT,LINRLX,1.)
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
 VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
 VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
 VARMAX(H1)=5.5E+05 ;VARMIN(H1)=3.34E+05
 VARMAX(CP)=1.0E+10 ;VARMIN(CP)=-1.0E+10
 VARMAX(COND)=1.0E+10 ;VARMIN(COND)=-1.0E+10
 VARMAX(EMU)=1.0E+10 ;VARMIN(EMU)=-1.0E+10
 VARMAX(T1)=1.0E+10 ;VARMIN(T1)=-1.0E+10
 VARMAX(VISL)=1.0E+10 ;VARMIN(VISL)=-1.0E+10
 VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=-1.0E+10
 VARMAX(WCRT)=1.0E+10 ;VARMIN(WCRT)=-1.0E+10
 VARMAX(VCRT)=1.0E+10 ;VARMIN(VCRT)=-1.0E+10
 VARMAX(UCRT)=1.0E+10 ;VARMIN(UCRT)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 NAMSAT =MOSG
 PARSOL = F
 ISG62 = 1
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=0.1158 ;EX(U1)=1.161E-03
 EX(W1)=0.02308 ;EX(H1)=5.252E+05
 EX(CP)=1723. ;EX(COND)=0.1298
 EX(EMU)=9.531E-04 ;EX(T1)=304.100006
 EX(VISL)=6.45E-07 ;EX(DEN1)=1489.
 EX(WCRT)=0.02281 ;EX(VCRT)=0.
 EX(UCRT)=1.148E-03
 ************************************************************
  Group 21. Print-out of Variables
 INIFLD = F ;SUBWGR = F
    * Y in OUTPUT argument list denotes:
    * 1-field 2-correction-eq. monitor 3-selective dumping      
    * 4-whole-field residual 5-spot-value table 6-residual table
 OUTPUT(P1,Y,N,Y,Y,Y,Y)
 OUTPUT(U1,Y,N,Y,Y,Y,Y)
 OUTPUT(W1,Y,N,Y,Y,Y,Y)
 OUTPUT(H1,Y,N,Y,Y,Y,Y)
 OUTPUT(CP,Y,N,Y,N,N,N)
 OUTPUT(COND,Y,N,Y,N,N,N)
 OUTPUT(EMU,Y,N,Y,N,N,N)
 OUTPUT(T1,Y,N,Y,N,N,N)
 OUTPUT(VISL,Y,N,Y,N,N,N)
 OUTPUT(DEN1,Y,N,Y,N,N,N)
 OUTPUT(WCRT,Y,N,Y,N,N,N)
 OUTPUT(VCRT,Y,N,Y,N,N,N)
 OUTPUT(UCRT,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 18 ;IYMON = 1 ;IZMON = 15
 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 120 ;NUMCLS = 5
 NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
 NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
 XZPR = F ;YZPR = F
 IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
 ISWPRF = 1 ;ISWPRL = 100000
 ITABL = 3 ;IPROF = 1
 ABSIZ =0.5 ;ORSIZ =0.4
 NTZPRF = 1 ;NCOLPF = 50
 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
 
 PATCH(XZ ,CONTUR, 1, 20, 1, 1, 1, 30, 1, 1)
 PLOT(XZ ,P1 ,0. ,20. )
 PLOT(XZ ,W1 ,0. ,20. )
 
 PATCH(CENTRE ,PROFIL, 1, 1, 1, 1, 1, 30, 1, 1)
 PLOT(CENTRE ,P1 ,0. ,0. )
 PLOT(CENTRE ,W1 ,0. ,0. )
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP