TALK=T;RUN(1,1)
  DISPLAY
   Turbulent buoyancy-driven air flow in a tall rectangular
   cavity of 5:1 aspect ratio with a Rayleigh number of 4.E10 is
   calculated using a buoyancy-extended low-Reynolds-number k-e
   and k-w models. The vertical walls are isothermal with a
   temperature difference of 45.8 degC; the horizontal walls are
   adiabatic.

   This problem has been studied experimentally by Cheesewright
   et al (1986) and numerically by Davidson (1990) and Peng &
   Davidson (1999). The primary objective is to compute the
   Nusselt number distribution along the hot and cold walls, and
   the overall Nusselt number. The calculation may be performed
   with or without the Boussinesq approximation, and the dynamic
   laminar viscosity is computed from Sutherland's formula. For
   the Boussinesq approximation, a uniform density and thermal
   conductivity are employed in the simulations.

   The main result is the average Nusselt number, and the table
   below compares the predicted values with the value deduced
   from the measurements:
                Data      LB     KW-R    KW-SST
       Nu,av     160      154      160     153
   All of the models produce good agreement with the measured
   value.

   ENDDIS
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(2D Turb. Buoyant Cavity Flow K-E:T219   )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 14
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 56
 XULAST =1.
 XFRAC(1)=1.294613E-04 ;XFRAC(12)=8.553162E-03
 XFRAC(23)=0.111097 ;XFRAC(34)=0.409077
 XFRAC(45)=0.493412 ;XFRAC(56)=0.5
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 56
 YVLAST =1.
 YFRAC(1)=6.473064E-04 ;YFRAC(12)=0.042766
 YFRAC(23)=0.555485 ;YFRAC(34)=2.045386
 YFRAC(45)=2.467059 ;YFRAC(56)=2.5
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 1
 ZWLAST =1.
 ZFRAC(1)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(1)=P1 ;NAME(3)=U1
 NAME(5)=V1 ;NAME(12)=KE
 NAME(13)=EP ;NAME(127)=C3EB
 NAME(128)=GENB ;NAME(129)=V1N
 NAME(130)=U1N ;NAME(131)=TNO
 NAME(132)=NUC ;NAME(133)=NUH
 NAME(134)=HTC ;NAME(135)=HTH
 NAME(136)=QC ;NAME(137)=QH
 NAME(138)=DWAL ;NAME(139)=LEN1
 NAME(140)=YPLS ;NAME(141)=ENUL
 NAME(142)=ENUT ;NAME(143)=KOND
 NAME(144)=DEN1 ;NAME(145)=TREY
 NAME(146)=REYN ;NAME(147)=LTLS
 NAME(148)=WDIS ;NAME(149)=EPKE
 NAME(150)=TEM1
    * 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(KE,Y,Y,N,N,N,N)
 SOLUTN(EP,Y,Y,N,N,N,N)
 SOLUTN(C3EB,Y,N,N,N,N,Y)
 SOLUTN(GENB,Y,N,N,N,N,Y)
 SOLUTN(V1N,Y,N,N,N,N,Y)
 SOLUTN(U1N,Y,N,N,N,N,Y)
 SOLUTN(TNO,Y,N,N,N,N,Y)
 SOLUTN(NUC,Y,N,N,N,N,Y)
 SOLUTN(NUH,Y,N,N,N,N,Y)
 SOLUTN(HTC,Y,N,N,N,N,Y)
 SOLUTN(HTH,Y,N,N,N,N,Y)
 SOLUTN(QC,Y,N,N,N,N,Y)
 SOLUTN(QH,Y,N,N,N,N,Y)
 SOLUTN(DWAL,Y,N,N,N,N,Y)
 SOLUTN(LEN1,Y,N,N,N,N,Y)
 SOLUTN(YPLS,Y,N,N,N,N,Y)
 SOLUTN(ENUL,Y,N,N,N,N,Y)
 SOLUTN(ENUT,Y,N,N,N,N,Y)
 SOLUTN(KOND,Y,N,N,N,N,Y)
 SOLUTN(DEN1,Y,N,N,N,N,Y)
 SOLUTN(TREY,Y,N,N,N,N,Y)
 SOLUTN(REYN,Y,N,N,N,N,Y)
 SOLUTN(LTLS,Y,Y,Y,N,N,Y)
 SOLUTN(WDIS,Y,N,N,N,N,N)
 SOLUTN(EPKE,Y,N,N,N,N,N)
 SOLUTN(TEM1,Y,Y,N,N,N,Y)
 DEN1 = 144
 VISL = 141
 VIST = 142
 LEN1 = 139
 ************************************************************
  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(KE,N,Y,Y,Y,Y,N)
 TERMS(EP,N,Y,Y,Y,Y,N)
 TERMS(LTLS,N,N,Y,N,Y,Y)
 TERMS(TEM1,N,Y,Y,Y,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWRH1 = T
 NEWENL = T
 NEWENT = T
 ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
 DENPCO = T
 ************************************************************
  Group 9. Properties used if PRPS is not
  stored, and where PRPS = -1.0 if it is!
 RHO1 = GRND10 ;TMP1 =0. ;EL1 = GRND4
 TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
 DVO1DT =1.070903 ;DRH1DP =0.
 RHO1A =0. ;RHO1B =2.828814E-03 ;RHO1C =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 EL1A =0. ;EL1B =0. ;EL1C =0.
 ENUL = GRND6 ;ENUT = GRND3
 ENULA =1.458E-06 ;ENULB =110.400002 ;ENULC =0. ;ENULD =0.
 ENULE =0. ;ENULF =0. ;ENULG =0.
 IENULA = 0 ;IENULB = 0
 ENUTA =0. ;ENUTB =0. ;ENUTC =0.
 IENUTA = 3
 PRNDTL(U1)=1. ;PRNDTL(V1)=1.
 PRNDTL(KE)=1. ;PRNDTL(EP)=1.
 PRNDTL(LTLS)=1. ;PRNDTL(TEM1)=-GRND4
 PRLH1A =3.5E-03 ;PRLH1B =7.5E-05 ;PRLH1C =0.
 PRT(U1)=1. ;PRT(V1)=1.
 PRT(KE)=1. ;PRT(EP)=1.314
 PRT(LTLS)=1. ;PRT(TEM1)=0.9
 CP1 =1008. ;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(V1)=1.0E-10 ;FIINIT(KE)=3.402741E-03
 FIINIT(EP)=1.368744E-03 ;FIINIT(C3EB)=1.0E-10
 FIINIT(GENB)=1.0E-10 ;FIINIT(V1N)=1.0E-10
 FIINIT(U1N)=1.0E-10 ;FIINIT(TNO)=1.0E-10
 FIINIT(NUC)=1.0E-10 ;FIINIT(NUH)=1.0E-10
 FIINIT(HTC)=1.0E-10 ;FIINIT(HTH)=1.0E-10
 FIINIT(QC)=1.0E-10 ;FIINIT(QH)=1.0E-10
 FIINIT(DWAL)=1.0E-10 ;FIINIT(LEN1)=1.0E-10
 FIINIT(YPLS)=1.0E-10 ;FIINIT(ENUL)=1.0E-10
 FIINIT(ENUT)=7.61339E-04 ;FIINIT(KOND)=0.028522
 FIINIT(DEN1)=1.070903 ;FIINIT(TREY)=1.0E-10
 FIINIT(REYN)=1.0E-10 ;FIINIT(LTLS)=1.0E-10
 FIINIT(WDIS)=0.1 ;FIINIT(EPKE)=1.0E-10
 FIINIT(TEM1)=330.100006
   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, 1)
 COVAL(KESOURCE,KE , GRND4 , GRND4 )
 COVAL(KESOURCE,EP , GRND4 , GRND4 )
 
 PATCH(HOT ,WWALL , 1, 1, 1, 56, 1, 1, 1, 1)
 COVAL(HOT ,V1 , GRND2 ,0. )
 COVAL(HOT ,KE ,1. ,0. )
 COVAL(HOT ,TEM1, GRND2 ,353. )
 
 PATCH(COLD ,EWALL , 56, 56, 1, 56, 1, 1, 1, 1)
 COVAL(COLD ,V1 , GRND2 ,0. )
 COVAL(COLD ,KE ,1. ,0. )
 COVAL(COLD ,TEM1, GRND2 ,307.200012 )
 
 PATCH(TOPWAL ,NWALL , 1, 56, 56, 56, 1, 1, 1, 1)
 COVAL(TOPWAL ,U1 , GRND2 ,0. )
 COVAL(TOPWAL ,KE ,1. ,0. )
 
 PATCH(BOTWAL ,SWALL , 1, 56, 1, 1, 1, 1, 1, 1)
 COVAL(BOTWAL ,U1 , GRND2 ,0. )
 COVAL(BOTWAL ,KE ,1. ,0. )
 
 PATCH(BUOY ,PHASEM, 1, 56, 1, 56, 1, 1, 1, 1)
 COVAL(BUOY ,V1 , FIXFLU , GRND2 )
 
 PATCH(REFP ,CELL , 28, 28, 28, 28, 1, 1, 1, 1)
 COVAL(REFP ,P1 ,1. ,0. )
 COVAL(REFP ,V1 ,0. ,0. )
 
 PATCH(KEBUOY ,PHASEM, 1, 56, 1, 56, 1, 1, 1, 1)
 COVAL(KEBUOY ,KE , GRND4 , GRND4 )
 COVAL(KEBUOY ,EP , GRND4 , GRND4 )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 BUOYA =0. ; BUOYB =-9.81 ; BUOYC =0.
 BUOYD =1.070903
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 8000 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-05
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=40 ;LITER(U1)=2
 LITER(V1)=2 ;LITER(KE)=2
 LITER(EP)=2 ;LITER(LTLS)=20
 LITER(TEM1)=50
 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
 ENDIT(V1)=1.0E-03 ;ENDIT(KE)=1.0E-03
 ENDIT(EP)=1.0E-03 ;ENDIT(LTLS)=1.0E-03
 ENDIT(TEM1)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(U1,FALSDT,0.024201)
 RELAX(V1,FALSDT,0.024201)
 RELAX(KE,FALSDT,0.024201)
 RELAX(EP,FALSDT,0.024201)
 RELAX(C3EB,LINRLX,1.)
 RELAX(GENB,LINRLX,1.)
 RELAX(V1N,LINRLX,1.)
 RELAX(U1N,LINRLX,1.)
 RELAX(TNO,LINRLX,1.)
 RELAX(NUC,LINRLX,1.)
 RELAX(NUH,LINRLX,1.)
 RELAX(HTC,LINRLX,1.)
 RELAX(HTH,LINRLX,1.)
 RELAX(QC,LINRLX,1.)
 RELAX(QH,LINRLX,1.)
 RELAX(DWAL,LINRLX,1.)
 RELAX(LEN1,LINRLX,1.)
 RELAX(YPLS,LINRLX,1.)
 RELAX(ENUL,LINRLX,1.)
 RELAX(ENUT,LINRLX,1.)
 RELAX(KOND,LINRLX,1.)
 RELAX(DEN1,LINRLX,1.)
 RELAX(TREY,LINRLX,1.)
 RELAX(REYN,LINRLX,1.)
 RELAX(LTLS,FALSDT,1.0E+09)
 RELAX(WDIS,LINRLX,1.)
 RELAX(EPKE,LINRLX,1.)
 RELAX(TEM1,FALSDT,1.210062)
 KELIN = 3
 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(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
 VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
 VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
 VARMAX(C3EB)=1.0E+10 ;VARMIN(C3EB)=-1.0E+10
 VARMAX(GENB)=1.0E+10 ;VARMIN(GENB)=-1.0E+10
 VARMAX(V1N)=1.0E+10 ;VARMIN(V1N)=-1.0E+10
 VARMAX(U1N)=1.0E+10 ;VARMIN(U1N)=-1.0E+10
 VARMAX(TNO)=1.0E+10 ;VARMIN(TNO)=-1.0E+10
 VARMAX(NUC)=1.0E+10 ;VARMIN(NUC)=-1.0E+10
 VARMAX(NUH)=1.0E+10 ;VARMIN(NUH)=-1.0E+10
 VARMAX(HTC)=1.0E+10 ;VARMIN(HTC)=-1.0E+10
 VARMAX(HTH)=1.0E+10 ;VARMIN(HTH)=-1.0E+10
 VARMAX(QC)=1.0E+10 ;VARMIN(QC)=-1.0E+10
 VARMAX(QH)=1.0E+10 ;VARMIN(QH)=-1.0E+10
 VARMAX(DWAL)=1.0E+10 ;VARMIN(DWAL)=-1.0E+10
 VARMAX(LEN1)=1.0E+10 ;VARMIN(LEN1)=-1.0E+10
 VARMAX(YPLS)=1.0E+10 ;VARMIN(YPLS)=-1.0E+10
 VARMAX(ENUL)=1.0E+10 ;VARMIN(ENUL)=-1.0E+10
 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
 VARMAX(KOND)=1.0E+10 ;VARMIN(KOND)=-1.0E+10
 VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=-1.0E+10
 VARMAX(TREY)=1.0E+10 ;VARMIN(TREY)=-1.0E+10
 VARMAX(REYN)=1.0E+10 ;VARMIN(REYN)=-1.0E+10
 VARMAX(LTLS)=1.0E+10 ;VARMIN(LTLS)=-1.0E+10
 VARMAX(WDIS)=1.0E+10 ;VARMIN(WDIS)=-1.0E+10
 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
 VARMAX(TEM1)=353. ;VARMIN(TEM1)=307.200012
 ************************************************************
  Group 19. Data transmitted to GROUND
 GENK = T
 PARSOL = F
 IENUTA = 3
 ISG62 = 1
 SPEDAT(SET,STORED,DWAL,C,=0.5*DXU!ZSLFIN)
 SPEDAT(SET,STORED,QH!HOT,C,=KOND*(353.-TEM1)/DWAL!ZSLFIN)
 SPEDAT(SET,STORED,QC!COLD,C,=KOND*(307.200012-TEM1)/DWAL!ZSLFIN)
 SPEDAT(SET,STORED,HTH!HOT,C,=QH/45.799999!ZSLFIN)
 SPEDAT(SET,STORED,HTC!COLD,C,=QC/45.799999!ZSLFIN)
 SPEDAT(SET,STORED,NUH!HOT,C,=((353.-TEM1)/DWAL)*2.5/45.799999!ZSL$)
 SPEDAT(SET,STORED,NUH!HOT,C,FIN)
 SPEDAT(SET,STORED,NUC!COLD,C,=ABS((307.200012-TEM1)/DWAL)*2.5/45.$)
 SPEDAT(SET,STORED,NUC!COLD,C,799999!ZSLFIN)
 SPEDAT(SET,MAKE,QHOT,C,=0.0)
 SPEDAT(SET,STORED,QHOT!HOT,C,=SUM(QH*AEAST)!ZSLFIN)
 SPEDAT(SET,PRINT,Q_HOT,C,=QHOT)
 SPEDAT(SET,MAKE,QCOLD,C,=0.0)
 SPEDAT(SET,STORED,QCOLD!COLD,C,=SUM(QC*AEAST)!ZSLFIN)
 SPEDAT(SET,PRINT,Q_COLD,C,=QCOLD)
 SPEDAT(SET,MAKE,NUAHOT,C,=0.0)
 SPEDAT(SET,STORED,NUAHOT!HOT,C,=QHOT/(45.799999*0.028522)!ZSLFIN)
 SPEDAT(SET,PRINT,NU_AV_HOT,C,=NUAHOT)
 SPEDAT(SET,MAKE,NUACOLD,C,=0.0)
 SPEDAT(SET,STORED,NUACOLD!COLD,C,=ABS(QCOLD)/(45.799999*0.028522)$)
 SPEDAT(SET,STORED,NUACOLD!COLD,C,!ZSLFIN)
 SPEDAT(SET,PRINT,NU_AV_COLD,C,=NUACOLD)
 SPEDAT(SET,STORED,U1N,C,=U1/1.844652!ZSLFIN)
 SPEDAT(SET,STORED,V1N,C,=V1/1.844652!ZSLFIN)
 SPEDAT(SET,STORED,TNO,C,=(TEM1-307.200012)/45.799999!ZSLFIN)
 SPEDAT(SET,KEBUOY,GCEB,R,-1.)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=0.2086 ;EX(U1)=8.015E-03
 EX(V1)=0.05177 ;EX(KE)=6.196E-04
 EX(EP)=8.658E-04 ;EX(C3EB)=0.606
 EX(GENB)=2.424E-05 ;EX(V1N)=0.02757
 EX(U1N)=4.267E-03 ;EX(TNO)=0.5034
 EX(NUC)=4.328 ;EX(NUH)=3.898
 EX(HTC)=0.04602 ;EX(HTH)=0.04668
 EX(QC)=2.108 ;EX(QH)=2.138
 EX(DWAL)=4.464E-03 ;EX(LEN1)=4.112E-03
 EX(YPLS)=0.01008 ;EX(ENUL)=1.858E-05
 EX(ENUT)=5.987E-05 ;EX(KOND)=0.02827
 EX(DEN1)=1.073 ;EX(TREY)=3.235
 EX(REYN)=27.629999 ;EX(LTLS)=4.122E-03
 EX(WDIS)=0.02636 ;EX(EPKE)=1.422E+06
 EX(TEM1)=330.299988
 ************************************************************
  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(V1,Y,N,Y,Y,Y,Y)
 OUTPUT(KE,Y,N,Y,Y,Y,Y)
 OUTPUT(EP,Y,N,Y,Y,Y,Y)
 OUTPUT(C3EB,Y,N,Y,N,N,N)
 OUTPUT(GENB,Y,N,Y,N,N,N)
 OUTPUT(V1N,Y,N,Y,N,N,N)
 OUTPUT(U1N,Y,N,Y,N,N,N)
 OUTPUT(TNO,Y,N,Y,N,N,N)
 OUTPUT(NUC,Y,N,Y,N,N,N)
 OUTPUT(NUH,Y,N,Y,N,N,N)
 OUTPUT(HTC,Y,N,Y,N,N,N)
 OUTPUT(HTH,Y,N,Y,N,N,N)
 OUTPUT(QC,Y,N,Y,N,N,N)
 OUTPUT(QH,Y,N,Y,N,N,N)
 OUTPUT(DWAL,Y,N,Y,N,N,N)
 OUTPUT(LEN1,Y,N,Y,N,N,N)
 OUTPUT(YPLS,Y,N,Y,N,N,N)
 OUTPUT(ENUL,Y,N,Y,N,N,N)
 OUTPUT(ENUT,Y,N,Y,N,N,N)
 OUTPUT(KOND,Y,N,Y,N,N,N)
 OUTPUT(DEN1,Y,N,Y,N,N,N)
 OUTPUT(TREY,Y,N,Y,N,N,N)
 OUTPUT(REYN,Y,N,Y,N,N,N)
 OUTPUT(LTLS,Y,N,Y,Y,Y,Y)
 OUTPUT(WDIS,Y,N,Y,N,N,N)
 OUTPUT(EPKE,Y,N,Y,N,N,N)
 OUTPUT(TEM1,Y,N,Y,Y,Y,Y)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 5 ;IYMON = 28 ;IZMON = 1
 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NXPRIN = 11 ;IXPRF = 1 ;IXPRL = 10000
 NYPRIN = 11 ;IYPRF = 1 ;IYPRL = 10000
 IPLTF = 1 ;IPLTL = -1 ;NPLT = 400
 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(NUHOT ,PROFIL, 1, 1, 1, 56, 1, 1, 1, 1)
 PLOT(NUHOT ,NUH ,0. ,0. )
 
 PATCH(NUCOLD ,PROFIL, 56, 56, 1, 56, 1, 1, 1, 1)
 PLOT(NUCOLD ,NUC ,0. ,0. )
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP