start_day=datenum(2021,11,29);
mypath='D:/AIRS/dATA/2021';
%mypath='H:/data/AIRS/AIRS3STD.006';
%datestr( start_day,'yyyy-mm-dd')
for i = 1:50
days=[];
mydate=datestr(start_day+8i,'yyyy.mm.dd');
for j =0:7
myc=(start_day+8i)-j;
day1=datestr( myc,'yyyy.mm.dd');
year1=datestr( myc,'yyyy');
info=dir([mypath '/' year1 '/AIRS.' day1 '.hdf']);
if ~isempty(info)
myfile=info.name;
myfilepath = [ mypath '/' year1 '/' myfile ];
days=[days; myfilepath];
end
end
if ~isempty(days)
datestr( start_day+8i,'yyyy-mm-dd')
C={};
days_size=length(days(:,1));
for j = 1:days_size
TotCH4_D = hdfread( days(j,:),'/descending/Data Fields/TotCH4_D', 'Index', {[1 1],[1 1],[180 360]});
C{j} = TotCH4_D;
end
D=zeros(180,360);
for j=1:180
for k=1:360
line=[];
for m=1:days_size
line=[line C{m}(j,k)];
end
line(find(line==-9999))=nan;
D(j,k)=nanmean(line);
end
end
%D(isnan(D))=-9999;
fid = fopen(['D:/AIRS/DATA/CH4_CO/TotCH4_D.' mydate '.dat' ],'w');
count = fwrite(fid, D, 'single');
fclose(fid);
%days
end
end