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