PHOTON USE
  p
  n4
 
 
 
  up z
  gr x 1
  msg contours of vfol at time = 2.0s
  con vfol x 1 sh;0.45 0.55 40
  msg velocity vectors
  vec x 1
  pause
  p
  n8
 
 
 
  up z
  gr x 1
  msg contours of vfol at time = 4.0s
  con vfol x 1 sh;0.45 0.55 40
  msg velocity vectors
  vec x 1
  pause
  p
  n12
 
 
 
  up z
  gr x 1
  msg contours of vfol at time = 6.0s
  con vfol x 1 sh;0.45 0.55 40
  msg velocity vectors
  vec x 1
  pause
  p
  n16
 
 
 
  up z
  gr x 1
  msg contours of vfol at time = 8.0s
  con vfol x 1 sh;0.45 0.55 40
  msg velocity vectors
  vec x 1
  pause
  p
  n20
 
 
 
  up z
  gr x 1
  msg contours of vfol at time = 10.0s
  con vfol x 1 sh;0.45 0.55 40
  msg velocity vectors
  vec x 1
  msg            -
  msg Press e to END
  enduse
 
    GROUP 1. Run title and other preliminaries
TEXT(FILLING OF A 2D DUCT BY HOL:P100
TITLE
  DISPLAY
  FREE SURFACE FLOW - Filling of a duct by H(eight) O(f) L(iguid)
  method
 
  2-dimensional (y-z), Cartesian, transient, elliptic simulation
 
  Vessel dimensions: 0.5 m x 0.5 m x 1m (height). Liquid enters
  at the bottom corner.
                 +-----------+
                 |           |
                 |           |      z ^
                 |           |        |
                 |           |        +---->
                 |           |             y
                 +-------  ^ +
                           |  liquid inflow
  enddis
 
REAL(VIN,TIMEINT);VIN=0.1
    GROUP 2. Transience; time-step specification
STEADY=F;LSTEP=50
TFRAC(1)=-LSTEP;TFRAC(2)=0.05/VIN
TIMEINT=TFRAC(2);LSTEP=20
    GROUP 3. X-direction grid specification
XULAST=0.5
    GROUP 4. Y-direction grid specification
GRDPWR(Y, 5,0.5,1.0)
    GROUP 5. Z-direction grid specification
GRDPWR(Z,10,1.0,1.0)
    GROUP 7. Variables stored, solved & named
STORE(DEN1,PRPS);SOLVE(VFOL)
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 (in differential equations) & devices
GALA=T;TERMS(VFOL,N,N,N,N,P,P)
    GROUP 9. Properties of the medium (or media)

    GROUP 11. Initialization of variable or porosity fields
FIINIT( P1 )=0.0;FIINIT( V1 )=0.0;FIINIT( W1 )=VIN
FIINIT(VFOL)=0.0;FIINIT(DEN1)=1.161
    GROUP 13. Boundary conditions and special sources
REAL(NCELLS);NCELLS=1
PATCH(IFLO,LOW,1,NX,NY+1-NCELLS,NY,1,1,1,LSTEP)
COVAL(IFLO,P1,FIXFLU,VIN*1000.5);COVAL(IFLO,W1,ONLYMS,VIN)
COVAL(IFLO,VFOL,ONLYMS,1.0/1000.5)
PATCH(GRAV,PHASEM,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(GRAV,W1,FIXFLU,-9.81)
PATCH(REFP,CELL,1,NX,NY,NY,NZ,NZ,1,LSTEP)
COVAL(REFP,P1,GRND,ZERO)
    GROUP 15. Termination of sweeps
LSWEEP=10
    GROUP 16. Termination of iterations
RESREF(P1)=1.E-8;RESREF(V1)=1.E-8;RESREF(W1)=1.E-8
    GROUP 17. Under-relaxation devices
RELAX(V1,FALSDT,0.001);RELAX(W1,FALSDT,0.001)
    GROUP 19. Data communicated by satellite to GROUND
IDISPB=1;IDISPC=lstep;IDISPA=4;CSG1=N;INIFLD=F
IPRPSA=67;IPRPSB=0;HOL=T
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
    GROUP 22. Spot-value print-out
TSTSWP=-1;NTPRIN=LSTEP/2
NPRMON=LSWEEP;IYMON=NY;IZMON=6
    GROUP 23. Field print-out and plot control
OUTPUT(P1,Y,Y,Y,Y,Y,Y);OUTPUT(V1,Y,Y,Y,Y,Y,Y)
OUTPUT(W1,Y,Y,Y,Y,Y,Y);OUTPUT(VFOL,Y,N,y,N,N,N)
OUTPUT(DEN1,Y,N,y,N,N,N);IHOLA=3