integer jstart,jend real correlated_sampling_noise_rms real uncorrelated_sampling_noise_rms common /filterpars/jstart,jend, + correlated_sampling_noise_rms, + uncorrelated_sampling_noise_rms integer max_filter_length integer max_covariance_length parameter (max_filter_length=100) parameter (max_covariance_length=101*100/2) double precision sum1 double precision sumP double precision sumPt double precision sumPa(max_filter_length) double precision sumPta(max_filter_length) double precision suma(max_filter_length) double precision sumaa(max_covariance_length) double precision sumaaInv(max_covariance_length) double precision sunit(max_covariance_length) common /filterstats/sum1,sumP,sumPt,sumPa,sumPta,suma, + sumaa,sumaaInv,sunit real g0P,g0Pt real gP(max_filter_length) real gPt(max_filter_length) common /filtervecs/g0P,gP,g0Pt,gPt integer idsum1 integer idsumP integer idsumPt integer idsumPa integer idsumPta integer idsuma integer idsumaa integer idsumaaInv integer idsunit integer idgP,idg0P integer idgPt,idg0Pt integer idPin,idPout,idPinout integer idPtin,idPtout,idPtinout integer idtin,idtout,idtinout common /filterids/idsum1,idsumP,idsumPt,idsumPa,idsumPta, + idsuma,idsumaa,idsumaaInv,idgP,idgPt, + idPin,idPout,idPinout, + idPtin,idPtout,idPtinout, + idtin,idtout,idtinout data idsum1,idsumP,idsumPt/10,11,12/ data idsuma,idsumPa,idsumPta/13,14,15/ data idsumaa,idsumaaInv,idsunit/20,21,22/ data idgP,idg0P,idgPt,idg0Pt/1,2,3,4/ data idPin,idPout,idPinout/100,101,102/ data idPtin,idPtout,idPtinout/200,201,202/ data idtin,idtout,idtinout/300,301,302/