implicit  real*8 (a-h,o-z)
 include 'parameter'
 dimension s(3,3), q(3,3), t(3,3), r(3,3), b(3), a(3), rot(mpla3,mpla3), bb(mpla3), aa(mpla3), rr(mpla3, mpla3), ss(mpla3, mpla3), tt(mpla3, mpla3), roti(mpla3, mpla3)
 character bpla(mpla)*2
 include 'constant'

 read (5,4) npla

 do i=1,npla {
    do j=1,3 {
       do k=1,3 {
          s(j,k)= .0
          q(j,k)= .0
          }
       }
    read (5,10) bpla(i), blat, blon, brad
    read (5,2) 
    read (5,2) 
    read (5,20) s(1,1), s(2,2), s(3,3), area1, trace1, s(1,2), s(1,3), s(2,3), sum1 
    read (5,2) 
    read (5,20) q(1,1), q(2,2), q(3,3), q(1,2), q(1,3), q(2,3)
    read (5,2) 
    q(2,1)= q(1,2)
    q(3,1)= q(1,3)
    q(3,2)= q(2,3)

    call sphcar(blat*co, blon*co, brad*co, b)

   call matvec(q, b, a, 3, 3)
   rms1= dot(b, a, 3)
    rms= sqrt(rms1/area1) *erad

   
  area9= area9 + area1
   rms9= rms9 + rms *area1
   
  write(6,14) bpla(i), blat, blon, brad, rms, area1
    }

 write(6,50) area9, rms9, rms9/area9
 50 format('  area', f8.4, '  rms9', f8.4, '  meanrms', f8.4)

  2 format()
  4 format(i4)
  6 format(i2)
 10 format(a2, f11.4, f11.4, f11.4)
 12 format(a2, f11.4, f11.4, f11.4, f11.8)
 14 format(a2, f9.2, ',', f9.2, ',', f8.3, ',', f7.1, ',', f9.4)
 20 format(9f12.8)
 23 format(f10.1, e25.16, i4, i4, '  ', a2, i1, ' ', a2, i1 )
 30 format(i4, i4, f11.4)
101 format(' ',a2,'  -',a2,'  ',f9.3,f9.3,f8.4)
103 format(3f10.1)

 end