c      main program  point1.fi
c      and symblv test program
c      parameter (nx=10,ny=10,nz=10)
       parameter (nx=20,ny=20,nz=20)
       parameter (n1=nx,n2=n1*ny,n3=n2*nz)
       parameter (x0=0.0,y0=0.0,z0=0.0)
       parameter (xl=9.0,yl=9.0,zl=9.0)
c      parameter (ico=5,rd=1.0)
       parameter (ico=7,rd=1.0)
c      parameter (x1=-14.0,y1=26.0,z1=0.0,h1=3.4)
c      parameter (x2=-14.0,y2=22.0,z2=0.0,h2=2.8)
c      parameter (x3=-18.0,y3=-1.0,z3=0.0,h3=2.0)
       parameter (x1= -3.0,y1=12.0,z1=0.0,h1=0.8)
       parameter (x2= -3.0,y2=10.0,z2=0.0,h2=0.7)
       parameter (x3= -4.0,y3= 7.0,z3=0.0,h3=0.5)
       parameter (nc1=38,nc2=44,nc3=6)
       dimension u(n3)
       character chr1*60,chr2*60,chr3*6
       data chr1/"Southward Turning from Northward IMF"/
       data chr2/"Bz=18.2nT Nsw=46/cc Vsw=412km/s t=342m(42m)"/
       data chr3/"X 8Re"/
c
       hx=xl/float(nx-1)
       hy=yl/float(ny-1)
       hz=zl/float(nz-1)
c
       i1=101
       write(6,*) i1
       do 10 k=1,nz
       do 10 j=1,ny
       do 10 i=1,nx
c      x=hx*float(i-1)
c      y=hy*float(j-1)
c      z=hz*float(k-1)
       x=0.5*hx*float(2*i-1-nx)
       y=0.5*hy*float(2*j-1-ny)
       z=0.5*hz*float(2*k-1-nz)
       r2=x*x+y*y+z*z
       r1=sqrt(r2)
       i1=i+n1*(j-1)+n2*(k-1)
       u(i1)=exp(-(r1/rd)**2)
       u(i1)=1.0/(1.0+r2/rd)
c      u(i1)=x/xl
   10  continue
c
       i1=101
       write(6,*) i1
c
       call initvrml
c
       call point3d(nx,ny,nz,x0,y0,z0,xl,yl,zl,ico,u)
c
c
       call symbvb
       call symblv(x1,y1,z1,h1,chr1,nc1)
       call symblv(x2,y2,z2,h2,chr2,nc2)
       call symblv(x3,y3,z3,h3,chr3,nc3)
       r=1.0
       g=0.0
       b=0.0
       y11=y1-10.0
       y22=y2-10.0
       y33=y3-10.0
       call symblvc(x1,y11,z1,h1,r,g,b,chr1,nc1)
       r=1.0
       g=1.0
       b=0.0
       call symblvc(x2,y22,z2,h2,r,g,b,chr2,nc2)
       r=0.0
       g=1.0
       b=1.0
       call symblvc(x3,y33,z3,h3,r,g,b,chr3,nc3)
c
       stop
       end