L(240)
    GROUP 1. Run title and other preliminaries
TEXT(2D Restricted Inlet/Outlet Channel
TITLE
  DISPLAY
  This run analyses the flow in a duct with restrictions at both
  ends. The walls are held at constant but unequal temperatures.
  Interesting parameters to vary include:- entry values of scalar
  quantities; the inlet and/or outlet area; and the wall
  temperatures.
  ENDDIS
  *** Introduce more Regions to allow correct Patch specification
    GROUP 4. Y-direction grid specification
NREGY=3
IREGY=1; GRDPWR(Y,2,0.4,1.0)
IREGY=2; GRDPWR(Y,1,0.2,1.0)
IREGY=3; GRDPWR(Y,2,0.4,1.0)
    GROUP 13. Boundary conditions and special sources
  ** Inlet
INLET(IN,LOW,#1,#1,#1,#1,#1,#1,1,1)
VALUE(IN,P1,5.0); VALUE(IN,W1,5.0)
VALUE(IN,V1,0.0); VALUE(IN,TEMP,9.0)
  ** Outlet
PATCH(OUTLET,HIGH,#1,#1,#1,#2,#NREGZ,#NREGZ,1,1)
COVAL(OUTLET,P1,FIXP,0.0)
COVAL(OUTLET,V1,ONLYMS,0.0); COVAL(OUTLET,W1,ONLYMS,0.0)
  ** North wall
WALL (NORTH,NORTH,#1,#1,#NREGY,#NREGY,1,NZ,1,1)
COVAL(NORTH,W1,1.0,0.0); COVAL(NORTH,TEMP,1.0,-1.0)
  ** High wall
WALL (HIGHWALL,HIGH,#1,#1,#NREGY,#NREGY,#NREGZ,#NREGZ,1,1)
COVAL(HIGHWALL,V1,1.0,0.0); COVAL(HIGHWALL,TEMP,1.0,-1.0)
  ** Inlet wall
WALL (LOWWALL,LOW,#1,#1,#2,#NREGY,#1,#1,1,1)
COVAL(LOWWALL,V1,1.0,0.0); COVAL(LOWWALL,TEMP,1.0,-20.0)
 
    GROUP 17. Under-relaxation devices
RELAX(V1,FALSDT,ZWLAST/5.0); RELAX(W1,FALSDT,ZWLAST/5.0)
 
    GROUP 22. Spot-value print-out
IYMON=1; IZMON=4
 
    GROUP 23. Field print-out and plot control
NPRINT=20; NPLT=2
PATCH(ZEQ3,PROFIL,1,1,1,NY,3,3,1,1)
PLOT(ZEQ3,W1,0.0,.0); PLOT(ZEQ3,TEMP,0.0,0.0)
PATCH(YEQ1,PROFIL,1,1,1,1,1,NZ,1,1)
PLOT(YEQ1,W1,0.0,0.0); PLOT(YEQ1,P1,0.0,0.0)
PLOT(YEQ1,TEMP,0.0,0.0)
PATCH(MAP,CONTUR,1,1,1,NY,1,NZ,1,1)
PLOT(MAP,W1,0.0,10.0); PLOT(MAP,P1,0.0,10.0)
PLOT(MAP,TEMP,0.0,10.0)