PHOTON USE
  p
 
 
 
  gr ou X 1
  VEC X 1 SH
  msg         WATER RECIRCULATION IN A STILLING POND
  msg                     Velocity vectors
  pause
  vec off;red
  con rho1 x 1 fil;0.1
  msg                     Water depth variation (*10^3 m))
  msg Press e to END
  ENDUSE
 
    GROUP 1. Run title and other preliminaries
TEXT( RECIRCULATION IN A STILLING POND:109
 
  DISPLAY
    This particular case deals with the viscous recirculating
    shallow   water flow under influence of bottom stresses.
    It shows the way of introducing the non-linear sources for 2D
    compressible flow model to include the effects of 3rd dimension
  ENDDIS

REAL(WIN,HIN,HEX)
HIN=0.2;WIN=1.
HEX=HIN
    GROUP 4. Y-direction grid specification
GRDPWR(Y,14,10.0,1.0)
    GROUP 5. Z-direction grid specification
GRDPWR(Z,24,20.0,1.0)
    GROUP 7. Variables stored, solved & named
SOLVE(P1,V1,W1);STORE(HPOR,RHO1,PRPS)
SOLUTN(P1,Y,Y,Y,N,N,N)
    GROUP 9. Properties of the medium (or media)
REAL(RHOR,CON);RHOR=1000.
ENUL=1.e-06;ENUT=100.*ENUL
RHO2=1.;RHO1=COMPRESS;DRH1DP=RHO1
CON=9.81/(2.*RHOR)
RHO1A=CON**-0.5;PRESS0=9.81*RHOR*HEX**2/2.
RHO1B=0.5;RHO1C=0.0
VARMIN(P1)=0.001-CON
    GROUP 11.
INIADD=F
FIINIT(W1)=0.;FIINIT(V1)=0.
PATCH(INITH,INIVAL,1,1,1,NY/2,NZ/4,NZ/4,1,1)
INIT(INITH,PRPS,0.0,198.0)
    GROUP 13. Boundary conditions and special sources
  ** Inlet
PATCH(INLET,LOW,1,1,1,4,1,1,1,1)
COVAL(INLET,P1,FIXFLU,RHOR*WIN*HIN)
COVAL(INLET,W1,ONLYMS,WIN)
  ** Outlet
PATCH(OUTLET,HIGH,1,1,1,8,NZ,NZ,1,1)
COVAL(OUTLET,P1,FIXP*RHOR,0.0)

 
  PLANTBEGIN
  ** Bottom stresses
PATCH(BOTSTRES,EAST,1,1,1,NY,1,NZ,1,1000)
    CO=0.003*1000.*SQRT(W1**2+V1**2)
COVAL(BOTSTRES,W1,GRND,0.0)
    CO=0.003*1000.*SQRT(W1**2+V1**2)
COVAL(BOTSTRES,V1,GRND,0.0)
  PLANTEND
 
    GROUP 15. Termination of sweeps
LSWEEP=150
RESREF(P1)=1.E-6;RESREF(V1)=1.E-6;RESREF(W1)=1.E-6
    GROUP 16. Termination of iterations
LITER(P1)=10;ENDIT(P1)=1.0E-04
    GROUP 17. Under-relaxation devices
RELAX(P1,LINRLX,0.5);RELAX(V1,FALSDT,100.)
RELAX(W1,FALSDT,100.)
    GROUP 18. Limits on variables or increments to them
VARMIN(W1)=-1.E11;VARMAX(W1)=5.0
    GROUP 19. Data communicated by satellite to GROUND
NAMSAT=MOSG
    GROUP 20. Preliminary print-out
    GROUP 22. Spot-value print-out
IYMON=5;IZMON=10
NPRMON=LSWEEP
    GROUP 23. Field print-out and plot control
NPRINT=LSWEEP;ITABL=2;NPLT=5
tstswp=-1
dmpstk=t
DISTIL=T
EX(P1)=3.940E+01; EX(V1)=2.297E-01; EX(W1)=3.304E-01
EX(PRPS)=9.792E-01; EX(RHO1)=2.137E+02;
EX(HPOR)=9.792E-01
 LIBREF=109
STOP