program rdmags_sp
      parameter(nmag=47)
      real tb(90),mt(90),mbol(90),mags(90,nmag)

      character*3  mcod(7)
      data mcod /'m23','m18','m13','m07','m03','p00','p02'/
      
      character cwhat*1,fn11*80 
c     _________________________________________________________

      write(*,*) 'ikind=?, 1: bb; 2:ss'
      read(*,*) ikind

      if(ikind.eq.1) cwhat='b'
      if(ikind.eq.2) cwhat='s'

      do nz=1,7
         fn11 = mcod(nz)//cwhat//cwhat//'_sp.mags'
         open (2,file=fn11,form='formatted',status='unknown') 
         read (2,*)
         do i=1,90
           read (2,59) tb(i), mt(i), mbol(i), (mags(i,k),k=1,nmag)
         enddo  
         close(2)
 59      format(e12.5, 49f10.5)     
      enddo
            
      stop
      end