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
  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>

  DISPLAY 
    This input  file  demonstrates  how  to  use  parametric
    analitics to generate the time varying 2D BFC grids.  It
    is  wholly  focused  on  calculation  of   grid   corner
    coordinates. Therefore, no other functions are supported
    by input data.
  ENDDIS
 
    PLANT information :
     * Data input groups used:  6
     * Ground groups planted : 19-2
     * Headings used  : MXYZ??
     * Functions used : None
     * Commands used  : None
 
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
 
   *  GROUP 1. Run title and other preliminaries.
TEXT( Gradually corrugated channel
 
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
dumpc(grid);readco(grid)
NAMSAT=MOSG

     PLANTBEGIN
    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 <<<<<<<<<<<<<<<<<<<<<
      PLANTEND
SPEDAT(SET,GXMONI,TRANSIENT,L,F) 
CSG1=PHI;CSG2=XYZ;IDISPA=1
SOLVE(MARK)
dmpstk=t
DISTIL=T
EX(MARK)=1.000E-10
 LIBREF=609