c      main program  lineb3.f
c      and symblv test program
       parameter (lasl=12,npo=3*lasl)
       parameter (x1=-14.0,y1=26.0,z1=0.0,h1=3.4)
       parameter (x2=-14.0,y2=22.0,z2=0.0,h2=2.8)
       parameter (x3=-18.0,y3=-1.0,z3=0.0,h3=2.0)
       parameter (n1=38,n2=44,n3=6)
       dimension po(npo)
       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
       r=1.0
       g=1.0
       b=1.0
c
       call initvrml
c
       do 100 jj=1,7
       r=1.0
       g=1.0
       b=1.0
       if(jj.eq.2) b=0.0
       if(jj.eq.3) g=0.0
       if(jj.eq.4) r=0.0
       if(jj.eq.5) b=0.0
       if(jj.eq.5) g=0.0
       if(jj.eq.6) b=0.0
       if(jj.eq.6) r=0.0
       if(jj.eq.7) g=0.0
       if(jj.eq.7) r=0.0
c
       do 40 ii=1,lasl
       lin=ii
       li1=3*(lin-1)+1
       li2=3*(lin-1)+2
       li3=3*lin
c      x=po(li1)
c      y=po(li2)
c      z=po(li3)
       x=2.1*float(ii-1)-0.0*float(jj)
       y=2.1*float(ii-1)-0.0*float(jj)
       z=2.1*float(ii-1)-2.2*float(jj)
       po(li1)=x
       po(li2)=y
       po(li3)=z
   40  continue
c
       call line3d(lasl,r,g,b,po)
  100  continue
c
c
       call symbvb
       call symblv(x1,y1,z1,h1,chr1,n1)
       call symblv(x2,y2,z2,h2,chr2,n2)
       call symblv(x3,y3,z3,h3,chr3,n3)
       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,n1)
       r=1.0
       g=1.0
       b=0.0
       call symblvc(x2,y22,z2,h2,r,g,b,chr2,n2)
       r=0.0
       g=1.0
       b=1.0
       call symblvc(x3,y33,z3,h3,r,g,b,chr3,n3)
c
       stop
       end