Changeset 16 for trunk/src/mode_sahelien/composite_uwnd700_EOF12_ER1030.m
- Timestamp:
- 01/06/09 10:53:18 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/mode_sahelien/composite_uwnd700_EOF12_ER1030.m
r15 r16 1 % 2 1 3 % initialisation 2 clear; close all; 4 clear; 5 close all; 3 6 4 7 … … 8 11 vwnd700=zeros(122,49,144,28); 9 12 10 siz=31+28+31+30+31+1; sizbis=31+28+31+30+31+30+31+31+30; 13 siz=31+28+31+30+31+1; 14 sizbis=31+28+31+30+31+30+31+31+30; 11 15 ncload('vwnd700.1979.nc'); vwnd700(:,:,:,1)=vwnd(siz:sizbis,1,:,:); clear vwnd; 12 16 ncload('vwnd700.1980.nc'); vwnd700(:,:,:,2)=vwnd(siz+1:sizbis+1,1,:,:); clear vwnd; … … 47 51 48 52 %%%%%%%%%%%%%%% FICHIER OLR 49 % ouverture / lecture du fichier 53 % ouverture / lecture du fichier 50 54 load eof12_rossby_1030.txt 51 55 ind_olr=reshape(eof12_rossby_1030,122,28); … … 54 58 % on calcule l ecart type sur la serie entiere de l OLR 55 59 ind_olr1=reshape(ind_olr, k,1); 56 ind_olr2=std(ind_olr1); clear ind_olr1;57 58 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 60 ind_olr2=std(ind_olr1); 61 clear ind_olr1; 62 59 63 % GEOPOTENTIEL 60 64 % intialisation des tableaux 61 vwnd700_filtre_compomax=NaN*ones(60,31,49,144); 65 vwnd700_filtre_compomax=NaN*ones(60,31,49,144); 62 66 vwnd700_filtre_compomin=NaN*ones(60,31,49,144); 63 67 64 compteurmax=0; 68 compteurmax=0; 65 69 compteurmin=0; 66 70 67 for a=1:28; % boucle sur les annees 68 clear delta; 71 % boucle sur les annees 72 for a=1:28; 73 clear delta; 69 74 delta=zeros(122,1); 70 for b=1:122; % boucle sur les jours : on prepare avant 75 % boucle sur les jours : on prepare avant 76 for b=1:122; 71 77 if b>1; 72 78 delta(b)=ind_olr(b,a)-ind_olr(b-1,a); 73 79 end; 74 80 end; 75 for b=1:122; % boucle sur les jours 76 if (b>20 && b<112); % 1ere condition (laisser une marge de 10) 77 % cas 1 : max deltas apres negetifs et deltas avant postifs 78 if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && ind_olr(b,a)>ind_olr2); 81 % boucle sur les jours 82 for b=1:122; 83 % 1ere condition (laisser une marge de 10) 84 if (b>20 && b<112); 85 % cas 1 : max deltas apres negetifs et deltas avant postifs 86 if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && ind_olr(b,a)>ind_olr2); 79 87 compteurmax=compteurmax+1; 80 88 vwnd700_filtre_compomax(compteurmax,1,:,:)=vwnd(b-20,a,:,:); … … 111 119 vwnd700_filtre_compomax(compteurmax,30,:,:)=vwnd(b+9,a,:,:); 112 120 vwnd700_filtre_compomax(compteurmax,31,:,:)=vwnd(b+10,a,:,:); 113 % cas 2 : min deltas apres positifs et deltas avant negatifs 121 % cas 2 : min deltas apres positifs et deltas avant negatifs 114 122 elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && ind_olr(b,a)<(-ind_olr2)); 115 123 compteurmin=compteurmin+1; … … 150 158 end; 151 159 end; 152 end; 160 end; 153 161 end; 154 162 … … 156 164 % on fait la somme des valeurs de vwnd700_filtreentiel sur le nbre d occurences 157 165 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 158 vwnd700_filtre_compomin1(:,:,:)=nanmean(vwnd700_filtre_compomin); clear vwnd700_filtre_compomin; 159 vwnd700_filtre_compomax1(:,:,:)=nanmean(vwnd700_filtre_compomax); clear vwnd700_filtre_compomax; 160 161 x=[-10:1:10]'; %' 162 163 figure(1); 166 vwnd700_filtre_compomin1(:,:,:)=nanmean(vwnd700_filtre_compomin); 167 clear vwnd700_filtre_compomin; 168 vwnd700_filtre_compomax1(:,:,:)=nanmean(vwnd700_filtre_compomax); 169 clear vwnd700_filtre_compomax; 170 171 x=[-10:1:10]'; 172 173 figure(1); 164 174 orient('landscape'); 165 175 lat=lat'; … … 217 227 else 218 228 set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 219 end; 229 end; 220 230 set (gca,'YTick',[-20:10:20],'YTickLabel',['20S';'10S';' 0 ';'10N';'20N';],'fontname','Arial','fontsize',6); 221 231 if e==4; 222 co=colorbar; 232 co=colorbar; 223 233 set (co,'ytick',[-5:1:5],'yticklabel',[-5:1:5],'fontname','Arial','fontsize',6); 224 234 end; 225 hold on; cartemonde1; 235 hold on; 236 cartemonde1; 226 237 set (gca,'tickDir','out'); 227 238 axis equal; … … 232 243 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 233 244 pos_vert=e/10; 234 set (gca,'position',[0.2 pos_vert 0.5 0.1]); 245 set (gca,'position',[0.2 pos_vert 0.5 0.1]); 235 246 end; 236 247 237 248 print -depsc2 composite_vwnd700_EOF12_ER1030_part1.eps; 238 249 239 figure(2); 250 figure(2); 240 251 orient('landscape'); 241 252 … … 291 302 else 292 303 set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 293 end; 304 end; 294 305 set (gca,'YTick',[-20:10:20],'YTickLabel',['20S';'10S';' 0 ';'10N';'20N';],'fontname','Arial','fontsize',6); 295 306 if e==4; 296 co=colorbar; 307 co=colorbar; 297 308 set (co,'ytick',[-5:1:5],'yticklabel',[-5:1:5],'fontname','Arial','fontsize',6); 298 309 end; 299 hold on; cartemonde1; 310 hold on; 311 cartemonde1; 300 312 set (gca,'tickDir','out'); 301 313 axis equal; … … 306 318 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 307 319 pos_vert=e/10; 308 set (gca,'position',[0.2 pos_vert 0.5 0.1]); 320 set (gca,'position',[0.2 pos_vert 0.5 0.1]); 309 321 end; 310 322 311 323 print -depsc2 composite_vwnd700_EOF12_ER1030_part2.eps; 312
Note: See TracChangeset
for help on using the changeset viewer.