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



  up z
  g ou ce X 1
  g ou ce x 1 y 1 5 z 9 12
  SET VEC REF 0.01;ve X 1 sh
  msg                    Velocity distribution:
  msg Press  to continue
  pause;ve off;red
  msg                   Temperature distribution :
  con temp x 1 fil;0.001
  g ou ce x 1 y 1 5 z 9 12 col 4
  msg Press e to END
  ENDUSE
  DISPLAY
   This case simulates the natural convection phenomena which occur
   when passive solids are placed in stably  stratified  fluids  of
   different  thermal conductivities.  Because of highly non-linear
   problem nature   the  number  of  sweeps  for  under-relaxations
   employed may not be enough for  completely  converged  solution.
   More details  of physical nature and mathematical statement plus
   the experimental and numeric data for comparisons can  be  found
   in  G.L.Quarini and K.H.  Winters ( Proc.  of 1st U.K.  National
   Conference on Heat Transfer, 1984, v.2, pp. 857-866)
  ENDDIS
  PLANTBEGIN
    VISL=VISC
    LAMPR(TEMP)=DIFS
PATCH(TMP,INIVAL,1,1,1,NY,1,NZ,1,1)
   VAL=ZGNZ/0.3
INIT(TMP,TEMP,0.0,GRND)
     Buoyancy forces
PATCH(BODY,VOLUME,1,1,1,NY,1,NZ,1,1)
    VAL=100.*5.4*TEMP
COVAL(BODY,W1,FIXFLU,GRND)
  PLANTEND
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT( NATURAL STIRRING IN A STABLY FLUID:124 )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 14
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = F
 NX = 1
 XULAST =1.
 XFRAC(1)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 10
 YVLAST =0.0875
 YFRAC(1)=0.121911 ;YFRAC(2)=0.229389
 YFRAC(3)=0.32015 ;YFRAC(4)=0.390239
 YFRAC(5)=0.428571 ;YFRAC(6)=0.479682
 YFRAC(7)=0.573133 ;YFRAC(8)=0.694147
 YFRAC(9)=0.837452 ;YFRAC(10)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 20
 ZWLAST =0.3
 ZFRAC(1)=0.087891 ;ZFRAC(2)=0.164063
 ZFRAC(3)=0.228516 ;ZFRAC(4)=0.28125
 ZFRAC(5)=0.322266 ;ZFRAC(6)=0.351563
 ZFRAC(7)=0.369141 ;ZFRAC(8)=0.375
 ZFRAC(9)=0.40625 ;ZFRAC(10)=0.5
 ZFRAC(11)=0.59375 ;ZFRAC(12)=0.625
 ZFRAC(13)=0.630859 ;ZFRAC(14)=0.648438
 ZFRAC(15)=0.677734 ;ZFRAC(16)=0.71875
 ZFRAC(17)=0.771484 ;ZFRAC(18)=0.835938
 ZFRAC(19)=0.912109 ;ZFRAC(20)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(1)=P1 ;NAME(5)=V1
 NAME(7)=W1 ;NAME(148)=DIFS
 NAME(149)=VISC ;NAME(150)=TEMP
    * 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(V1,Y,Y,N,N,N,Y)
 SOLUTN(W1,Y,Y,N,N,N,Y)
 SOLUTN(DIFS,Y,N,N,N,N,Y)
 SOLUTN(VISC,Y,N,N,N,N,Y)
 SOLUTN(TEMP,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(V1,Y,Y,Y,Y,Y,Y)
 TERMS(W1,Y,Y,Y,Y,Y,Y)
 TERMS(TEMP,N,Y,Y,Y,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 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 =1. ;TMP1 =0. ;EL1 =0.
 TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
 DVO1DT =0. ;DRH1DP =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 ENUL = GRND ;ENUT =0.
 ENULA =0. ;ENULB =0. ;ENULC =0. ;ENULD =0.
 ENULE =0. ;ENULF =0. ;ENULG =0.
 IENULA = 0 ;IENULB = 0
 PRNDTL(V1)=1. ;PRNDTL(W1)=1.
 PRNDTL(TEMP)=-GRND
 PRT(V1)=1. ;PRT(W1)=1.
 PRT(TEMP)=1.
 CP1 =1. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1.0E-10 ;FIINIT(V1)=1.0E-10
 FIINIT(W1)=1.0E-02 ;FIINIT(DIFS)=1.0E-10
 FIINIT(VISC)=1.0E-10 ;FIINIT(TEMP)=1.0E-10
 
 PATCH(VISCF ,INIVAL, 1, 1, 1, 10, 1, 20, 1, 1)
 INIT(VISCF ,VISC,0. ,5.4 )
 
 PATCH(VISCS ,INIVAL, 1, 1, 1, 5, 9, 12, 1, 1)
 INIT(VISCS ,VISC,0. ,1.0E+10 )
 
 PATCH(KOVRF ,INIVAL, 1, 1, 1, 10, 1, 20, 1, 1)
 INIT(KOVRF ,DIFS,0. ,1. )
 
 PATCH(KOVRM ,INIVAL, 1, 1, 1, 5, 9, 12, 1, 1)
 INIT(KOVRM ,DIFS,0. ,600. )
 
 PATCH(TMP ,INIVAL, 1, 1, 1, 10, 1, 20, 1, 1)
 INIT(TMP ,TEMP,0. , GRND )
 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(REP ,VOLUME, 1, 1, 1, 1, 20, 20, 1, 1)
 COVAL(REP ,P1 ,1. ,0. )
 COVAL(REP ,V1 ,0. ,0. )
 COVAL(REP ,W1 ,0. ,0. )
 COVAL(REP ,TEMP,0. , SAME )
 
 PATCH(COLD ,LWALL , 1, 1, 1, 10, 1, 1, 1, 1)
 COVAL(COLD ,V1 ,1. ,0. )
 COVAL(COLD ,TEMP,0.185185 ,0. )
 
 PATCH(HOT ,HWALL , 1, 1, 1, 10, 20, 20, 1, 1)
 COVAL(HOT ,V1 ,1. ,0. )
 COVAL(HOT ,TEMP,0.185185 ,1. )
 
 PATCH(RIGHT ,NWALL , 1, 1, 10, 10, 1, 20, 1, 1)
 COVAL(RIGHT ,W1 ,1. ,0. )
 
 PATCH(LEFT ,SWALL , 1, 1, 1, 1, 1, 20, 1, 1)
 COVAL(LEFT ,W1 ,1. ,0. )
 
 PATCH(ZEROW ,CELL , 1, 1, 1, 5, 8, 12, 1, 1)
 COVAL(ZEROW ,W1 , FIXVAL ,0. )
 
 PATCH(ZEROV ,CELL , 1, 1, 1, 5, 9, 12, 1, 1)
 COVAL(ZEROV ,V1 , FIXVAL ,0. )
 
 PATCH(BODY ,VOLUME, 1, 1, 1, 10, 1, 20, 1, 1)
 COVAL(BODY ,W1 , FIXFLU , GRND )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 220 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-05
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(V1)=10
 LITER(W1)=10 ;LITER(TEMP)=20
 ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03
 ENDIT(W1)=1.0E-03 ;ENDIT(TEMP)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,0.5)
 RELAX(V1,FALSDT,2.0E-04)
 RELAX(W1,FALSDT,1.0E-04)
 RELAX(DIFS,LINRLX,1.)
 RELAX(VISC,LINRLX,1.)
 RELAX(TEMP,FALSDT,0.1)
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
 VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
 VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
 VARMAX(DIFS)=1.0E+10 ;VARMIN(DIFS)=-1.0E+10
 VARMAX(VISC)=1.0E+10 ;VARMIN(VISC)=-1.0E+10
 VARMAX(TEMP)=1.0E+10 ;VARMIN(TEMP)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 NAMSAT =MOSG
 PARSOL = F
 ISG62 = 1
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 ************************************************************
  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,N,Y,Y,Y,Y)
 OUTPUT(V1,Y,N,Y,Y,Y,Y)
 OUTPUT(W1,Y,N,Y,Y,Y,Y)
 OUTPUT(DIFS,Y,N,Y,N,N,N)
 OUTPUT(VISC,Y,N,Y,N,N,N)
 OUTPUT(TEMP,Y,N,Y,Y,Y,Y)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 1 ;IYMON = 6 ;IZMON = 8
 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 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
   No PATCHes yet used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP