GROUP 1. Run title and other preliminaries
 
  DISPLAY
    This library case is the first of a two-stage calculation.
    It calculates steady flow in a section of a circular tube.
    Library case 487 picks up the resulting flow field and solves
    the equations for the time-dependent concentration fields which
    result from the sudden injection of a dye.
  ENDDIS
 
TEXT(Test Of Whole Field Solver        
TITLE
    GROUP 2. Transience; time-step specification
TFIRST=0.0
    GROUP 3. X-direction grid specification
CARTES=F;GRDPWR(X,2,0.01,1.0)
    GROUP 4. Y-direction grid specification
GRDPWR(Y,5,0.2,1.0)
    GROUP 5. Z-direction grid specification
GRDPWR(Z,5,1.0,1.0)
    GROUP 7. Variables stored, solved & named
SOLVE(P1,U1,V1,W1);STORE(C1)
    GROUP 8. Terms (in differential equations) & devices
TERMS(C1,Y,Y,N,Y,Y,N)
    GROUP 9. Properties of the medium (or media)
ENUT=1.E-3
    GROUP 13. Boundary conditions and special sources
REAL(VELIN);VELIN=20.0
    Inlet 1
PATCH(INLET,HIGH,1,NX,1,2,1,1,1,1)
COVAL(INLET,P1,FIXFLU,VELIN);COVAL(INLET,U1,ONLYMS,0.0)
COVAL(INLET,V1,ONLYMS,0.0);COVAL(INLET,W1,ONLYMS,VELIN)
    Outlet Boundary
PATCH(OUTLET,HIGH,1,NX,1,NY,NZ,NZ,1,1)
COVAL(OUTLET,P1,FIXP,0.0)
    GROUP 15. Termination of sweeps
FSWEEP=1;LSWEEP=30
    GROUP 16. Termination of iterations
REAL(MDOT);MDOT=RHO1*VELIN
RESREF(P1)=MDOT*1.E-12;RESREF(C1)=MDOT*1.E-12
RESREF(U1)=VELIN*MDOT*1.E-12;RESREF(V1)=VELIN*MDOT*1.E-12
RESREF(W1)=VELIN*MDOT*1.E-12
LITER(P1)=25;LITER(C1)=-30;LITER(U1)=2;LITER(V1)=2;LITER(W1)=2
    GROUP 17. Under-relaxation devices
REAL(RELAXN);RELAXN=ZWLAST/VELIN/NZ;RELAX(P1,LINRLX,1.0)
RELAX(U1,FALSDT,RELAXN*5.);RELAX(V1,FALSDT,RELAXN*5.)
RELAX(W1,FALSDT,RELAXN*5.)
    GROUP 21. Print-out of variables
OUTPUT(C1,Y,Y,Y,Y,Y,Y)
    GROUP 22. Spot-value print-out
IXMON=1;IYMON=3;IZMON=3;NPRMON=1;TSTSWP=1
    GROUP 23. Field print-out and plot control
ITABL=3;IPLTL=LSWEEP;IPLTF=FSWEEP;NPRINT=LSWEEP
IZPRF=1;IZPRL=NZ;NUMCLS=8;NPLT=1
    GROUP 24. Dumps for restarts