GROUP 1. Run title and other preliminaries
TEXT( Flow in stirred 3D vessel
REAL(PROP);PROP=50.;RG(7)=PROP
  GROUP 2. Transience; time-step specification
STEADY=F;GRDPWR(T,8,8.,1.0)
REAL(REVL)
     Number of revolutions, 1
REVL=2.
     Angular velocity, 1/sec
RG(10)=REVL/TLAST
     Origin of rotation axes
     XIC  , m
RG(11)=7.5
     YIC  , m
RG(12)=7.5
 
    GROUP 3. X-direction grid specification
GRDPWR(X,15,15.,1.0)
    GROUP 4. Y-direction grid specification
GRDPWR(Y,15,15.,1.0)
GRDPWR(Z,15,15.,1.0)
    GROUP 7. Variables stored, solved & named
SOLVE(P1,U1,V1,W1)
STORE(MARK)
    GROUP 8. Terms (in differential equations) & devices
    GROUP 9. Properties of the medium (or media)
RHO1=1000.; ENUL=1.0
    GROUP 11. Initialization of variable or porosity fields
FIINIT(P1)=0.0;FIINIT(U1)=0.0;FIINIT(V1)=0.0;FIINIT(W1)=0.0
iniadd=f;FIINIT(MARK)=0.0
 
PATCH(INIBLADE,INIVAL,8,8,4,12,6,10,6,10)
INIT( INIBLADE,MARK, 0.0, PROP)
 
    GROUP 13. Boundary conditions and special sources
 
PATCH(PRESREL,CELL,NX,NX,NY,NY,NZ,NZ,1,LSTEP)
COVAL(PRESREL,P1,1000.,0.0)

    GROUP 15. Termination of sweeps
LSWEEP=50
    GROUP 17. Under-relaxation devices
relax(p1,linrlx,0.5)
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
        Automatic local under-relaxation
RG(99)=ENUL
 
  PLANTBEGIN
PATCH(SS050U,EAST,1,NX,1,NY,1,NZ,1,lstep)
   CO=1.e10
   VAL= RG(10)*(YG2D-RG(12))
COVAL(SS050U,U1,GRND,GRND)
 
PATCH(SS050V,NORTH,1,NX,1,NY,1,NZ,1,lstep)
   CO=1.e10
   VAL=-RG(10)*(XG2D-RG(11))
COVAL(SS050V,V1,GRND,GRND)
 
PATCH(RELAX,PHASEM,1,NX,1,NY,1,NZ,1,1)
   CO=SQRT(U1**2+W1**2+V1**2)/$
              AMIN1(DXU2D*1,AMIN1(DYV2D*1,DZ*1))+$
              RG(99)/AMIN1(DXU2D*1,AMIN1(DYV2D*1,DZ*1))**2
COVAL(RELAX,U1,GRND,SAME)
   CO=SQRT(U1**2+W1**2+V1**2)/$
              AMIN1(DXU2D*1,AMIN1(DYV2D*1,DZ*1))+$
              RG(99)/AMIN1(DXU2D*1,AMIN1(DYV2D*1,DZ*1))**2
COVAL(RELAX,V1,GRND,SAME)
   CO=SQRT(U1**2+W1**2+V1**2)/$
              AMIN1(DXU2D*1,AMIN1(DYV2D*1,DZ*1))+$
              RG(99)/AMIN1(DXU2D*1,AMIN1(DYV2D*1,DZ*1))**2
COVAL(RELAX,W1,GRND,SAME)
STORE(TFAL);OUTPUT(TFAL,Y,Y,Y,Y,Y,Y)
   TFAL=1/(SQRT(U1**2+W1**2+V1**2)/$
               AMIN1(DXU2D*1,AMIN1(DYV2D*1,DZ*1))+$
               RG(99)/AMIN1(DXU2D*1,AMIN1(DYV2D*1,DZ*1))**2)
STORE(CONT,ISTP)
 
   ISTP=FLOAT(ISTEP)*$
           AMAX1(AMIN1(1.0,5.-FLOAT(ISTEP)),0.0)
 
   CONT=ISTP+AMAX1(0.0,FLOAT(ISTEP)-4.)
 
   MARK=RG(7)*AMAX1(0.0,1.-CONT)
  REGION(8,8,4,7,6,10)
 
   MARK=RG(7)*AMAX1(0.0,1.-CONT)
  REGION(8,8,9,12,6,10)
 
   MARK=RG(7)*AMAX1(0.0,2.-CONT)
  REGION(4,7,8,8,6,10)
 
   MARK=RG(7)*AMAX1(0.0,2.-CONT)
  REGION(9,12,8,8,6,10)
 
   MARK=RG(7)*AMAX1(-AMIN1(0.0,1.-CONT),0.0)
  REGION(8,8,4,7,6,10)
 
   MARK=RG(7)*AMAX1(-AMIN1(0.0,1.-CONT),0.0)
  REGION(8,8,9,12,6,10)
 
   MARK=MARK+2.*RG(7)*AMIN1(0.0,2.-CONT)
  REGION(8,8,4,7,6,10)
 
   MARK=MARK+2.*RG(7)*AMIN1(0.0,2.-CONT)
  REGION(8,8,9,12,6,10)
 
   MARK=MARK-RG(7)*AMIN1(0.0,2.-CONT)
  REGION(4,7,8,8,6,10)
 
   MARK=MARK-RG(7)*AMIN1(0.0,2.-CONT)
  REGION(9,12,8,8,6,10)
 
   MARK=MARK-2.*RG(7)*AMIN1(0.0,3.-CONT)
  REGION(8,8,4,7,6,10)
 
   MARK=MARK-2.*RG(7)*AMIN1(0.0,3.-CONT)
  REGION(8,8,9,12,6,10)
 
   MARK=MARK+2.*RG(7)*AMIN1(0.0,3.-CONT)
  REGION(4,7,8,8,6,10)
 
   MARK=MARK+2.*RG(7)*AMIN1(0.0,3.-CONT)
  REGION(9,12,8,8,6,10)
  PLANTEND

NAMSAT=MOSG
 
    GROUP 22. Spot-value print-out
NPRMON=LSWEEP;IXMON=NX/2+1;IYMON=NY/2+1;IZMON=NZ/2+1
    GROUP 23. Print-out & plot control
TSTSWP=-1
idispa=1; csg1=p; selref=t; resfac=1.e-2
  PHOTON USE
  p
  p1;
 
  con mark z 8 sh
  25 25.1 0;int 1
  1
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec y 8 sh;vec z 8 sh
  msg at 1 sec.
  pause
  p
  p2;
 
  con mark z 8 sh
  25 25.1 0;int 1
  1
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec y 8 sh;vec z 8 sh
  msg at 2 sec.
  pause
  p
  p3;
 
  con mark z 8 sh
  25 25.1 0;int 1
  1
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec y 8 sh;vec z 8 sh
  msg at 3 sec.
  pause
  p
  p4;
 
  con mark z 8 sh
  25 25.1 0;int 1
  1
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec y 8 sh;vec z 8 sh
  msg at 4 sec.
  pause
  p
  p5;
 
  con mark z 8 sh
  25 25.1 0;int 1
  1
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec y 8 sh;vec z 8 sh
  msg at 5 sec.
  pause
  p
  p6;
 
  con mark z 8 sh
  25 25.1 0;int 1
  1
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec y 8 sh;vec z 8 sh
  msg at 6 sec.
  pause
  p
  p7;
 
  con mark z 8 sh
  25 25.1 0;int 1
  1
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec y 8 sh;vec z 8 sh
  msg at 7 sec.
  pause
  p
  p8;
 
  con mark z 8 sh
  25 25.1 0;int 1
  1
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec y 8 sh;vec z 8 sh
  msg at 8 sec.
  ENDUSE
dmpstk=t
DISTIL=T
EX(P1)=1.010E+02; EX(U1)=1.093E-01; EX(V1)=1.072E-01
EX(W1)=2.242E-02; EX(CONT)=4.000E+00
EX(TFAL)=8.665E-01; EX(MARK)=6.667E-01
 LIBREF=116
STOP