Changeset 16 for trunk/src/mode_sahelien/composite_olr_phase5.m
- Timestamp:
- 01/06/09 10:53:18 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/mode_sahelien/composite_olr_phase5.m
r15 r16 1 % 2 1 3 % initialisation 2 clear; 4 clear; 3 5 4 6 cd('..'); … … 25 27 olr1=zeros(siz,21,144); 26 28 olr1(:,:,1:72)=olr(:,:,73:144); 27 olr1(:,:,73:144)=olr(:,:,1:72); clear olr; 29 olr1(:,:,73:144)=olr(:,:,1:72); 30 clear olr; 28 31 29 32 % on permture le temps et l'espce … … 37 40 vecteur=zeros(siz,1); 38 41 for c=1:siz; 39 vecteur(c)=olr1(a,b,c); 42 vecteur(c)=olr1(a,b,c); 40 43 end 41 44 std_c=nanmean(vecteur); 42 45 for c=1:siz; 43 olr1(a,b,c)=olr1(a,b,c)-std_c; 44 end 46 olr1(a,b,c)=olr1(a,b,c)-std_c; 47 end 45 48 end 46 end; 47 48 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 49 end; 50 49 51 % on met sous la forme longitude latitude jour 50 51 52 52 53 % on selectionne les mois d ete (attention aux annees bixestiles!) … … 151 152 clear olr; 152 153 olr2=reshape(olr1,3416,21,144); 153 olr3(:,:)=mean(olr2); clear olr2;154 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 154 olr3(:,:)=mean(olr2); 155 clear olr2; 155 156 156 157 load phase5.txt; 157 158 phase5=reshape(phase5,122,28,3); 158 159 olr_filtre_compomax=zeros(90,21,144); 159 compteurmax=0; 160 161 for a=1:28; % boucle sur les annees 162 for b=1:122; % boucle sur les jours 163 if (b>11 && b<112); % 1ere condition (laisser une marge de 10) 164 if phase5(b,a,1)==1; 160 compteurmax=0; 161 162 % boucle sur les annees 163 for a=1:28; 164 % boucle sur les jours 165 for b=1:122; 166 % 1ere condition (laisser une marge de 10) 167 if (b>11 && b<112); 168 if phase5(b,a,1)==1; 165 169 compteurmax=compteurmax+1; 166 170 olr_filtre_compomax(compteurmax,:,:)=olr1(b,a,:,:); 167 171 end; 168 172 end; 169 end; 173 end; 170 174 end; 171 175 … … 175 179 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 176 180 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 177 olr_filtre_compomax1(:,:)=mean(olr_filtre_compomax); clear olr_filtre_compomax;178 179 180 x=[-10:1:10]'; %'181 182 lat=lat(21:41)'; 181 olr_filtre_compomax1(:,:)=mean(olr_filtre_compomax); 182 clear olr_filtre_compomax; 183 184 x=[-10:1:10]'; 185 186 lat=lat(21:41)'; 183 187 lon=lon-180; 184 delta=olr_filtre_compomax1-olr3; clear olr3; 188 delta=olr_filtre_compomax1-olr3; 189 clear olr3; 185 190 for b=1:21; 186 191 for c=1:144; … … 198 203 set (gca,'XTick',[-180:60:180],'XTickLabel',[' 180';'120W';' 60W';' 0';' 60E';'120E';' 180'],'fontsize',6); 199 204 set (gca,'YTick',[-10:10:40],'YTickLabel',[' ';' 0';'10N';'20N';'30N';' '],'fontname','Arial','fontsize',6); 200 co=colorbar; 205 co=colorbar; 201 206 set(co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6); 202 hold on; cartemonde1; 207 hold on; 208 cartemonde1; 203 209 set (gca,'tickDir','out'); 204 210 axis equal; 205 211 axis([-180 180 -10 40]); 206 207 208 209 210 211 212
Note: See TracChangeset
for help on using the changeset viewer.