```
TALK=T;RUN( 1, 1)
DISPLAY

WATER RECIRCULATION IN A STILLING POND

Steady free surface turbulent flow of water recirculating under
influence of bottom stresses and internal blockage in an open
pond.

The task is to simulate the velocity distribution along with free
water surface  elevation.

PLANT is used to introduce the momentum sinks caused by bottom
stresses.

ENDDIS

************************************************************
Q1 created by VDI menu, Version 3.3, Date 24/07/00
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
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
************************************************************
IRUNN   =       1 ;LIBREF =     302
************************************************************
Group 1. Run Title
TEXT( recirculation in a stilling pond       )
************************************************************
Group 2. Transience
************************************************************
Groups 3, 4, 5  Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,1,14,24)
* Set overall domain extent:
*        xulast  yvlast  zwlast
name
XSI= 1.000000E+00; YSI= 1.000000E+01; ZSI= 2.000000E+01
RSET(D,CHAM    )
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS  =    T
* Non-default variable names
NAME(149) =RHO1 ; NAME(150) =HPOR
* Solved variables list
SOLVE(P1  ,V1  ,W1  )
* Stored variables list
STORE(HPOR,RHO1)
SOLUTN(P1  ,Y,Y,Y,N,N,N)

************************************************************
Group 8. Terms & Devices
NEWRH1  =    T
************************************************************
Group 9. Properties
REAL(WIN,HIN)
HIN=0.5;WIN=1.
REAL(RHOR,CON);RHOR=1000.
ENUL=1.e-06;ENUT=100.*ENUL
RHO2=1.;RHO1=GRND3;DRH1DP=RHO1
CON=9.81/(2.*RHOR)
RHO1A=CON**-0.5;PRESS0=9.81*RHOR*HIN**2/2.
RHO1B=0.5;RHO1C=0.0
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************

PATCH(SHOWBOT,EAST,1,1,1,NY,1,NZ,1,1000)
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 11.Initialise Var/Porosity Fields
FIINIT(V1  ) =  0.000000E+00 ;FIINIT(W1  ) =  0.000000E+00
FIINIT(HPOR) =  1.000000E+00
No PATCHes used for this Group

************************************************************
Group 12. Convection and diffusion adjustments
No PATCHes used for this Group
************************************************************
Group 13. Boundary & Special Sources
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP  =     150
RESREF(P1  ) =  1.000000E-06 ;RESREF(V1  ) =  1.000000E-06
RESREF(W1  ) =  1.000000E-06
RESFAC  = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
LITER (P1  ) =   10
ENDIT (P1  ) =  1.000000E-04
************************************************************
Group 17. Relaxation
RELAX(P1  ,LINRLX, 5.000000E-01)
RELAX(V1  ,FALSDT, 1.000000E+01)
RELAX(W1  ,FALSDT, 1.000000E+01)
RELAX(RHO1,LINRLX, 3.000000E-01)
************************************************************
Group 18. Limits
VARMAX(P1  ) = 1.000000E+10 ;VARMIN(P1  ) =-3.905000E-03
VARMAX(V1  ) = 1.000000E+06 ;VARMIN(V1  ) =-1.000000E+06
VARMAX(W1  ) = 5.000000E+00 ;VARMIN(W1  ) =-1.000000E+11
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD  =    T  ;USEGRX =    T
ASAP    =    T
namsat=mosg
************************************************************
Group 20. Preliminary Printout
ECHO    =    T
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON   =       1 ;IYMON  =       5 ;IZMON  =      10
NPRMON  =     150
NPRMNT  =       1
TSTSWP  =      -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT  =     150
NPLT    =       5
ISWPRF  =       1 ;ISWPRL =  100000
ITABL   =       2
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
NOWIPE  =    T

GVIEW(P,-7.182436E-01,-6.541613E-01,-2.370635E-01)
GVIEW(UP,-6.853593E-01,7.239238E-01,7.884788E-02)

> DOM,    SIZE,        1.000000E+00, 1.000000E+01, 2.000000E+01
> DOM,    MONIT,       5.000000E-01, 3.214290E+00, 7.916670E+00
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> GRID,   RSET_Y_1,      4, 1.000000E+00
> GRID,   RSET_Y_2,      3, 1.000000E+00
> GRID,   RSET_Y_3,      1, 1.000000E+00
> GRID,   RSET_Y_4,      6, 1.000000E+00
> GRID,   RSET_Z_1,      5, 1.000000E+00
> GRID,   RSET_Z_2,      1, 1.000000E+00
> GRID,   RSET_Z_3,     18, 1.000000E+00
> DOM,    RELAX,       5.000000E-01

> OBJ1,   NAME,        INITH
> OBJ1,   POSITION,    0.000000E+00, 0.000000E+00, 5.000000E+00
> OBJ1,   SIZE,        1.000000E+00, 5.000000E+00, 0.000000E+00
> OBJ1,   CLIPART,     cube11
> OBJ1,   ROTATION24,        1
> OBJ1,   TYPE,        PLATE
> OBJ1,   POROSITY,     -1.000000E+00
> OBJ1,   SIDE,        BOTH

> OBJ2,   NAME,        INLET
> OBJ2,   POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ2,   SIZE,        1.000000E+00, 2.857143E+00, 0.000000E+00
> OBJ2,   CLIPART,     cube3t
> OBJ2,   ROTATION24,        1
> OBJ2,   TYPE,        INLET
> OBJ2,   DENSITY,       2.000000E+02
> OBJ2,   VELOCITY,      0.000000E+00, 0.000000E+00, 1.000000E+00
> OBJ2,   TEMPERATURE,   0.000000E+00

> OBJ3,   NAME,        OUTLET
> OBJ3,   POSITION,    0.000000E+00, 0.000000E+00, 2.000000E+01
> OBJ3,   SIZE,        1.000000E+00, 5.714286E+00, 0.000000E+00
> OBJ3,   CLIPART,     cube12t
> OBJ3,   ROTATION24,        1
> OBJ3,   TYPE,        OUTLET
> OBJ3,   PRESSURE,      0.000000E+00
> OBJ3,   TEMPERATURE,  -1.026000E+04
> OBJ3,   COEFFICIENT,   1.000000E+00

> OBJ4,   NAME,        RESTRIC
> OBJ4,   POSITION,    0.000000E+00, 0.000000E+00, 4.166667E+00
> OBJ4,   SIZE,        1.000000E+00, 5.000000E+00, 8.333330E-01
> OBJ4,   CLIPART,     cube14
> OBJ4,   ROTATION24,        1
> OBJ4,   TYPE,        BLOCKAGE
> OBJ4,   MATERIAL,      198
STOP
distil=t
EX(P1  )=   1.081E+02
EX(V1  )=   1.435E-01
EX(W1  )=   1.840E-01
EX(PRPS)=   5.104E+00
EX(RHO1)=   5.111E+02
EX(HPOR)=   9.792E-01

```