for icur=6:6 switch icur case 1 repr='/raid2/king/data_1x1_92-03/forcing_ECMWF/'; list_files=dir([repr '*19*']); case 2 repr='/raid2/king/data_1x1_92-03/forcing_ECMWF/'; list_files=dir([repr '*20*']); case 3 repr='/raid2/king/data_1x1_92-03/bulk_ECMWF/'; list_files=dir([repr '*19*']); case 4 repr='/raid2/king/data_1x1_92-03/bulk_ECMWF/'; list_files=dir([repr '*20*']); case 5 repr='/net/ross/raid2/king/ERA40_ECMWF_gau6_78-06/'; list_files=dir([repr '*20*']); case 6 repr='/net/ross/raid2/king/ERA40_ECMWF_gau6_78-06/'; list_files=dir([repr '*19*']); end stat_files=zeros(length(list_files),3); for fcur=1:length(list_files) file_cur=deblank(list_files(fcur).name); fid=fopen([repr file_cur],'r','b'); tmp4=list_files(fcur).bytes/4; tmp1=fread(fid,tmp4,'float32'); stat_files(fcur,1)=sum(~isfinite(tmp1)); stat_files(fcur,2)=sum(tmp1~=0); stat_files(fcur,3)=sum(tmp1==0); if stat_files(fcur,1)~=0 fprintf(['\n ' file_cur ' ' num2str(stat_files(fcur,:)) '\n']); else fprintf('x'); end fclose(fid); clear tmp1; end eval(['list_files' num2str(icur) '=list_files;']); eval(['stat_files' num2str(icur) '=stat_files;']); end%for icur=1:4