L(322)
    GROUP 1. Run title and other preliminaries
TEXT(Central Obstacle In Laminar Pipe  
TITLE
  DISPLAY
  This run performs the same analysis as the preceeding case
  except that an area restriction is added at half length.
 
                     hot pipe wall
               /////////////////////////////////
               ---------------------------------
           --->          - - - - -
  uniform flow          -    ||     -  - - - - ->  exit
  at inlet --->   - -  -     ||obstacle
             __________ _ ___||____________axis____
 
  ENDDIS
 
    4. Obstacle--- at the half length of the pipe
 
   *  Write objects
RSET(B,OBSTC,0,0,0.25,1.,0.006,0,12,0)
 
    GROUP 4. Y-direction grid specification
NREGY=2
IREGY=1; GRDPWR(Y,6,0.006,1.0)
IREGY=2; GRDPWR(Y,4,0.004,1.0)
 
    GROUP 5. Z-direction grid specification
NREGZ=2
IREGZ=1; GRDPWR(Z,5,0.25,1.0)
IREGZ=2; GRDPWR(Z,5,0.25,1.0)
 
PATCH(OBSTC,HIGH,#1,#1,#1,#1,#1,#1,#1,#NREGT)
COVAL(OBSTC,W1,FIXVAL,0.0)
 
    GROUP 23. Field print-out and plot control
NYPRIN=1;NZPRIN=1;NTPRIN=1
PATCH(TIMEPLOT,CONTUR,1,1,1,NY,1,NZ,1,LSTEP)
PLOT (TIMEPLOT,W1,0.0,10);PLOT(TIMEPLOT,H1,0.0,10)
PLOT (TIMEPLOT,A,0.0,10) ;PLOT(TIMEPLOT,B,0.0,10.0)
PLOT (TIMEPLOT,C,0.0,10)