domaine_global_def; domaine;

nb_days=365;
nb_days=13;

ustress=zeros(jpi,jpj,nb_days); 
vstress=zeros(jpi,jpj,nb_days);
qnet=zeros(jpi,jpj,nb_days); 
empmr=zeros(jpi,jpj,nb_days);

ycur=2004;
%fcur=1;

switch fcur;
case 1; sol_cur=['.ncep.flux.' num2str(ycur) '.daily.bin'];
case 2; sol_cur=['.ecmwf.flux.' num2str(ycur) '.daily.bin'];
case 3; sol_cur=['.ncep_noicemask.flux.' num2str(ycur) '.daily.bin'];
case 4; sol_cur=['.ecmwf_noicemask.flux.' num2str(ycur) '.daily.bin'];
case 5; sol_cur=['.ncep.LY.' num2str(ycur) '.daily.bin'];
case 6; sol_cur=['.ecmwf.LY.' num2str(ycur) '.daily.bin'];
case 7; sol_cur=['.core.LY.' num2str(ycur) '.daily.bin'];
case 8; sol_cur=['.ncep.coare.' num2str(ycur) '.daily.bin'];
case 9; sol_cur=['.ncep.LY4ncep.' num2str(ycur) '.daily.bin'];
end;

rep_in='./';

fid_ustress=fopen([rep_in '/ustress' sol_cur ],'r','b'); 
fid_vstress=fopen([rep_in '/vstress' sol_cur ],'r','b');
fid_qnet=fopen([rep_in '/qnet' sol_cur ],'r','b'); 
fid_empmr=fopen([rep_in '/empmr' sol_cur ],'r','b');

ustress=ustress+reshape(fread(fid_ustress,jpi*jpj*nb_days,'float32'),jpi,jpj,nb_days);
vstress=vstress+reshape(fread(fid_vstress,jpi*jpj*nb_days,'float32'),jpi,jpj,nb_days);
qnet=qnet+reshape(fread(fid_qnet,jpi*jpj*nb_days,'float32'),jpi,jpj,nb_days);
empmr=empmr+reshape(fread(fid_empmr,jpi*jpj*nb_days,'float32'),jpi,jpj,nb_days);

fclose('all');


ustress_mean=nanmean(ustress,3); vstress_mean=nanmean(vstress,3);
qnet_mean=nanmean(qnet,3); empmr_mean=nanmean(empmr,3);


figure;
subplot(2,2,1); pcolor(lon2D_t,lat2D_t,ustress_mean); shading flat; caxis([-1 1]*0.2); colorbar;
subplot(2,2,2); pcolor(lon2D_t,lat2D_t,vstress_mean); shading flat; caxis([-1 1]*0.1); colorbar;
subplot(2,2,3); pcolor(lon2D_t,lat2D_t,-qnet_mean); shading flat; caxis([-1 1]*200); colorbar;
subplot(2,2,4); pcolor(lon2D_t,lat2D_t,-empmr_mean); shading flat; caxis([-1 1]*1e-4); colorbar;

pause(2);