Changeset 16 for trunk/src/mode_sahelien/composite_hgt925_EOF12_ER30100.m
- Timestamp:
- 01/06/09 10:53:18 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/mode_sahelien/composite_hgt925_EOF12_ER30100.m
r15 r16 1 % 1 2 % initialisation 2 clear; close all;3 3 clear; 4 close all; 4 5 5 6 % etape 1 :on selectionne les mois d ete (attention aux annees bixestiles!) … … 8 9 hgt925=zeros(122,49,144,28); 9 10 10 siz=31+28+31+30+31+1; sizbis=31+28+31+30+31+30+31+31+30; 11 siz=31+28+31+30+31+1; 12 sizbis=31+28+31+30+31+30+31+31+30; 11 13 ncload('hgt925.1979.nc'); hgt925(:,:,:,1)=hgt(siz:sizbis,1,:,:); clear hgt; 12 14 ncload('hgt925.1980.nc'); hgt925(:,:,:,2)=hgt(siz+1:sizbis+1,1,:,:); clear hgt; … … 38 40 ncload('hgt925.2006.nc'); hgt925(:,:,:,28)=hgt(siz:sizbis,1,:,:);clear hgt; 39 41 40 41 42 clear siz1 siz1bis; 42 43 … … 46 47 47 48 %%%%%%%%%%%%%%% FICHIER OLR 48 % ouverture / lecture du fichier 49 % ouverture / lecture du fichier 49 50 load eof12_rossby_30100.txt 50 51 ind_olr=reshape(eof12_rossby_30100,122,28); … … 53 54 % on calcule l ecart type sur la serie entiere de l OLR 54 55 ind_olr1=reshape(ind_olr, k,1); 55 ind_olr2=std(ind_olr1); clear ind_olr1;56 57 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 56 ind_olr2=std(ind_olr1); 57 clear ind_olr1; 58 58 59 % GEOPOTENTIEL 59 60 % intialisation des tableaux 60 hgt925_filtre_compomax=NaN*ones(60,31,49,144); 61 hgt925_filtre_compomax=NaN*ones(60,31,49,144); 61 62 hgt925_filtre_compomin=NaN*ones(60,31,49,144); 62 63 63 compteurmax=0; 64 compteurmax=0; 64 65 compteurmin=0; 65 66 66 for a=1:28; % boucle sur les annees 67 clear delta; 67 % boucle sur les annees 68 for a=1:28; 69 clear delta; 68 70 delta=zeros(122,1); 69 for b=1:122; % boucle sur les jours : on prepare avant 71 % boucle sur les jours : on prepare avant 72 for b=1:122; 70 73 if b>1; 71 74 delta(b)=ind_olr(b,a)-ind_olr(b-1,a); 72 75 end; 73 76 end; 74 for b=1:122; % boucle sur les jours 75 if (b>20 && b<112); % 1ere condition (laisser une marge de 10) 76 % cas 1 : max deltas apres negetifs et deltas avant postifs 77 if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && ind_olr(b,a)>ind_olr2); 77 % boucle sur les jours 78 for b=1:122; 79 % 1ere condition (laisser une marge de 10) 80 if (b>20 && b<112); 81 % cas 1 : max deltas apres negetifs et deltas avant postifs 82 if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && ind_olr(b,a)>ind_olr2); 78 83 compteurmax=compteurmax+1; 79 84 hgt925_filtre_compomax(compteurmax,1,:,:)=hgt(b-20,a,:,:); … … 110 115 hgt925_filtre_compomax(compteurmax,30,:,:)=hgt(b+9,a,:,:); 111 116 hgt925_filtre_compomax(compteurmax,31,:,:)=hgt(b+10,a,:,:); 112 % cas 2 : min deltas apres positifs et deltas avant negatifs 117 % cas 2 : min deltas apres positifs et deltas avant negatifs 113 118 elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && ind_olr(b,a)<(-ind_olr2)); 114 119 compteurmin=compteurmin+1; … … 149 154 end; 150 155 end; 151 end; 156 end; 152 157 end; 153 158 … … 155 160 % on fait la somme des valeurs de hgt925_filtreentiel sur le nbre d occurences 156 161 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 157 hgt925_filtre_compomin1(:,:,:)=nanmean(hgt925_filtre_compomin); clear hgt925_filtre_compomin; 158 hgt925_filtre_compomax1(:,:,:)=nanmean(hgt925_filtre_compomax); clear hgt925_filtre_compomax; 159 160 x=[-10:1:10]'; %' 161 162 figure(1); 162 hgt925_filtre_compomin1(:,:,:)=nanmean(hgt925_filtre_compomin); 163 clear hgt925_filtre_compomin; 164 hgt925_filtre_compomax1(:,:,:)=nanmean(hgt925_filtre_compomax); 165 clear hgt925_filtre_compomax; 166 167 x=[-10:1:10]'; 168 169 figure(1); 163 170 orient('landscape'); 164 171 lat=lat'; … … 216 223 else 217 224 set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 218 end; 225 end; 219 226 set (gca,'YTick',[-20:10:20],'YTickLabel',['20S';'10S';' 0 ';'10N';'20N';],'fontname','Arial','fontsize',6); 220 227 if e==4; 221 co=colorbar; 228 co=colorbar; 222 229 set (co,'ytick',[-20:5:20],'yticklabel',[-20:5:20],'fontname','Arial','fontsize',6); 223 230 end; 224 hold on; cartemonde1; 231 hold on; 232 cartemonde1; 225 233 set (gca,'tickDir','out'); 226 234 axis equal; … … 231 239 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 232 240 pos_vert=e/10; 233 set (gca,'position',[0.2 pos_vert 0.5 0.1]); 241 set (gca,'position',[0.2 pos_vert 0.5 0.1]); 234 242 end; 235 243 236 244 print -depsc2 composite_hgt925_EOF12_ER30100_part1.eps; 237 245 238 figure(2); 246 figure(2); 239 247 orient('landscape'); 240 248 … … 289 297 else 290 298 set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 291 end; 299 end; 292 300 set (gca,'YTick',[-20:10:20],'YTickLabel',['20S';'10S';' 0 ';'10N';'20N';],'fontname','Arial','fontsize',6); 293 301 if e==4; 294 co=colorbar; 302 co=colorbar; 295 303 set (co,'ytick',[-20:5:20],'yticklabel',[-20:5:20],'fontname','Arial','fontsize',6); 296 304 end; 297 hold on; cartemonde1; 305 hold on; 306 cartemonde1; 298 307 set (gca,'tickDir','out'); 299 308 axis equal; … … 304 313 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 305 314 pos_vert=e/10; 306 set (gca,'position',[0.2 pos_vert 0.5 0.1]); 315 set (gca,'position',[0.2 pos_vert 0.5 0.1]); 307 316 end; 308 317
Note: See TracChangeset
for help on using the changeset viewer.