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