TALK=T;RUN(1,1)
  DISPLAY

  Simulation of laminar steady flow in z direction around 2 spheres
  in a uniform stream of atmospheric air, treated as incompressible.

  The basic grid is cartesian and rather coarse (18*18*31), but the
  flow follows the curved surface smoothly because the PARSOL (i.e.
  cut-cell) technique  is used.

  Provision is made for:
  1. Solving for only one quarter of the domain; this is allowed
  by reason of symmetry and desirable for economy and accuracy.
  2. Refining the grid near the spheres to increase accuracy more.
  3. Varying the Reynolds number.

  The variables computed are:
  * p1, the pressure,
  * u1, the x-direction velocity,
  * v1, the y-direction velocity, and
  * w1, the z-direction velocity.

  This file contains a macro of commands which cause the Viewer
  (when the macro button is pressed) to display contours
  automatically.

  [Case derived from core linrary case 805 by dbs 11.09.08]
  ENDDIS
   VRV USE
  * Start of frame
  VARIABLE  Pressure; CON ON
  msg Pressure contours
  PAUSE
  * Start of frame
  VARIABLE  Velocity; VEC ON
  msg Pressure contours superimposed by velocity vectors
  PAUSE
  * Start of frame
  * Setting object visibility and painting status
  OBJECT SHOW TYPE BLOCKAGE
  OBJECT PAINT TYPE BLOCKAGE ON
  VARIABLE  Velocity; VEC OFF
  CONTOUR OFF
  msg Velocity contours on the surface of the sphere
  PAUSE
  * Start of frame
  * Setting object visibility and painting status
  OBJECT SHOW TYPE BLOCKAGE
  OBJECT PAINT TYPE BLOCKAGE ON
  VARIABLE  Pressure; VEC OFF
  msg Pressure contours on the surface of the sphere
  PAUSE
  ENDUSE
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(2 spheres Re=40. quarter=T finegrid=T   )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 0
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 18
 XULAST =2.
 XFRAC(1)=0.0625 ;XFRAC(2)=0.125
 XFRAC(3)=0.1875 ;XFRAC(4)=0.25
 XFRAC(5)=0.3 ;XFRAC(6)=0.35
 XFRAC(7)=0.4 ;XFRAC(8)=0.45
 XFRAC(9)=0.5 ;XFRAC(10)=0.55
 XFRAC(11)=0.6 ;XFRAC(12)=0.65
 XFRAC(13)=0.7 ;XFRAC(14)=0.76
 XFRAC(15)=0.82 ;XFRAC(16)=0.88
 XFRAC(17)=0.94 ;XFRAC(18)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 18
 YVLAST =2.
 YFRAC(1)=0.0625 ;YFRAC(2)=0.125
 YFRAC(3)=0.1875 ;YFRAC(4)=0.25
 YFRAC(5)=0.3 ;YFRAC(6)=0.35
 YFRAC(7)=0.4 ;YFRAC(8)=0.45
 YFRAC(9)=0.5 ;YFRAC(10)=0.55
 YFRAC(11)=0.6 ;YFRAC(12)=0.65
 YFRAC(13)=0.7 ;YFRAC(14)=0.76
 YFRAC(15)=0.82 ;YFRAC(16)=0.88
 YFRAC(17)=0.94 ;YFRAC(18)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 31
 ZWLAST =10.
 ZFRAC(1)=0.0364 ;ZFRAC(2)=0.0728
 ZFRAC(3)=0.1092 ;ZFRAC(4)=0.1456
 ZFRAC(5)=0.182 ;ZFRAC(6)=0.202
 ZFRAC(7)=0.222 ;ZFRAC(8)=0.255333
 ZFRAC(9)=0.288667 ;ZFRAC(10)=0.322
 ZFRAC(11)=0.355333 ;ZFRAC(12)=0.388667
 ZFRAC(13)=0.422 ;ZFRAC(14)=0.442
 ZFRAC(15)=0.462 ;ZFRAC(16)=0.4988
 ZFRAC(17)=0.5356 ;ZFRAC(18)=0.5724
 ZFRAC(19)=0.6092 ;ZFRAC(20)=0.646
 ZFRAC(21)=0.666 ;ZFRAC(22)=0.699333
 ZFRAC(23)=0.732667 ;ZFRAC(24)=0.766
 ZFRAC(25)=0.786 ;ZFRAC(26)=0.821667
 ZFRAC(27)=0.857333 ;ZFRAC(28)=0.893
 ZFRAC(29)=0.928667 ;ZFRAC(30)=0.964333
 ZFRAC(31)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(1)=P1 ;NAME(3)=U1
 NAME(5)=V1 ;NAME(7)=W1
 NAME(148)=STGP ;NAME(149)=VLSQ
 NAME(150)=PRPS
    * 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,Y,N,N,Y)
 SOLUTN(V1,Y,Y,Y,N,N,Y)
 SOLUTN(W1,Y,Y,Y,N,N,Y)
 SOLUTN(STGP,Y,N,N,N,N,N)
 SOLUTN(VLSQ,Y,N,N,N,N,Y)
 SOLUTN(PRPS,Y,N,N,N,N,Y)
 PRPS = 150
 ************************************************************
  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)
 TERMS(W1,Y,Y,Y,Y,Y,Y)
 DIFCUT =0.5 ;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. ;TMP1 =0.
 EL1 =0.
 TSURR =0. ;TEMP0 =0.
 PRESS0 =0.
 DVO1DT =0. ;DRH1DP =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 ENUL =0.05 ;ENUT =0.
 PRNDTL(U1)=1. ;PRNDTL(V1)=1.
 PRNDTL(W1)=1.
 PRT(U1)=1. ;PRT(V1)=1.
 PRT(W1)=1.
 CP1 =1005. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1.0E-10 ;FIINIT(U1)=1.0E-10
 FIINIT(V1)=1.0E-10 ;FIINIT(W1)=1.
 FIINIT(STGP)=1.0E-10 ;FIINIT(VLSQ)=1.0E-10
 FIINIT(PRPS)=-1.
 
 PATCH(^OB5 ,INIVAL, 5, 0, 0, 0, 0, 0, 1, 1)
 INIT(^OB5 ,PRPS,0. ,198. )
 
 PATCH(^OB6 ,INIVAL, 6, 0, 0, 0, 0, 0, 1, 1)
 INIT(^OB6 ,PRPS,0. ,198. )
 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(IMBL3&2 ,VOLUME, 3, 16, 3, 16, 3, 29, 1, 1)
 COVAL(IMBL3&2 ,W1 ,0. ,0. )
 
 PATCH(IMBL4&3 ,VOLUME, 3, 16, 3, 16, 4, 28, 1, 1)
 COVAL(IMBL4&3 ,W1 ,0. ,0. )
 
 PATCH(IMBL5&4 ,VOLUME, 3, 16, 3, 16, 5, 27, 1, 1)
 COVAL(IMBL5&4 ,W1 ,0. ,0. )
 
 PATCH(IMBL810 ,VOLUME, 3, 16, 3, 16, 8, 21, 1, 1)
 COVAL(IMBL810 ,W1 ,0. ,0. )
 
 PATCH(IMBL815 ,VOLUME, 5, 14, 5, 14, 8, 16, 1, 1)
 COVAL(IMBL815 ,W1 ,0. ,0. )
 
 PATCH(^OB1 ,CELL , 1, 0, 0, 0, 0, 0, 1, 1)
 
 PATCH(^OB2 ,CELL , 2, 0, 0, 0, 0, 0, 1, 1)
 
 PATCH(^OB3 ,LOW , 3, 0, 0, 0, 0, 0, 1, 1)
 COVAL(^OB3 ,P1 , FIXFLU ,1. )
 COVAL(^OB3 ,U1 ,0. ,0. )
 COVAL(^OB3 ,V1 ,0. ,0. )
 COVAL(^OB3 ,W1 ,0. ,1. )
 
 PATCH(^OB4 ,HIGH , 4, 0, 0, 0, 0, 0, 1, 1)
 COVAL(^OB4 ,P1 ,1000. ,0. )
 COVAL(^OB4 ,U1 ,0. ,0. )
 COVAL(^OB4 ,V1 ,0. ,0. )
 COVAL(^OB4 ,W1 ,0. ,0. )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 30 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-03
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=250 ;LITER(U1)=10
 LITER(V1)=10 ;LITER(W1)=10
 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
 ENDIT(V1)=1.0E-03 ;ENDIT(W1)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(U1,FALSDT,1.)
 RELAX(V1,FALSDT,1.)
 RELAX(W1,FALSDT,1.)
 RELAX(STGP,LINRLX,1.)
 RELAX(VLSQ,LINRLX,1.)
 RELAX(PRPS,LINRLX,1.)
 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(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
 VARMAX(STGP)=1.0E+10 ;VARMIN(STGP)=-1.0E+10
 VARMAX(VLSQ)=1.0E+10 ;VARMIN(VLSQ)=-1.0E+10
 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 ISG50 = 1
 ISG52 = 2
 ISG62 = 0
 SPEDAT(SET,STORED,STGP,C,=P1+0.5*RHO1*VLSQ!SWPFIN)
 SPEDAT(SET,OBJNAM,^OB1,C,G1)
 SPEDAT(SET,OBJTYP,^OB1,C,FINE_GRID_VOL)
 SPEDAT(SET,^OB1,FINEX,I,4)
 SPEDAT(SET,^OB1,FINEY,I,4)
 SPEDAT(SET,^OB1,FINEZ,I,4)
 SPEDAT(SET,OBJNAM,^OB2,C,G2)
 SPEDAT(SET,OBJTYP,^OB2,C,FINE_GRID_VOL)
 SPEDAT(SET,^OB2,FINEX,I,4)
 SPEDAT(SET,^OB2,FINEY,I,4)
 SPEDAT(SET,^OB2,FINEZ,I,4)
 SPEDAT(SET,OBJNAM,^OB3,C,AIRIN)
 SPEDAT(SET,OBJTYP,^OB3,C,INLET)
 SPEDAT(SET,OBJNAM,^OB4,C,AIROUT)
 SPEDAT(SET,OBJTYP,^OB4,C,OUTLET)
 SPEDAT(SET,ARATIO,^OB4,R,1.)
 SPEDAT(SET,OBJNAM,^OB5,C,FIRSTSPH)
 SPEDAT(SET,OBJTYP,^OB5,C,BLOCKAGE)
 SPEDAT(SET,FIRSTSPH,MATERIAL,R,198.)
 SPEDAT(SET,OBJNAM,^OB6,C,SECONDSP)
 SPEDAT(SET,OBJTYP,^OB6,C,BLOCKAGE)
 SPEDAT(SET,SECONDSP,MATERIAL,R,198.)
 SPEDAT(SET,FACETDAT,NUMOBJ,I,6)
 SPEDAT(SET,MATERIAL,198,L,T)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=0.1298 ;EX(U1)=0.0324
 EX(V1)=0.03335 ;EX(W1)=1.012
 EX(STGP)=0. ;EX(VLSQ)=0.
 EX(PRPS)=3.354
 ************************************************************
  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,Y,Y,N,Y,Y)
 OUTPUT(V1,Y,N,Y,Y,Y,Y)
 OUTPUT(W1,Y,N,Y,Y,Y,Y)
 OUTPUT(STGP,Y,N,Y,N,N,N)
 OUTPUT(VLSQ,Y,N,Y,N,N,N)
 OUTPUT(PRPS,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 9 ;IYMON = 9 ;IZMON = 17
 NPRMON = 100000 ;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
 NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
 XZPR = F ;YZPR = F
 IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
 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