L(170)
    GROUP 1. Run title and other preliminaries
TEXT(Channel Flow  K-L Model           
TITLE
  DISPLAY
  This is like case 170; the Prandtl one-equation model is used with
  the mixing length specified according to Nikuradse's formula.
  ENDDIS
 
  GXLEN, GXENUT, GXKESO, GXGENK, GXSQUR and GXWFUN are used for
  this case.
 
REAL(TKEIN)
 
    GROUP 7. Variables stored,named and solved
TURMOD(KLMODL);SOLVE(TEMP)
EL1=NIKURPIP
WD23=WD2**3;EL1A=0.14*WD2;EL1B=-0.08/WD2;EL1C=-0.06/WD23
 
    GROUP 8. Terms (in differential equations) & devices
TERMS(TEMP,N,Y,Y,Y,Y,Y)
 
    GROUP 11. Initialization of variable or porosity fields
   **TKEIN=0.25*WIN*WIN*FRIC  WHERE FRIC=0.018 AT REYNO=1.E5
TKEIN=WIN*WIN*0.018*0.25;FIINIT(KE)=TKEIN
 
    GROUP 13. Boundary conditions and special sources
   **NORTH WALL AT IY=NY
COVAL(WFUN,KE,LOGLAW,LOGLAW)
   **Inlet Boundary
VALUE(IN,KE,TKEIN)
 
    GROUP 16. Termination of iterations
RESREF(KE)=1.E-5;RESREF(TEMP)=1.E-5
 
    GROUP 18. Limits on variables or increments to them
VARMIN(W1)=-1.E+10
VARMIN(KE)=1.E-10
 
    GROUP 19. Data communicated by SATELLITE to GROUND
  ** Select strain-rate for turbulence production term
GENK=F
 
    GROUP 21. Print-out of variables
OUTPUT(KE,Y,Y,Y,N,Y,Y)
 
    GROUP 23. Field print-out and plot control
PLOT(IZEQNZ,KE,0,20.0)