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