TALK=T;RUN(1,1)
  photon use
  p
  p1


  msg the 40*40*625 computational grid
  msg the 625-cell direction is normal to the screen
  gr z 1;  pause;  gr off;  red
  gr ou z 1
  msg Upstream smoke injection region and smoke contours
  gr z 1 x 18 23 y 18 23
  con h1 z 1 fi;0.1
  pause
  p
  p2


  msg Smoke contours 1/5 of distance along tunnel
  con h1 z 1 fi;0.1
  pause
  p
  p3


  msg Smoke contours 2/5 of distance along tunnel
  con h1 z 1 fi;0.1
  pause
  p
  p4


  msg Smoke contours 3/5 of distance along tunnel
  con h1 z 1 fi;0.1
  pause
  p
  p5


  msg Smoke contours 4/5 of distance along tunnel
  con h1 z 1 fi;0.1
  pause
  p
  phi


  msg Smoke contours 5/5 of distance along tunnel
  con h1 z 1 fi;0.1
  pause;con off;red
  msg contours of longitudinal velocity at tunnel exit
  con w1 z 1 fi;0.1
  pause
  con off
  red
  msg velocity vectors at the tunnel exit
  set vec ref 0.1
  vec z 1 sh
  enduse
  DISPLAY
  This case shows how the parabolic option can be enployed with
  advantage for fine-grid calculations on machines of modest size.

  The problem concerns the way in which a smoke plume rises as
  it is carried within a wind-tunnel having an imposed shear-flow
  velocity profile at inlet.

  The cross-stream grid is 40 * 40, and there are 625 z-direction
  intervals, so that one million cells feature in the calculation.

  The calculation is rather lengthy, but a tolerably good answer
  can be obtained by halving the grid fineness in every direction
  and reducing the number of iterations.

  Also, since the flow is symmetrical, only one half needs to be
  simulated; but both halves are simulated so as to demonstrate that
  a truly symmetrical solution is indeed arrived at.
  ENDDIS
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(Buoyant Plume in tunnel, 1.0E+06 Cells  )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 14
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 40
 XULAST =1.5
 XFRAC(1)=0.025 ;XFRAC(9)=0.225
 XFRAC(17)=0.425 ;XFRAC(25)=0.625
 XFRAC(33)=0.825
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 40
 YVLAST =1.
 YFRAC(1)=0.025 ;YFRAC(9)=0.225
 YFRAC(17)=0.425 ;YFRAC(25)=0.625
 YFRAC(33)=0.825
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = T
 NZ = 625
 ZWLAST =15.625
   Method of pairs used for grid setting.
 ZFRAC(2)=1.6E-03
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(1)=P1 ;NAME(3)=U1
 NAME(5)=V1 ;NAME(7)=W1
 NAME(14)=H1
    * 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,N,N,N,Y)
 SOLUTN(U1,Y,Y,N,N,N,Y)
 SOLUTN(V1,Y,Y,N,N,N,Y)
 SOLUTN(W1,Y,Y,N,N,N,Y)
 SOLUTN(H1,Y,Y,N,N,N,Y)
 ************************************************************
  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(V1,Y,Y,Y,Y,Y,Y)
 TERMS(W1,Y,Y,Y,Y,Y,Y)
 TERMS(H1,N,Y,Y,Y,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWRH1 = 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 = GRND4 ;TMP1 = GRND2
 EL1 =0.
 TSURR =0. ;TEMP0 =0.
 PRESS0 =0.
 DVO1DT =0. ;DRH1DP =0.
 RHO1A =1. ;RHO1B =-0.02
 RHO1C =0.
 TMP1A =0. ;TMP1B =0.
 TMP1C =0.
 TMP2A =0. ;TMP2B =0.
 TMP2C =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 ENUL =1.0E-05 ;ENUT =0.
 PRNDTL(U1)=1. ;PRNDTL(V1)=1.
 PRNDTL(W1)=1. ;PRNDTL(H1)=0.7
 PRT(U1)=1. ;PRT(V1)=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)=0. ;FIINIT(U1)=0.
 FIINIT(V1)=0. ;FIINIT(W1)=1.0E-10
 FIINIT(H1)=0.
   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(INLET1 ,LOW , 1, 40, 1, 8, 1, 1, 1, 1)
 COVAL(INLET1 ,P1 , FIXFLU ,0.8 )
 COVAL(INLET1 ,W1 ,0. ,0.8 )
 COVAL(INLET1 ,H1 ,0. ,0. )
 
 PATCH(INLET2 ,LOW , 1, 40, 9, 16, 1, 1, 1, 1)
 COVAL(INLET2 ,P1 , FIXFLU ,0.9 )
 COVAL(INLET2 ,W1 ,0. ,0.9 )
 
 PATCH(INLET3 ,LOW , 1, 40, 17, 24, 1, 1, 1, 1)
 COVAL(INLET3 ,P1 , FIXFLU ,1. )
 COVAL(INLET3 ,W1 ,0. ,1. )
 
 PATCH(INLET4 ,LOW , 1, 40, 25, 32, 1, 1, 1, 1)
 COVAL(INLET4 ,P1 , FIXFLU ,1.1 )
 COVAL(INLET4 ,W1 ,0. ,1.1 )
 
 PATCH(INLET5 ,LOW , 1, 40, 33, 40, 1, 1, 1, 1)
 COVAL(INLET5 ,P1 , FIXFLU ,1.2 )
 COVAL(INLET5 ,W1 ,0. ,1.2 )
 
 PATCH(SMOKE ,LOW , 18, 23, 18, 23, 1, 1, 1, 1)
 COVAL(SMOKE ,P1 , FIXFLU ,1.0E-02 )
 COVAL(SMOKE ,W1 ,0. ,1.0E-02 )
 COVAL(SMOKE ,H1 ,0. ,100. )
 
 PATCH(BUOYANCY,PHASEM, 1, 40, 1, 40, 1, 625, 1, 1)
 COVAL(BUOYANCY,V1 , FIXFLU ,-9.81 )
 
 PATCH(BOTTOM ,SWALL , 1, 40, 1, 1, 1, 625, 1, 1)
 COVAL(BOTTOM ,U1 ,1. ,0. )
 COVAL(BOTTOM ,W1 ,1. ,0. )
 
 PATCH(TOP ,NWALL , 1, 40, 40, 40, 1, 625, 1, 1)
 COVAL(TOP ,U1 ,1. ,0. )
 COVAL(TOP ,W1 ,1. ,0. )
 
 PATCH(WEST ,WWALL , 1, 1, 1, 40, 1, 625, 1, 1)
 COVAL(WEST ,V1 ,1. ,0. )
 COVAL(WEST ,W1 ,1. ,0. )
 
 PATCH(EAST ,EWALL , 40, 40, 1, 40, 1, 625, 1, 1)
 COVAL(EAST ,V1 ,1. ,0. )
 COVAL(EAST ,W1 ,1. ,0. )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 IPARAB = 0
 AZPH =0. ;PBAR =0.
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 1 ;ISWC1 = 1
 LITHYD = 100 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-04
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=100 ;LITER(U1)=10
 LITER(V1)=10 ;LITER(W1)=10
 LITER(H1)=10
 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-06
 ENDIT(V1)=1.0E-06 ;ENDIT(W1)=1.0E-03
 ENDIT(H1)=1.0E-06
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(U1,FALSDT,0.1)
 RELAX(V1,FALSDT,0.05)
 RELAX(W1,FALSDT,0.1)
 RELAX(H1,FALSDT,0.1)
 OVRRLX =1.7
 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(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
 VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
 VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 PARSOL = F
 IDISPB = 1 ;IDISPC = 625 ;IDISPD = 125
 ISG62 = 1
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  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,Y,Y,Y,Y,Y)
 OUTPUT(U1,Y,Y,Y,Y,Y,Y)
 OUTPUT(V1,Y,Y,Y,Y,Y,Y)
 OUTPUT(W1,Y,Y,Y,Y,Y,Y)
 OUTPUT(H1,Y,Y,Y,Y,Y,Y)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 20 ;IYMON = 20 ;IZMON = 1
 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -5
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NXPRIN = 8 ;IXPRF = 1 ;IXPRL = 10000
 NYPRIN = 8 ;IYPRF = 1 ;IYPRL = 10000
 NZPRIN = 625 ;IZPRF = 1 ;IZPRL = 10000
 IPLTF = 1 ;IPLTL = -1 ;NPLT = 100
 ISWPRF = 1 ;ISWPRL = 100000
 ITABL = 3 ;IPROF = 1
 ABSIZ =0.5 ;ORSIZ =0.4
 NTZPRF = 1 ;NCOLPF = 50
 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
   No PATCHes yet used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP