GROUP 1. Run title and other preliminaries
TEXT(Hydrostatic Equil In Annular Cav. 
TITLE
  DISPLAY
  This case tests the correctness of the gravitational terms in the
  momentum equations when the grid is polar, coarse and non-uniform.
  Since there is no inlet or outlet, and the density is uniform, all
  velocity components should equal zero (insofar as round-off error
  permits).
 
  The momentum-source terms are computed in subroutine GXBUOY.
 
 
  ENDDIS
    GROUP 3. X-direction grid specification
CARTES=F; GRDPWR(X,5,3.1416,1.5)
    GROUP 4. Y-direction grid specification
GRDPWR(Y,-5,0.1,2.0); RINNER=0.01
    GROUP 7. Variables stored, solved & named
#solvel
    GROUP 9. Properties of the medium (or media)
ENUL=1.0E-03
    GROUP 13. Boundary conditions and special sources
    ** Gravitational force
#gravity
gravdir=2
#density
   ** Reference pressure
PATCH(REFP,CELL,NX,NX,NY,NY,1,1,1,1); COVAL(REFP,P1,FIXVAL,0.0)
COVAL(REFP,U1,ONLYMS,0.0); COVAL(REFP,V1,ONLYMS,0.0)
    GROUP 15. Termination of sweeps
LSWEEP=200
RESFAC=1.E-6
    GROUP 17. Under-relaxation devices
RELAX(U1,FALSDT,0.1);RELAX(V1,FALSDT,0.1)
    GROUP 22. Spot-value print-out
IXMON=3;IYMON=4;tstswp=-1