combine

 file{1}='atest80-512.fits'
 file{2}='atest81-512.fits'
 file{3}='atest84-512.fits'
 file{4}='atest86-512.fits'
 file{5}='atest87-512.fits'
 file{6}='atest88-512.fits'
 file{7}='atest89-512.fits'
 file{8}='atest90-512.fits'
 file{9}='atest91-512.fits'
 file{10}='atest92-512.fits'
 file{11}='atest93-512.fits'
 file{12}='atest95-512.fits'
 file{13}='atest97-512.fits'
 file{14}='atest98-512.fits'
 file{15}='atest99-512.fits'

0001 % file{1}='atest80-512.fits'
0002 % file{2}='atest81-512.fits'
0003 % file{3}='atest84-512.fits'
0004 % file{4}='atest86-512.fits'
0005 % file{5}='atest87-512.fits'
0006 % file{6}='atest88-512.fits'
0007 % file{7}='atest89-512.fits'
0008 % file{8}='atest90-512.fits'
0009 % file{9}='atest91-512.fits'
0010 % file{10}='atest92-512.fits'
0011 % file{11}='atest93-512.fits'
0012 % file{12}='atest95-512.fits'
0013 % file{13}='atest97-512.fits'
0014 % file{14}='atest98-512.fits'
0015 % file{15}='atest99-512.fits'
0016 
0017 
0018 numscans=24;
0019 nside = 256
0020 file{1}='atest80-256.fits'
0021 file{2}='atest81-256.fits'
0022 file{3}='atest84-256.fits'
0023 file{4}='atest86-256.fits'
0024 file{5}='atest87-256.fits'
0025 file{6}='atest88-256.fits'
0026 file{7}='atest89-256.fits'
0027 file{8}='atest90-256.fits'
0028 file{9}='atest91-256.fits'
0029 file{10}='atest92-256.fits'
0030 file{11}='atest93-256.fits'
0031 file{12}='atest95-256.fits'
0032 file{13}='atest97-256.fits'
0033 file{14}='atest98-256.fits'
0034 file{15}='atest99-256.fits'
0035 file{16}='atest101-256.fits'
0036 file{17}='atest102-256.fits'
0037 file{18}='atest103-256.fits'
0038 file{19}='atest104-256.fits'
0039 file{20}='atest105-256.fits'
0040 file{21}='atest106-256.fits'
0041 file{22}='atest107-256.fits'
0042 file{23}='atest108-256.fits'
0043 file{24}='atest109-256.fits'
0044 file{25}='atest110-256.fits'
0045 file{26}='atest111-256.fits'
0046 file{27}='atest112-256.fits'
0047 file{28}='atest113-256.fits'
0048 file{29}='atest114-256.fits'
0049 
0050 
0051 
0052 
0053 
0054 unix(['rm combine.fits']);
0055 
0056 
0057 sum=[];
0058 f=[];
0059 g=[];
0060 
0061 h = gov.fnal.eag.healpix.PixTools;
0062 npix= h.Nside2Npix(nside);
0063 sum=zeros(npix,6);
0064 
0065 for k=1:numscans
0066     for i=1:6
0067         f(:,i) = cbass_read_bin_col(file{k},2,i); 
0068     end
0069     %zero the additions for certain conditions
0070     indI=find(f(:,1)>=8);
0071     f(indI,6)=0;
0072     indI=find(f(:,6)<=6);
0073     f(indI,6)=0;
0074     for i=1:5
0075         g(:,i) = f(:,i).*f(:,6);
0076     end
0077     g(:,6)=f(:,6);
0078     
0079     sum=sum+g;
0080 end
0081 %
0082 % for i=1:6
0083 %     f(:,i) = cbass_read_bin_col(file{2},2,i);
0084 % end
0085 %
0086 % for i=1:5
0087 %     g(:,i) = f(:,i).*f(:,6);
0088 % end
0089 % g(:,6)=f(:,6);
0090 % sum = sum+g;
0091 %
0092 % for i=1:6
0093 %     f(:,i) = cbass_read_bin_col(file{3},2,i);
0094 % end
0095 %
0096 % for i=1:5
0097 %     g(:,i) = f(:,i).*f(:,6);
0098 % end
0099 % g(:,6)=f(:,6);
0100 % sum = sum+g;
0101 %
0102 % for i=1:6
0103 %     f(:,i) = cbass_read_bin_col(file{4},2,i);
0104 % end
0105 %
0106 % for i=1:5
0107 %     g(:,i) = f(:,i).*f(:,6);
0108 % end
0109 % g(:,6)=f(:,6);
0110 % sum = sum+g;
0111 
0112 
0113 for i=1:5
0114     meansum(:,i)=sum(:,i)./sum(:,6);
0115 end
0116 meansum(:,6)=sum(:,6);
0117 meansum(:,1)=meansum(:,1)-5;
0118 %meansum(find(isnan(meansum))) = -1;
0119 
0120  cbass_write_image('combine.fits',meansum,nside,npix,'test','test2')
0121