TALK=T;RUN(1,1)
  DISPLAY
  The case considered is 2d laminar incompressible flow over a
  thin fence of height H located in a planar channel with a
  blockage ratio S/H=0.75. The flow Reynolds number is 82.5 based
  on fence height S and inlet bulk velocity Uin. This situation
  has been studied experimentally and numerically by Carvalho et
  al (Appl.Math.Modelling, Vol.11, p23, (1987) ). The boundary
  conditions correspond to an inlet flow of fully-developed
  parabolic velocity profile located at 8 fence heights upstream
  of the fence, and an outlet condition of fixed pressure at 15
  fence heights downstream of the fence, and no-slip conditions
  at the walls. In practice the fence thickness is t/H=0.133,
  but in the calculations it is taken as zero.

  The present calculation is performed in all six planes as a
  consistency test. The default calculation is made with the
  linear-upwind scheme (LUS) for momentum, and the option exists
  to make the calculation with the HYBRID scheme. A mesh of 24
  vertical by 48 horizontal cells has been used in the
  calculations.
  ENDDIS
  PHOTON USE
  P


   0.20443E+04 0.15633E+04 CR
  GR OU Z 1;GR Z 1 COL 3;GR X 13 13 Y 1 12 COL 15
  STREAM 2D Z 1 X 4 35 Y 1 M
  -.199E-2 .642E-3 10
  set order
  1
  4
  5
  2
  6
  7
  3
  MAG GR 2
   0.16213E+04 0.17326E+04 CR
  msg Streamlines
  msg press 
  pause
  stream off;grid off 2;redr
  msg vectors
  vec z 1 sh
  msg press  and then  to end
  pause
  ENDUSE
  AUTOPLOT USE
  file
  phi 5

  d 1 u1 y 1;d 1 u1 y m;div x 7.5e-3 1 2;shift x -8 1 2
  col3 1;colf 2;level y 0;scale x 0 10;scale y -.3 .25;
  msg horizontal velocity along bottom & top walls behind fence
  msg separation & reattachment points when U1 passes through zero
  msg blue line - U1 along bottom wall
  msg red  line - U1 along top wall
  msg press  to continue
  msg press  to end
  ENDUSE
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(2DYX LAMINAR FLOW OVER A FENCE: N108    )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 14
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 48
 XULAST =0.1725
 XFRAC(1)=0.034486 ;XFRAC(2)=0.06835
 XFRAC(3)=0.101542 ;XFRAC(4)=0.134004
 XFRAC(5)=0.165662 ;XFRAC(6)=0.196426
 XFRAC(7)=0.226177 ;XFRAC(8)=0.254755
 XFRAC(9)=0.281925 ;XFRAC(10)=0.307314
 XFRAC(11)=0.330192 ;XFRAC(12)=0.347826
 XFRAC(13)=0.354009 ;XFRAC(14)=0.363049
 XFRAC(15)=0.373615 ;XFRAC(16)=0.38531
 XFRAC(17)=0.397925 ;XFRAC(18)=0.411325
 XFRAC(19)=0.425415 ;XFRAC(20)=0.440123
 XFRAC(21)=0.455395 ;XFRAC(22)=0.471185
 XFRAC(23)=0.487457 ;XFRAC(24)=0.504179
 XFRAC(25)=0.521325 ;XFRAC(26)=0.538872
 XFRAC(27)=0.556798 ;XFRAC(28)=0.575088
 XFRAC(29)=0.593723 ;XFRAC(30)=0.612691
 XFRAC(31)=0.631977 ;XFRAC(32)=0.651571
 XFRAC(33)=0.671461 ;XFRAC(34)=0.691637
 XFRAC(35)=0.71209 ;XFRAC(36)=0.732812
 XFRAC(37)=0.753794 ;XFRAC(38)=0.77503
 XFRAC(39)=0.796513 ;XFRAC(40)=0.818235
 XFRAC(41)=0.840192 ;XFRAC(42)=0.862376
 XFRAC(43)=0.884784 ;XFRAC(44)=0.90741
 XFRAC(45)=0.930249 ;XFRAC(46)=0.953296
 XFRAC(47)=0.976548 ;XFRAC(48)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 24
 YVLAST =1.0E-02
 YFRAC(1)=0.074361 ;YFRAC(2)=0.14738
 YFRAC(3)=0.218951 ;YFRAC(4)=0.288946
 YFRAC(5)=0.357209 ;YFRAC(6)=0.423544
 YFRAC(7)=0.487694 ;YFRAC(8)=0.549315
 YFRAC(9)=0.607902 ;YFRAC(10)=0.662647
 YFRAC(11)=0.711977 ;YFRAC(12)=0.75
 YFRAC(13)=0.776738 ;YFRAC(14)=0.802756
 YFRAC(15)=0.828004 ;YFRAC(16)=0.852422
 YFRAC(17)=0.87594 ;YFRAC(18)=0.898468
 YFRAC(19)=0.919894 ;YFRAC(20)=0.940065
 YFRAC(21)=0.958765 ;YFRAC(22)=0.975659
 YFRAC(23)=0.990114 ;YFRAC(24)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 1
 ZWLAST =1.
 ZFRAC(1)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(1)=P1 ;NAME(3)=U1
 NAME(5)=V1 ;NAME(150)=EPOR
    * Y in SOLUTN argument list denotes:
    * 1-stored 2-solved 3-whole-field
    * 4-point-by-point 5-explicit 6-harmonic averaging 
 SOLUTN(P1,Y,Y,Y,N,N,Y)
 SOLUTN(U1,Y,Y,N,N,N,Y)
 SOLUTN(V1,Y,Y,N,N,N,Y)
 SOLUTN(EPOR,Y,N,N,N,N,N)
 EPOR = 150 ;HPOR = 0 ;NPOR = 0 ;VPOR = 0
 ************************************************************
  Group 8. Terms & Devices
    * Y in TERMS argument list denotes:
    * 1-built-in source 2-convection 3-diffusion 4-transient
    * 5-first phase variable 6-interphase transport         
 TERMS(P1,Y,Y,Y,N,Y,Y)
 TERMS(U1,Y,Y,Y,Y,Y,Y)
 TERMS(V1,Y,Y,Y,Y,Y,Y)
 DIFCUT =0. ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
 ************************************************************
  Group 9. Properties used if PRPS is not
  stored, and where PRPS = -1.0 if it is!
 RHO1 =1.2 ;TMP1 =0.
 EL1 =0.
 TSURR =0. ;TEMP0 =0.
 PRESS0 =0.
 DVO1DT =0. ;DRH1DP =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 ENUL =1.8E-05 ;ENUT =0.
 PRNDTL(U1)=1. ;PRNDTL(V1)=1.
 PRT(U1)=1. ;PRT(V1)=1.
 CP1 =1. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1.0E-10 ;FIINIT(U1)=0.198
 FIINIT(V1)=0.0198 ;FIINIT(EPOR)=1.
 
 PATCH(FENCE ,INIVAL, 12, 12, 1, 12, 1, 1, 1, 1)
 INIT(FENCE ,U1 ,0. ,0. )
 INIT(FENCE ,EPOR,0. ,0. )
 INIADD = F
 FSWEEP = 1
 NAMFI =CHAM
 ************************************************************
  Group 12. Patchwise adjustment of terms
  Patches for this group are printed with those
  for Group 13.
  Their names begin either with GP12 or &
 ************************************************************
  Group 13. Boundary & Special Sources
 
 PATCH(HOCS ,CELL , 0, 0, 0, 0, 0, 0, 1, 1)
 COVAL(HOCS ,U1 , FIXFLU , GRND1 )
 COVAL(HOCS ,V1 , FIXFLU , GRND1 )
 
 PATCH(OUTLET ,EAST , 48, 48, 1, 24, 1, 1, 1, 1)
 COVAL(OUTLET ,P1 ,1000. ,0. )
 COVAL(OUTLET ,U1 ,0. ,0. )
 COVAL(OUTLET ,V1 ,0. ,0. )
 
 PATCH(TOP ,NWALL , 1, 48, 24, 24, 1, 1, 1, 1)
 COVAL(TOP ,U1 ,1. ,0. )
 
 PATCH(BOT ,SWALL , 1, 48, 1, 1, 1, 1, 1, 1)
 COVAL(BOT ,U1 ,1. ,0. )
 
 PATCH(IN1 ,WEST , 1, 1, 1, 1, 1, 1, 1, 1)
 COVAL(IN1 ,P1 , FIXFLU ,0.051033 )
 COVAL(IN1 ,U1 ,0. ,0.042528 )
 COVAL(IN1 ,V1 ,0. ,0. )
 
 PATCH(IN2 ,WEST , 1, 1, 2, 2, 1, 1, 1, 1)
 COVAL(IN2 ,P1 , FIXFLU ,0.140533 )
 COVAL(IN2 ,U1 ,0. ,0.117111 )
 COVAL(IN2 ,V1 ,0. ,0. )
 
 PATCH(IN3 ,WEST , 1, 1, 3, 3, 1, 1, 1, 1)
 COVAL(IN3 ,P1 , FIXFLU ,0.213292 )
 COVAL(IN3 ,U1 ,0. ,0.177743 )
 COVAL(IN3 ,V1 ,0. ,0. )
 
 PATCH(IN4 ,WEST , 1, 1, 4, 4, 1, 1, 1, 1)
 COVAL(IN4 ,P1 , FIXFLU ,0.270092 )
 COVAL(IN4 ,U1 ,0. ,0.225077 )
 COVAL(IN4 ,V1 ,0. ,0. )
 
 PATCH(IN5 ,WEST , 1, 1, 5, 5, 1, 1, 1, 1)
 COVAL(IN5 ,P1 , FIXFLU ,0.311777 )
 COVAL(IN5 ,U1 ,0. ,0.259814 )
 COVAL(IN5 ,V1 ,0. ,0. )
 
 PATCH(IN6 ,WEST , 1, 1, 6, 6, 1, 1, 1, 1)
 COVAL(IN6 ,P1 , FIXFLU ,0.339268 )
 COVAL(IN6 ,U1 ,0. ,0.282723 )
 COVAL(IN6 ,V1 ,0. ,0. )
 
 PATCH(IN7 ,WEST , 1, 1, 7, 7, 1, 1, 1, 1)
 COVAL(IN7 ,P1 , FIXFLU ,0.353592 )
 COVAL(IN7 ,U1 ,0. ,0.29466 )
 COVAL(IN7 ,V1 ,0. ,0. )
 
 PATCH(IN8 ,WEST , 1, 1, 8, 8, 1, 1, 1, 1)
 COVAL(IN8 ,P1 , FIXFLU ,0.355912 )
 COVAL(IN8 ,U1 ,0. ,0.296593 )
 COVAL(IN8 ,V1 ,0. ,0. )
 
 PATCH(IN9 ,WEST , 1, 1, 9, 9, 1, 1, 1, 1)
 COVAL(IN9 ,P1 , FIXFLU ,0.347591 )
 COVAL(IN9 ,U1 ,0. ,0.289659 )
 COVAL(IN9 ,V1 ,0. ,0. )
 
 PATCH(IN10 ,WEST , 1, 1, 10, 10, 1, 1, 1, 1)
 COVAL(IN10 ,P1 , FIXFLU ,0.330313 )
 COVAL(IN10 ,U1 ,0. ,0.275261 )
 COVAL(IN10 ,V1 ,0. ,0. )
 
 PATCH(IN11 ,WEST , 1, 1, 11, 11, 1, 1, 1, 1)
 COVAL(IN11 ,P1 , FIXFLU ,0.306382 )
 COVAL(IN11 ,U1 ,0. ,0.255318 )
 COVAL(IN11 ,V1 ,0. ,0. )
 
 PATCH(IN12 ,WEST , 1, 1, 12, 12, 1, 1, 1, 1)
 COVAL(IN12 ,P1 , FIXFLU ,0.280336 )
 COVAL(IN12 ,U1 ,0. ,0.233613 )
 COVAL(IN12 ,V1 ,0. ,0. )
 
 PATCH(IN13 ,WEST , 1, 1, 13, 13, 1, 1, 1, 1)
 COVAL(IN13 ,P1 , FIXFLU ,0.257516 )
 COVAL(IN13 ,U1 ,0. ,0.214596 )
 COVAL(IN13 ,V1 ,0. ,0. )
 
 PATCH(IN14 ,WEST , 1, 1, 14, 14, 1, 1, 1, 1)
 COVAL(IN14 ,P1 , FIXFLU ,0.236716 )
 COVAL(IN14 ,U1 ,0. ,0.197264 )
 COVAL(IN14 ,V1 ,0. ,0. )
 
 PATCH(IN15 ,WEST , 1, 1, 15, 15, 1, 1, 1, 1)
 COVAL(IN15 ,P1 , FIXFLU ,0.214604 )
 COVAL(IN15 ,U1 ,0. ,0.178836 )
 COVAL(IN15 ,V1 ,0. ,0. )
 
 PATCH(IN16 ,WEST , 1, 1, 16, 16, 1, 1, 1, 1)
 COVAL(IN16 ,P1 , FIXFLU ,0.191394 )
 COVAL(IN16 ,U1 ,0. ,0.159495 )
 COVAL(IN16 ,V1 ,0. ,0. )
 
 PATCH(IN17 ,WEST , 1, 1, 17, 17, 1, 1, 1, 1)
 COVAL(IN17 ,P1 , FIXFLU ,0.167326 )
 COVAL(IN17 ,U1 ,0. ,0.139438 )
 COVAL(IN17 ,V1 ,0. ,0. )
 
 PATCH(IN18 ,WEST , 1, 1, 18, 18, 1, 1, 1, 1)
 COVAL(IN18 ,P1 , FIXFLU ,0.142664 )
 COVAL(IN18 ,U1 ,0. ,0.118887 )
 COVAL(IN18 ,V1 ,0. ,0. )
 
 PATCH(IN19 ,WEST , 1, 1, 19, 19, 1, 1, 1, 1)
 COVAL(IN19 ,P1 , FIXFLU ,0.117713 )
 COVAL(IN19 ,U1 ,0. ,0.098094 )
 COVAL(IN19 ,V1 ,0. ,0. )
 
 PATCH(IN20 ,WEST , 1, 1, 20, 20, 1, 1, 1, 1)
 COVAL(IN20 ,P1 , FIXFLU ,0.092832 )
 COVAL(IN20 ,U1 ,0. ,0.07736 )
 COVAL(IN20 ,V1 ,0. ,0. )
 
 PATCH(IN21 ,WEST , 1, 1, 21, 21, 1, 1, 1, 1)
 COVAL(IN21 ,P1 , FIXFLU ,0.068466 )
 COVAL(IN21 ,U1 ,0. ,0.057055 )
 COVAL(IN21 ,V1 ,0. ,0. )
 
 PATCH(IN22 ,WEST , 1, 1, 22, 22, 1, 1, 1, 1)
 COVAL(IN22 ,P1 , FIXFLU ,0.04521 )
 COVAL(IN22 ,U1 ,0. ,0.037675 )
 COVAL(IN22 ,V1 ,0. ,0. )
 
 PATCH(IN23 ,WEST , 1, 1, 23, 23, 1, 1, 1, 1)
 COVAL(IN23 ,P1 , FIXFLU ,0.02398 )
 COVAL(IN23 ,U1 ,0. ,0.019983 )
 COVAL(IN23 ,V1 ,0. ,0. )
 
 PATCH(IN24 ,WEST , 1, 1, 24, 24, 1, 1, 1, 1)
 COVAL(IN24 ,P1 , FIXFLU ,7.011771E-03 )
 COVAL(IN24 ,U1 ,0. ,5.843142E-03 )
 COVAL(IN24 ,V1 ,0. ,0. )
 
 PATCH(FENCE-EW,EWALL , 12, 12, 1, 12, 1, 1, 1, 1)
 COVAL(FENCE-EW,V1 ,1. ,0. )
 
 PATCH(FENCE-WW,WWALL , 13, 13, 1, 12, 1, 1, 1, 1)
 COVAL(FENCE-WW,V1 ,1. ,0. )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 1600 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-04
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(U1)=10
 LITER(V1)=10
 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
 ENDIT(V1)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(U1,FALSDT,1.089015E-03)
 RELAX(V1,FALSDT,1.089015E-03)
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
 VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
 VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
 VARMAX(EPOR)=1.0E+10 ;VARMIN(EPOR)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 PARSOL = F
 ISG62 = 1
 SPEDAT(SET,SCHEME,INLCS003,I,1)
 SPEDAT(SET,SCHEME,INLCS005,I,1)
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  Group 21. Print-out of Variables
 INIFLD = F ;SUBWGR = F
    * Y in OUTPUT argument list denotes:
    * 1-field 2-correction-eq. monitor 3-selective dumping      
    * 4-whole-field residual 5-spot-value table 6-residual table
 OUTPUT(P1,Y,N,Y,Y,Y,Y)
 OUTPUT(U1,Y,N,Y,Y,Y,Y)
 OUTPUT(V1,Y,N,Y,Y,Y,Y)
 OUTPUT(EPOR,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 24 ;IYMON = 10 ;IZMON = 1
 NPRMON = 100 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
 NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
 IPLTF = 1 ;IPLTL = 1600 ;NPLT = 20
 ISWPRF = 1 ;ISWPRL = 100000
 ITABL = 3 ;IPROF = 1
 ABSIZ =0.5 ;ORSIZ =0.4
 NTZPRF = 1 ;NCOLPF = 50
 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
   No PATCHes yet used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP