c----------------------------------------------
c      defusem1.f
c    
c      subrtn file : pk1subrtn
c       
c----------------------------------------------
c      main program  defusem1.f
c      defusem1.f    for DEF & USE program
c----( nx: x point, ny: y point )
       parameter (nx=10,ny=10)
       parameter (pai=3.1415926,rd=1.0)
c----( p: x,y,z no jun ni set   u: r,g,b no jun ni set )
       dimension p(3*nx*ny)
       dimension u(3*nx*ny)
c----( x1,y1,z1: heikou ido  x2,y2,z2: kaiten ido )
       parameter (x1=0,y1=0,z1=0,x2=0,y2=0.1,z2=0)
c----( r=kaitendo  n=kaisu )
       parameter (r=0.5708,n=100)
c
       call initvrml
c
c################################( test data )
       dd=2.0*pai/float(nx)
       ny1=ny-1
c
       do 10  jj=1,ny1
       do 20  i=1,nx
       do 20  j=1,2
       j1=j+jj-1
       i1=3*(i-1)+1+3*nx*(j1-1)
       i2=i1+1
       i3=i1+2
       p(i1)=i
       p(i3)=j1
       p(i2)=rd*sin(float(i)*dd)
       u(i1)=1.0-j1*0.1
       u(i2)=j1*0.1
       u(i3)=i*0.1
   20  continue
   10  continue
c#############################################
c
       call defusem(nx,ny,u,p,x1,y1,z1,x2,y2,z2,r,n)
c
       stop
       end

