PHOTON USE
    x
    x1


    msg(  PLANT generated grid
    gr z 1
    ENDUSE

  DISPLAY

    PLANT is used to create body-fittted-coordinate grids
    defined by formulae.
    

    Example: 2D periodically broken channel
  ENDDIS


TEXT( 2D periodically broken channel

    GROUP 6. Body-fitted coordinates or grid distortion
nx=150;ny=20;nz=1
bfc=t;lsweep=1
dumpc(grid);readco(grid)
CSG1=PHI;CSG2=XYZ;IDISPA=1
store(mark)
NAMSAT=MOSG
    
    PLANTBEGIN
real(twopi)
twopi=2.0*3.14157
real(a,alfa)
a=0.75;alfa=twopi/12;RG(1)=a;RG(2)=alfa
    XC=12*1.5714/FLOAT(NX)*FLOAT(I-1)
    YC=8/:TWOPI:*RG(1)/RG(2)*$
               (SIN(RG(2))*SIN(XC)+$
               SIN(3*RG(2))*SIN(3*XC)/9+$
               SIN(5*RG(2))*SIN(5*XC)/25+$
               SIN(7*RG(2))*SIN(7*XC)/49+$
               SIN(9*RG(2))*SIN(9*XC)/81+$
               SIN(11*RG(2))*SIN(11*XC)/121+$
               SIN(13*RG(2))*SIN(13*XC)/169+$
               SIN(15*RG(2))*SIN(15*XC)/225+$
               SIN(17*RG(2))*SIN(17*XC)/289+$
               SIN(19*RG(2))*SIN(19*XC)/361)+$
                2.*FLOAT(J-1)/FLOAT(NY)
    ZC=1.*FLOAT(K-1)/FLOAT(NZ)
    PLANTEND
dmpstk=t
DISTIL=T
EX(MARK)=1.000E-10
 LIBREF=512
STOP