TALK=T;RUN(1,1)
  DISPLAY
  The case considered is 2d transient turbulent fluid flow and
  heat transfer of liquid steel (melt) in a ladle. This case
  has been studied numerically by Chakraborty and Sahai (Metall.
  Trans. B, Vol.23B, p135, 1992).

  The cylindrical ladle is 3.35m in diameter, and the melt depth
  is 4m, and the holding period before release of the melt into
  the mold is 20 minutes.

  The melt is taken to be stagnant and isothermal at the start
  of the holding period. The free surface is presumed flat, and
  heat loss rates are prescribed at both walls and the free
  surface. Convection currents are set up in the melt due to
  buoyancy forces resulting from the cooling of the melt with
  time. After about 5 or 10 minutes, it is expected that a
  clockwise overall recirculation pattern will be established
  within the melt.
  ENDDIS
  PHOTON USE
    P



    UP Z
    VEC X 1 SH
    PAU;CL
    CON TMP1 X 1 FI;.01
  ENDUSE
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(2DTrans Buoyant Flow Of Melt In A Ladle )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 14
 ************************************************************
  Group 2. Time dependence
 STEADY = F
    * Set overall time and no. of steps
 TFIRST =0. ;TLAST =1200.
 FSTEP = 1 ;LSTEP = 300
 TFRAC(1)=3.333333E-03 ;TFRAC(76)=0.253333
 TFRAC(151)=0.503333 ;TFRAC(226)=0.753333
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = F
 NX = 1
 XULAST =0.1
 XFRAC(1)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 12
 YVLAST =1.675
 YFRAC(1)=0.061703 ;YFRAC(2)=0.151932
 YFRAC(3)=0.257375 ;YFRAC(4)=0.374099
 YFRAC(5)=0.5 ;YFRAC(6)=0.590797
 YFRAC(7)=0.677148 ;YFRAC(8)=0.758443
 YFRAC(9)=0.833812 ;YFRAC(10)=0.901897
 YFRAC(11)=0.960158 ;YFRAC(12)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 16
 ZWLAST =4.
 ZFRAC(1)=0.033493 ;ZFRAC(2)=0.082469
 ZFRAC(3)=0.139705 ;ZFRAC(4)=0.203063
 ZFRAC(5)=0.271403 ;ZFRAC(6)=0.343993
 ZFRAC(7)=0.42032 ;ZFRAC(8)=0.5
 ZFRAC(9)=0.57968 ;ZFRAC(10)=0.656007
 ZFRAC(11)=0.728597 ;ZFRAC(12)=0.796937
 ZFRAC(13)=0.860295 ;ZFRAC(14)=0.917531
 ZFRAC(15)=0.966507 ;ZFRAC(16)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(1)=P1 ;NAME(5)=V1
 NAME(7)=W1 ;NAME(12)=KE
 NAME(13)=EP ;NAME(14)=H1
 NAME(148)=ENUT ;NAME(149)=TMP1
 NAME(150)=RHO1
    * 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(V1,Y,Y,Y,N,N,N)
 SOLUTN(W1,Y,Y,Y,N,N,N)
 SOLUTN(KE,Y,Y,N,N,N,N)
 SOLUTN(EP,Y,Y,N,N,N,N)
 SOLUTN(H1,Y,Y,Y,N,N,N)
 SOLUTN(ENUT,Y,N,N,N,N,Y)
 SOLUTN(TMP1,Y,N,N,N,N,Y)
 SOLUTN(RHO1,Y,N,N,N,N,Y)
 DEN1 = 150
 VIST = 148
 TEMP1 = 149
 ************************************************************
  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(KE,N,Y,Y,Y,Y,N)
 TERMS(EP,N,Y,Y,Y,Y,N)
 TERMS(H1,Y,Y,Y,Y,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWRH1 = T
 NEWENT = T
 HUNIT =1.
 ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
 ************************************************************
  Group 9. Properties used if PRPS is not
  stored, and where PRPS = -1.0 if it is!
 RHO1 = GRND1 ;TMP1 = GRND2 ;EL1 = GRND4
 TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
 DVO1DT =7000. ;DRH1DP =0.
 RHO1A =7000. ;RHO1B =-1.866667E-03 ;RHO1C =0.
 TMP1A =1580. ;TMP1B =1.333333E-03 ;TMP1C =0.
 TMP2A =0. ;TMP2B =0. ;TMP2C =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 EL1A =0. ;EL1B =0. ;EL1C =0.
 ENUL =7.142857E-07 ;ENUT = GRND3
 ENUTA =0. ;ENUTB =0. ;ENUTC =0.
 IENUTA = 0
 PRNDTL(V1)=1. ;PRNDTL(W1)=1.
 PRNDTL(KE)=1. ;PRNDTL(EP)=1.
 PRNDTL(H1)=0.091463
 PRT(V1)=1. ;PRT(W1)=1.
 PRT(KE)=1. ;PRT(EP)=1.314
 PRT(H1)=0.9
 CP1 =750. ;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-10 ;FIINIT(KE)=1.0E-10
 FIINIT(EP)=1.0E-10 ;FIINIT(H1)=0.
 FIINIT(ENUT)=1.0E-10 ;FIINIT(TMP1)=1580.
 FIINIT(RHO1)=7000.
   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(KESOURCE,PHASEM, 0, 0, 0, 0, 0, 0, 1, 300)
 COVAL(KESOURCE,KE , GRND4 , GRND4 )
 COVAL(KESOURCE,EP , GRND4 , GRND4 )
 
 PATCH(QTOP ,HIGH , 1, 1, 1, 12, 16, 16, 1, 300)
 COVAL(QTOP ,H1 , FIXFLU ,-1.0E+05 )
 
 PATCH(BOTTOM ,LWALL , 1, 1, 1, 12, 1, 1, 1, 300)
 COVAL(BOTTOM ,V1 , GRND2 ,0. )
 COVAL(BOTTOM ,KE , GRND2 , GRND2 )
 COVAL(BOTTOM ,EP , GRND2 , GRND2 )
 
 PATCH(QBOT ,LOW , 1, 1, 1, 12, 1, 1, 1, 300)
 COVAL(QBOT ,H1 , FIXFLU ,-1.25E+04 )
 
 PATCH(SIDEW ,NWALL , 1, 1, 12, 12, 1, 16, 1, 300)
 COVAL(SIDEW ,W1 , GRND2 ,0. )
 COVAL(SIDEW ,KE , GRND2 , GRND2 )
 COVAL(SIDEW ,EP , GRND2 , GRND2 )
 
 PATCH(QSIDE ,NORTH , 1, 1, 12, 12, 1, 16, 1, 300)
 COVAL(QSIDE ,H1 , FIXFLU ,-1.25E+04 )
 
 PATCH(BUOYANCY,PHASEM, 1, 1, 1, 12, 1, 16, 1, 300)
 COVAL(BUOYANCY,W1 , FIXFLU , GRND2 )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 BUOYA =0. ; BUOYB =0. ; BUOYC =-9.81
 BUOYD =7000.
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 30 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = F
 RESREF(P1)=1.0E-12 ;RESREF(V1)=1.0E-08
 RESREF(W1)=1.0E-08 ;RESREF(KE)=1.0E-08
 RESREF(EP)=1.0E-08 ;RESREF(H1)=1.0E-08
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=100 ;LITER(V1)=10
 LITER(W1)=10 ;LITER(KE)=20
 LITER(EP)=20 ;LITER(H1)=100
 ENDIT(P1)=GRND1 ;ENDIT(V1)=1.0E-03
 ENDIT(W1)=1.0E-03 ;ENDIT(KE)=1.0E-03
 ENDIT(EP)=1.0E-03 ;ENDIT(H1)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(V1,FALSDT,1.)
 RELAX(W1,FALSDT,1.)
 RELAX(KE,LINRLX,0.5)
 RELAX(EP,LINRLX,0.5)
 RELAX(H1,FALSDT,1.0E+09)
 RELAX(ENUT,LINRLX,1.)
 RELAX(TMP1,LINRLX,1.)
 RELAX(RHO1,LINRLX,1.)
 KELIN = 0
 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(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
 VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
 VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10
 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
 VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=-1.0E+10
 VARMAX(RHO1)=1.0E+10 ;VARMIN(RHO1)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 GENK = T
 PARSOL = F
 ISG62 = 1
 SPEDAT(SET,GXMONI,TRANSIENT,L,F)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=313.5 ;EX(V1)=0.
 EX(W1)=0.0219 ;EX(KE)=0.
 EX(EP)=0. ;EX(H1)=7389.
 EX(ENUT)=0. ;EX(TMP1)=1570.
 EX(RHO1)=7014.
 ************************************************************
  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(KE,Y,N,Y,Y,Y,Y)
 OUTPUT(EP,Y,N,Y,Y,Y,Y)
 OUTPUT(H1,Y,N,Y,Y,Y,Y)
 OUTPUT(ENUT,Y,N,Y,N,N,N)
 OUTPUT(TMP1,Y,Y,Y,Y,Y,Y)
 OUTPUT(RHO1,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 1 ;IYMON = 11 ;IZMON = 15
 NPRMON = 10000 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NTPRIN = 75 ;ISTPRF = 1 ;ISTPRL = 100000
 NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
 NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
 XZPR = F ;YZPR = F
 IPLTF = 1 ;IPLTL = -1 ;NPLT = 5
 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
 IDISPA = 75 ;IDISPB = 0 ;IDISPC = 0
 CSG1    ='G'
STOP