GROUP 1. Run title and other preliminaries
TEXT(1D Hydrostatic; X-Direction      
TITLE
  DISPLAY
  This run analyses the situation in which a fluid is at rest
  within a container of non-uniform cross section area.
  Porosities are used to represent a linear variation of the area
  with the x-distance. A gravitational acceleration of unity acts
  along the x-direction. The results demonstrates accord with
  the laws of hydrostatics.
    Interesting variations include:change of the fluid density;
  and changes to the shape of the container.
  ENDDIS
 
    GROUP 3. X-direction grid specification
  ** X-direction grid distances vary as IX**2.0
GRDPWR(X,5,1.0,2.0)
 
    GROUP 7. Variables stored, solved & named
SOLVE(P1,U1)
  ** Provide storage for porosities on the east faces.
STORE(EPOR)
 
    GROUP 11. Initialization of variable or porosity fields
  ** Set varying-cross section by means of non-uniform
     porosities (varying linearly with IX).
INIADD=T     
PATCH(SECTION,LINVLX,1,NX,1,1,1,1,1,1);INIT(SECTION,EPOR,1.0,0.1)
 
    GROUP 13. Boundary conditions and special sources
  ** Fix external pressure
PATCH(ANCHOR,CELL,1,1,1,1,1,1,1,1)
COVAL(ANCHOR,P1,FIXP,0.02);COVAL(ANCHOR,U1,ONLYMS,0.0)
  ** Add gravitational force
PATCH(GRAVITY,PHASEM,1,NX,1,1,1,1,1,1);COVAL(GRAVITY,U1,FIXFLU,1.0)
 
    GROUP 15. Termination of sweeps
LSWEEP=20
 
    GROUP 17. Under-relaxation devices
RELAX(U1,FALSDT,1.0)
 
    GROUP 23. Field print-out and plot control
PATCH(LONGPLOT,PROFIL,1,NX,1,1,1,1,1,1)
PLOT(LONGPLOT,P1,0.0,0.0);PLOT(LONGPLOT,U1,0.0,0.0)