```
GROUP 1. Run title and other preliminaries
TEXT(Viscous Heating In Couette Flow
TITLE
DISPLAY
|             h   h
|         h    u  | moving wall
fixed  |      h    u     |
wall   |   h    u        |
| h   u           |
|  u              |
u -------> y      |ywall

The flow is laminar, with uniform properties. The non-linearity
of the temperature profile is due to viscous dissipation.

The solution should be:  u/uwall = y/ywall;
(h-hzero)/(hwall-hzero) = y/ywall + H*(y/ywall)*(1 - y/ywall)
where:  H = 0.5 * Prandtl No * uwall**2 /(hwall -hzero)

In this case, H = 1, which implies zero heat flux at the moving
wall. The numerical calculations conform completely with the
analytical solution.

ENDDIS
REAL(HCONST)
HCONST=1.0
GROUP 4. Y-direction grid specification
NY=40;GRDPWR(Y,NY,YVLAST,1.0)
GROUP 7. Variables stored, solved & named
SOLVE(U1,H1)

GROUP 8. Terms (in differential equations) & devices
TERMS(H1,Y,N,Y,N,Y,N);TERMS(U1,Y,N,Y,N,Y,N)

GROUP 9. Properties of the medium (or media)
ENUL=1.0;PRNDTL(H1)=2.0*HCONST

GROUP 13. Boundary conditions and special sources
WALL (INNER,SOUTH,1,1,1,1,1,1,1,1)
COVAL(INNER,U1,1.,0.0);COVAL(INNER,H1,1.0,0.0)
WALL (OUTER,NORTH,1,1,NY,NY,1,1,1,1)
COVAL(OUTER,U1,1.0,1.0);COVAL(OUTER,H1,1.0,1.0)

GROUP 15. Termination of sweeps
LSWEEP=3

GROUP 16. Termination of iterations
RESREF(U1)=1.E-10;RESREF(H1)=1.E-10
GROUP 22. Spot-value print-out
;IYMON=20
GROUP 23. Field print-out and plot control
PATCH(PROFUH,PROFIL,1,1,1,NY,1,1,1,1)
PLOT(PROFUH,U1,0.0,0.0);PLOT(PROFUH,H1,0.0,0.0)
ORSIZ=0.5;ITABL=2
**END OF LIBRARY CASE 228
```