TALK=T;RUN(1,1)
  DISPLAY

  This In-Form case uses "IMAT" In-Form option for simulation
  flow in stirred 3D vessel.

  The following statements make the vertical and horizontal
  velocities proportional to the corresponding distances from
  the paddle axis.

  (SOURCE of U1 at EI is 1.E10*(VEL*(YIC-YG)-U1) with IMAT>=100!LIN$
E)
  (SOURCE of V1 at NI is 1.E10*(VEL*(XG-XIC)-V1) with IMAT>=100!LIN$
E)

  These momentum sources are active only for cell having IMAT
  (the property index, as set by PRPS) values greater or equal
  to 100 .

  The following In-Form formulae describe density and viscosity:

  (PROPERTY of RHO1 is DENS)
  (PROPERTY VISL is VISC)

  The next statement sets PRPS to the IMAT value of water
  throughout the domain; but the value is changed below for
  the cells ocupied by the paddle

  (STORED of PRPS is 67.)

  Cells occupied by paddle material are marked by setting the
  value of PRPS in them to 100.

  (STORED of PRPS at TIME=1 is :MAT:)
  (STORED of PRPS at TIME=2 is :MAT:)
  (STORED of PRPS at TIME=3 is :MAT:)
  (STORED of PRPS at TIME=4 is :MAT:)
  (STORED of PRPS at TIME=5 is :MAT:)
  (STORED of PRPS at TIME=6 is :MAT:)
  (STORED of PRPS at TIME=7 is :MAT:)
  (STORED of PRPS at TIME=8 is :MAT:)

  The Q1 contains PHOTON USE commands
  ENDDIS
  PHOTON USE
  p
  p1;

  set prop off ; vi 1 1 5

  surf prps z 89
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec z 8 sh; vec y 1 sh; vec x 1 sh
  msg at 1 sec.
  * con v1 y 5 fi;0.001;*dump picb1
  pause
  p
  p2;

  set prop off ; vi 1 1 5
  surf prps z 89
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec z 8 sh; vec y 1 sh; vec x 1 sh
  msg at 2 sec.
  * con v1 y 5 fi;0.001;*dump picb2
  pause
  p
  p3;

  set prop off ; vi 1 1 5
  surf prps z 89
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec z 8 sh; vec y 1 sh; vec x 1 sh
  msg at 3 sec.
  * con v1 y 5 fi;0.001;*dump picb3
  pause
  p
  p4;

  set prop off ; vi 1 1 5
  surf prps z 89
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec z 8 sh; vec y 1 sh; vec x 1 sh
  msg at 4 sec.
  * con v1 y 5 fi;0.001;*dump picb4
  pause
  p
  p5;

  set prop off ; vi 1 1 5
  surf prps z 89
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec z 8 sh; vec y 1 sh; vec x 1 sh
  msg at 5 sec.
  * con v1 y 5 fi;0.001;*dump picb5
  pause
  p
  p6;

  set prop off ; vi 1 1 5
  surf prps z 89
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec z 8 sh; vec y 1 sh; vec x 1 sh
  msg at 6 sec.
  * con v1 y 5 fi;0.001;*dump picb6
  pause
  p
  p7;

  set prop off ; vi 1 1 5
  surf prps z 89
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec z 8 sh; vec y 1 sh; vec x 1 sh
  msg at 7 sec.
  * con v1 y 5 fi;0.001;*dump picb7
  pause
  p
  p8;

  set prop off ; vi 1 1 5
  surf prps z 89
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec z 8 sh; vec y 1 sh; vec x 1 sh
  msg at 8 sec.
  * con v1 y 5 fi;0.001;*dump picb8
  ENDUSE
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(Flow in stirred 3D vessel: 757          )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 757
 ************************************************************
  Group 2. Time dependence
 STEADY = F
    * Set overall time and no. of steps
 TFIRST =0. ;TLAST =8.
 FSTEP = 1 ;LSTEP = 8
 TFRAC(1)=0.125 ;TFRAC(2)=0.25
 TFRAC(3)=0.375 ;TFRAC(4)=0.5
 TFRAC(5)=0.625 ;TFRAC(6)=0.75
 TFRAC(7)=0.875 ;TFRAC(8)=1.
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 15
 XULAST =15.
 XFRAC(1)=0.066667 ;XFRAC(2)=0.133333
 XFRAC(3)=0.2 ;XFRAC(4)=0.266667
 XFRAC(5)=0.333333 ;XFRAC(6)=0.4
 XFRAC(7)=0.466667 ;XFRAC(8)=0.533333
 XFRAC(9)=0.6 ;XFRAC(10)=0.666667
 XFRAC(11)=0.733333 ;XFRAC(12)=0.8
 XFRAC(13)=0.866667 ;XFRAC(14)=0.933333
 XFRAC(15)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 15
 YVLAST =15.
 YFRAC(1)=0.066667 ;YFRAC(2)=0.133333
 YFRAC(3)=0.2 ;YFRAC(4)=0.266667
 YFRAC(5)=0.333333 ;YFRAC(6)=0.4
 YFRAC(7)=0.466667 ;YFRAC(8)=0.533333
 YFRAC(9)=0.6 ;YFRAC(10)=0.666667
 YFRAC(11)=0.733333 ;YFRAC(12)=0.8
 YFRAC(13)=0.866667 ;YFRAC(14)=0.933333
 YFRAC(15)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 15
 ZWLAST =15.
 ZFRAC(1)=0.066667 ;ZFRAC(2)=0.133333
 ZFRAC(3)=0.2 ;ZFRAC(4)=0.266667
 ZFRAC(5)=0.333333 ;ZFRAC(6)=0.4
 ZFRAC(7)=0.466667 ;ZFRAC(8)=0.533333
 ZFRAC(9)=0.6 ;ZFRAC(10)=0.666667
 ZFRAC(11)=0.733333 ;ZFRAC(12)=0.8
 ZFRAC(13)=0.866667 ;ZFRAC(14)=0.933333
 ZFRAC(15)=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(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,N,N,N,Y)
 SOLUTN(U1,Y,Y,N,N,N,Y)
 SOLUTN(V1,Y,Y,N,N,N,Y)
 SOLUTN(W1,Y,Y,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 =1000. ;TMP1 =0.
 EL1 =0.
 TSURR =0. ;TEMP0 =0.
 PRESS0 =0.
 DVO1DT =0. ;DRH1DP =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 ENUL =1. ;ENUT =0.
 PRNDTL(U1)=1. ;PRNDTL(V1)=1.
 PRNDTL(W1)=1.
 PRT(U1)=1. ;PRT(V1)=1.
 PRT(W1)=1.
 CP1 =1. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=0. ;FIINIT(U1)=0.
 FIINIT(V1)=0. ;FIINIT(W1)=0.
 FIINIT(PRPS)=67.
   No PATCHes yet used for this Group
 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(PRESREL ,CELL , 15, 15, 15, 15, 15, 15, 1, 8)
 COVAL(PRESREL ,P1 ,1000. ,0. )
 
 PATCH(EI ,EAST , 1, 15, 1, 15, 1, 15, 1, 8)
 COVAL(EI ,U1 ,In-Form:source - see Grp 19)
 
 PATCH(NI ,NORTH , 1, 15, 1, 15, 1, 15, 1, 8)
 COVAL(NI ,V1 ,In-Form:source - see Grp 19)
 
 PATCH(TIME=1 ,CELL , 8, 8, 4, 12, 6, 10, 1, 1)
 
 PATCH(TIME=2 ,CELL , 4, 12, 8, 8, 6, 10, 2, 2)
 
 PATCH(TIME=3 ,CELL , 8, 8, 4, 12, 6, 10, 3, 3)
 
 PATCH(TIME=4 ,CELL , 4, 12, 8, 8, 6, 10, 4, 4)
 
 PATCH(TIME=5 ,CELL , 8, 8, 4, 12, 6, 10, 5, 5)
 
 PATCH(TIME=6 ,CELL , 4, 12, 8, 8, 6, 10, 6, 6)
 
 PATCH(TIME=7 ,CELL , 8, 8, 4, 12, 6, 10, 7, 7)
 
 PATCH(TIME=8 ,CELL , 4, 12, 8, 8, 6, 10, 8, 8)
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 50 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-02
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;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(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(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 PARSOL = F
 ISG62 = 1
 SPEDAT(SET,PROPERTY,RHO1,C,=1000.)
 SPEDAT(SET,PROPERTY,ENUL,C,=1.)
 SPEDAT(SET,SOURCE,U1!EI,C,=1.E10*(2.*(7.5-YG)-U1)!IMAT>=100!LINE)
 SPEDAT(SET,SOURCE,V1!NI,C,=1.E10*(2.*(XG-7.5)-V1)!IMAT>=100!LINE)
 SPEDAT(SET,MATERIAL,67,L,T)
 SPEDAT(SET,STORED,PRPS,C,=67.!ZSLFIN)
 SPEDAT(SET,MATERIAL,100,L,T)
 SPEDAT(SET,STORED,PRPS!TIME=1,C,=100.!ZSLFIN)
 SPEDAT(SET,STORED,PRPS!TIME=2,C,=100.!ZSLFIN)
 SPEDAT(SET,STORED,PRPS!TIME=3,C,=100.!ZSLFIN)
 SPEDAT(SET,STORED,PRPS!TIME=4,C,=100.!ZSLFIN)
 SPEDAT(SET,STORED,PRPS!TIME=5,C,=100.!ZSLFIN)
 SPEDAT(SET,STORED,PRPS!TIME=6,C,=100.!ZSLFIN)
 SPEDAT(SET,STORED,PRPS!TIME=7,C,=100.!ZSLFIN)
 SPEDAT(SET,STORED,PRPS!TIME=8,C,=100.!ZSLFIN)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=4077. ;EX(U1)=1.238
 EX(V1)=1.336 ;EX(W1)=0.8445
 EX(PRPS)=67.440002
 ************************************************************
  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(W1,Y,N,Y,Y,Y,Y)
 OUTPUT(PRPS,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 8 ;IYMON = 8 ;IZMON = 8
 NPRMON = 50 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NTPRIN = 100000 ;ISTPRF = 1 ;ISTPRL = 100000
 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
 IDISPA = 1 ;IDISPB = 0 ;IDISPC = 0
 CSG1    ='P'
STOP