GROUP 1. Run title and other preliminaries
    use for error report number 171
 
TEXT(Filling Of A 2d Duct - HOL Method 
TITLE
REAL(VIN); VIN=0.1
    GROUP 2. Transience; time-step specification
STEADY=T
    GROUP 3. X-direction grid specification
GRDPWR(X,2,0.2,1.0)
    GROUP 5. Z-direction grid specification
GRDPWR(Z,3,0.3,1.0)
    GROUP 7. Variables stored, solved & named
SOLUTN(P1,Y,Y,Y,N,N,N); SOLUTN(U1,Y,Y,N,N,N,Y)
SOLUTN(W1,Y,Y,N,N,N,Y)
    GROUP 11. Initialization of variable or porosity fields
FIINIT(P1)=0.0; FIINIT(U1)=0.0; FIINIT(W1)=VIN
    GROUP 13. Boundary conditions and special sources
PATCH(IFLO,LOW,1,NX,1,NY,1,1,1,LSTEP)
COVAL(IFLO,P1,FIXFLU,VIN); COVAL(IFLO,W1,ONLYMS,VIN)
PATCH(GRAV,PHASEM,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(GRAV,W1,FIXFLU,-9.81)
PATCH(REFP,CELL,1,NX,1,NY,NZ,NZ,1,LSTEP)
COVAL(REFP,P1, FIXP ,0.00)
    GROUP 15. Termination of sweeps
LSWEEP=20
    GROUP 16. Termination of iterations
RESREF(P1)=1.E-8; RESREF(U1)=1.E-8; RESREF(W1)=1.E-8
    GROUP 17. Under-relaxation devices
RELAX(U1,FALSDT,1.E-1); RELAX(W1,FALSDT,1.E-1)

    GROUP 22. Spot-value print-out
TSTSWP=1; NTPRIN=LSTEP/5; NPRMON=LSWEEP; IYMON=NY; IZMON=6
    GROUP 23. Field print-out and plot control
OUTPUT(P1,Y,Y,Y,Y,Y,Y); OUTPUT(U1,Y,Y,Y,Y,Y,Y)
OUTPUT(W1,Y,Y,Y,Y,Y,Y)
    GROUP 24. Dumps for restarts