TALK=T;RUN( 1, 1)
  DISPLAY

  WATER POURING INTO A BUND

  FREE SURFACE EXAMPLE - Water rise in a channel with
  obstruction: 3D, elliptic, Scalar-Equation-Model.

  The case represents a water flow within an open channel.
  The channel has an immersed solid obstruction. The model
  predicts the elevation of the water surface over the
  blockage and the gravity-driven sideway 'spills'.

  ENDDIS

 ************************************************************
   Q1 created by VDI menu, Version 3.3, Date 24/07/00
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
  Echo DISPLAY / USE settings
  PHOTON USE
  p
  n4



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 0.08s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n8



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 0.16s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n12



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 0.24s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n16



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 0.32s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n20



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 0.40s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n24



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 0.48s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n28



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 0.56s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n32



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 0.64s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n36



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 0.72s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n40



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 0.8s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n44



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 0.88s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n48



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 0.96s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n52



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 1.04s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n56



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 1.12s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n60



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 1.2s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n64



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 1.28s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n68



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 1.36s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n72



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 1.44s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n76



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 1.52s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n80



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 1.6s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n84



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 1.68s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n88



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 1.76s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n92



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 1.84s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n96



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 1.92s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  pause
  p
  n100



  up z
  gr x 1
  use patgeo
  msg contours of den1 at time = 2.00s
  con den1 z 1 fi;.001
  con den1 x 1 fi;.001
  con den1 x 4 fi;.001
  con den1 x 8 fi;.001
  msg velocity vectors
  vec x 1;vec z 1
  msg            -
  msg Press e to END
  enduse

 ************************************************************
 IRUNN   =       1 ;LIBREF =     106
 ************************************************************
  Group 1. Run Title
 TEXT(Water pouring into a bund)
 ************************************************************
  Group 2. Transience
 STEADY=F
    * Set overall time and no. of steps
 RSET(U,0.000000E+00,2.000000E+00,100)
    * Modify regions
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,10,10,10)
    * Set overall domain extent:
    *        xulast  yvlast  zwlast
            name
 XSI= 1.000000E+00; YSI= 1.000000E+00; ZSI= 3.125000E-01
 RSET(D,CHAM    )
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 NAME(147) =SURN ; NAME(148) =VFOL
 NAME(149) =PRPS ; NAME(150) =DEN1
    * Solved variables list
 SOLVE(P1  ,U1  ,V1  ,W1  ,SURN,VFOL)
    * Stored variables list
 STORE(DEN1,PRPS)
    * Additional solver options
 SOLUTN(P1  ,Y,Y,Y,N,N,N)
 SOLUTN(V1  ,Y,Y,N,N,N,N)
 SOLUTN(W1  ,Y,Y,N,N,N,N)

 ************************************************************
  Group 8. Terms & Devices
 TERMS (SURN,N,N,N,N,N,Y)
 TERMS (VFOL,N,N,N,N,Y,Y)
 GALA    =    T
 NEWRH1  =    T
 NEWENL  =    T
 RUPLIM  = 6.000000E-01 ;RLOLIM = 4.000000E-01
 ************************************************************
  Group 9. Properties
 DVO1DT  = 3.330000E-03
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(P1  ) =  0.000000E+00 ;FIINIT(V1  ) =  0.000000E+00
 FIINIT(W1  ) =  0.000000E+00 ;FIINIT(SURN) =  0.000000E+00
 FIINIT(PRPS) =  0.000000E+00 ;FIINIT(DEN1) =  1.189000E+00
   No PATCHes used for this Group


 INIADD  =    F
 ************************************************************
  Group 12. Convection and diffusion adjustments
   No PATCHes used for this Group
 ************************************************************
  Group 13. Boundary & Special Sources

 PATCH (BUOYANCY,PHASEM,0,0,0,0,0,0,1,100)
 COVAL (BUOYANCY,U1  , FIXFLU      , GRND1       )
 COVAL (BUOYANCY,V1  , FIXFLU      , GRND1       )
 COVAL (BUOYANCY,W1  , FIXFLU      , GRND1       )

 BUOYA   = 0.000000E+00 ; BUOYB = 0.000000E+00
 BUOYC   =-9.810000E+00
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =      40
 RESREF(P1  ) =  1.000000E-06 ;RESREF(V1  ) =  1.000000E-06
 RESREF(W1  ) =  1.000000E-06
 RESFAC  = 1.000000E-03
 ************************************************************
  Group 16. Terminate Iterations
 LITER (SURN) =    1
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 5.000000E-01)
 RELAX(U1  ,FALSDT, 3.000000E-02)
 RELAX(V1  ,FALSDT, 3.000000E-02)
 RELAX(W1  ,FALSDT, 3.000000E-02)
 RELAX(SURN,FALSDT, 1.000000E+01)
 RELAX(VFOL,FALSDT, 1.000000E+00)
 RELAX(DEN1,LINRLX, 4.000000E-01)
 ************************************************************
  Group 18. Limits
 VARMAX(U1  ) = 1.000000E+06 ;VARMIN(U1  ) =-1.000000E+06
 VARMAX(V1  ) = 1.000000E+06 ;VARMIN(V1  ) =-1.000000E+06
 VARMAX(W1  ) = 1.000000E+06 ;VARMIN(W1  ) =-1.000000E+06
 VARMAX(SURN) = 1.000000E+00 ;VARMIN(SURN) = 0.000000E+00
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 USEGRD  =    T  ;USEGRX =    T
 SURF    =    T
 ASAP    =    T
 IPRPSA  =      67
 CSG1    ='N'
 SPEDAT(SET,GXMONI,TRANSIENT,L,F)
 ************************************************************
  Group 20. Preliminary Printout
 ECHO    =    T
 distil=t
 EX(P1  )=   4.275E+02  
 EX(U1  )=   3.480E-01  
 EX(V1  )=   7.489E-01  
 EX(W1  )=   3.694E-01  
 EX(SURN)=   4.914E-01  
 EX(VFOL)=   4.885E-01  
 EX(PRPS)=   3.839E+01  
 EX(DEN1)=   5.117E+02  
 ************************************************************
  Group 21. Print-out of Variables
 OUTPUT(P1  ,Y,Y,Y,Y,Y,Y)
 OUTPUT(U1  ,Y,Y,Y,N,Y,Y)
 OUTPUT(V1  ,Y,Y,Y,Y,Y,Y)
 OUTPUT(W1  ,Y,Y,Y,Y,Y,Y)
 OUTPUT(SURN,Y,Y,Y,Y,Y,Y)
 OUTPUT(VFOL,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =       5 ;IYMON  =       5 ;IZMON  =       5
 NPRMON  =  100000
 NPRMNT  =       1
 TSTSWP  =      -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT  =  100000
 NTPRIN  =      10 ;ISTPRF =       1 ;ISTPRL =   10000
 ISWPRF  =       1 ;ISWPRL =  100000
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 NOWIPE  =    T
 IDISPA  =       2 ;IDISPB =       0 ;IDISPC =       0
 CSG1    ='N'

 GVIEW(P,-9.841644E-01,0.000000E+00,1.772584E-01)
 GVIEW(UP,1.772584E-01,0.000000E+00,9.841644E-01)

> DOM,    SIZE,        1.000000E+00, 1.000000E+00, 3.125000E-01
> DOM,    MONIT,       4.500000E-01, 4.500000E-01, 1.406250E-01
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> DOM,    RELAX,       1.000000E+00

> OBJ1,   NAME,        SOUTHOUT
> OBJ1,   POSITION,    0.000000E+00, 0.000000E+00, 9.375000E-02
> OBJ1,   SIZE,        1.000000E+00, 0.000000E+00, 2.187500E-01
> OBJ1,   CLIPART,     cube13
> OBJ1,   ROTATION24,        1
> OBJ1,   GRID,        2.000000E+00
> OBJ1,   TYPE,        PLATE
> OBJ1,   SURF_TEMP,     0.000000E+00,-1.026000E+04

> OBJ2,   NAME,        HIGHOUT
> OBJ2,   POSITION,    0.000000E+00, 0.000000E+00, 3.125000E-01
> OBJ2,   SIZE,        1.000000E+00, 1.000000E+00, 0.000000E+00
> OBJ2,   CLIPART,     cube12t
> OBJ2,   ROTATION24,        1
> OBJ2,   TYPE,        OUTLET
> OBJ2,   PRESSURE,      0.000000E+00
> OBJ2,   TEMPERATURE,  -1.026000E+04
> OBJ2,   COEFFICIENT,   1.000000E+00
> OBJ2,   TIME_LIMITS,   ALWAYS_ACTIVE

> OBJ3,   NAME,        NORTHOUT
> OBJ3,   POSITION,    8.000000E-01, 1.000000E+00, 0.000000E+00
> OBJ3,   SIZE,        2.000000E-01, 0.000000E+00, 3.125000E-01
> OBJ3,   CLIPART,     cube13
> OBJ3,   ROTATION24,        1
> OBJ3,   TYPE,        PLATE
> OBJ3,   SURF_TEMP,     0.000000E+00,-1.026000E+04

> OBJ4,   NAME,        HEAVYINL
> OBJ4,   POSITION,    5.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ4,   SIZE,        5.000000E-01, 0.000000E+00, 9.375000E-02
> OBJ4,   CLIPART,     cube13
> OBJ4,   ROTATION24,        1
> OBJ4,   GRID,        2.000000E+00
> OBJ4,   TYPE,        PLATE
> OBJ4,   SURF_TEMP,     0.000000E+00, 0.000000E+00
> OBJ4,   VELOCITY,      0.000000E+00, 2.000000E+00, 0.000000E+00

> OBJ5,   NAME,        BLOCK
> OBJ5,   POSITION,    0.000000E+00, 5.000000E-01, 0.000000E+00
> OBJ5,   SIZE,        1.000000E+00, 9.999999E-02, 9.450000E-02
> OBJ5,   CLIPART,     cube4
> OBJ5,   ROTATION24,        1
> OBJ5,   GRID,        2.000000E+00
> OBJ5,   TYPE,        BLOCKAGE
> OBJ5,   MATERIAL,      100
> OBJ5,   TIME_LIMITS,   ALWAYS_ACTIVE
> OBJ5,   FIXED_TMP,     0.000000E+00, 0.000000E+00

> OBJ6,   NAME,        EASTOUT
> OBJ6,   POSITION,    1.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ6,   SIZE,        0.000000E+00, 1.000000E+00, 3.125000E-01
> OBJ6,   CLIPART,     cube13
> OBJ6,   ROTATION24,        1
> OBJ6,   GRID,        2.000000E+00
> OBJ6,   TYPE,        PLATE
> OBJ6,   FIXED_TMP,     0.000000E+00, 0.000000E+00

> OBJ7,   NAME,        B7
> OBJ7,   POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ7,   SIZE,        5.000000E-01, 0.000000E+00, 9.375000E-02
> OBJ7,   CLIPART,     cube3t
> OBJ7,   ROTATION24,        1
> OBJ7,   GRID,        2.000000E+00
> OBJ7,   TYPE,        INLET
> OBJ7,   FLUID,         1.000000E+00
> OBJ7,   VELOCITY,      0.000000E+00, 2.000000E+00, 0.000000E+00
> OBJ7,   TEMPERATURE,   0.000000E+00
> OBJ7,   TIME_LIMITS,   ALWAYS_ACTIVE

> OBJ8,   NAME,        B8
> OBJ8,   POSITION,    2.000000E-01, 1.000000E+00, 0.000000E+00
> OBJ8,   SIZE,        6.000000E-01, 0.000000E+00, 6.000000E-02
> OBJ8,   CLIPART,     cube12t
> OBJ8,   ROTATION24,        1
> OBJ8,   GRID,        2.000000E+00
> OBJ8,   TYPE,        OUTLET
> OBJ8,   PRESSURE,      0.000000E+00
> OBJ8,   TEMPERATURE,  -1.026000E+04
> OBJ8,   COEFFICIENT,   1.000000E+00
> OBJ8,   TIME_LIMITS,   ALWAYS_ACTIVE

> OBJ9,   NAME,        B10
> OBJ9,   POSITION,    0.000000E+00, 1.000000E+00, 0.000000E+00
> OBJ9,   SIZE,        2.000000E-01, 0.000000E+00, 3.125000E-01
> OBJ9,   CLIPART,     cube13
> OBJ9,   ROTATION24,        1
> OBJ9,   GRID,        2.000000E+00
> OBJ9,   TYPE,        PLATE
> OBJ9,   SURF_TEMP,     0.000000E+00,-1.026000E+04

> OBJ10,  NAME,        B12
> OBJ10,  POSITION,    2.000000E-01, 1.000000E+00, 6.000000E-02
> OBJ10,  SIZE,        6.000000E-01, 0.000000E+00, 2.525000E-01
> OBJ10,  CLIPART,     cube13
> OBJ10,  ROTATION24,        1
> OBJ10,  GRID,        2.000000E+00
> OBJ10,  TYPE,        PLATE
> OBJ10,  SURF_TEMP,     0.000000E+00,-1.026000E+04
STOP