atmos-reanalysis/jra-25/monthly.m


for yr=2013
        cyr=num2str(yr);
ndays=[31 28 31 30 31 30 31 31 30 31 30 31 ];
if(mod(yr,4)==0);ndays(2)=29;end
        [yr ndays]
t1=dir(['/net/nares/raid11/ecco-shared/atmos-reanalysis/jra-25/J*' cyr]);
for i=1:length(t1)
  fn=['/bay/scratch2/JRA25/' t1(i).name]
  fid=fopen(fn,'r','b');
  fo=['/bay/scratch2/JRA25/' t1(i).name '_monav'];
  fout=fopen(fo,'w','b');
  for nmo=1:12
          nrecs=4*ndays(nmo);
          JRA25g=fread(fid,[320*160 nrecs],'float32');
          JRA25gm=mean(JRA25g');
          fwrite(fout,JRA25gm,'float32');
  end
  fclose(fid);
  fclose(fout);
end  %  for length(t1);
end  % for yr


1	gforget	1.1
2			for yr=2013
3			cyr=num2str(yr);
4			ndays=[31 28 31 30 31 30 31 31 30 31 30 31 ];
5			if(mod(yr,4)==0);ndays(2)=29;end
6			[yr ndays]
7			t1=dir(['/net/nares/raid11/ecco-shared/atmos-reanalysis/jra-25/J*' cyr]);
8			for i=1:length(t1)
9			fn=['/bay/scratch2/JRA25/' t1(i).name]
10			fid=fopen(fn,'r','b');
11			fo=['/bay/scratch2/JRA25/' t1(i).name '_monav'];
12			fout=fopen(fo,'w','b');
13			for nmo=1:12
14			nrecs=4*ndays(nmo);
15			JRA25g=fread(fid,[320*160 nrecs],'float32');
16			JRA25gm=mean(JRA25g');
17			fwrite(fout,JRA25gm,'float32');
18			end
19			fclose(fid);
20			fclose(fout);
21			end % for length(t1);
22			end % for yr
23
24