nav=74 name='qqq' isig=[7,8,9,10,11,18+lindgen(6)] nsig=n_elements(isig) pars=dblarr(nav,320,50,24) openr,1,'avlist' for iav=0,nav-1 do begin readf,1,name print,name par=readx(name) nmodes=n_elements(par)/24 for j=0,nmodes-1 do pars(iav,par(0,j),par(1,j),*)=par(*,j) end close,1 count=nint(reform(rebin(pars(*,*,*,2)<1,1,320,50,1)*nav)) av=dblarr(320,50,24) for i=0,23 do av(*,*,i)=reform(rebin(pars(*,*,*,i),1,320,50,1))*nav/count for i=0,nsig-1 do av(*,*,isig(i))=sqrt(reform(rebin(pars(*,*,*,isig(i))^2,1,320,50,1))*nav)/count openw,1,'av' openw,2,'av75' for l=0,319 do for n=0,49 do begin if (count(l,n) gt 0) then printf,1,l,n,av(l,n,2:23),format='(2i4,22g14.6)' if (count(l,n) ge 0.75*nav) then printf,2,l,n,av(l,n,2:23),format='(2i4,22g14.6)' end close,1 close,2 relerr=av(*,*,7)/av(*,*,2) wall=where(count ne 0) w75=where(count ge 0.75*nav) sigall=sqrt(1/total(1/relerr(wall)^2)) sig75=sqrt(1/total(1/relerr(w75)^2)) print,sigall,sig75 print,1/sigall,1/sig75 end