TALK=F;RUN( 1, 1);VDU=VGAMOUSE

TEXT( Library case Y606: Gradually corrugated channel

  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    PLANT information :
     * Data input groups used:  6
     * Ground groups planted : 19-2
     * Headings used  : MXYZ??
     * Functions used : None
     * Commands used  : None
  <<<<<<<<<<<<<<<<<<<<<<<<<  Comment ends <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<

REAL(LENGTH,WIDTH,RE,PR,PI,UIN)
PI=3.14157;LENGTH=6*PI;WIDTH=2.;RG(1)=LENGTH;RG(2)=WIDTH
steady=f
grdpwr(t,8,4800,1.)
   *  GROUP 3. X-direction grid specification.
NX=150
   *  GROUP 4. Y-direction grid specification.
NY=20
bfc=t

NAMSAT=MOSG
    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=ZWLAST*FLOAT(K-1)/FLOAT(NZ)
  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    The coefficients of formula for the Y coordinate of cell corners are
    made linear function of the time moments.
  <<<<<<<<<<<<<<<<<<<<<<<<<  Comment ends <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<

CSG1=PHI;CSG2=XYZ;IDISPA=1
STORE(MARK)
  PHOTON USE
   x;x1;;
   gr z 1
  msg            Grid at t=600 sec
  msg  Press Enter to continue
   pause
   x;x2;;
   gr z 1
  msg            Grid at t= 1200 sec
  msg  Press Enter to continue
   pause
   x;x3;;
   gr z 1
  msg            Grid at t= 1800 sec
  msg  Press Enter to continue
   pause
   x;x4;;
   gr z 1
  msg            Grid at t=2400 sec
  msg  Press Enter to continue
   pause
   x;x5;;
   gr z 1
  msg            Grid at t= 3000 sec
  msg  Press Enter to continue
   pause
   x;x6;;
   gr z 1
  msg            Grid at t= 3600 sec
  msg  Press Enter to continue
   pause
   x;x7;;
   gr z 1
  msg            Grid at t= 4200 sec
  msg  Press Enter to continue
   pause
   x;x8;;
   gr z 1
  msg            Grid at t= 4800 sec
  msg  Press Enter to continue
  msg Press Enter to end
  ENDUSE
STOP