PHOTON USE
   p
   p1
 
   x1
 
   gr ou z 1
   vec z 1 sh
  msg            Velocity distribution at t=600 sec
  msg  Press Enter to continue
   pause
   p
   p2
 
   x2
 
   gr ou z 1
   vec z 1 sh
  msg            Velocity distribution at t= 1200 sec
  msg  Press Enter to continue
   pause
   p
   p3
 
   x3
 
   gr ou z 1
   vec z 1 sh
  msg            Velocity distribution at t= 1800 sec
  msg  Press Enter to continue
   pause
   p
   p4
 
   x4
 
   gr ou z 1
   vec z 1 sh
  msg            Velocity distribution at t=2400 sec
  msg  Press Enter to continue
   pause
   p
   p5
 
   x5
 
   gr ou z 1
   vec z 1 sh
  msg            Velocity distribution at t= 3000 sec
  msg  Press Enter to continue
   pause
   p
   p6
 
   x6
 
   gr ou z 1
   vec z 1 sh
  msg            Velocity distribution at t= 3600 sec
  msg  Press Enter to continue
   pause
   p
   p7
 
   x7
 
   gr ou z 1
   vec z 1 sh
  msg            Velocity distribution at t= 4200 sec
  msg  Press Enter to continue
   pause
   p
   p8
 
   x8
 
   gr ou z 1
   vec z 1 sh
  msg            Velocity distribution at t= 4800 sec
  msg  Press Enter to continue
  msg Press Enter to end
  ENDUSE
 
   *  GROUP 1. Run title and other preliminaries.
 
  DISPLAY
   Flow in a gradually corrugated channel
  ENDDIS
 
TEXT( Flow in gradually corrugated channel
 
steady=f
grdpwr(t,8,4800,1.)
   *  GROUP 3. X-direction grid specification.
NX=150
GRDPWR(X,NX,LENGTH,1.0)
   *  GROUP 4. Y-direction grid specification.
NY=20
GRDPWR(Y,NY,1,1.0)
bfc=t
nonort=t
NAMSAT=MOSG

   PLANTBEGIN
REAL(LENGTH,WIDTH,RE,PR,PI,UIN)
PI=3.14157;LENGTH=6*PI;WIDTH=2.;RG(1)=LENGTH;RG(2)=WIDTH
UIN=1.0;RE=800.0;PR = 1.0
    XC=RG(1)/FLOAT(NX)*FLOAT(I-1)
    YC=1/:PI:+$
               0.5*((TIM-600)/4200)*SIN(XC-1.6)-$
               2/3.14157*($
                   ((TIM-600)/4200)*COS(2*(XC-1.6))/3+$
                   ((TIM-600)/4200)*COS(4*(XC-1.6))/15+$
                   ((TIM-600)/4200)*COS(6*(XC-1.6))/35+$
                   ((TIM-600)/4200)*COS(8*(XC-1.6))/63+$
                   ((TIM-600)/4200)*COS(10*(XC-1.6))/99+$
                   ((TIM-600)/4200)*COS(12*(XC-1.6))/143)+$
                    RG(2)*FLOAT(J-1)/FLOAT(NY)
    ZC=1.*FLOAT(K-1)/FLOAT(NZ)
   PLANTEND
 
   *  GROUP 7. Variables stored, solved & named.
SOLVE(P1); SOLUTN(P1,Y,Y,Y,N,N,N)
SOLVE(U1,V1)
   *  GROUP 8. Terms (in differential equations) & devices.
  TERMS(H1,N,Y,Y,N,Y,N)
RHO1= 1.000E+00   ;  ENUL= 1./RE
   *  GROUP 11. Initialization of variable or porosity fields.
FIINIT(P1)= 1.300E-04;FIINIT(U1)= 1.000E+00
   *  GROUP 13. Boundary conditions and special sources.
PATCH(INLET,WEST,1,1,1,NY,1,1,1,lstep)
COVAL(INLET,P1,FIXFLU,UIN*RHO1)
COVAL(INLET,U1,ONLYMS,UIN)
     ** Fixed pressure outlet
PATCH(OUTLET,EAST,NX,NX,1,NY,1,1,1,lstep)
COVAL(OUTLET,P1,1000.*fixp,0.0)
COVAL(OUTLET,U1,onlyms,0.0)
COVAL(OUTLET,V1,onlyms,0.0)
     ** HOT top wall
WALL(TOPWALL,NORTH,1,NX,NY,NY,1,1,1,lstep)
     ** COLD bottom wall
WALL(BOTWALL,SOUTH,1,NX,1,1,1,1,1,lstep)
 
   *  GROUP 15. Termination of sweeps.
LSWEEP=100
RESREF(P1)= 1.000E-06;RESREF(U1)= 1.000E-06
RESREF(V1)= 1.000E-06
   *  GROUP 16. Termination of iterations.
   *  GROUP 17. Under-relaxation devices.
RELAX(P1,LINRLX, 0.5)
REAL(SCALEL,SCALEU);SCALEL= 1.0;SCALEU=1.00
RELAX(U1,FALSDT, .30*SCALEL/SCALEU)
RELAX(V1,FALSDT, .30*SCALEL/SCALEU)
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
   *  GROUP 21. Print-out of variables.
OUTPUT(P1,Y,Y,Y,Y,Y,Y);OUTPUT(U1,Y,Y,Y,Y,Y,Y)
OUTPUT(V1,Y,Y,Y,Y,Y,Y)
   *  GROUP 22. Spot-value print-out.
IXMON=7;IYMON=8
   *  GROUP 23. Field print-out and plot control.
tstswp=-1;CSG1=PHI;CSG2=XYZ;IDISPA=1
  bfcsym=t

tstswp=-1
dmpstk=t
DISTIL=T
EX(P1)=4.093E-01; EX(U1)=9.360E-01; EX(V1)=1.790E-01
EX(WCRT)=1.000E-10; EX(VCRT)=1.752E-01; EX(UCRT)=9.345E-01
 LIBREF=515
STOP