photon use
  p
  s1
 
 
 
  con temp z 1 fi;0.02
  msg One sector only of a baffled vessel stirred by a paddle
  msg temperature contours and velocity vectors at time step 1
  gr ou z 1
  msg the paddle and baffle are on the right
  vec z 1;pause
  p
  s3
 
 
 
  msg temperature contours and velocity vectors at time step 3
  gr ou z 1
  con temp z 1 fi;0.02
  vec z 1;pause
  p
  s5
 
 
 
  msg temperature contours and velocity vectors at time step 5
  gr ou z 1
  con temp z 1 fi;0.02
  vec z 1;pause
  clear
  p
  s7
 
 
 
  msg temperature contours and velocity vectors at time step 7
  gr ou z 1
  con temp z 1 fi;0.02
  vec z 1;pause
  p
  s9
 
 
 
  msg temperature contours and velocity vectors at time step 9
  gr ou z 1
  con temp z 1 fi;0.02
  vec z 1;pause
  p
  s11
 
 
 
  con temp z 1 fi;0.02
  gr ou z 1
  msg One sector only of a baffled vessel stirred by a paddle
  msg temperature contours and velocity vectors at time step 11
  gr ou z 1
  vec z 1;pause
  p
  s13
 
 
 
  msg temperature contours and velocity vectors at time step 13
  gr ou z 1
  con temp z 1 fi;0.02
  vec z 1;pause
  p
  s15
 
 
 
  msg temperature contours and velocity vectors at time step 15
  gr ou z 1
  con temp z 1 fi;0.02
  vec z 1;pause
  clear
  p
  s17
 
 
 
  msg temperature contours and velocity vectors at time step 17
  gr ou z 1
  con temp z 1 fi;0.02
  vec z 1;pause
  p
  s19
 
 
 
  msg temperature contours and velocity vectors at time step 19
  gr ou z 1
  con temp z 1 fi;0.02
  vec z 1;pause
  ext
 
 
 
 
  view z
  msg last time step
  gr ou z 1
  vec z 1;pause sh
 
  msg pressure contours
  con p1 fi
  z 1
  1 20 1 11
  0.01
  pause
  msg press ENTER for the outer half
 
  con p1 fi
  z 1
  1 20 12 22
  0.01
 
  clear
  con off
  red
  msg temperature contours
  con temp fi
  z 1
  1 20 1 11
  0.01
  pause
  msg press ENTER for the outer half
 
  con temp fi
  z 1
  1 20 12 22
  0.01
  vec z 1
  enduse
 
    GROUP 1. Run title and other preliminaries
TEXT(Rotating Inner And Fixed Outer Grid 
TITLE
REAL(XLENGTH,YLENGTH,ZLENGTH,OMEGA,RELVEL,JETVEL,DELTAT,DTF,V1IN)
REAL(TIN,TOUT)
INTEGER(NYNOM,IXSHFT,ISHIFT)
INTEGER(IXPAD,IYPAD,IXBAF,IYBAF,IXJET,IYJET)
    **** data input ****
TIN=1.0;TOUT=0.0
RELVEL=1.0;V1IN=RELVEL*0.1;JETVEL=RELVEL*2.0
XLENGTH=1.0;YLENGTH=1.0;ZLENGTH=1.0
OMEGA=RELVEL*2.0/YLENGTH
DELTAT=XLENGTH/(OMEGA*NX)
DTF=DELTAT
NYNOM=10;NX=10;NZ=1
NYNOM=20;NX=20
NPHI=20
NY=NYNOM+2
IXBAF=NX;IYBAF=NY-NYNOM/3
IXPAD=NX;IYPAD=NYNOM/3
IXJET=IXPAD;IYJET=IYPAD
MESG(NX=:NX:; NYNOM=:NYNOM:; NZ=:NZ:
MESGB(XLENGTH=:XLENGTH:; YLENGTH=:YLENGTH:; ZLENGTH=:ZLENGTH:
mesgb(angular velocity is :omega: radians per second
mesgb(       time step is :deltat: seconds
DELAY(500)
ISHIFT=2-IXSHFT*NY
mesgm(ishift is the offset between linked points, =2 - ixshft*ny
ISHIFT
    GROUP 2. Time dependence
STEADY=F;LSTEP=NX;TFRAC(1)=-LSTEP;TFRAC(2)=DELTAT
    GROUP 3. X-direction grid specification
   **Domain is XLENGTH m long in x-direction, with equal intervals
CARTES=F;GRDPWR(X,NX,XLENGTH,1.0)
    GROUP 4. Y-direction grid specification
YVLAST=YLENGTH
YFRAC(1)=-NYNOM/2
YFRAC(2)=1/NYNOM
YFRAC(3)=2
YFRAC(4)=0.1*YFRAC(2)
YFRAC(5)=-YFRAC(1)
YFRAC(6)=YFRAC(2)
    GROUP 5. Z-direction grid specification
   **Domain is ZLENGTH m long in z-direction, with equal intervals
GRDPWR(Z,NZ,ZLENGTH,1.0)
    GROUP 7. Variables stored, solved & named
   **Choose first-phase enthalpy (H1) as dependent variable
     and activate the whole-field elliptic solver
SOLUTN(H1,Y,Y,Y,N,N,N);NAME(H1)=TEMP
SOLVE(P1,U1,V1);SOLUTN(P1,Y,Y,Y,N,N,N)
    GROUP 8. Terms (in differential equations) & devices
   ** cut out built-in source term
TERMS(TEMP,N,Y,Y,Y,Y,Y)
    GROUP 9. Properties of the medium (or media)
ENUT=0.002*RELVEL*YLENGTH
PRT(TEMP)=0.7;PRNDTL(TEMP)=0.7
    GROUP 11. Initialization of variable or porosity fields
IURINI=-1
FIINIT(TEMP)=0.0
PATCH(LOWER,INIVAL,1,NX,1,NYNOM/2+1,1,1,1,1)
COVAL(LOWER,TEMP,0.0,TIN)
    GROUP 13. Boundary conditions and special sources
XCYCLE=T
   ** inflow boundary
PATCH(INFLOW,SOUTH,1,NX,1,1,1,1,1,LSTEP)
COVAL(INFLOW,V1,ONLYMS,V1IN)
COVAL(INFLOW,TEMP,ONLYMS,TIN)
COVAL(INFLOW,P1,FIXFLU,V1IN)
IURVAL=-1;IURPRN=-1
   ** outlet boundary
PATCH(OUTFLOW,CELL,1,NX,NY,NY,1,1,1,LSTEP)
COVAL(OUTFLOW,TEMP,ONLYMS,SAME)
COVAL(OUTFLOW,P1,1.E-3,0.0)
PATCH(+TE,NORTH,1,NX,NYNOM/2+1,NYNOM/2+1,1,1,1,LSTEP)
COVAL(+TE,TEMP,0.0,ISHIFT)
COVAL(+TE,P1,fixp,ISHIFT)
COVAL(+TE,U1,-OMEGA,ISHIFT)
COVAL(+TE,V1,0.0,ISHIFT)
PATCH(VRELAX,CELL,1,NX,NYNOM/2+1,NYNOM/2+1,1,1,1,LSTEP)
COVAL(VRELAX,V1,1.E3,SAME)
PATCH(+TW, SOUTH,1,NX,NYNOM/2+2,NYNOM/2+2,1,1,1,LSTEP)
COVAL(+TW,TEMP,fixval,-ISHIFT)
COVAL(+TW,P1,fixp,-ISHIFT)
COVAL(+TW,U1,OMEGA,-ISHIFT)
PATCH(GP12DFN1,NORTH,1,NX,NYNOM/2,NYNOM/2,1,1,1,LSTEP)
COVAL(GP12DFN1,TEMP,0.5,0.0)
PATCH(GP12DFN2,NORTH,1,NX,NYNOM/2+2,NYNOM/2+2,1,1,1,LSTEP)
COVAL(GP12DFN2,TEMP,0.5,0.0)
COVAL(GP12DFN2,U1,0.5,0.0)
PATCH(ROTA,PHASEM,1,NX,1,NYNOM/2+1,1,1,1,LSTEP)
COVAL(ROTA,U1,FIXFLU,GRND)
COVAL(ROTA,V1,FIXFLU,GRND)
ANGVEL=OMEGA
PATCH(PADDLE,EAST,IXPAD,IXPAD,1,IYPAD,1,1,1,LSTEP)
COVAL(PADDLE,U1,FIXVAL,0.0)
PATCH(BAFFLE,EAST,IXBAF,IXBAF,IYBAF,NY,1,1,1,LSTEP)
COVAL(BAFFLE,U1,FIXVAL,0.0)
PATCH(JET,NORTH,IXJET,IXJET,IYJET,IYJET,1,1,1,LSTEP)
COVAL(JET,V1,FIXVAL,JETVEL)
    GROUP 15. Termination of sweeps
LSWEEP=50;RESREF(P1)=1.E-10
    GROUP 16. Termination of iterations
LITER(TEMP)=10;LITER(P1)=20
    GROUP 17. Under-relaxation devices
RELAX(P1,LINRLX,0.5);RELAX(V1,FALSDT,DTF);RELAX(U1,FALSDT,DTF)
    GROUP 21. Print-out of variables
   **Print fields of temperature
OUTPUT(TEMP,Y,Y,Y,Y,Y,Y)
    GROUP 22. Spot-value print-out
IYMON=NY/2+1;IZMON=NZ/2+1;ixmon=nx-1;itabl=1
    GROUP 23. Field print-out and plot control
IXPRF=NYNOM/2-1;IXPRL=NYNOM/2+3
NXPRIN=1;NYPRIN=1
    GROUP 24. Dumps for restarts
CSG1=STIR;IDISPA=1;IDISPB=1;IDISPC=LSTEP
UWATCH=T;NPLT=1;TSTSWP=-1;SELREF=T;RESFAC=1.E-1
NYPRIN=NY/5;NXPRIN=NX/5
IXPRF=1;IXPRL=NX
RINNER=YVLAST/10
LITER(U1)=10;LITER(V1)=10