subroutine bandread(vi,arad,pcol,dgn,wgt,rdmol,qrot)

      implicit double precision (o-z, a-h)
      include "parameters.dat"
      integer*2 di
      character name(4)*4

      dimension nso(4),nfund(4),qrot(4),vis(12,3,4)
      dimension strs(12,3,4),arads(12,3,4),rmass(3,4),rdmol(4),iso(3,4)
     & ,pcols(12,3,4),dgns(12,3,4),wgt(4)
      dimension str(12,4),arad(12,4),pcol(12,4),vi(12,4),dgn(12,4)

      data nso/3,2,1,1/, nmol/4/

      open(unit=60,file='bandparams.dat',status='old')
      do im = 1, nmol
         do is = 1, nso(im)
            read(unit=60,fmt=901) name(im),nfund(im),iso(is,im)
     &       ,rmass(is,im),rdmol(im),qrot(im)
            do nf = 1, nfund(im)
               read(unit=60,fmt=902) vis(nf,is,im),strs(nf,is,im)
     &          ,arads(nf,is,im),pcols(nf,is,im),di
               dgns(nf,is,im) = di
            end do
         end do
      end do
      close(unit=60)

      do im = 1, nmol 
         wgt(im) = rmass(1,im)
         do nf = 1, nfund(im)
            vi(nf,im) = vis(nf,1,im)
            str(nf,im) = strs(nf,1,im)
            arad(nf,im) = arads(nf,1,im)
            pcol(nf,im) = pcols(nf,1,im)
            dgn(nf,im) = dgns(nf,1,im)
         end do
      end do

      return

 901  format(1x,a4,1x,i2,1x,i4,1x,0pf7.4,1x,0pf3.1,1x,0pf3.1)
 902  format(1x,0pf7.2,1x,1pe10.3,1x,0pf6.2,1x,1pe10.3,1x,i2)

      end