L(456)
TEXT(Variable Porosity, 2D             
TITLE
  In this case a new method of surface-wave emulation is activated,
  namely, varying the porosity of the cells according to the
  pressure in the cells. The porosity variations are to be thought
  of as adjustments to the heights of the cells needed to move to
  satisfy hydrostatics. Of all the methods, this is the most
  general.
NREGZ=1; IREGZ=1; GRDPWR(Z,1,2.5,1.0)
  * Add storage of POROSITIES since CONPOR creates them at the end
STORE(EPOR,NPOR)
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)
VARMIN(P1)=-0.9