TALK=T;RUN(1,1)
  DISPLAY
  ----------------------------------------------------------
  This case concerns simulation of a stratified flow through
  Oresund Channel. Originally, problem had been set by Urban
  Svensson (CFE). The library case is its simplified version.

  The  main  objective  of  this  case is to demonstrate new
  features of CCM specially  developed for simulation of  so
  called shallow water flows, i.e flows with character  size
  in  one  direction  much  less  than  these in two others.
  'Shallow' direction is Z axis. (See MBFGE entry in PHENC).

  There are two options:
   LSG10 when TRUE activates simulation of shallow flows  by
         truncated 3D Navier-Stokes equations;
   LSG8  when TRUE  activates same as LSG10  treatment; plus
         pressure-correction equation is solved 2D, with WC1
         defined by integration of continuity equation.
  ----------------------------------------------------------
  ENDDIS
  PHOTON USE
   p ; ; ; ; ;

   set vec av off
   msg Velocity Vectors:
   vi z
   vec k 1 sh
   pause
   cl
   vec k m sh
  ENDUSE
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT( CCM: Bay flow (3D,K-E                  )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 0
 ************************************************************
  Group 2. Time dependence
 STEADY = F
    * Set overall time and no. of steps
 TFIRST =0. ;TLAST =6000.
 FSTEP = 1 ;LSTEP = 10
 TFRAC(1)=0.1 ;TFRAC(2)=0.2
 TFRAC(3)=0.3 ;TFRAC(4)=0.4
 TFRAC(5)=0.5 ;TFRAC(6)=0.6
 TFRAC(7)=0.7 ;TFRAC(8)=0.8
 TFRAC(9)=0.9 ;TFRAC(10)=1.
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 14
 XULAST =1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 27
 YVLAST =1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 8
 ZWLAST =1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 BFC = T ;NONORT = F
 NCRT = -1
 RSTGEO = F ;SAVGEO = F
 UUP = F ;VUP = F ;WUP = F
 NGEOM =CHAM
 NAMXYZ =CHAM
 ANGMIN =20.
 DOMAIN(   1,  15,   1,  28,   1,   9)
    * 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
 CCM = T
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(1)=P1 ;NAME(3)=U1
 NAME(5)=V1 ;NAME(7)=W1
 NAME(12)=KE ;NAME(13)=EP
 NAME(14)=SALT ;NAME(16)=UC1
 NAME(17)=VC1 ;NAME(18)=WC1
 NAME(132)=PRPS ;NAME(133)=EPKE
 NAME(134)=LEN1 ;NAME(135)=RHO1
 NAME(136)=GEN1 ;NAME(137)=ENUT
 NAME(138)=DWDZ ;NAME(139)=DWDY
 NAME(140)=DWDX ;NAME(141)=DVDZ
 NAME(142)=DVDY ;NAME(143)=DVDX
 NAME(144)=DUDZ ;NAME(145)=DUDY
 NAME(146)=DUDX ;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(V1,Y,Y,N,N,N,Y)
 SOLUTN(W1,Y,Y,N,N,N,Y)
 SOLUTN(KE,Y,Y,Y,N,N,N)
 SOLUTN(EP,Y,Y,Y,N,N,N)
 SOLUTN(SALT,Y,Y,Y,N,N,Y)
 SOLUTN(UC1,Y,Y,Y,N,N,Y)
 SOLUTN(VC1,Y,Y,Y,N,N,Y)
 SOLUTN(WC1,Y,Y,Y,N,N,Y)
 SOLUTN(PRPS,Y,N,N,N,N,N)
 SOLUTN(EPKE,Y,N,N,N,N,N)
 SOLUTN(LEN1,Y,N,N,N,N,Y)
 SOLUTN(RHO1,Y,N,N,N,N,Y)
 SOLUTN(GEN1,Y,N,N,N,N,Y)
 SOLUTN(ENUT,Y,N,N,N,N,Y)
 SOLUTN(DWDZ,Y,N,N,N,N,Y)
 SOLUTN(DWDY,Y,N,N,N,N,Y)
 SOLUTN(DWDX,Y,N,N,N,N,Y)
 SOLUTN(DVDZ,Y,N,N,N,N,Y)
 SOLUTN(DVDY,Y,N,N,N,N,Y)
 SOLUTN(DVDX,Y,N,N,N,N,Y)
 SOLUTN(DUDZ,Y,N,N,N,N,Y)
 SOLUTN(DUDY,Y,N,N,N,N,Y)
 SOLUTN(DUDX,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 = 135
 VIST = 137
 LEN1 = 134
 PRPS = 132
 ************************************************************
  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,N,N,N,N,N,N)
 TERMS(V1,N,N,N,N,N,N)
 TERMS(W1,N,N,N,N,N,N)
 TERMS(KE,N,Y,Y,Y,Y,N)
 TERMS(EP,N,Y,Y,Y,Y,N)
 TERMS(SALT,N,Y,Y,Y,Y,Y)
 TERMS(UC1,N,Y,Y,Y,Y,Y)
 TERMS(VC1,N,Y,Y,Y,N,Y)
 TERMS(WC1,N,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 =0. ;EL1 = GRND4
 TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
 DVO1DT =0. ;DRH1DP =0.
 RHO1A =1000. ;RHO1B =0.8 ;RHO1C =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 EL1A =0. ;EL1B =0. ;EL1C =0.
 ENUL =1.0E-06 ;ENUT = GRND3
 ENUTA =0. ;ENUTB =0. ;ENUTC =0.
 IENUTA = 0
 PRNDTL(KE)=1. ;PRNDTL(EP)=1.
 PRNDTL(SALT)=1000. ;PRNDTL(UC1)=1.
 PRNDTL(VC1)=1. ;PRNDTL(WC1)=1.
 PRT(KE)=1. ;PRT(EP)=1.314
 PRT(SALT)=2. ;PRT(UC1)=1.
 PRT(VC1)=1. ;PRT(WC1)=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(V1)=1.0E-10 ;FIINIT(W1)=1.0E-10
 FIINIT(KE)=1.0E-04 ;FIINIT(EP)=1.0E-06
 FIINIT(SALT)=10. ;FIINIT(UC1)=1.0E-06
 FIINIT(VC1)=1.0E-06 ;FIINIT(WC1)=0.
 FIINIT(PRPS)=-1. ;FIINIT(EPKE)=1.0E-10
 FIINIT(LEN1)=1.0E-10 ;FIINIT(RHO1)=1.0E-10
 FIINIT(GEN1)=1.0E-10 ;FIINIT(ENUT)=1.0E-03
 FIINIT(DWDZ)=1.0E-10 ;FIINIT(DWDY)=1.0E-10
 FIINIT(DWDX)=1.0E-10 ;FIINIT(DVDZ)=1.0E-10
 FIINIT(DVDY)=1.0E-10 ;FIINIT(DVDX)=1.0E-10
 FIINIT(DUDZ)=1.0E-10 ;FIINIT(DUDY)=1.0E-10
 FIINIT(DUDX)=1.0E-10 ;FIINIT(WCRT)=1.0E-10
 FIINIT(VCRT)=1.0E-10 ;FIINIT(UCRT)=1.0E-10
 
 PATCH(CMP1 ,INIVAL, 7, 7, 10, 14, 1, 8, 1, 1)
 INIT(CMP1 ,PRPS,0. ,199. )
 
 PATCH(CMP2 ,INIVAL, 1, 4, 26, 27, 1, 8, 1, 1)
 INIT(CMP2 ,PRPS,0. ,199. )
 
 PATCH(CMP3 ,INIVAL, 1, 3, 25, 25, 1, 8, 1, 1)
 INIT(CMP3 ,PRPS,0. ,199. )
 
 PATCH(CMP4 ,INIVAL, 1, 2, 24, 24, 1, 8, 1, 1)
 INIT(CMP4 ,PRPS,0. ,199. )
 
 PATCH(CMP5 ,INIVAL, 1, 1, 23, 23, 1, 8, 1, 1)
 INIT(CMP5 ,PRPS,0. ,199. )
 
 PATCH(CMP6 ,INIVAL, 8, 8, 21, 21, 1, 8, 1, 1)
 INIT(CMP6 ,PRPS,0. ,199. )
 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, 10)
 COVAL(KESOURCE,KE , GRND4 , GRND4 )
 COVAL(KESOURCE,EP , GRND4 , GRND4 )
 
 PATCH(NORR ,NORTH , 5, 14, 27, 27, 5, 8, 1, 10)
 COVAL(NORR ,P1 ,10. ,0. )
 COVAL(NORR ,KE ,0. ,1.0E-04 )
 COVAL(NORR ,EP ,0. ,1.0E-06 )
 COVAL(NORR ,SALT,0. , SAME )
 
 PATCH(NORRL ,NORTH , 5, 14, 27, 27, 1, 4, 1, 10)
 COVAL(NORRL ,P1 ,1.0E-20 ,3.0E+22 )
 COVAL(NORRL ,KE ,0. ,1.0E-04 )
 COVAL(NORRL ,EP ,0. ,1.0E-06 )
 COVAL(NORRL ,SALT,0. ,30. )
 
 PATCH(SYD ,SOUTH , 5, 14, 1, 1, 1, 8, 1, 10)
 COVAL(SYD ,P1 ,1000. ,1000. )
 COVAL(SYD ,KE ,0. ,1.0E-04 )
 COVAL(SYD ,EP ,0. ,1.0E-06 )
 COVAL(SYD ,SALT,0. ,10. )
 
 PATCH(BOTFRI ,LWALL , 1, 14, 1, 27, 1, 1, 1, 10)
 COVAL(BOTFRI ,KE , GRND2 , GRND2 )
 COVAL(BOTFRI ,EP , GRND2 , GRND2 )
 COVAL(BOTFRI ,UC1 , GRND2 ,0. )
 COVAL(BOTFRI ,VC1 , GRND2 ,0. )
 COVAL(BOTFRI ,WC1 , GRND2 ,0. )
 
 PATCH(KEBUOY ,PHASEM, 1, 14, 1, 27, 1, 8, 1, 10)
 COVAL(KEBUOY ,KE , GRND4 , GRND4 )
 COVAL(KEBUOY ,EP , GRND4 , GRND4 )
 XCYCLE = F
 CORIOL =1.0E-04
 EGWF = T
 WALLCO = GRND2
 BUOYA =0. ; BUOYB =0. ; BUOYC =-9.81
 BUOYE =1008.
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 50 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-06
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(U1)=10
 LITER(V1)=10 ;LITER(W1)=10
 LITER(KE)=10 ;LITER(EP)=10
 LITER(SALT)=20 ;LITER(UC1)=5
 LITER(VC1)=5 ;LITER(WC1)=5
 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
 ENDIT(V1)=1.0E-03 ;ENDIT(W1)=1.0E-03
 ENDIT(KE)=1.0E-03 ;ENDIT(EP)=1.0E-03
 ENDIT(SALT)=1.0E-03 ;ENDIT(UC1)=1.0E-03
 ENDIT(VC1)=1.0E-03 ;ENDIT(WC1)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,0.5)
 RELAX(U1,FALSDT,1.)
 RELAX(V1,FALSDT,1.)
 RELAX(W1,FALSDT,1.)
 RELAX(KE,LINRLX,0.25)
 RELAX(EP,LINRLX,0.25)
 RELAX(SALT,FALSDT,1000.)
 RELAX(UC1,FALSDT,1000.)
 RELAX(VC1,FALSDT,1000.)
 RELAX(WC1,FALSDT,1.)
 RELAX(PRPS,LINRLX,1.)
 RELAX(EPKE,LINRLX,1.)
 RELAX(LEN1,LINRLX,1.)
 RELAX(RHO1,LINRLX,1.)
 RELAX(GEN1,LINRLX,1.)
 RELAX(ENUT,LINRLX,1.)
 RELAX(DWDZ,LINRLX,1.)
 RELAX(DWDY,LINRLX,1.)
 RELAX(DWDX,LINRLX,1.)
 RELAX(DVDZ,LINRLX,1.)
 RELAX(DVDY,LINRLX,1.)
 RELAX(DVDX,LINRLX,1.)
 RELAX(DUDZ,LINRLX,1.)
 RELAX(DUDY,LINRLX,1.)
 RELAX(DUDX,LINRLX,1.)
 RELAX(WCRT,LINRLX,1.)
 RELAX(VCRT,LINRLX,1.)
 RELAX(UCRT,LINRLX,1.)
 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(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(SALT)=30. ;VARMIN(SALT)=10.
 VARMAX(UC1)=1.0E+10 ;VARMIN(UC1)=-1.0E+10
 VARMAX(VC1)=1.0E+10 ;VARMIN(VC1)=-1.0E+10
 VARMAX(WC1)=1.0E+10 ;VARMIN(WC1)=-1.0E+10
 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
 VARMAX(LEN1)=1.0E+10 ;VARMIN(LEN1)=-1.0E+10
 VARMAX(RHO1)=1.0E+10 ;VARMIN(RHO1)=-1.0E+10
 VARMAX(GEN1)=1.0E+10 ;VARMIN(GEN1)=-1.0E+10
 VARMAX(ENUT)=1. ;VARMIN(ENUT)=-1.0E+10
 VARMAX(DWDZ)=1.0E+10 ;VARMIN(DWDZ)=-1.0E+10
 VARMAX(DWDY)=1.0E+10 ;VARMIN(DWDY)=-1.0E+10
 VARMAX(DWDX)=1.0E+10 ;VARMIN(DWDX)=-1.0E+10
 VARMAX(DVDZ)=1.0E+10 ;VARMIN(DVDZ)=-1.0E+10
 VARMAX(DVDY)=1.0E+10 ;VARMIN(DVDY)=-1.0E+10
 VARMAX(DVDX)=1.0E+10 ;VARMIN(DVDX)=-1.0E+10
 VARMAX(DUDZ)=1.0E+10 ;VARMIN(DUDZ)=-1.0E+10
 VARMAX(DUDY)=1.0E+10 ;VARMIN(DUDY)=-1.0E+10
 VARMAX(DUDX)=1.0E+10 ;VARMIN(DUDX)=-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
 GENK = T
 PARSOL = F
 LSG4 = T
 LSG10 = T
 ISG62 = 1
 SPEDAT(SET,GXMONI,TRANSIENT,L,F)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 SPEDAT(SET,MATERIAL,199,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(U1,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(SALT,Y,N,Y,Y,Y,Y)
 OUTPUT(UC1,Y,N,Y,Y,Y,Y)
 OUTPUT(VC1,Y,N,Y,Y,Y,Y)
 OUTPUT(WC1,Y,N,Y,Y,Y,Y)
 OUTPUT(PRPS,Y,N,Y,N,N,N)
 OUTPUT(EPKE,Y,N,Y,N,N,N)
 OUTPUT(LEN1,Y,N,Y,N,N,N)
 OUTPUT(RHO1,Y,N,Y,N,N,N)
 OUTPUT(GEN1,Y,N,Y,N,N,N)
 OUTPUT(ENUT,Y,N,N,Y,Y,Y)
 OUTPUT(DWDZ,Y,N,Y,N,N,N)
 OUTPUT(DWDY,Y,N,Y,N,N,N)
 OUTPUT(DWDX,Y,N,Y,N,N,N)
 OUTPUT(DVDZ,Y,N,Y,N,N,N)
 OUTPUT(DVDY,Y,N,Y,N,N,N)
 OUTPUT(DVDX,Y,N,Y,N,N,N)
 OUTPUT(DUDZ,Y,N,Y,N,N,N)
 OUTPUT(DUDY,Y,N,Y,N,N,N)
 OUTPUT(DUDX,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 = 7 ;IYMON = 5 ;IZMON = 5
 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NTPRIN = 100000 ;ISTPRF = 1 ;ISTPRL = 100000
 NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
 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
 IDISPA = 10 ;IDISPB = 0 ;IDISPC = 0
 CSG1    ='URB'
STOP