GROUP 1. Run title and other preliminaries
TEXT(HYD. EQ. In YZ Tank With Blockages
TITLE
             Density of the liquid ---------- DLIQ
REAL(DLIQ);DLIQ=1.0
 
    GROUP 2. Transience; time-step specification
STEADY=F;GRDPWR(T,1,0.1,1.0)
 
    GROUP 4. Y-direction grid specification
GRDPWR(Y,6,6.0,1.0)
 
    GROUP 5. Z-direction grid specification
GRDPWR(Z,7,8.0,1.0)
 
    GROUP 7. Variables stored, solved & named
SOLVE(P1,V1,W1);SOLUTN(P1,Y,Y,Y,N,N,N)
 
    GROUP 9. Properties of the medium (or media)
RHO1=DLIQ;ENUL=1.0E-06
 
    GROUP 11. Initialization of variable or porosity fields
PATCH(PINIT,LINVLZ,1,1,1,NY,1,NZ,1,1);INIT(PINIT,P1,9.81,0.0)
CONPOR(0.5,HIGH,1,1,1,1,4,4);CONPOR(0.5,HIGH,1,1,6,6,4,4)
CONPOR(0.75,VOLUME,1,1,1,1,4,4);CONPOR(0.75,VOLUME,1,1,6,6,4,4)
CONPOR(0.5,VOLUME,1,1,1,1,5,5);CONPOR(0.5,VOLUME,1,1,6,6,5,5)
 
    GROUP 13. Boundary conditions and special sources
   ** Reference pressure at the top of the tank (low-y)
PATCH(REFP,CELL,1,1,1,NY,1,1,1,LSTEP);COVAL(REFP,P1,FIXVAL,0.0)
COVAL(REFP,V1,ONLYMS,0.0);COVAL(REFP,W1,ONLYMS,0.0)
   ** Gravity
PATCH(GRAV,PHASEM,1,1,1,NY,1,NZ,1,LSTEP)
COVAL(GRAV,W1,FIXFLU,9.81)
 
    GROUP 15. Termination of sweeps
LSWEEP=5
 
    GROUP 17. Under-relaxation devices
RELAX(V1,FALSDT,0.1);RELAX(W1,FALSDT,0.1)
SPEDAT(SET,GXMONI,TRANSIENT,L,F) 
    GROUP 21. Print-out of variables
OUTPUT(HPOR,Y,N,N,N,N,N);OUTPUT(VPOR,Y,N,N,N,N,N)
    GROUP 22. Spot-value print-out
IXMON=5;IYMON=5