c----------------------------------------------
c      mesh1.f
c    
c      subrtn file : pk1subrtn
c       
c----------------------------------------------
c      main program 
c      mesh1.f     for mesh program << 1 color >>
c----( ln: point su, k1: belt su )
       parameter (ln=100,ln3=ln*3,ny=20,ln3y=ln3*ny,k1=2)
       parameter (r=0.0,g=0.7,b=0.7,pai=3.1415926,rd=10.0)  
c----( po: x,y,z no jun ni set )
       dimension po(ln3y)
c
       call initvrml
c
       dd=2.0*pai/float(ln)
c
       do 30  k=1,k1
c
       do 20  i=1,ln
       do 20  j=1,2
       i1=3*(i-1)+1+3*ln*(j-1)
       i2=i1+1
       i3=i2+1
       po(i1)=i
c      po(i2)=rd*cos(float(i+ny*(j-1))*dd)
       po(i2)=ny*(j-1)+ny*(k-1)
       po(i3)=0.0
   20  continue
c
       call face1(ln,po,r,g,b)
c
   30  continue
c
       stop
       end
