gael/bulkMatlab/plot_bulk.m


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);

1	gforget	1.1
2			domaine_global_def; domaine;
3
4			nb_days=365;
5			nb_days=13;
6
7			ustress=zeros(jpi,jpj,nb_days);
8			vstress=zeros(jpi,jpj,nb_days);
9			qnet=zeros(jpi,jpj,nb_days);
10			empmr=zeros(jpi,jpj,nb_days);
11
12			ycur=2004;
13			%fcur=1;
14
15			switch fcur;
16			case 1; sol_cur=['.ncep.flux.' num2str(ycur) '.daily.bin'];
17			case 2; sol_cur=['.ecmwf.flux.' num2str(ycur) '.daily.bin'];
18			case 3; sol_cur=['.ncep_noicemask.flux.' num2str(ycur) '.daily.bin'];
19			case 4; sol_cur=['.ecmwf_noicemask.flux.' num2str(ycur) '.daily.bin'];
20			case 5; sol_cur=['.ncep.LY.' num2str(ycur) '.daily.bin'];
21			case 6; sol_cur=['.ecmwf.LY.' num2str(ycur) '.daily.bin'];
22			case 7; sol_cur=['.core.LY.' num2str(ycur) '.daily.bin'];
23			case 8; sol_cur=['.ncep.coare.' num2str(ycur) '.daily.bin'];
24			case 9; sol_cur=['.ncep.LY4ncep.' num2str(ycur) '.daily.bin'];
25			end;
26
27			rep_in='./';
28
29			fid_ustress=fopen([rep_in '/ustress' sol_cur ],'r','b');
30			fid_vstress=fopen([rep_in '/vstress' sol_cur ],'r','b');
31			fid_qnet=fopen([rep_in '/qnet' sol_cur ],'r','b');
32			fid_empmr=fopen([rep_in '/empmr' sol_cur ],'r','b');
33
34			ustress=ustress+reshape(fread(fid_ustress,jpijpjnb_days,'float32'),jpi,jpj,nb_days);
35			vstress=vstress+reshape(fread(fid_vstress,jpijpjnb_days,'float32'),jpi,jpj,nb_days);
36			qnet=qnet+reshape(fread(fid_qnet,jpijpjnb_days,'float32'),jpi,jpj,nb_days);
37			empmr=empmr+reshape(fread(fid_empmr,jpijpjnb_days,'float32'),jpi,jpj,nb_days);
38
39			fclose('all');
40
41
42			ustress_mean=nanmean(ustress,3); vstress_mean=nanmean(vstress,3);
43			qnet_mean=nanmean(qnet,3); empmr_mean=nanmean(empmr,3);
44
45
46			figure;
47			subplot(2,2,1); pcolor(lon2D_t,lat2D_t,ustress_mean); shading flat; caxis([-1 1]*0.2); colorbar;
48			subplot(2,2,2); pcolor(lon2D_t,lat2D_t,vstress_mean); shading flat; caxis([-1 1]*0.1); colorbar;
49			subplot(2,2,3); pcolor(lon2D_t,lat2D_t,-qnet_mean); shading flat; caxis([-1 1]*200); colorbar;
50			subplot(2,2,4); pcolor(lon2D_t,lat2D_t,-empmr_mean); shading flat; caxis([-1 1]*1e-4); colorbar;
51
52			pause(2);
53