L(457)
TEXT(Variable Porosity, 3D             
TITLE
  NZ is set back to 5 making the simulation 3D. The top layer of
  cells are stretched & contracted by porosity adjustments so as
  to fit the free-surface waves.
    To read this run type L(458)
NREGZ=3
IREGZ=1; GRDPWR(Z,1,0.5,1.0)
IREGZ=2; GRDPWR(Z,2,1.0,1.0)
IREGZ=3; GRDPWR(Z,2,1.0,1.0)
STORE(HPOR)
RHO1=1.0;DRH1DP=1.0;IPORIA=1;PORIA=1.0;PORIB=1.0
  * Re-define Patches with new limits
INLET(UPSTRM,WEST,#1,#1,#1,#NREGY,#1,#NREGZ,1,1)
PATCH(DWNSTRM,EAST,#NREGX,#NREGX,#1,#NREGY,#1,#NREGZ,1,1)
 
OUTPUT(VPOR,Y,N,N,N,N,N)
RELAX(W1,FALSDT,1.0)
LSWEEP=5;VARMIN(P1)=-0.9