n0=0 n1=0 n2=0 n3=0 n4=0 n5=0 n10=0 nu=0 sumf=0 sums=0 sumfovers=0 1 read(5,5,end=100) f,sig 5 Format(12x,2f14.2) n0=n0+1 if(f.gt.sig) n1=n1+1 if(f.gt.2*sig) n2=n2+1 if(f.gt.3*sig) n3=n3+1 if(f.gt.4*sig) n4=n4+1 if(f.gt.5*sig) n5=n5+1 if(f.gt.10*sig) n10=n10+1 if(sig.eq.0.) nu=nu+1 sumf=sumf+f sums=sums+sig if (sig.gt.0) sumfovers=sumfovers+f/sig goto 1 100 write(6,*)' / (excluding onobserved):' write(6,*) sumfovers/(n0-nu),sumf/sums c write(6,*)' / (including onobserved):' c write(6,*) sumfovers/(n0),sumf/sums c write(6,*) 'nu, sumfovers, sumf, sums:',nu ,sumfovers, sumf, sums write(6,*) ' Total Unobs >1s >2s >3s >4s >5s >10s' write(6,6) n0,nu,n1,n2,n3,n4,n5,n10 write(6,7) float(n0)/n0,float(nu)/n0,float(n1)/n0,float(n2)/n0, . float(n3)/n0,float(n4)/n0,float(n5)/n0,float(n10)/n0 6 FORMAT(8I8) 7 FORMAT(8F8.4) end