DISPLAY
  The problem is similar to that of library case 230; but the
  Saffman-Spalding k-W model is used, the variable being VOSQ,
  which stands for vorticity-fluctuations squared, and  which
  is related to KE and EPS as (KE/EPS) squared.
  ENDDIS
    GROUP 1. Run title and other preliminaries
TEXT(X-Spatial Decay Of Turbu, K-Vosq 
TITLE
REAL(XLEN,UIN,TKEIN,VOIN);XLEN=1.0;UIN=1.0; TKEIN=1.0; VOIN=1.0
    GROUP 3. X-direction grid specification
GRDPWR(X,40,XLEN,2.0)
    GROUP 7. Variables stored, solved & named
STORE(U1,VIST)
    ** The next statement activates the Safmann-Spalding
    ** two-equation turbulence model
TURMOD(KWMODL)
    GROUP 8. Terms (in differential equations) & devices
    ** cut out diffusion of KE and VOSQ
TERMS(KE,N,Y,N,N,N,N); TERMS(VOSQ,N,Y,N,N,N,N)
    GROUP 11. Initialization of variable or porosity fields
FIINIT(U1)=UIN; FIINIT(KE)=TKEIN; FIINIT(VOSQ)=VOIN
    GROUP 13. Boundary conditions and special sources
    ** Inlet Boundary
PATCH(INLET,WEST,1,1,1,1,1,1,1,1)
COVAL(INLET,KE,1.E5,TKEIN); COVAL(INLET,VOSQ,1.E5,VOIN)
    GROUP 16. Termination of iterations
LSWEEP=10; RESREF(KE)=1.E-7; RESREF(VOSQ)=1.E-7
    GROUP 21. Print-out of variables
OUTPUT(U1,N,N,N,N,N,N); OUTPUT(VIST,N,N,N,N,N,N)
    GROUP 22. Monitor print-out
IXMON=NX; ITABL=1; TSTSWP=LSWEEP/2
    GROUP 23. Field print-out and plot control
IPROF=1; NXPRIN=2; ORSIZ=0.4
PATCH(PROFILE,PROFIL,1,NX,1,1,1,1,1,1)
PLOT(PROFILE,KE,0.0,TKEIN); PLOT(PROFILE,VOSQ,0.0,2.0*VOIN)
PLOT(PROFILE,VIST,0.0,0.0)