GROUP 1. Run title and other preliminaries
TEXT(Forced Vortex In A Pipe - Inviscid
TITLE
REAL(OMEGA,PLEN,PRAD,WIN)
PRAD=0.5;PLEN=1.0;WIN=30.0;OMEGA=500.0
    GROUP 3. X-direction grid specification
CARTES=F
NX=3;GRDPWR(X,NX,0.5,1.0)
    GROUP 4. Y-direction grid specification
NY=5;GRDPWR(Y,NY,PRAD,1.0)
    GROUP 5. Z-direction grid specification
NZ=3;GRDPWR(Z,NZ,PLEN,1.0)
    GROUP 7. Variables stored, solved & named
SOLVE(P1,U1,V1,W1);STORE(DEN1)
    GROUP 9. Properties of the medium (or media)
ENUL=1.0;ENUT=0.0;DIFCUT=0.0
    GROUP 11. Initialization of variable or porosity fields
IURINI=-1;FIINIT(W1)=WIN;FIINIT(U1)=OMEGA;INIADD=F
FIINIT(DEN1)=RHO1;INIFLD=T;FIINIT(V1)=0.0;FIINIT(P1)=3.E3
    GROUP 13. Boundary conditions and special sources
XCYCLE=T
IURVAL=-1
PATCH(INLET,LOW,1,NX,1,NY,1,1,1,1)
COVAL(INLET,P1,FIXFLU,RHO1*WIN)
COVAL(INLET,U1,ONLYMS,OMEGA)
COVAL(INLET,W1,ONLYMS,WIN)
PATCH(EQUIRADI,HIGH,1,NX,1,NY,NZ,NZ,1,1)
COVAL(EQUIRADI,P1,0.001,0.0)
    GROUP 15. Termination of sweeps
LSWEEP=100;TSTSWP=10
RESREF(P1)=1E-12;RESREF(U1)=1E-12
RESREF(V1)=1E-12;RESREF(W1)=1E-12
    GROUP 16. Termination of iterations
LITER(P1)=200;LITER(U1)=5
LITER(V1)=5;LITER(W1)=5
    GROUP 17. Under-relaxation devices
RELAX(P1,LINRLX,1.0)
RELAX(U1,FALSDT,1.0E-2)
RELAX(V1,FALSDT,1.0E-20)
RELAX(W1,FALSDT,1.0E-1)
    GROUP 21. Print-out of variables
NXPRIN=1;NYPRIN=1;NZPRIN=1
IURPRN=-1
    GROUP 22. Spot-value print-out
IYMON=NY/2;IZMON=NZ-1;ITABL=3;NPLT=1
    GROUP 23. Field print-out and plot control
NPRINT=LSWEEP
ISG21=3 ! ensure at least 3 sweeps