source: trunk/src/mode_sahelien/composite_wind925_sahel_1090.m @ 37

Last change on this file since 37 was 26, checked in by pinsard, 15 years ago

add semi-colon after end everywhere

File size: 21.2 KB
Line 
1%
2
3% initialisation
4close all
5clear;
6
7% etape 1 :on selectionne les mois d ete (attention aux annees bixestiles!)
8%          on compile sous forme d'un seul fichier
9
10uwnd700=zeros(122,49,144,22);
11
12siz=31+28+31+30+31+1;
13sizbis=31+28+31+30+31+30+31+31+30;
14ncload('uwnd700.1979.nc'); uwnd700(:,:,:,1)=uwnd(siz:sizbis,1,:,:); clear uwnd;
15ncload('uwnd700.1980.nc'); uwnd700(:,:,:,2)=uwnd(siz+1:sizbis+1,1,:,:); clear uwnd;
16ncload('uwnd700.1981.nc'); uwnd700(:,:,:,3)=uwnd(siz:sizbis,1,:,:);clear uwnd;
17ncload('uwnd700.1982.nc'); uwnd700(:,:,:,4)=uwnd(siz:sizbis,1,:,:);clear uwnd;
18ncload('uwnd700.1983.nc'); uwnd700(:,:,:,5)=uwnd(siz:sizbis,1,:,:);clear uwnd;
19ncload('uwnd700.1984.nc'); uwnd700(:,:,:,6)=uwnd(siz+1:sizbis+1,1,:,:);clear uwnd;
20ncload('uwnd700.1985.nc'); uwnd700(:,:,:,7)=uwnd(siz:sizbis,1,:,:);clear uwnd;
21ncload('uwnd700.1986.nc'); uwnd700(:,:,:,8)=uwnd(siz:sizbis,1,:,:);clear uwnd;
22ncload('uwnd700.1987.nc'); uwnd700(:,:,:,9)=uwnd(siz:sizbis,1,:,:);clear uwnd;
23ncload('uwnd700.1988.nc'); uwnd700(:,:,:,10)=uwnd(siz+1:sizbis+1,1,:,:);clear uwnd;
24ncload('uwnd700.1989.nc'); uwnd700(:,:,:,11)=uwnd(siz:sizbis,1,:,:);clear uwnd;
25ncload('uwnd700.1990.nc'); uwnd700(:,:,:,12)=uwnd(siz:sizbis,1,:,:);clear uwnd;
26ncload('uwnd700.1991.nc'); uwnd700(:,:,:,13)=uwnd(siz:sizbis,1,:,:);clear uwnd;
27ncload('uwnd700.1992.nc'); uwnd700(:,:,:,14)=uwnd(siz+1:sizbis+1,1,:,:);clear uwnd;
28ncload('uwnd700.1993.nc'); uwnd700(:,:,:,15)=uwnd(siz:sizbis,1,:,:);clear uwnd;
29ncload('uwnd700.1994.nc'); uwnd700(:,:,:,16)=uwnd(siz:sizbis,1,:,:);clear uwnd;
30ncload('uwnd700.1995.nc'); uwnd700(:,:,:,17)=uwnd(siz:sizbis,1,:,:);clear uwnd;
31ncload('uwnd700.1996.nc'); uwnd700(:,:,:,18)=uwnd(siz+1:sizbis+1,1,:,:);clear uwnd;
32ncload('uwnd700.1997.nc'); uwnd700(:,:,:,19)=uwnd(siz:sizbis,1,:,:);clear uwnd;
33ncload('uwnd700.1998.nc'); uwnd700(:,:,:,20)=uwnd(siz:sizbis,1,:,:);clear uwnd;
34ncload('uwnd700.1999.nc'); uwnd700(:,:,:,21)=uwnd(siz:sizbis,1,:,:);clear uwnd;
35ncload('uwnd700.2000.nc'); uwnd700(:,:,:,22)=uwnd(siz+1:sizbis+1,1,:,:);clear uwnd;
36clear siz1 siz1bis;
37
38uwnd=permute(uwnd700,[1 4 2 3]);
39clear uwnd700;
40
41
42%%%%%%%%%%%%%%% FICHIER OLR
43% ouverture / lecture du fichier
44status=load('ind_olr_10_90_eof23.txt');
45ind_olr_10_90_eof23=status;
46clear status;
47ind_olr=reshape(ind_olr_10_90_eof23,122,28);
48clear ind_olr_10_90_eof23;
49
50% on selectionne les annees de 1979 a 2000
51ind_olr1=ind_olr(:,1:22);
52clear ind_olr;
53k=122*22;
54ind_olr=ind_olr1;
55clear ind_olr1;
56
57% on calcule l ecart type sur la serie entiere de l OLR
58ind_olr1=reshape(ind_olr, k,1);
59ind_olr2=std(ind_olr1);
60clear ind_olr1;
61
62
63% intialisation des tableaux
64uwnd700_filtre_compomax=NaN*ones(60,31,49,144);
65uwnd700_filtre_compomin=NaN*ones(60,31,49,144);
66
67compteurmax=0;
68compteurmin=0;
69
70% boucle sur les annees
71for a=1:22;
72   clear delta;
73   delta=zeros(122,1);
74   % boucle sur les jours : on prepare avant
75   for b=1:122;
76        if b>1;
77            delta(b)=ind_olr(b,a)-ind_olr(b-1,a);
78        end;
79   end;
80   % boucle sur les jours
81   for b=1:122;
82      % 1ere condition (laisser une marge de 10)
83      if (b>20 && b<112);
84% cas 1 : max deltas apres negetifs et deltas avant postifs
85         if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && ind_olr(b,a)>ind_olr2);
86            compteurmax=compteurmax+1;
87            uwnd700_filtre_compomax(compteurmax,1,:,:)=uwnd(b-20,a,:,:);
88            uwnd700_filtre_compomax(compteurmax,2,:,:)=uwnd(b-19,a,:,:);
89            uwnd700_filtre_compomax(compteurmax,3,:,:)=uwnd(b-18,a,:,:);
90            uwnd700_filtre_compomax(compteurmax,4,:,:)=uwnd(b-17,a,:,:);
91            uwnd700_filtre_compomax(compteurmax,5,:,:)=uwnd(b-16,a,:,:);
92            uwnd700_filtre_compomax(compteurmax,6,:,:)=uwnd(b-15,a,:,:);
93            uwnd700_filtre_compomax(compteurmax,7,:,:)=uwnd(b-14,a,:,:);
94            uwnd700_filtre_compomax(compteurmax,8,:,:)=uwnd(b-13,a,:,:);
95            uwnd700_filtre_compomax(compteurmax,9,:,:)=uwnd(b-12,a,:,:);
96            uwnd700_filtre_compomax(compteurmax,10,:,:)=uwnd(b-11,a,:,:);
97
98            uwnd700_filtre_compomax(compteurmax,11,:,:)=uwnd(b-10,a,:,:);
99            uwnd700_filtre_compomax(compteurmax,12,:,:)=uwnd(b-9,a,:,:);
100            uwnd700_filtre_compomax(compteurmax,13,:,:)=uwnd(b-8,a,:,:);
101            uwnd700_filtre_compomax(compteurmax,14,:,:)=uwnd(b-7,a,:,:);
102            uwnd700_filtre_compomax(compteurmax,15,:,:)=uwnd(b-6,a,:,:);
103            uwnd700_filtre_compomax(compteurmax,16,:,:)=uwnd(b-5,a,:,:);
104            uwnd700_filtre_compomax(compteurmax,17,:,:)=uwnd(b-4,a,:,:);
105            uwnd700_filtre_compomax(compteurmax,18,:,:)=uwnd(b-3,a,:,:);
106            uwnd700_filtre_compomax(compteurmax,19,:,:)=uwnd(b-2,a,:,:);
107            uwnd700_filtre_compomax(compteurmax,20,:,:)=uwnd(b-1,a,:,:);
108
109            uwnd700_filtre_compomax(compteurmax,21,:,:)=uwnd(b,a,:,:);
110            uwnd700_filtre_compomax(compteurmax,22,:,:)=uwnd(b+1,a,:,:);
111            uwnd700_filtre_compomax(compteurmax,23,:,:)=uwnd(b+2,a,:,:);
112            uwnd700_filtre_compomax(compteurmax,24,:,:)=uwnd(b+3,a,:,:);
113            uwnd700_filtre_compomax(compteurmax,25,:,:)=uwnd(b+4,a,:,:);
114            uwnd700_filtre_compomax(compteurmax,26,:,:)=uwnd(b+5,a,:,:);
115            uwnd700_filtre_compomax(compteurmax,27,:,:)=uwnd(b+6,a,:,:);
116            uwnd700_filtre_compomax(compteurmax,28,:,:)=uwnd(b+7,a,:,:);
117            uwnd700_filtre_compomax(compteurmax,29,:,:)=uwnd(b+8,a,:,:);
118            uwnd700_filtre_compomax(compteurmax,30,:,:)=uwnd(b+9,a,:,:);
119            uwnd700_filtre_compomax(compteurmax,31,:,:)=uwnd(b+10,a,:,:);
120% cas 2 : min deltas apres positifs et deltas avant negatifs
121         elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && ind_olr(b,a)<(-ind_olr2));
122            compteurmin=compteurmin+1;
123
124            uwnd700_filtre_compomin(compteurmin,1,:,:)=uwnd(b-20,a,:,:);
125            uwnd700_filtre_compomin(compteurmin,2,:,:)=uwnd(b-19,a,:,:);
126            uwnd700_filtre_compomin(compteurmin,3,:,:)=uwnd(b-18,a,:,:);
127            uwnd700_filtre_compomin(compteurmin,4,:,:)=uwnd(b-17,a,:,:);
128            uwnd700_filtre_compomin(compteurmin,5,:,:)=uwnd(b-16,a,:,:);
129            uwnd700_filtre_compomin(compteurmin,6,:,:)=uwnd(b-15,a,:,:);
130            uwnd700_filtre_compomin(compteurmin,7,:,:)=uwnd(b-14,a,:,:);
131            uwnd700_filtre_compomin(compteurmin,8,:,:)=uwnd(b-13,a,:,:);
132            uwnd700_filtre_compomin(compteurmin,9,:,:)=uwnd(b-12,a,:,:);
133            uwnd700_filtre_compomin(compteurmin,10,:,:)=uwnd(b-11,a,:,:);
134
135            uwnd700_filtre_compomin(compteurmin,11,:,:)=uwnd(b-10,a,:,:);
136            uwnd700_filtre_compomin(compteurmin,12,:,:)=uwnd(b-9,a,:,:);
137            uwnd700_filtre_compomin(compteurmin,13,:,:)=uwnd(b-8,a,:,:);
138            uwnd700_filtre_compomin(compteurmin,14,:,:)=uwnd(b-7,a,:,:);
139            uwnd700_filtre_compomin(compteurmin,15,:,:)=uwnd(b-6,a,:,:);
140            uwnd700_filtre_compomin(compteurmin,16,:,:)=uwnd(b-5,a,:,:);
141            uwnd700_filtre_compomin(compteurmin,17,:,:)=uwnd(b-4,a,:,:);
142            uwnd700_filtre_compomin(compteurmin,18,:,:)=uwnd(b-3,a,:,:);
143            uwnd700_filtre_compomin(compteurmin,19,:,:)=uwnd(b-2,a,:,:);
144            uwnd700_filtre_compomin(compteurmin,20,:,:)=uwnd(b-1,a,:,:);
145
146            uwnd700_filtre_compomin(compteurmin,21,:,:)=uwnd(b,a,:,:);
147            uwnd700_filtre_compomin(compteurmin,22,:,:)=uwnd(b+1,a,:,:);
148            uwnd700_filtre_compomin(compteurmin,23,:,:)=uwnd(b+2,a,:,:);
149            uwnd700_filtre_compomin(compteurmin,24,:,:)=uwnd(b+3,a,:,:);
150            uwnd700_filtre_compomin(compteurmin,25,:,:)=uwnd(b+4,a,:,:);
151            uwnd700_filtre_compomin(compteurmin,26,:,:)=uwnd(b+5,a,:,:);
152            uwnd700_filtre_compomin(compteurmin,27,:,:)=uwnd(b+6,a,:,:);
153            uwnd700_filtre_compomin(compteurmin,28,:,:)=uwnd(b+7,a,:,:);
154            uwnd700_filtre_compomin(compteurmin,29,:,:)=uwnd(b+8,a,:,:);
155            uwnd700_filtre_compomin(compteurmin,30,:,:)=uwnd(b+9,a,:,:);
156            uwnd700_filtre_compomin(compteurmin,31,:,:)=uwnd(b+10,a,:,:);
157         end;
158      end;
159   end;
160end;
161
162
163% on fait la somme des valeurs de uwnd700_filtreentiel sur le nbre d occurences
164% on obtent un tableau de la forme (jours-avant-apres,lon,lat)
165uwnd700_filtre_compomin1(:,:,:)=nanmean(uwnd700_filtre_compomin);
166clear uwnd700_filtre_compomin;
167uwnd700_filtre_compomax1(:,:,:)=nanmean(uwnd700_filtre_compomax);
168clear uwnd700_filtre_compomax;
169
170% etape 1 :on selectionne les mois d ete (attention aux annees bixestiles!)
171%          on compile sous forme d'un seul fichier
172
173vwnd700=zeros(122,49,144,22);
174
175siz=31+28+31+30+31+1;
176sizbis=31+28+31+30+31+30+31+31+30;
177ncload('vwnd700.1979.nc'); vwnd700(:,:,:,1)=vwnd(siz:sizbis,1,:,:); clear vwnd;
178ncload('vwnd700.1980.nc'); vwnd700(:,:,:,2)=vwnd(siz+1:sizbis+1,1,:,:); clear vwnd;
179ncload('vwnd700.1981.nc'); vwnd700(:,:,:,3)=vwnd(siz:sizbis,1,:,:);clear vwnd;
180ncload('vwnd700.1982.nc'); vwnd700(:,:,:,4)=vwnd(siz:sizbis,1,:,:);clear vwnd;
181ncload('vwnd700.1983.nc'); vwnd700(:,:,:,5)=vwnd(siz:sizbis,1,:,:);clear vwnd;
182ncload('vwnd700.1984.nc'); vwnd700(:,:,:,6)=vwnd(siz+1:sizbis+1,1,:,:);clear vwnd;
183ncload('vwnd700.1985.nc'); vwnd700(:,:,:,7)=vwnd(siz:sizbis,1,:,:);clear vwnd;
184ncload('vwnd700.1986.nc'); vwnd700(:,:,:,8)=vwnd(siz:sizbis,1,:,:);clear vwnd;
185ncload('vwnd700.1987.nc'); vwnd700(:,:,:,9)=vwnd(siz:sizbis,1,:,:);clear vwnd;
186ncload('vwnd700.1988.nc'); vwnd700(:,:,:,10)=vwnd(siz+1:sizbis+1,1,:,:);clear vwnd;
187ncload('vwnd700.1989.nc'); vwnd700(:,:,:,11)=vwnd(siz:sizbis,1,:,:);clear vwnd;
188ncload('vwnd700.1990.nc'); vwnd700(:,:,:,12)=vwnd(siz:sizbis,1,:,:);clear vwnd;
189ncload('vwnd700.1991.nc'); vwnd700(:,:,:,13)=vwnd(siz:sizbis,1,:,:);clear vwnd;
190ncload('vwnd700.1992.nc'); vwnd700(:,:,:,14)=vwnd(siz+1:sizbis+1,1,:,:);clear vwnd;
191ncload('vwnd700.1993.nc'); vwnd700(:,:,:,15)=vwnd(siz:sizbis,1,:,:);clear vwnd;
192ncload('vwnd700.1994.nc'); vwnd700(:,:,:,16)=vwnd(siz:sizbis,1,:,:);clear vwnd;
193ncload('vwnd700.1995.nc'); vwnd700(:,:,:,17)=vwnd(siz:sizbis,1,:,:);clear vwnd;
194ncload('vwnd700.1996.nc'); vwnd700(:,:,:,18)=vwnd(siz+1:sizbis+1,1,:,:);clear vwnd;
195ncload('vwnd700.1997.nc'); vwnd700(:,:,:,19)=vwnd(siz:sizbis,1,:,:);clear vwnd;
196ncload('vwnd700.1998.nc'); vwnd700(:,:,:,20)=vwnd(siz:sizbis,1,:,:);clear vwnd;
197ncload('vwnd700.1999.nc'); vwnd700(:,:,:,21)=vwnd(siz:sizbis,1,:,:);clear vwnd;
198ncload('vwnd700.2000.nc'); vwnd700(:,:,:,22)=vwnd(siz+1:sizbis+1,1,:,:);clear vwnd;
199clear siz1 siz1bis;
200
201vwnd=permute(vwnd700,[1 4 2 3]);
202clear vwnd700;
203
204%%%%%%%%%%%%%% FICHIER OLR
205% ouverture / lecture du fichier
206fid=fopen('sahel_pb.dat','r');
207v=fread(fid,3538,'float');
208fclose(fid);
209
210% on met sous la forme annee mois
211ind_olr=reshape(v,122,29);
212clear v;
213
214% on selectionne les annees de 1979 a 2000
215ind_olr1=ind_olr(:,1:22);
216clear ind_olr;
217k=122*22;
218ind_olr=ind_olr1;
219clear ind_olr1;
220
221% on calcule l ecart type sur la serie entiere de l OLR
222ind_olr1=reshape(ind_olr, k,1);
223ind_olr2=std(ind_olr1);
224clear ind_olr1;
225
226
227% intialisation des tableaux
228vwnd700_filtre_compomax=NaN*ones(60,31,49,144);
229vwnd700_filtre_compomin=NaN*ones(60,31,49,144);
230
231compteurmax=0;
232compteurmin=0;
233
234% boucle sur les annees
235for a=1:22;
236   clear delta;
237   delta=zeros(122,1);
238   % boucle sur les jours : on prepare avant
239   for b=1:122;
240        if b>1;
241            delta(b)=ind_olr(b,a)-ind_olr(b-1,a);
242        end;
243   end;
244   % boucle sur les jours
245   for b=1:122;
246      % 1ere condition (laisser une marge de 10)
247      if (b>20 && b<112);
248% cas 1 : max deltas apres negetifs et deltas avant postifs
249         if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && ind_olr(b,a)>ind_olr2);
250            compteurmax=compteurmax+1;
251            vwnd700_filtre_compomax(compteurmax,1,:,:)=vwnd(b-20,a,:,:);
252            vwnd700_filtre_compomax(compteurmax,2,:,:)=vwnd(b-19,a,:,:);
253            vwnd700_filtre_compomax(compteurmax,3,:,:)=vwnd(b-18,a,:,:);
254            vwnd700_filtre_compomax(compteurmax,4,:,:)=vwnd(b-17,a,:,:);
255            vwnd700_filtre_compomax(compteurmax,5,:,:)=vwnd(b-16,a,:,:);
256            vwnd700_filtre_compomax(compteurmax,6,:,:)=vwnd(b-15,a,:,:);
257            vwnd700_filtre_compomax(compteurmax,7,:,:)=vwnd(b-14,a,:,:);
258            vwnd700_filtre_compomax(compteurmax,8,:,:)=vwnd(b-13,a,:,:);
259            vwnd700_filtre_compomax(compteurmax,9,:,:)=vwnd(b-12,a,:,:);
260            vwnd700_filtre_compomax(compteurmax,10,:,:)=vwnd(b-11,a,:,:);
261
262            vwnd700_filtre_compomax(compteurmax,11,:,:)=vwnd(b-10,a,:,:);
263            vwnd700_filtre_compomax(compteurmax,12,:,:)=vwnd(b-9,a,:,:);
264            vwnd700_filtre_compomax(compteurmax,13,:,:)=vwnd(b-8,a,:,:);
265            vwnd700_filtre_compomax(compteurmax,14,:,:)=vwnd(b-7,a,:,:);
266            vwnd700_filtre_compomax(compteurmax,15,:,:)=vwnd(b-6,a,:,:);
267            vwnd700_filtre_compomax(compteurmax,16,:,:)=vwnd(b-5,a,:,:);
268            vwnd700_filtre_compomax(compteurmax,17,:,:)=vwnd(b-4,a,:,:);
269            vwnd700_filtre_compomax(compteurmax,18,:,:)=vwnd(b-3,a,:,:);
270            vwnd700_filtre_compomax(compteurmax,19,:,:)=vwnd(b-2,a,:,:);
271            vwnd700_filtre_compomax(compteurmax,20,:,:)=vwnd(b-1,a,:,:);
272
273            vwnd700_filtre_compomax(compteurmax,21,:,:)=vwnd(b,a,:,:);
274            vwnd700_filtre_compomax(compteurmax,22,:,:)=vwnd(b+1,a,:,:);
275            vwnd700_filtre_compomax(compteurmax,23,:,:)=vwnd(b+2,a,:,:);
276            vwnd700_filtre_compomax(compteurmax,24,:,:)=vwnd(b+3,a,:,:);
277            vwnd700_filtre_compomax(compteurmax,25,:,:)=vwnd(b+4,a,:,:);
278            vwnd700_filtre_compomax(compteurmax,26,:,:)=vwnd(b+5,a,:,:);
279            vwnd700_filtre_compomax(compteurmax,27,:,:)=vwnd(b+6,a,:,:);
280            vwnd700_filtre_compomax(compteurmax,28,:,:)=vwnd(b+7,a,:,:);
281            vwnd700_filtre_compomax(compteurmax,29,:,:)=vwnd(b+8,a,:,:);
282            vwnd700_filtre_compomax(compteurmax,30,:,:)=vwnd(b+9,a,:,:);
283            vwnd700_filtre_compomax(compteurmax,31,:,:)=vwnd(b+10,a,:,:);
284% cas 2 : min deltas apres positifs et deltas avant negatifs
285         elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && ind_olr(b,a)<(-ind_olr2));
286            compteurmin=compteurmin+1;
287
288            vwnd700_filtre_compomin(compteurmin,1,:,:)=vwnd(b-20,a,:,:);
289            vwnd700_filtre_compomin(compteurmin,2,:,:)=vwnd(b-19,a,:,:);
290            vwnd700_filtre_compomin(compteurmin,3,:,:)=vwnd(b-18,a,:,:);
291            vwnd700_filtre_compomin(compteurmin,4,:,:)=vwnd(b-17,a,:,:);
292            vwnd700_filtre_compomin(compteurmin,5,:,:)=vwnd(b-16,a,:,:);
293            vwnd700_filtre_compomin(compteurmin,6,:,:)=vwnd(b-15,a,:,:);
294            vwnd700_filtre_compomin(compteurmin,7,:,:)=vwnd(b-14,a,:,:);
295            vwnd700_filtre_compomin(compteurmin,8,:,:)=vwnd(b-13,a,:,:);
296            vwnd700_filtre_compomin(compteurmin,9,:,:)=vwnd(b-12,a,:,:);
297            vwnd700_filtre_compomin(compteurmin,10,:,:)=vwnd(b-11,a,:,:);
298
299            vwnd700_filtre_compomin(compteurmin,11,:,:)=vwnd(b-10,a,:,:);
300            vwnd700_filtre_compomin(compteurmin,12,:,:)=vwnd(b-9,a,:,:);
301            vwnd700_filtre_compomin(compteurmin,13,:,:)=vwnd(b-8,a,:,:);
302            vwnd700_filtre_compomin(compteurmin,14,:,:)=vwnd(b-7,a,:,:);
303            vwnd700_filtre_compomin(compteurmin,15,:,:)=vwnd(b-6,a,:,:);
304            vwnd700_filtre_compomin(compteurmin,16,:,:)=vwnd(b-5,a,:,:);
305            vwnd700_filtre_compomin(compteurmin,17,:,:)=vwnd(b-4,a,:,:);
306            vwnd700_filtre_compomin(compteurmin,18,:,:)=vwnd(b-3,a,:,:);
307            vwnd700_filtre_compomin(compteurmin,19,:,:)=vwnd(b-2,a,:,:);
308            vwnd700_filtre_compomin(compteurmin,20,:,:)=vwnd(b-1,a,:,:);
309
310            vwnd700_filtre_compomin(compteurmin,21,:,:)=vwnd(b,a,:,:);
311            vwnd700_filtre_compomin(compteurmin,22,:,:)=vwnd(b+1,a,:,:);
312            vwnd700_filtre_compomin(compteurmin,23,:,:)=vwnd(b+2,a,:,:);
313            vwnd700_filtre_compomin(compteurmin,24,:,:)=vwnd(b+3,a,:,:);
314            vwnd700_filtre_compomin(compteurmin,25,:,:)=vwnd(b+4,a,:,:);
315            vwnd700_filtre_compomin(compteurmin,26,:,:)=vwnd(b+5,a,:,:);
316            vwnd700_filtre_compomin(compteurmin,27,:,:)=vwnd(b+6,a,:,:);
317            vwnd700_filtre_compomin(compteurmin,28,:,:)=vwnd(b+7,a,:,:);
318            vwnd700_filtre_compomin(compteurmin,29,:,:)=vwnd(b+8,a,:,:);
319            vwnd700_filtre_compomin(compteurmin,30,:,:)=vwnd(b+9,a,:,:);
320            vwnd700_filtre_compomin(compteurmin,31,:,:)=vwnd(b+10,a,:,:);
321         end;
322      end;
323   end;
324end;
325
326
327% on fait la somme des valeurs de vwnd700_filtreentiel sur le nbre d occurences
328% on obtent un tableau de la forme (jours-avant-apres,lon,lat)
329vwnd700_filtre_compomin1(:,:,:)=nanmean(vwnd700_filtre_compomin);
330clear vwnd700_filtre_compomin;
331vwnd700_filtre_compomax1(:,:,:)=nanmean(vwnd700_filtre_compomax);
332clear vwnd700_filtre_compomax;
333
334figure(1)
335orient('landscape')
336status=load('pression.mat','pression');
337pression=status.pression;
338clear status;
339
340for e=1:7;
341       a=7-(e-1)+1;
342       f=2*(a-1)+1;
343
344       clear delta_uwnd700_filtre_min;
345       clear delta_vwnd700_filtre_min;
346       delta_uwnd700_filtre_min=zeros(49,144);
347       delta_vwnd700_filtre_min=zeros(49,144);
348       for b=1:49;
349          for c=1:144;
350             delta_uwnd700_filtre_min(b,c)=uwnd700_filtre_compomin1(f,b,c)-uwnd700_filtre_compomax1(f,b,c);
351             delta_vwnd700_filtre_min(b,c)=vwnd700_filtre_compomin1(f,b,c)-vwnd700_filtre_compomax1(f,b,c);
352
353         end;
354       end;
355       subplot(7,1,8-e);
356
357       delta_vwnd700_filtre=zeros(49,144);
358       delta_uwnd700_filtre=zeros(49,144);
359       delta_vwnd700_filtre(:,1:72)=delta_vwnd700_filtre_min(:,73:144);
360       delta_vwnd700_filtre(:,73:144)=delta_vwnd700_filtre_min(:,1:72);
361       delta_uwnd700_filtre(:,1:72)=delta_uwnd700_filtre_min(:,73:144);
362       delta_uwnd700_filtre(:,73:144)=delta_uwnd700_filtre_min(:,1:72);
363       xi=[-180:3.75:176.25];
364       xi=xi';
365       yi=[-87.1591:3.708895:87.1591];
366       clear carteu cartev;
367       carteu=interp2(lon-180,lat,delta_uwnd700_filtre,xi,yi);
368       cartev=interp2(lon-180,lat,delta_vwnd700_filtre,xi,yi);
369       clear coco ;
370       coco(:,:)=pression(e,:,:);
371       hgt=interp2(lon-180,lat,coco,xi,yi);
372       palette
373       aa=contourf(xi,yi,hgt,[-20:2:20],'LineStyle','none');
374       hold on;
375       aa=quiver(xi,yi,carteu,cartev,'k');
376       if e==1;
377       set (gca,'XTick',[-180:60:180],'XTickLabel',[' 180';'120W';' 60W';'   0';' 60E';'120E';' 180'],'fontname','Arial','fontsize',6);
378       else
379      set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6);
380       end;
381       set (gca,'YTick',[-10:10:40],'YTickLabel',['   ';'  0';'10N';'20N';'30N';'   '],'fontname','Arial','fontsize',6);
382       cartemonde1;
383       if e==4;
384          co=colorbar;
385          set (co,'xtick',[-40:20:40],'xticklabel',[-40:20:40],'fontname','Arial','fontsize',6);
386       end;
387       set (gca,'tickDir','out');
388       axis equal
389       axis([-180 180 -10 40]);
390
391       if e==7;
392           title('Geopotential and winds at 700 Hpa - Sahelian mode','fontsize',8) ;
393       end;
394       text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10);
395       pos_vert=e/10;
396       set (gca,'position',[0.2 pos_vert 0.5 0.1])
397end;
398
399print('-depsc2','composite_sahel_1_u700_hgt.eps');
400
401figure(2)
402orient('landscape')
403status=load('pression.mat','pression');
404pression=status.pression;
405clear status;
406
407for e=1:7;
408       a=7-(e-1)+1+7;
409       f=2*(a-1)+1;
410
411       clear delta_uwnd700_filtre_min;
412       clear delta_vwnd700_filtre_min;
413       delta_uwnd700_filtre_min=zeros(49,144);
414       delta_vwnd700_filtre_min=zeros(49,144);
415       for b=1:49;
416          for c=1:144;
417             delta_uwnd700_filtre_min(b,c)=uwnd700_filtre_compomin1(f,b,c)-uwnd700_filtre_compomax1(f,b,c);
418             delta_vwnd700_filtre_min(b,c)=vwnd700_filtre_compomin1(f,b,c)-vwnd700_filtre_compomax1(f,b,c);
419
420         end;
421       end;
422       subplot(7,1,8-e);
423
424       delta_vwnd700_filtre=zeros(49,144);
425       delta_uwnd700_filtre=zeros(49,144);
426       delta_vwnd700_filtre(:,1:72)=delta_vwnd700_filtre_min(:,73:144);
427       delta_vwnd700_filtre(:,73:144)=delta_vwnd700_filtre_min(:,1:72);
428       delta_uwnd700_filtre(:,1:72)=delta_uwnd700_filtre_min(:,73:144);
429       delta_uwnd700_filtre(:,73:144)=delta_uwnd700_filtre_min(:,1:72);
430       xi=[-180:3.75:176.25];
431       xi=xi';
432       yi=[-87.1591:3.708895:87.1591];
433       clear carteu cartev;
434       carteu=interp2(lon-180,lat,delta_uwnd700_filtre,xi,yi);
435       cartev=interp2(lon-180,lat,delta_vwnd700_filtre,xi,yi);
436       clear coco ;
437       coco(:,:)=pression(e+7,:,:);
438       hgt=interp2(lon-180,lat,coco,xi,yi);
439       palette
440       aa=contourf(xi,yi,hgt,[-20:2:20],'LineStyle','none');
441       hold on;
442       aa=quiver(xi,yi,carteu,cartev,'k');
443       if e==1;
444       set (gca,'XTick',[-180:60:180],'XTickLabel',[' 180';'120W';' 60W';'   0';' 60E';'120E';' 180'],'fontname','Arial','fontsize',6);
445       else
446      set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6);
447       end;
448       set (gca,'YTick',[-10:10:40],'YTickLabel',['   ';'  0';'10N';'20N';'30N';'   '],'fontname','Arial','fontsize',6);
449       cartemonde1;
450       if e==4;
451          co=colorbar;
452          set (co,'xtick',[-40:20:40],'xticklabel',[-40:20:40],'fontname','Arial','fontsize',6);
453       end;
454       set (gca,'tickDir','out');
455       axis equal
456       axis([-180 180 -10 40]);
457
458       if e==7;
459           title('Geopotential and winds at 700 Hpa - Sahelian mode','fontsize',8) ;
460       end;
461       text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10);
462       pos_vert=e/10;
463       set (gca,'position',[0.2 pos_vert 0.5 0.1])
464end;
465
466print('-depsc2','composite_sahel_2_u700_hgt.eps');
Note: See TracBrowser for help on using the repository browser.