function [res snrs]=read_res(fich) resA=zeros(6,2,30,10); cdts=zeros(6,2,30,10); fid = fopen('missing_qsub.txt','w'); for mint=1:3600 [snr pen t wv]=ind2sub([6 2 30 1],mint); m=int2str(mint); try load(strcat(fich,m)); resA(snr,pen,t,wv)=SRRA(Aind,hatA); cdts(snr,pen,t,wv)=niveau_bruit_db; catch a_str=['qsub -cwd -l h_rss=10G -t ' int2str(mint) ' gretsi_fig3.sh'] fprintf(fid, '%s\n', a_str); resA(snr,pen,t)=0; end end save(['res-' fich],'resA','cdts','Ts'); fclose(fid); end