source: trunk/src/mode_sahelien/composite_hgt500_eof234_1025_3j.m @ 26

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

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