1      subroutine argos_diana_report(sgmnam,xrms,isel)
2c
3      implicit none
4c
5#include "argos_diana_common.fh"
6#include "global.fh"
7#include "msgids.fh"
8c
9      character*16 sgmnam(nsa)
10      integer isel(nsa)
11      real*8 xrms(nsa)
12c
13      integer i,n
14      real*8 arms,bfact
15c
16c     report results
17c
18      bfact=128.0d0*atan(1.0d0)*atan(1.0d0)/3.0d0
19c
20      if(me.eq.0) then
21      write(lfnrms,1000)
22 1000 format('analysis')
23      endif
24      if(np.gt.1) then
25      call ga_dgop(mag_d12,xrms,nsa,'+')
26      endif
27      if(me.eq.0) then
28      write(lfnrms,1001) (sgmnam(i),i,isel(i),
29     + sqrt(xrms(i)/dble(ndata)),bfact*xrms(i)/dble(ndata),i=1,nsa)
30 1001 format(a16,i10,i5,2f12.6)
31      arms=0.0d0
32      n=0
33      if(me.eq.0) write(lfnrms,1000)
34      do 1 i=1,nsa
35      if(isel(i).ne.0) then
36      arms=arms+xrms(i)
37      n=n+1
38      endif
39      if(i.eq.nsa) goto 2
40      if(sgmnam(i)(11:16).ne.sgmnam(i+1)(11:16)) goto 2
41      goto 1
42    2 continue
43      if(ndata*n.gt.0) then
44      write(lfnrms,1002) sgmnam(i)(1:5),sgmnam(i)(11:16),
45     + sqrt(arms/dble(ndata*n)),bfact*arms/dble(ndata*n)
46 1002 format(a5,a6,2f12.6)
47      endif
48      arms=0.0d0
49      n=0
50    1 continue
51      close(unit=lfnrms)
52      endif
53      return
54      end
55c $Id$
56