source: trunk/src/mode_sahelien/composite_pluie_ird_sahel.m @ 23

Last change on this file since 23 was 23, checked in by pinsard, 16 years ago

replace fcommands save, load and print by function calls
File size: 19.8 KB
Line 
1%
2
3%%%%%%%%%%%%%%% FICHIER pluie IRD
4% ouverture / lecture du fichier
5clear;
6close all;
7
8% on selectionne les mois d ete (attention aux annees bixestiles!)
9
10siz1=31+28+31+30+31+1;
11siz1bis=31+28+31+30+31+30+31+31+30;
12pluie=zeros(41,17,122,12);
13
14clear v fid; fid=fopen('pluvio.79.haut','r','b'); siz=41*17*365;
15v=fread(fid,siz,'float'); fclose(fid);
16v=reshape(v,41,17,365);
17clear a b c;
18for a=1:41;
19   for b=1:17;
20       for c=1:365;
21          if v(a,b,c)>990;
22              v(a,b,c)=NaN;
23          elseif v(a,b,c)<0;
24              v(a,b,c)=NaN;
25          end;
26       end;
27   end;
28end;
29pluie(:,:,:,1)=v(:,:,siz1:siz1bis); % 1979
30
31clear v fid; fid=fopen('pluvio.80.haut','r','b'); siz=41*17*366;
32v=fread(fid,siz,'float'); fclose(fid);
33v=reshape(v,41,17,366);
34clear a b c;
35for a=1:41;
36   for b=1:17;
37       for c=1:366;
38          if v(a,b,c)>990;
39              v(a,b,c)=NaN;
40          elseif v(a,b,c)<0;
41              v(a,b,c)=NaN;
42          end;
43       end;
44   end;
45end;
46pluie(:,:,:,2)=v(:,:,siz1+1:siz1bis+1); % 1980
47
48clear v fid; fid=fopen('pluvio.81.haut','r','b'); siz=41*17*365;
49v=fread(fid,siz,'float'); fclose(fid);
50v=reshape(v,41,17,365);
51clear a b c;
52for a=1:41;
53   for b=1:17;
54       for c=1:365;
55          if v(a,b,c)>990;
56              v(a,b,c)=NaN;
57          elseif v(a,b,c)<0;
58              v(a,b,c)=NaN;
59          end;
60       end;
61   end;
62end;
63pluie(:,:,:,3)=v(:,:,siz1:siz1bis); % 1981
64
65clear v fid; fid=fopen('pluvio.82.haut','r','b'); siz=41*17*365;
66v=fread(fid,siz,'float'); fclose(fid);
67v=reshape(v,41,17,365);
68clear a b c;
69for a=1:41;
70   for b=1:17;
71       for c=1:365;
72          if v(a,b,c)>990;
73              v(a,b,c)=NaN;
74          elseif v(a,b,c)<0;
75              v(a,b,c)=NaN;
76          end;
77       end;
78   end;
79end;
80pluie(:,:,:,4)=v(:,:,siz1:siz1bis); % 1982
81
82clear v fid; fid=fopen('pluvio.83.haut','r','b'); siz=41*17*365;
83v=fread(fid,siz,'float'); fclose(fid);
84v=reshape(v,41,17,365);
85clear a b c;
86for a=1:41;
87   for b=1:17;
88       for c=1:365;
89          if v(a,b,c)>990;
90              v(a,b,c)=NaN;
91          elseif v(a,b,c)<0;
92              v(a,b,c)=NaN;
93          end;
94       end;
95   end;
96end;
97
98pluie(:,:,:,5)=v(:,:,siz1:siz1bis); % 1983
99
100clear v fid; fid=fopen('pluvio.84.haut','r','b'); siz=41*17*366;
101v=fread(fid,siz,'float'); fclose(fid);
102v=reshape(v,41,17,366);
103clear a b c;
104for a=1:41;
105   for b=1:17;
106       for c=1:366;
107          if v(a,b,c)>990;
108              v(a,b,c)=NaN;
109          elseif v(a,b,c)<0;
110              v(a,b,c)=NaN;
111          end;
112       end;
113   end;
114end;
115pluie(:,:,:,6)=v(:,:,siz1+1:siz1bis+1); % 1984
116
117clear v fid; fid=fopen('pluvio.85.haut','r','b'); siz=41*17*365;
118v=fread(fid,siz,'float'); fclose(fid);
119v=reshape(v,41,17,365);
120clear a b c;
121for a=1:41;
122   for b=1:17;
123       for c=1:365;
124          if v(a,b,c)>990;
125              v(a,b,c)=NaN;
126          elseif v(a,b,c)<0;
127              v(a,b,c)=NaN;
128          end;
129       end;
130   end;
131end;
132pluie(:,:,:,7)=v(:,:,siz1:siz1bis); % 1985
133
134clear v fid; fid=fopen('pluvio.86.haut','r','b'); siz=41*17*365;
135v=fread(fid,siz,'float'); fclose(fid);
136v=reshape(v,41,17,365);
137clear a b c;
138for a=1:41;
139   for b=1:17;
140       for c=1:365;
141          if v(a,b,c)>990;
142              v(a,b,c)=NaN;
143          elseif v(a,b,c)<0;
144              v(a,b,c)=NaN;
145          end;
146       end;
147   end;
148end;
149pluie(:,:,:,8)=v(:,:,siz1:siz1bis); % 1986
150
151clear v fid; fid=fopen('pluvio.87.haut','r','b'); siz=41*17*365;
152v=fread(fid,siz,'float'); fclose(fid);
153v=reshape(v,41,17,365);
154clear a b c;
155for a=1:41;
156   for b=1:17;
157       for c=1:365;
158          if v(a,b,c)>990;
159              v(a,b,c)=NaN;
160          elseif v(a,b,c)<0;
161              v(a,b,c)=NaN;
162          end;
163       end;
164   end;
165end;
166pluie(:,:,:,9)=v(:,:,siz1:siz1bis); % 1987
167
168clear v fid; fid=fopen('pluvio.88.haut','r','b'); siz=41*17*366;
169v=fread(fid,siz,'float'); fclose(fid);
170v=reshape(v,41,17,366);
171clear a b c;
172for a=1:41;
173   for b=1:17;
174       for c=1:366;
175          if v(a,b,c)>990;
176              v(a,b,c)=NaN;
177          elseif v(a,b,c)<0;
178              v(a,b,c)=NaN;
179          end;
180       end;
181   end;
182end;
183pluie(:,:,:,10)=v(:,:,siz1:siz1bis); % 1988
184
185clear v fid; fid=fopen('pluvio.89.haut','r','b'); siz=41*17*365;
186v=fread(fid,siz,'float'); fclose(fid);
187v=reshape(v,41,17,365);
188clear a b c;
189for a=1:41;
190   for b=1:17;
191       for c=1:365;
192          if v(a,b,c)>990;
193              v(a,b,c)=NaN;
194          elseif v(a,b,c)<0;
195              v(a,b,c)=NaN;
196          end;
197       end;
198   end;
199end;
200pluie(:,:,:,11)=v(:,:,siz1:siz1bis); % 1989
201
202clear v fid; fid=fopen('pluvio.90.haut','r','b'); siz=41*17*365;
203v=fread(fid,siz,'float'); fclose(fid);
204v=reshape(v,41,17,365);
205clear a b c;
206for a=1:41;
207   for b=1:17;
208       for c=1:365;
209          if v(a,b,c)>990;
210              v(a,b,c)=NaN;
211          elseif v(a,b,c)<0;
212              v(a,b,c)=NaN;
213          end;
214       end;
215   end;
216end;
217pluie(:,:,:,12)=v(:,:,siz1:siz1bis); % 1990
218
219clear siz sizbis;
220clear v;
221
222status=save('pluie.mat','pluie;');
223clear
224
225status=load('pluie.mat');
226% pluie dimension : lon lat jours annee (41 17 122 12) => jour annee lat lon  (17 41 122 12)
227
228status=load('eof24_1030.txt');
229
230% on met sous la forme annee mois
231ind_olr=reshape(eof24_1030,122,28);
232clear v;
233
234% on selectionne les annees de 1979 a 2000
235ind_olr1=ind_olr(:,1:12);
236clear ind_olr;
237k=122*12;
238ind_olr=ind_olr1;
239clear ind_olr1;
240
241% on calcule l ecart type sur la serie entiere de l OLR
242ind_olr1=reshape(ind_olr, k,1);
243ind_olr2=std(ind_olr1);
244clear ind_olr1;
245
246% intialisation des tableaux
247pluie_filtre_compomax=NaN*ones(30,31,17,41);
248pluie_filtre_compomin=NaN*ones(30,31,17,41);
249pluie=permute(pluie,[3 4 2 1]);
250compteurmax=0;
251compteurmin=0;
252
253% boucle sur les annees
254for a=1:12;
255   clear delta;
256   delta=zeros(122,1);
257   % boucle sur les jours : on prepare avant
258   for b=1:122;
259        if b>1;
260            delta(b)=ind_olr(b,a)-ind_olr(b-1,a);
261        end;
262   end;
263   % boucle sur les jours
264   for b=1:122;
265      % 1ere condition (laisser une marge de 10)
266      if (b>20 && b<112);
267% cas 1 : max deltas apres negetifs et deltas avant postifs
268         if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && ind_olr(b,a)>ind_olr2);
269            compteurmax=compteurmax+1;
270            pluie_filtre_compomax(compteurmax,1,:,:)=pluie(b-20,a,:,:);
271            pluie_filtre_compomax(compteurmax,2,:,:)=pluie(b-19,a,:,:);
272            pluie_filtre_compomax(compteurmax,3,:,:)=pluie(b-18,a,:,:);
273            pluie_filtre_compomax(compteurmax,4,:,:)=pluie(b-17,a,:,:);
274            pluie_filtre_compomax(compteurmax,5,:,:)=pluie(b-16,a,:,:);
275            pluie_filtre_compomax(compteurmax,6,:,:)=pluie(b-15,a,:,:);
276            pluie_filtre_compomax(compteurmax,7,:,:)=pluie(b-14,a,:,:);
277            pluie_filtre_compomax(compteurmax,8,:,:)=pluie(b-13,a,:,:);
278            pluie_filtre_compomax(compteurmax,9,:,:)=pluie(b-12,a,:,:);
279            pluie_filtre_compomax(compteurmax,10,:,:)=pluie(b-11,a,:,:);
280            pluie_filtre_compomax(compteurmax,11,:,:)=pluie(b-10,a,:,:);
281            pluie_filtre_compomax(compteurmax,12,:,:)=pluie(b-9,a,:,:);
282            pluie_filtre_compomax(compteurmax,13,:,:)=pluie(b-8,a,:,:);
283            pluie_filtre_compomax(compteurmax,14,:,:)=pluie(b-7,a,:,:);
284            pluie_filtre_compomax(compteurmax,15,:,:)=pluie(b-6,a,:,:);
285            pluie_filtre_compomax(compteurmax,16,:,:)=pluie(b-5,a,:,:);
286            pluie_filtre_compomax(compteurmax,17,:,:)=pluie(b-4,a,:,:);
287            pluie_filtre_compomax(compteurmax,18,:,:)=pluie(b-3,a,:,:);
288            pluie_filtre_compomax(compteurmax,19,:,:)=pluie(b-2,a,:,:);
289            pluie_filtre_compomax(compteurmax,20,:,:)=pluie(b-1,a,:,:);
290            pluie_filtre_compomax(compteurmax,21,:,:)=pluie(b,a,:,:);
291            pluie_filtre_compomax(compteurmax,22,:,:)=pluie(b+1,a,:,:);
292            pluie_filtre_compomax(compteurmax,23,:,:)=pluie(b+2,a,:,:);
293            pluie_filtre_compomax(compteurmax,24,:,:)=pluie(b+3,a,:,:);
294            pluie_filtre_compomax(compteurmax,25,:,:)=pluie(b+4,a,:,:);
295            pluie_filtre_compomax(compteurmax,26,:,:)=pluie(b+5,a,:,:);
296            pluie_filtre_compomax(compteurmax,27,:,:)=pluie(b+6,a,:,:);
297            pluie_filtre_compomax(compteurmax,28,:,:)=pluie(b+7,a,:,:);
298            pluie_filtre_compomax(compteurmax,29,:,:)=pluie(b+8,a,:,:);
299            pluie_filtre_compomax(compteurmax,30,:,:)=pluie(b+9,a,:,:);
300            pluie_filtre_compomax(compteurmax,31,:,:)=pluie(b+10,a,:,:);
301% cas 2 : min deltas apres positifs et deltas avant negatifs
302         elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && ind_olr(b,a)<(-ind_olr2));
303            compteurmin=compteurmin+1;
304
305            pluie_filtre_compomin(compteurmin,1,:,:)=pluie(b-20,a,:,:);
306            pluie_filtre_compomin(compteurmin,2,:,:)=pluie(b-19,a,:,:);
307            pluie_filtre_compomin(compteurmin,3,:,:)=pluie(b-18,a,:,:);
308            pluie_filtre_compomin(compteurmin,4,:,:)=pluie(b-17,a,:,:);
309            pluie_filtre_compomin(compteurmin,5,:,:)=pluie(b-16,a,:,:);
310            pluie_filtre_compomin(compteurmin,6,:,:)=pluie(b-15,a,:,:);
311            pluie_filtre_compomin(compteurmin,7,:,:)=pluie(b-14,a,:,:);
312            pluie_filtre_compomin(compteurmin,8,:,:)=pluie(b-13,a,:,:);
313            pluie_filtre_compomin(compteurmin,9,:,:)=pluie(b-12,a,:,:);
314            pluie_filtre_compomin(compteurmin,10,:,:)=pluie(b-11,a,:,:);
315            pluie_filtre_compomin(compteurmin,11,:,:)=pluie(b-10,a,:,:);
316            pluie_filtre_compomin(compteurmin,12,:,:)=pluie(b-9,a,:,:);
317            pluie_filtre_compomin(compteurmin,13,:,:)=pluie(b-8,a,:,:);
318            pluie_filtre_compomin(compteurmin,14,:,:)=pluie(b-7,a,:,:);
319            pluie_filtre_compomin(compteurmin,15,:,:)=pluie(b-6,a,:,:);
320            pluie_filtre_compomin(compteurmin,16,:,:)=pluie(b-5,a,:,:);
321            pluie_filtre_compomin(compteurmin,17,:,:)=pluie(b-4,a,:,:);
322            pluie_filtre_compomin(compteurmin,18,:,:)=pluie(b-3,a,:,:);
323            pluie_filtre_compomin(compteurmin,19,:,:)=pluie(b-2,a,:,:);
324            pluie_filtre_compomin(compteurmin,20,:,:)=pluie(b-1,a,:,:);
325            pluie_filtre_compomin(compteurmin,21,:,:)=pluie(b,a,:,:);
326            pluie_filtre_compomin(compteurmin,22,:,:)=pluie(b+1,a,:,:);
327            pluie_filtre_compomin(compteurmin,23,:,:)=pluie(b+2,a,:,:);
328            pluie_filtre_compomin(compteurmin,24,:,:)=pluie(b+3,a,:,:);
329            pluie_filtre_compomin(compteurmin,25,:,:)=pluie(b+4,a,:,:);
330            pluie_filtre_compomin(compteurmin,26,:,:)=pluie(b+5,a,:,:);
331            pluie_filtre_compomin(compteurmin,27,:,:)=pluie(b+6,a,:,:);
332            pluie_filtre_compomin(compteurmin,28,:,:)=pluie(b+7,a,:,:);
333            pluie_filtre_compomin(compteurmin,29,:,:)=pluie(b+8,a,:,:);
334            pluie_filtre_compomin(compteurmin,30,:,:)=pluie(b+9,a,:,:);
335            pluie_filtre_compomin(compteurmin,31,:,:)=pluie(b+10,a,:,:);
336         end;
337      end;
338   end;
339end;
340
341% on fait la somme des valeurs de pluie_filtreentiel sur le nbre d occurences
342% on obtent un tableau de la forme (jours-avant-apres,lon,lat)
343pluie_filtre_compomin1(:,:,:)=nanmean(pluie_filtre_compomin);
344clear pluie_filtre_compomin;
345pluie_filtre_compomax1(:,:,:)=nanmean(pluie_filtre_compomax);
346clear pluie_filtre_compomax;
347
348figure(1);
349orient('landscape');
350figure(2);
351orient('landscape');
352figure(3);
353orient('landscape');
354
355for e=1:7;
356       a=7-(e-1)+1;
357       f=2*(a-1)+1;
358        clear pluie_min pluie_max;
359       pluie_min=zeros(17,41);
360       pluie_max=zeros(17,41);
361
362       pluie_min(:,:)=pluie_filtre_compomin1(f,:,:);
363       pluie_max(:,:)=pluie_filtre_compomax1(f,:,:);
364
365       for b=1:17;
366          for c=1:41;
367             if pluie_min(b,c)<0;
368                  pluie_min(b,c)=0;
369             elseif pluie_min(b,c)>10;
370                  pluie_min(b,c)=10;
371             end;
372             if pluie_max(b,c)<0;
373                  pluie_max(b,c)=0;
374             elseif pluie_max(b,c)>10;
375                  pluie_max(b,c)=10;
376             end;
377          end;
378       end
379       pluie_min(1,1)=-0.0001;
380       pluie_min(1,2)=10.0001;
381       pluie_max(1,1)=-0.0001;
382       pluie_max(1,2)=10.0001;
383       lat=[-10:2.5:30];
384       lon=[-70:2.5:30]';
385       figure(1);
386       subplot(7,1,8-e);
387       palette
388       aa=contourf(lon,lat,pluie_min,[0:0.5:10],'LineStyle','none');
389       if e==1;
390         set (gca,'XTick',[-180:60:180],'XTickLabel',[' 180';'120W';' 60W';'   0';' 60E';'120E';' 180'],'fontname','Arial','fontsize',6);
391       else
392         set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6);
393       end;
394       set (gca,'YTick',[-10:10:40],'YTickLabel',['   ';'  0';'10N';'20N';'30N';'   '],'fontname','Arial','fontsize',6);
395       if e==4;
396          co=colorbar;
397          set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6);
398       end;
399       hold on;
400       cartemonde1;
401       set (gca,'tickDir','out');
402       axis equal;
403       axis([-180 180 -10 40]);
404       if e==7;
405           title('Pluie IRD : composite min (mode sahelien)','fontsize',8);
406       end;
407       text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10);
408       pos_vert=e/10;
409       set (gca,'position',[0.2 pos_vert 0.5 0.1]);
410%         ***************************
411       figure(2);
412       subplot(7,1,8-e);
413       palette
414       aa=contourf(lon,lat,pluie_max,[0:0.5:10],'LineStyle','none');
415       if e==1;
416         set (gca,'XTick',[-180:60:180],'XTickLabel',[' 180';'120W';' 60W';'   0';' 60E';'120E';' 180'],'fontname','Arial','fontsize',6);
417       else
418         set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6);
419       end;
420       set (gca,'YTick',[-10:10:40],'YTickLabel',['   ';'  0';'10N';'20N';'30N';'   '],'fontname','Arial','fontsize',6);
421       if e==4;
422          co=colorbar;
423          set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6);
424       end;
425       hold on;
426       cartemonde1;
427       set (gca,'tickDir','out');
428       axis equal;
429       axis([-180 180 -10 40]);
430       if e==7;
431           title('Pluies IRD : composite max (mode Sahelien)','fontsize',8);
432       end;
433       text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10);
434       pos_vert=e/10;
435       set (gca,'position',[0.2 pos_vert 0.5 0.1]);
436%********************************************
437
438       figure(3);
439       subplot(7,1,8-e);
440       palette
441       diff=pluie_min-pluie_max;
442       diff(1,1)=-5.01;
443       diff(1,2)=5.01;
444       aa=contourf(lon,lat,diff,[-5:0.5:5],'LineStyle','none');
445       if e==1;
446         set (gca,'XTick',[-180:60:180],'XTickLabel',[' 180';'120W';' 60W';'   0';' 60E';'120E';' 180'],'fontname','Arial','fontsize',6);
447       else
448         set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6);
449       end;
450       set (gca,'YTick',[-10:10:40],'YTickLabel',['   ';'  0';'10N';'20N';'30N';'   '],'fontname','Arial','fontsize',6);
451       if e==4;
452          co=colorbar;
453          set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6);
454       end;
455       hold on;
456       cartemonde1;
457       set (gca,'tickDir','out');
458       axis equal;
459       axis([-180 180 -10 40]);
460       if e==7;
461           title('Pluie IRD : differences composite min - max (mode Sahelien)','fontsize',8);
462       end;
463       text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10);
464       pos_vert=e/10;
465       set (gca,'position',[0.2 pos_vert 0.5 0.1]);
466
467end;
468figure(1);
469status=print('-depsc2','composite_pluie_min_mode_sahel1.eps');
470figure(2);
471status=print('-depsc2','composite_pluie_max_mode_sahel1.eps');
472figure(3);
473status=print('-depsc2','composite_pluie_diff_mode_sahel1.eps');
474close all
475
476figure(1);
477orient('landscape');
478figure(2);
479orient('landscape');
480figure(3);
481orient('landscape');
482
483for e=1:7;
484       a=7-(e-1)+7;
485       f=2*(a-1)+1;
486        clear pluie_min pluie_max;
487       pluie_min=zeros(17,41);
488       pluie_max=zeros(17,41);
489
490       pluie_min(:,:)=pluie_filtre_compomin1(f,:,:);
491       pluie_max(:,:)=pluie_filtre_compomax1(f,:,:);
492
493       for b=1:17;
494          for c=1:41;
495             if pluie_min(b,c)<0;
496                  pluie_min(b,c)=0;
497             elseif pluie_min(b,c)>10;
498                  pluie_min(b,c)=10;
499             end;
500             if pluie_max(b,c)<0;
501                  pluie_max(b,c)=0;
502             elseif pluie_max(b,c)>10;
503                  pluie_max(b,c)=10;
504             end;
505          end;
506       end
507       pluie_min(1,1)=-0.0001;
508       pluie_min(1,2)=10.0001;
509       pluie_max(1,1)=-0.0001;
510       pluie_max(1,2)=10.0001;
511       lat=[-10:2.5:30];
512       lon=[-70:2.5:30]';
513       figure(1);
514       subplot(7,1,8-e);
515       palette
516       aa=contourf(lon,lat,pluie_min,[0:0.5:10],'LineStyle','none');
517       if e==1;
518         set (gca,'XTick',[-180:60:180],'XTickLabel',[' 180';'120W';' 60W';'   0';' 60E';'120E';' 180'],'fontname','Arial','fontsize',6);
519       else
520         set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6);
521       end;
522       set (gca,'YTick',[-10:10:40],'YTickLabel',['   ';'  0';'10N';'20N';'30N';'   '],'fontname','Arial','fontsize',6);
523       if e==4;
524          co=colorbar;
525          set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6);
526       end;
527       hold on;
528       cartemonde1;
529       set (gca,'tickDir','out');
530       axis equal;
531       axis([-180 180 -10 40]);
532       if e==7;
533           title('Pluie IRD : composite min (mode Sahelien)','fontsize',8);
534       end;
535       text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10);
536       pos_vert=e/10;
537       set (gca,'position',[0.2 pos_vert 0.5 0.1]);
538%         ***************************
539       figure(2);
540       subplot(7,1,8-e);
541       palette
542       aa=contourf(lon,lat,pluie_max,[0:0.5:10],'LineStyle','none');
543       if e==1;
544         set (gca,'XTick',[-180:60:180],'XTickLabel',[' 180';'120W';' 60W';'   0';' 60E';'120E';' 180'],'fontname','Arial','fontsize',6);
545       else
546         set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6);
547       end;
548       set (gca,'YTick',[-10:10:40],'YTickLabel',['   ';'  0';'10N';'20N';'30N';'   '],'fontname','Arial','fontsize',6);
549       if e==4;
550          co=colorbar;
551          set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6);
552       end;
553       hold on;
554       cartemonde1;
555       set (gca,'tickDir','out');
556       axis equal;
557       axis([-180 180 -10 40]);
558       if e==7;
559           title('Pluie IRD : composite max (mode Sahelien)','fontsize',8);
560       end;
561       text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10);
562       pos_vert=e/10;
563       set (gca,'position',[0.2 pos_vert 0.5 0.1]);
564%********************************************
565
566       figure(3);
567       subplot(7,1,8-e);
568       palette
569       diff=pluie_min-pluie_max;
570       diff(1,1)=-5.01;
571       diff(1,2)=5.01;
572       aa=contourf(lon,lat,diff,[-5:0.5:5],'LineStyle','none');
573       if e==1;
574         set (gca,'XTick',[-180:60:180],'XTickLabel',[' 180';'120W';' 60W';'   0';' 60E';'120E';' 180'],'fontname','Arial','fontsize',6);
575       else
576         set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6);
577       end;
578       set (gca,'YTick',[-10:10:40],'YTickLabel',['   ';'  0';'10N';'20N';'30N';'   '],'fontname','Arial','fontsize',6);
579       if e==4;
580          co=colorbar;
581          set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6);
582       end;
583       hold on;
584       cartemonde1;
585       set (gca,'tickDir','out');
586       axis equal;
587       axis([-180 180 -10 40]);
588       if e==7;
589           title('Pluie IRD : diffrences composites min -max (mode Sahelien)','fontsize',8);
590       end;
591       text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10);
592       pos_vert=e/10;
593       set (gca,'position',[0.2 pos_vert 0.5 0.1]);
594
595end;
596figure(1);
597status=print('-depsc2','composite_pluie_min_mode_sahel2.eps');
598figure(2);
599status=print('-depsc2','composite_pluie_max_mode_sahel2.eps');
600figure(3);
601status=print('-depsc2','composite_pluie_diff_mode_sahel2.eps');
Note: See TracBrowser for help on using the repository browser.