- Timestamp:
- 01/06/09 12:58:01 (15 years ago)
- Location:
- trunk/src
- Files:
-
- 146 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/SIMULS_IRCAAM/anomalies_olr_AFR_1979_2000.m
r16 r18 16 16 17 17 olr=zeros(122,30,64,128); 18 19 18 20 19 ncload('rlut_d.AfNQIVIV.197106-09.nc','rlut'); -
trunk/src/SIMULS_IRCAAM/carte_eof_ircaam_1030_as.m
r17 r18 100 100 set (gca,'tickDir','out'); 101 101 102 103 104 102 clear carte1 x y a b carte coco; 105 103 carte1=C(:,357); -
trunk/src/SIMULS_IRCAAM/carte_eof_ircaam_1030_ctl.m
r17 r18 102 102 set (gca,'tickDir','out'); 103 103 104 105 106 104 clear carte1 x y a b carte coco; 107 105 carte1=C(:,357); -
trunk/src/SIMULS_IRCAAM/carte_eof_ircaam_1030_tr.m
r17 r18 100 100 set (gca,'tickDir','out'); 101 101 102 103 104 102 clear carte1 x y a b carte coco; 105 103 carte1=C(:,357); -
trunk/src/SIMULS_IRCAAM/carte_eof_noaa_1_8.m
r17 r18 141 141 set (gca,'tickDir','out'); 142 142 143 144 145 146 %********************147 143 subplot(3,2,4); 148 144 vaude=0; -
trunk/src/SIMULS_IRCAAM/carte_eof_noaa_8_28.m
r17 r18 141 141 set (gca,'tickDir','out'); 142 142 143 144 145 146 %********************147 143 subplot(3,2,4); 148 144 vaude=0 -
trunk/src/SIMULS_IRCAAM/carte_eof_olr_noaa_jas.m
r17 r18 142 142 set (gca,'tickDir','out'); 143 143 144 145 146 147 144 subplot(3,2,4); 148 145 vaude=0; -
trunk/src/SIMULS_IRCAAM/carte_eof_olr_noaa_jas_20.m
r17 r18 141 141 set (gca,'tickDir','out'); 142 142 143 144 145 146 %********************147 143 subplot(3,2,4); 148 144 vaude=0; -
trunk/src/SIMULS_IRCAAM/climato_olr_noaa.m
r16 r18 73 73 siz28=siz27+365; siz28bis=siz28+a; % 2000 74 74 75 76 75 olr(:,:,:,1)=olr1(:,:,siz1:siz1bis); % 1979 77 76 olr(:,:,:,2)=olr1(:,:,siz2:siz2bis); % 1980 … … 102 101 olr(:,:,:,27)=olr1(:,:,siz27:siz27bis); % 2005 103 102 olr(:,:,:,28)=olr1(:,:,siz28:siz28bis); % 2006 104 105 103 106 104 clear siz1 siz1bis; -
trunk/src/SIMULS_IRCAAM/composite_olr_AFR_eof1_1030.m
r17 r18 222 222 end 223 223 224 225 224 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 226 225 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/SIMULS_IRCAAM/composite_olr_AFR_eof1_30100.m
r17 r18 221 221 end 222 222 223 224 223 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 225 224 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/SIMULS_IRCAAM/composite_olr_AFR_eof2_1030.m
r17 r18 209 209 end; 210 210 211 212 211 % on fait le test de Student 213 212 nb_value=min(compteurmin,compteurmax); -
trunk/src/SIMULS_IRCAAM/composite_olr_AFR_eof2_30100.m
r17 r18 221 221 end 222 222 223 224 223 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 225 224 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/SIMULS_IRCAAM/composite_olr_AS_eof1_1030.m
r17 r18 222 222 end 223 223 224 225 224 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 226 225 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/SIMULS_IRCAAM/composite_olr_AS_eof1_30100.m
r17 r18 222 222 end 223 223 224 225 224 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 226 225 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/SIMULS_IRCAAM/composite_olr_AS_eof2_1030.m
r17 r18 222 222 end 223 223 224 225 224 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 226 225 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/SIMULS_IRCAAM/composite_olr_AS_eof2_30100.m
r17 r18 222 222 end 223 223 224 225 224 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 226 225 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/SIMULS_IRCAAM/composite_olr_CTL_eof1_1030.m
r17 r18 224 224 end 225 225 226 227 226 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 228 227 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/SIMULS_IRCAAM/composite_olr_CTL_eof1_30100.m
r17 r18 93 93 ncload('rlut_d.CtIV.200006-09.nc','rlut','lat','lon'); 94 94 olr(:,29,:,:)=rlut; clear rlut; 95 96 95 97 96 cd('..'); … … 224 223 end 225 224 226 227 225 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 228 226 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/SIMULS_IRCAAM/composite_olr_CTL_eof2_1030.m
r17 r18 223 223 end 224 224 225 226 225 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 227 226 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/SIMULS_IRCAAM/composite_olr_CTL_eof2_30100.m
r17 r18 223 223 end 224 224 225 226 225 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 227 226 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/SIMULS_IRCAAM/composite_olr_TR_eof1_1030.m
r17 r18 222 222 end 223 223 224 225 226 224 olr_filtre_compomin1(:,:,:)=nanmean(olr_filtre_compomin); 227 225 clear olr_filtre_compomin; -
trunk/src/SIMULS_IRCAAM/composite_olr_TR_eof2_1030.m
r17 r18 230 230 end 231 231 232 233 234 232 olr_filtre_compomin1(:,:,:)=nanmean(olr_filtre_compomin); 235 233 clear olr_filtre_compomin; -
trunk/src/SIMULS_IRCAAM/composite_olr_TR_eof2_30100.m
r17 r18 230 230 end 231 231 232 233 234 232 olr_filtre_compomin1(:,:,:)=nanmean(olr_filtre_compomin); 235 233 clear olr_filtre_compomin; -
trunk/src/SIMULS_IRCAAM/correlations_eofs_olr_jas.m
r17 r18 86 86 siz27=siz26+365; siz27bis=siz27+a; % 2000 87 87 siz28=siz27+365; siz28bis=siz28+a; % 2000 88 89 88 90 89 olr(:,:,:,1)=olr1(:,:,siz1:siz1bis); % 1979 … … 117 116 olr(:,:,:,28)=olr1(:,:,siz28:siz28bis); % 2006 118 117 119 120 118 clear siz1 siz1bis; 121 119 clear siz2 siz2bis; … … 157 155 load eof1_noaa_annuel.txt 158 156 159 160 157 correlations_eof1=zeros(37,144); 161 158 correlations_eof2=zeros(37,144); -
trunk/src/SIMULS_IRCAAM/correlations_eofs_olr_jas_20.m
r17 r18 87 87 siz27=siz26+365; siz27bis=siz27+a; % 2000 88 88 siz28=siz27+365; siz28bis=siz28+a; % 2000 89 90 89 91 90 olr(:,:,:,1)=olr1(:,:,siz1:siz1bis); % 1979 … … 118 117 olr(:,:,:,28)=olr1(:,:,siz28:siz28bis); % 2006 119 118 120 121 119 clear siz1 siz1bis; 122 120 clear siz2 siz2bis; … … 158 156 load eof1_noaa_annuel_20.txt 159 157 160 161 158 correlations_eof1=zeros(37,144); 162 159 correlations_eof2=zeros(37,144); -
trunk/src/SIMULS_IRCAAM/palette.m
r17 r18 24 24 palette2(1,:) = [0.9 0 0.9]; % Mauve 25 25 26 27 26 colormap(palette2); -
trunk/src/SIMULS_IRCAAM/palettecomplet.m
r17 r18 25 25 palette2(1,:) = [0.9 0 0.9]; % Mauve 26 26 27 28 27 colormap(palette2); -
trunk/src/SIMULS_IRCAAM/prepare_pluie_AFR_annuel.m
r16 r18 16 16 17 17 pluie=zeros(122,30,64,128); 18 19 18 20 19 ncload('pr_d.AfNQIVIV.197106-09.nc','pr'); … … 134 133 end 135 134 136 137 138 135 % calcul de la matrice de covariance 139 136 vaudou=zeros(360,360); -
trunk/src/SIMULS_IRCAAM/time_serie_afr_1030.m
r17 r18 8 8 PC1=olr*C; 9 9 10 11 10 %%%%%%% cronique de la VP1 12 11 time_serie_olr1=PC1(:,360); … … 15 14 dt_olr1=std(time_serie_olr1); 16 15 dt_olr2=std(time_serie_olr2); 17 18 16 19 17 for a=1:3538; … … 60 58 axis([1 122 -3 3]); 61 59 62 63 60 set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); 64 61 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); -
trunk/src/SIMULS_IRCAAM/time_serie_afr_30100.m
r17 r18 13 13 dt_olr1=std(time_serie_olr1); 14 14 dt_olr2=std(time_serie_olr2); 15 16 15 17 16 for a=1:3538; … … 58 57 axis([1 122 -3 3]); 59 58 60 61 59 set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); 62 60 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); -
trunk/src/SIMULS_IRCAAM/time_serie_as_1030.m
r17 r18 9 9 PC1=olr*C; 10 10 11 12 11 %%%%%%% cronique de la VP1 13 12 time_serie_olr1=PC1(:,360); … … 16 15 dt_olr1=std(time_serie_olr1); 17 16 dt_olr2=std(time_serie_olr2); 18 19 17 20 18 for a=1:3538; … … 61 59 axis([1 122 -3 3]); 62 60 63 64 61 set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); 65 62 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); -
trunk/src/SIMULS_IRCAAM/time_serie_as_30100.m
r17 r18 9 9 PC1=olr*C; 10 10 11 12 11 %%%%%%% cronique de la VP1 13 12 time_serie_olr1=PC1(:,360); … … 16 15 dt_olr1=std(time_serie_olr1); 17 16 dt_olr2=std(time_serie_olr2); 18 19 17 20 18 for a=1:3538; … … 61 59 axis([1 122 -3 3]); 62 60 63 64 61 set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); 65 62 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); -
trunk/src/SIMULS_IRCAAM/time_serie_ctl_1030.m
r17 r18 8 8 PC1=olr*C; 9 9 10 11 10 %%%%%%% cronique de la VP1 12 11 time_serie_olr1=PC1(:,360); … … 15 14 dt_olr1=std(time_serie_olr1); 16 15 dt_olr2=std(time_serie_olr2); 17 18 16 19 17 for a=1:3538; … … 60 58 axis([1 122 -3 3]); 61 59 62 63 60 set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); 64 61 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); -
trunk/src/SIMULS_IRCAAM/time_serie_ctl_30100.m
r17 r18 8 8 PC1=olr*C; 9 9 10 11 10 %%%%%%% cronique de la VP1 12 11 time_serie_olr1=PC1(:,360); … … 15 14 dt_olr1=std(time_serie_olr1); 16 15 dt_olr2=std(time_serie_olr2); 17 18 16 19 17 for a=1:3538; … … 60 58 axis([1 122 -3 3]); 61 59 62 63 60 set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); 64 61 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); -
trunk/src/SIMULS_IRCAAM/time_serie_noaa_annuel.m
r17 r18 8 8 olr=olr'; 9 9 PC1=olr*C; 10 11 10 12 11 %%%%%%% cronique de la VP1 -
trunk/src/SIMULS_IRCAAM/time_serie_noaa_annuel_20.m
r17 r18 6 6 olr=olr'; 7 7 PC1=olr*C; 8 9 8 10 9 %%%%%%% cronique de la VP1 -
trunk/src/SIMULS_IRCAAM/time_serie_tr_1030.m
r17 r18 7 7 PC1=olr*C; 8 8 9 10 9 %%%%%%% cronique de la VP1 11 10 time_serie_olr1=PC1(:,360); … … 14 13 dt_olr1=std(time_serie_olr1); 15 14 dt_olr2=std(time_serie_olr2); 16 17 15 18 16 for a=1:3538; … … 59 57 axis([1 122 -3 3]); 60 58 61 62 59 set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); 63 60 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); -
trunk/src/SIMULS_IRCAAM/time_serie_tr_30100.m
r17 r18 8 8 PC1=olr*C; 9 9 10 11 10 %%%%%%% cronique de la VP1 12 11 time_serie_olr1=PC1(:,360); … … 15 14 dt_olr1=std(time_serie_olr1); 16 15 dt_olr2=std(time_serie_olr2); 17 18 16 19 17 for a=1:3538; … … 60 58 axis([1 122 -3 3]); 61 59 62 63 60 set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); 64 61 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); -
trunk/src/cartemonde.m
r17 r18 53 53 load monde.mat 54 54 55 56 55 % trace fond de carte continents 57 56 [nll,ncc]=size(npconti); -
trunk/src/cartemonde1.m
r17 r18 52 52 load monde.mat 53 53 54 55 54 % trace fond de carte continents 56 55 [nll,ncc]=size(npconti); -
trunk/src/mode_sahelien/composite_hgt500_eof12_ER30100.m
r17 r18 194 194 end; 195 195 196 197 196 % on fait la somme des valeurs de hgt500_filtreentiel sur le nbre d occurences 198 197 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/mode_sahelien/composite_hgt500_eof1_1030.m
r17 r18 208 208 end 209 209 210 211 212 213 210 % on fait la somme des valeurs de hgt500_filtreentiel sur le nbre d occurences 214 211 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/mode_sahelien/composite_hgt500_eof1_30100.m
r17 r18 209 209 end 210 210 211 212 213 214 211 % on fait la somme des valeurs de hgt500_filtreentiel sur le nbre d occurences 215 212 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/mode_sahelien/composite_hgt500_eof23_1030.m
r17 r18 57 57 end 58 58 59 60 59 hgt1=zeros(122,28,49,144); 61 60 hgt1(:,:,:,1:72)=hgt500(:,:,:,73:144); … … 82 81 end; 83 82 clear hgt2 hgt4; 84 85 83 86 84 %%%%%%%%%%%%%%% FICHIER OLR … … 196 194 end; 197 195 198 199 196 % on fait la somme des valeurs de hgt500_filtreentiel sur le nbre d occurences 200 197 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/mode_sahelien/composite_hgt500_eof23_2590.m
r17 r18 46 46 clear hgt925; 47 47 48 49 48 %%%%%%%%%%%%%%% FICHIER OLR 50 49 % ouverture / lecture du fichier … … 158 157 end; 159 158 160 161 159 % on fait la somme des valeurs de hgt925_filtreentiel sur le nbre d occurences 162 160 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/mode_sahelien/composite_hgt500_eof23_30100.m
r17 r18 83 83 clear hgt2 hgt4; 84 84 85 86 85 %%%%%%%%%%%%%%% FICHIER OLR 87 86 load eof23_30100.txt; … … 210 209 end 211 210 212 213 214 215 211 % on fait la somme des valeurs de hgt500_filtreentiel sur le nbre d occurences 216 212 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/mode_sahelien/composite_hgt700_EOF12_ER1030.m
r17 r18 46 46 clear hgt700; 47 47 48 49 48 %%%%%%%%%%%%%%% FICHIER OLR 50 49 % ouverture / lecture du fichier … … 157 156 end; 158 157 end; 159 160 158 161 159 % on fait la somme des valeurs de hgt700_filtreentiel sur le nbre d occurences -
trunk/src/mode_sahelien/composite_hgt700_EOF12_ER30100_dat.m
r17 r18 41 41 ncload('hgt700.2006.nc'); hgt700(:,:,:,28)=hgt(siz:sizbis,1,:,:);clear hgt; 42 42 43 44 43 clear siz1 siz1bis; 45 44 46 45 hgt=permute(hgt700,[1 4 2 3]); 47 46 clear hgt700; 48 49 47 50 48 %%%%%%%%%%%%%%% FICHIER OLR … … 159 157 end; 160 158 end; 161 162 159 163 160 % on fait la somme des valeurs de hgt700_filtreentiel sur le nbre d occurences -
trunk/src/mode_sahelien/composite_hgt925_EOF12_ER1030.m
r17 r18 46 46 clear hgt925; 47 47 48 49 48 %%%%%%%%%%%%%%% FICHIER OLR 50 49 % ouverture / lecture du fichier … … 157 156 end; 158 157 end; 159 160 158 161 159 % on fait la somme des valeurs de hgt925_filtreentiel sur le nbre d occurences -
trunk/src/mode_sahelien/composite_hgt925_EOF12_ER30100.m
r17 r18 45 45 clear hgt925; 46 46 47 48 47 %%%%%%%%%%%%%%% FICHIER OLR 49 48 % ouverture / lecture du fichier … … 157 156 end; 158 157 159 160 158 % on fait la somme des valeurs de hgt925_filtreentiel sur le nbre d occurences 161 159 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 216 214 end 217 215 218 219 delta_hgt925=flipud(delta_hgt925); 216 delta_hgt925=flipud(delta_hgt925); 220 217 aa=contourf(lon-180,lat,delta_hgt925,[-20:2:20],'LineStyle','none'); 221 218 if e==1; -
trunk/src/mode_sahelien/composite_hgt925_EOF12_ER30100_dat.m
r17 r18 1 1 % 2 2 3 % initialisation 3 4 clear; … … 45 46 clear hgt925; 46 47 47 48 48 %%%%%%%%%%%%%%% FICHIER OLR 49 49 % ouverture / lecture du fichier … … 157 157 end; 158 158 end; 159 160 159 161 160 % on fait la somme des valeurs de hgt925_filtreentiel sur le nbre d occurences -
trunk/src/mode_sahelien/composite_hgt925_EOF1_30100.m
r17 r18 46 46 clear hgt925; 47 47 48 49 48 %%%%%%%%%%%%%%% FICHIER OLR 50 49 % ouverture / lecture du fichier … … 158 157 end; 159 158 160 161 159 % on fait la somme des valeurs de hgt925_filtreentiel sur le nbre d occurences 162 160 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 283 281 284 282 print -depsc2 composite_hgt925_EOF1_30100_part2.eps; 285 -
trunk/src/mode_sahelien/composite_hgt925_EOF23_1030.m
r17 r18 46 46 clear hgt925; 47 47 48 49 48 %%%%%%%%%%%%%%% FICHIER OLR 50 49 % ouverture / lecture du fichier … … 157 156 end; 158 157 end; 159 160 158 161 159 % on fait la somme des valeurs de hgt925_filtreentiel sur le nbre d occurences … … 294 292 delta_hgt925=flipud(delta_hgt925); 295 293 296 297 294 aa=contourf(lon-180,lat,delta_hgt925,[-20:2:20],'LineStyle','none'); 298 295 if e==1; -
trunk/src/mode_sahelien/composite_hgt925_EOF23_30100.m
r17 r18 41 41 ncload('hgt925.2006.nc'); hgt925(:,:,:,28)=hgt(siz:sizbis,1,:,:);clear hgt; 42 42 43 44 43 clear siz1 siz1bis; 45 44 … … 158 157 end; 159 158 160 161 159 % on fait la somme des valeurs de hgt925_filtreentiel sur le nbre d occurences 162 160 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/mode_sahelien/composite_hgt925_EOF24_1030.m
r17 r18 41 41 ncload('hgt925.2006.nc'); hgt925(:,:,:,28)=hgt(siz:sizbis,1,:,:);clear hgt; 42 42 43 44 43 clear siz1 siz1bis; 45 44 … … 157 156 end; 158 157 end; 159 160 158 161 159 % on fait la somme des valeurs de hgt925_filtreentiel sur le nbre d occurences … … 294 292 delta_hgt925=flipud(delta_hgt925); 295 293 296 297 294 aa=contourf(lon-180,lat,delta_hgt925,[-20:2:20],'LineStyle','none'); 298 295 if e==1; -
trunk/src/mode_sahelien/composite_hgt925_eof23_2590.m
r17 r18 46 46 clear hgt925; 47 47 48 49 48 %%%%%%%%%%%%%%% FICHIER OLR 50 49 % ouverture / lecture du fichier … … 158 157 end; 159 158 160 161 159 % on fait la somme des valeurs de hgt925_filtreentiel sur le nbre d occurences 162 160 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 232 230 lat=lat'; 233 231 234 235 232 for e=1:7; 236 233 a=7-(e-1)+1+7; -
trunk/src/mode_sahelien/composite_hgt925_eof23_2590_bis.m
r17 r18 46 46 clear hgt925; 47 47 48 49 48 %%%%%%%%%%%%%%% FICHIER OLR 50 49 % ouverture / lecture du fichier … … 158 157 end; 159 158 160 161 159 % on fait la somme des valeurs de hgt925_filtreentiel sur le nbre d occurences 162 160 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 280 278 281 279 print -depsc2 compo_hgt925_eof2590_part2_lat_moy.eps; 282 -
trunk/src/mode_sahelien/composite_hgt925_rossby.m
r17 r18 149 149 end; 150 150 151 152 151 % on fait la somme des valeurs de hgt700_filtreentiel sur le nbre d occurences 153 152 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 216 215 set (gca,'position',[0.2 pos_vert 0.5 0.1]); 217 216 end; 218 219 217 220 218 %%%%%%%%%%%%%%%%%%%% -
trunk/src/mode_sahelien/composite_hgt925_sahel.m
r17 r18 38 38 hgt=permute(hgt700,[1 4 2 3]); 39 39 clear hgt700; 40 41 40 42 41 %%%%%%%%%%%%%%% FICHIER OLR … … 162 161 end; 163 162 164 165 163 % on fait la somme des valeurs de hgt700_filtreentiel sur le nbre d occurences 166 164 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 235 233 lat=lat'; 236 234 237 238 235 for e=1:7; 239 236 a=7-(e-1)+1+7; -
trunk/src/mode_sahelien/composite_olr_EOF12_ER1030.m
r17 r18 51 51 52 52 % on met sous la forme longitude latitude jour 53 54 53 % on selectionne les mois d ete (attention aux annees bixestiles!) 55 54 … … 280 279 end 281 280 282 283 281 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 284 282 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/mode_sahelien/composite_olr_EOF12_ER1030_centre.m
r17 r18 278 278 end 279 279 280 281 280 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 282 281 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/mode_sahelien/composite_olr_EOF12_ER30100.m
r17 r18 114 114 olr(:,:,:,28)=olr1(:,:,siz28:siz28bis); % 2006 115 115 116 117 116 clear siz1 siz1bis; 118 117 clear siz2 siz2bis; … … 277 276 end 278 277 279 280 278 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 281 279 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 369 367 print -depsc2 composite_olr_EOF12_ER30100_part1.eps; 370 368 371 372 369 figure(2); 373 370 orient('landscape'); -
trunk/src/mode_sahelien/composite_olr_EOF12_ER30100_dat.m
r17 r18 51 51 52 52 % on met sous la forme longitude latitude jour 53 54 53 55 54 % on selectionne les mois d ete (attention aux annees bixestiles!) … … 273 272 end 274 273 end 275 276 274 277 275 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences -
trunk/src/mode_sahelien/composite_olr_EOF1_1030.m
r17 r18 85 85 siz28=siz27+365; siz28bis=siz28+a; % 2000 86 86 87 88 87 olr(:,:,:,1)=olr1(:,:,siz1:siz1bis); % 1979 89 88 olr(:,:,:,2)=olr1(:,:,siz2:siz2bis); % 1980 … … 268 267 end 269 268 270 271 269 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 272 270 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/mode_sahelien/composite_olr_EOF1_30100.m
r17 r18 85 85 siz28=siz27+365; siz28bis=siz28+a; % 2000 86 86 87 88 87 olr(:,:,:,1)=olr1(:,:,siz1:siz1bis); % 1979 89 88 olr(:,:,:,2)=olr1(:,:,siz2:siz2bis); % 1980 … … 115 114 olr(:,:,:,28)=olr1(:,:,siz28:siz28bis); % 2006 116 115 117 118 116 clear siz1 siz1bis; 119 117 clear siz2 siz2bis; … … 269 267 end 270 268 271 272 269 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 273 270 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/mode_sahelien/composite_olr_EOF23_1030.m
r17 r18 85 85 siz28=siz27+365; siz28bis=siz28+a; % 2000 86 86 87 88 87 olr(:,:,:,1)=olr1(:,:,siz1:siz1bis); % 1979 89 88 olr(:,:,:,2)=olr1(:,:,siz2:siz2bis); % 1980 … … 115 114 olr(:,:,:,28)=olr1(:,:,siz28:siz28bis); % 2006 116 115 117 118 116 clear siz1 siz1bis; 119 117 clear siz2 siz2bis; … … 277 275 end 278 276 279 280 277 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 281 278 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/mode_sahelien/composite_olr_EOF23_30100.m
r17 r18 85 85 siz28=siz27+365; siz28bis=siz28+a; % 2000 86 86 87 88 87 olr(:,:,:,1)=olr1(:,:,siz1:siz1bis); % 1979 89 88 olr(:,:,:,2)=olr1(:,:,siz2:siz2bis); % 1980 … … 115 114 olr(:,:,:,28)=olr1(:,:,siz28:siz28bis); % 2006 116 115 117 118 116 clear siz1 siz1bis; 119 117 clear siz2 siz2bis; … … 269 267 end 270 268 271 272 269 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 273 270 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/mode_sahelien/composite_olr_EOF24_1030.m
r17 r18 85 85 siz28=siz27+365; siz28bis=siz28+a; % 2000 86 86 87 88 87 olr(:,:,:,1)=olr1(:,:,siz1:siz1bis); % 1979 89 88 olr(:,:,:,2)=olr1(:,:,siz2:siz2bis); % 1980 … … 115 114 olr(:,:,:,28)=olr1(:,:,siz28:siz28bis); % 2006 116 115 117 118 116 clear siz1 siz1bis; 119 117 clear siz2 siz2bis; … … 159 157 olr_moy=mean(olrbis); 160 158 161 162 159 for a=1:k; 163 160 olrbis(a)=(olrbis(a)-olr_moy)/olr_std; … … 279 276 end; 280 277 281 282 278 % on fait la somme des valeurs de olr_filtre sur le nbre d occurences 283 279 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 344 340 axis([-180 180 -20 40]); 345 341 346 347 342 if e==7; 348 title('OLR composite EOF234 1030j ITCZ 7.5-12.5 ','fontsize',8) ;343 title('OLR composite EOF234 1030j ITCZ 7.5-12.5 ','fontsize',8) ; 349 344 end 350 345 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 409 404 axis([-180 180 -20 40]); 410 405 411 412 406 if e==7; 413 407 title('OLR composite EOF234 1030j ITCZ 7.5-12.5 ','fontsize',8) ; -
trunk/src/mode_sahelien/composite_olr_EOF24_1030_passbande.m
r17 r18 85 85 siz28=siz27+365; siz28bis=siz28+a; % 2000 86 86 87 88 87 olr(:,:,:,1)=olr1(:,:,siz1:siz1bis); % 1979 89 88 olr(:,:,:,2)=olr1(:,:,siz2:siz2bis); % 1980 … … 115 114 olr(:,:,:,28)=olr1(:,:,siz28:siz28bis); % 2006 116 115 117 118 116 clear siz1 siz1bis; 119 117 clear siz2 siz2bis; … … 169 167 olr_moy=mean(olrbis); 170 168 171 172 169 for a=1:k; 173 170 olrbis(a)=(olrbis(a)-olr_moy)/olr_std; … … 289 286 end; 290 287 291 292 288 % on fait la somme des valeurs de olr_filtre sur le nbre d occurences 293 289 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 354 350 axis([-180 180 -20 40]); 355 351 356 357 352 if e==7; 358 353 title('OLR composite EOF234 1030j ITCZ 7.5-12.5 ','fontsize',8) ; … … 419 414 axis([-180 180 -20 40]); 420 415 421 422 416 if e==7; 423 417 title('OLR composite EOF234 1030j ITCZ 7.5-12.5 ','fontsize',8) ; -
trunk/src/mode_sahelien/composite_olr_ER1030.m
r17 r18 85 85 siz28=siz27+365; siz28bis=siz28+a; % 2006 86 86 87 88 89 90 87 olr(:,:,:,1)=olr1(:,:,siz1:siz1bis); % 1979 91 88 olr(:,:,:,2)=olr1(:,:,siz2:siz2bis); % 1980 … … 271 268 end 272 269 273 274 270 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 275 271 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/mode_sahelien/composite_olr_ER1030_filtre.m
r17 r18 154 154 end; 155 155 156 157 158 156 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 159 157 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 181 179 fin=-30+24*2.5; 182 180 183 184 181 clear fin; 185 182 fin=-10+16*2.5; … … 194 191 clabel(cs,h,[-30 -20 -10],'color','k','rotation',0,'fontsize',5); 195 192 196 197 193 cartemonde1; 198 194 -
trunk/src/mode_sahelien/composite_olr_eof234_1030_filtre.m
r17 r18 1 1 % 2 3 2 4 3 %%%%%%%%%%%%%%% FICHIER OLR … … 22 21 23 22 load eof24_1030.txt 24 25 23 26 24 % on met sous la forme annee mois … … 69 67 compteurmin_eof234_1030=compteurmin; 70 68 71 72 69 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 73 70 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 97 94 fin=-30+24*2.5; 98 95 99 100 96 clear fin; 101 97 fin=-10+16*2.5; -
trunk/src/mode_sahelien/composite_olr_eof23_1030_filtre.m
r17 r18 21 21 22 22 load eof23_1030.txt 23 24 23 25 24 % on met sous la forme annee mois … … 68 67 compteurmax_eof23_1030=compteurmax; 69 68 70 71 69 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 72 70 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 96 94 fin=-30+24*2.5; 97 95 98 99 96 clear fin; 100 97 fin=-10+16*2.5; -
trunk/src/mode_sahelien/composite_olr_phase1.m
r17 r18 86 86 siz27=siz26+365; siz27bis=siz27+a; % 2005 87 87 siz28=siz27+365; siz28bis=siz28+a; % 2006 88 89 90 91 88 92 89 olr(:,:,:,1)=olr1(:,:,siz1:siz1bis); % 1979 … … 182 179 clear olr_filtre_compomax; 183 180 184 185 181 x=[-10:1:10]'; 186 182 -
trunk/src/mode_sahelien/composite_olr_phase2.m
r17 r18 86 86 siz27=siz26+365; siz27bis=siz27+a; % 2005 87 87 siz28=siz27+365; siz28bis=siz28+a; % 2006 88 89 90 91 88 92 89 olr(:,:,:,1)=olr1(:,:,siz1:siz1bis); % 1979 … … 181 178 clear olr_filtre_compomax; 182 179 183 184 180 x=[-10:1:10]'; 185 181 -
trunk/src/mode_sahelien/composite_olr_phase3.m
r17 r18 86 86 siz27=siz26+365; siz27bis=siz27+a; % 2005 87 87 siz28=siz27+365; siz28bis=siz28+a; % 2006 88 89 90 91 88 92 89 olr(:,:,:,1)=olr1(:,:,siz1:siz1bis); % 1979 … … 182 179 clear olr_filtre_compomax; 183 180 184 185 181 x=[-10:1:10]'; 186 182 -
trunk/src/mode_sahelien/composite_olr_phase4.m
r17 r18 86 86 siz27=siz26+365; siz27bis=siz27+a; % 2005 87 87 siz28=siz27+365; siz28bis=siz28+a; % 2006 88 89 90 91 88 92 89 olr(:,:,:,1)=olr1(:,:,siz1:siz1bis); % 1979 … … 181 178 clear olr_filtre_compomax; 182 179 183 184 180 x=[-10:1:10]'; 185 181 -
trunk/src/mode_sahelien/composite_olr_phase5.m
r17 r18 86 86 siz27=siz26+365; siz27bis=siz27+a; % 2005 87 87 siz28=siz27+365; siz28bis=siz28+a; % 2006 88 89 90 91 88 92 89 olr(:,:,:,1)=olr1(:,:,siz1:siz1bis); % 1979 -
trunk/src/mode_sahelien/composite_olr_phase6.m
r17 r18 86 86 siz27=siz26+365; siz27bis=siz27+a; % 2005 87 87 siz28=siz27+365; siz28bis=siz28+a; % 2006 88 89 90 91 88 92 89 olr(:,:,:,1)=olr1(:,:,siz1:siz1bis); % 1979 -
trunk/src/mode_sahelien/composite_olr_phase7.m
r17 r18 86 86 siz27=siz26+365; siz27bis=siz27+a; % 2005 87 87 siz28=siz27+365; siz28bis=siz28+a; % 2006 88 89 90 91 88 92 89 olr(:,:,:,1)=olr1(:,:,siz1:siz1bis); % 1979 … … 181 178 clear olr_filtre_compomax; 182 179 183 184 180 x=[-10:1:10]'; 185 181 -
trunk/src/mode_sahelien/composite_olr_phase8.m
r17 r18 86 86 siz27=siz26+365; siz27bis=siz27+a; % 2005 87 87 siz28=siz27+365; siz28bis=siz28+a; % 2006 88 89 90 91 88 92 89 olr(:,:,:,1)=olr1(:,:,siz1:siz1bis); % 1979 … … 182 179 clear olr_filtre_compomax; 183 180 184 185 181 x=[-10:1:10]'; 186 182 -
trunk/src/mode_sahelien/composite_olr_sahel_1025.m
r17 r18 85 85 siz28=siz27+365; siz28bis=siz28+a; % 2000 86 86 87 88 87 olr(:,:,:,1)=olr1(:,:,siz1:siz1bis); % 1979 89 88 olr(:,:,:,2)=olr1(:,:,siz2:siz2bis); % 1980 … … 115 114 olr(:,:,:,28)=olr1(:,:,siz28:siz28bis); % 2006 116 115 117 118 116 clear siz1 siz1bis; 119 117 clear siz2 siz2bis; … … 270 268 end 271 269 272 273 270 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 274 271 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 344 341 print -depsc2 composite_olr_jjas_sahel_10_25_eof234_part1.eps; 345 342 346 347 343 figure(2); 348 344 orient('landscape'); -
trunk/src/mode_sahelien/composite_olr_sahel_pb_3j.m
r17 r18 86 86 siz28=siz27+365; siz28bis=siz28+a; % 2000 87 87 88 89 88 olr(:,:,:,1)=olr1(:,:,siz1:siz1bis); % 1979 90 89 olr(:,:,:,2)=olr1(:,:,siz2:siz2bis); % 1980 … … 115 114 olr(:,:,:,27)=olr1(:,:,siz27:siz27bis); % 2005 116 115 olr(:,:,:,28)=olr1(:,:,siz28:siz28bis); % 2006 117 118 116 119 117 clear siz1 siz1bis; … … 292 290 olr_filtre_compomax(3,compteurmax,31,:,:)=olr(b+11,a,:,:); 293 291 294 295 292 % cas 2 : min deltas apres positifs et deltas avant negatifs 296 293 elseif (delta(b+1)>0 && delta(b)<0 && ind_olr(b,a)<-1.50); … … 393 390 olr_filtre_compomin(3,compteurmin,31,:,:)=olr(b+10,a,:,:); 394 391 395 396 392 end; 397 393 end; … … 451 447 set (gca,'YTick',[-10:10:40],'YTickLabel',['10S';' 0';'10N';'20N';'30N';'40N'],'fontsize',6); 452 448 453 454 449 if e==4; 455 450 co=colorbar; -
trunk/src/mode_sahelien/composite_pluie_eof234_1030.m
r17 r18 104 104 end; 105 105 end; 106 107 106 108 107 pluie(:,:,:,5)=v(:,:,siz1:siz1bis); % 1983 … … 283 282 end; 284 283 285 286 284 % on fait la somme des valeurs de pluie_filtreentiel sur le nbre d occurences 287 285 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/mode_sahelien/composite_pluie_eof23_1030.m
r17 r18 98 98 end; 99 99 end; 100 101 100 102 101 pluie(:,:,:,5)=v(:,:,siz1:siz1bis); % 1983 … … 276 275 end; 277 276 278 279 277 % on fait la somme des valeurs de pluie_filtreentiel sur le nbre d occurences 280 278 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/mode_sahelien/composite_pluie_ird_rossby.m
r17 r18 96 96 end; 97 97 end; 98 99 98 100 99 pluie(:,:,:,5)=v(:,:,siz1:siz1bis); % 1983 … … 334 333 end; 335 334 336 337 335 % on fait la somme des valeurs de pluie_filtreentiel sur le nbre d occurences 338 336 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 459 457 pos_vert=e/10; 460 458 set (gca,'position',[0.2 pos_vert 0.5 0.1]); 461 462 463 459 464 460 end; -
trunk/src/mode_sahelien/composite_pluie_ird_rossby1030.m
r17 r18 96 96 end; 97 97 end; 98 99 98 100 99 pluie(:,:,:,5)=v(:,:,siz1:siz1bis); % 1983 … … 336 335 end; 337 336 338 339 337 % on fait la somme des valeurs de pluie_filtreentiel sur le nbre d occurences 340 338 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 461 459 pos_vert=e/10; 462 460 set (gca,'position',[0.2 pos_vert 0.5 0.1]); 463 464 465 461 466 462 end; -
trunk/src/mode_sahelien/composite_pluie_ird_rossby1030b.m
r17 r18 95 95 end; 96 96 end; 97 98 97 99 98 pluie(:,:,:,5)=v(:,:,siz1:siz1bis); % 1983 … … 274 273 end; 275 274 276 277 275 % on fait la somme des valeurs de pluie_filtreentiel sur le nbre d occurences 278 276 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/mode_sahelien/composite_pluie_ird_rossby30100b.m
r17 r18 95 95 end; 96 96 end; 97 98 97 99 98 pluie(:,:,:,5)=v(:,:,siz1:siz1bis); % 1983 … … 274 273 end; 275 274 276 277 275 % on fait la somme des valeurs de pluie_filtreentiel sur le nbre d occurences 278 276 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/mode_sahelien/composite_pluie_ird_sahel.m
r17 r18 96 96 end; 97 97 end; 98 99 98 100 99 pluie(:,:,:,5)=v(:,:,siz1:siz1bis); % 1983 … … 245 244 ind_olr2=std(ind_olr1); 246 245 clear ind_olr1; 247 248 249 246 250 247 % intialisation des tableaux … … 343 340 end; 344 341 345 346 342 % on fait la somme des valeurs de pluie_filtreentiel sur le nbre d occurences 347 343 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 470 466 set (gca,'position',[0.2 pos_vert 0.5 0.1]); 471 467 472 473 474 468 end; 475 469 figure(1); … … 480 474 print -depsc2 composite_pluie_diff_mode_sahel1.eps; 481 475 close all 482 483 %*******************************484 485 476 486 477 figure(1); … … 603 594 set (gca,'position',[0.2 pos_vert 0.5 0.1]); 604 595 605 606 607 596 end; 608 597 figure(1); -
trunk/src/mode_sahelien/composite_uwnd700_EOF12_ER1030.m
r17 r18 4 4 clear; 5 5 close all; 6 7 6 8 7 % etape 1 :on selectionne les mois d ete (attention aux annees bixestiles!) … … 49 48 clear vwnd700; 50 49 51 52 50 %%%%%%%%%%%%%%% FICHIER OLR 53 51 % ouverture / lecture du fichier … … 161 159 end; 162 160 163 164 161 % on fait la somme des valeurs de vwnd700_filtreentiel sur le nbre d occurences 165 162 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/mode_sahelien/composite_uwnd700_EOF12_ER30100_dat.m
r17 r18 4 4 clear; 5 5 close all; 6 7 6 8 7 % etape 1 :on selectionne les mois d ete (attention aux annees bixestiles!) … … 49 48 clear uwnd700; 50 49 51 52 50 %%%%%%%%%%%%%%% FICHIER OLR 53 51 % ouverture / lecture du fichier … … 161 159 end; 162 160 end; 163 164 161 165 162 % on fait la somme des valeurs de uwnd700_filtreentiel sur le nbre d occurences -
trunk/src/mode_sahelien/composite_uwnd700_EOF24_1030.m
r17 r18 4 4 clear; 5 5 close all; 6 7 6 8 7 % etape 1 :on selectionne les mois d ete (attention aux annees bixestiles!) … … 49 48 clear uwnd700; 50 49 51 52 50 %%%%%%%%%%%%%%% FICHIER OLR 53 51 % ouverture / lecture du fichier … … 160 158 end; 161 159 end; 162 163 160 164 161 % on fait la somme des valeurs de uwnd700_filtreentiel sur le nbre d occurences … … 297 294 delta_uwnd700=flipud(delta_uwnd700); 298 295 299 300 296 aa=contourf(lon-180,lat,delta_uwnd700,[-5:0.5:5],'LineStyle','none'); 301 297 if e==1; -
trunk/src/mode_sahelien/composite_uwnd925_EOF12_ER1030.m
r17 r18 48 48 clear vwnd925; 49 49 50 51 50 %%%%%%%%%%%%%%% FICHIER OLR 52 51 % ouverture / lecture du fichier … … 160 159 end; 161 160 162 163 161 % on fait la somme des valeurs de vwnd925_filtreentiel sur le nbre d occurences 164 162 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/mode_sahelien/composite_uwnd925_EOF12_ER30100_dat.m
r17 r18 49 49 clear uwnd925; 50 50 51 52 51 %%%%%%%%%%%%%%% FICHIER OLR 53 52 % ouverture / lecture du fichier … … 161 160 end; 162 161 end; 163 164 162 165 163 % on fait la somme des valeurs de uwnd925_filtreentiel sur le nbre d occurences -
trunk/src/mode_sahelien/composite_uwnd925_EOF23_1030.m
r17 r18 4 4 clear; 5 5 close all; 6 7 6 8 7 % etape 1 :on selectionne les mois d ete (attention aux annees bixestiles!) … … 49 48 clear uwnd925; 50 49 51 52 50 %%%%%%%%%%%%%%% FICHIER OLR 53 51 % ouverture / lecture du fichier … … 160 158 end; 161 159 end; 162 163 160 164 161 % on fait la somme des valeurs de uwnd925_filtreentiel sur le nbre d occurences … … 297 294 delta_uwnd925=flipud(delta_uwnd925); 298 295 299 300 296 aa=contourf(lon-180,lat,delta_uwnd925,[-5:0.5:5],'LineStyle','none'); 301 297 if e==1; -
trunk/src/mode_sahelien/composite_uwnd925_EOF24_1030.m
r17 r18 48 48 clear vwnd925; 49 49 50 51 50 %%%%%%%%%%%%%%% FICHIER OLR 52 51 % ouverture / lecture du fichier … … 159 158 end; 160 159 end; 161 162 160 163 161 % on fait la somme des valeurs de vwnd925_filtreentiel sur le nbre d occurences … … 296 294 delta_vwnd925=flipud(delta_vwnd925); 297 295 298 299 296 aa=contourf(lon-180,lat,delta_vwnd925,[-5:0.5:5],'LineStyle','none'); 300 297 if e==1; -
trunk/src/mode_sahelien/composite_vwnd700_EOF12_ER30100_dat.m
r17 r18 43 43 vwnd700=vwnd700-187.65; 44 44 45 46 45 clear siz1 siz1bis; 47 46 48 47 vwnd=permute(vwnd700,[1 4 2 3]); 49 48 clear vwnd700; 50 51 49 52 50 %%%%%%%%%%%%%%% FICHIER OLR … … 161 159 end; 162 160 end; 163 164 161 165 162 % on fait la somme des valeurs de vwnd700_filtreentiel sur le nbre d occurences -
trunk/src/mode_sahelien/composite_vwnd925_EOF12_ER30100_dat.m
r17 r18 42 42 vwnd925=vwnd925-187.65; 43 43 44 45 44 clear siz1 siz1bis; 46 45 47 46 vwnd=permute(vwnd925,[1 4 2 3]); 48 47 clear vwnd925; 49 50 48 51 49 %%%%%%%%%%%%%%% FICHIER OLR … … 160 158 end; 161 159 end; 162 163 160 164 161 % on fait la somme des valeurs de vwnd925_filtreentiel sur le nbre d occurences -
trunk/src/mode_sahelien/composite_vwnd925_EOF23_1030.m
r17 r18 47 47 clear vwnd925; 48 48 49 50 49 %%%%%%%%%%%%%%% FICHIER OLR 51 50 % ouverture / lecture du fichier … … 158 157 end; 159 158 end; 160 161 159 162 160 % on fait la somme des valeurs de vwnd925_filtreentiel sur le nbre d occurences … … 295 293 delta_vwnd925=flipud(delta_vwnd925); 296 294 297 298 295 aa=contourf(lon-180,lat,delta_vwnd925,[-5:0.5:5],'LineStyle','none'); 299 296 if e==1; -
trunk/src/mode_sahelien/composite_wind200_rossby.m
r17 r18 148 148 end; 149 149 150 151 150 % on fait la somme des valeurs de uwnd200_filtreentiel sur le nbre d occurences 152 151 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 299 298 end; 300 299 301 302 300 % on fait la somme des valeurs de vwnd200_filtreentiel sur le nbre d occurences 303 301 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 306 304 vwnd200_filtre_compomax1(:,:,:)=nanmean(vwnd200_filtre_compomax); 307 305 clear vwnd200_filtre_compomax; 308 309 306 310 307 figure(1) … … 371 368 end; 372 369 373 374 370 print -depsc2 composite_rossby_eof2a_u200_hgt.eps; 375 371 … … 436 432 end; 437 433 438 439 434 print -depsc2 composite_rossby_eof2b_u200_hgt.eps; -
trunk/src/mode_sahelien/composite_wind200_sahel.m
r17 r18 38 38 uwnd=permute(uwnd200,[1 4 2 3]); 39 39 clear uwnd200; 40 41 40 42 41 %%%%%%%%%%%%%%% FICHIER OLR … … 161 160 end; 162 161 163 164 162 % on fait la somme des valeurs de uwnd200_filtreentiel sur le nbre d occurences 165 163 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 168 166 uwnd200_filtre_compomax1(:,:,:)=nanmean(uwnd200_filtre_compomax); 169 167 clear uwnd200_filtre_compomax; 170 171 168 172 169 % etape 1 :on selectionne les mois d ete (attention aux annees bixestiles!) … … 225 222 ind_olr2=std(ind_olr1); 226 223 clear ind_olr1; 227 228 224 229 225 % intialisation des tableaux … … 325 321 end; 326 322 end; 327 328 323 329 324 % on fait la somme des valeurs de vwnd200_filtreentiel sur le nbre d occurences … … 397 392 end; 398 393 399 400 394 print -depsc2 composite_sahel_1_u200_hgt.eps; 401 395 -
trunk/src/mode_sahelien/composite_wind500_rossby_europe.m
r17 r18 148 148 end; 149 149 150 151 150 % on fait la somme des valeurs de uwnd500_filtreentiel sur le nbre d occurences 152 151 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 298 297 end; 299 298 end; 300 301 299 302 300 % on fait la somme des valeurs de vwnd500_filtreentiel sur le nbre d occurences … … 370 368 end; 371 369 372 373 370 print -depsc2 composite_rossby_eof2a_u500_hgt.eps; 374 375 371 376 372 for e=1:7; … … 433 429 end; 434 430 435 436 431 print -depsc2 composite_rossby_eof2b_u500_hgt.eps; -
trunk/src/mode_sahelien/composite_wind925_EOF12_ER1030.m
r17 r18 41 41 ncload('uwnd925.2006.nc'); uwnd925(:,:,:,28)=uwnd(siz:sizbis,1,:,:);clear uwnd; 42 42 43 44 43 clear siz1 siz1bis; 45 44 46 45 uwnd=permute(uwnd925,[1 4 2 3]); 47 46 clear uwnd925; 48 49 47 50 48 %%%%%%%%%%%%%%% FICHIER OLR … … 60 58 ind_olr2=std(ind_olr1); 61 59 clear ind_olr1; 62 63 60 64 61 % intialisation des tableaux … … 160 157 end; 161 158 end; 162 163 159 164 160 % on fait la somme des valeurs de uwnd925_filtreentiel sur le nbre d occurences … … 222 218 clear ind_olr1; 223 219 224 225 220 % intialisation des tableaux 226 221 vwnd925_filtre_compomax=NaN*ones(60,31,49,144); … … 322 317 end; 323 318 324 325 319 % on fait la somme des valeurs de vwnd925_filtreentiel sur le nbre d occurences 326 320 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 383 377 end; 384 378 385 386 379 print -depsc2 composite_wind925_EOF12_ER1030_part1.eps; 387 380 … … 438 431 end; 439 432 440 441 433 print -depsc2 composite_wind925_EOF12_ER1030_part2.eps; 442 434 -
trunk/src/mode_sahelien/composite_wind925_EOF12_ER30100.m
r17 r18 41 41 ncload('uwnd925.2006.nc'); uwnd925(:,:,:,28)=uwnd(siz:sizbis,1,:,:);clear uwnd; 42 42 43 44 43 clear siz1 siz1bis; 45 44 46 45 uwnd=permute(uwnd925,[1 4 2 3]); 47 46 clear uwnd925; 48 49 47 50 48 %%%%%%%%%%%%%%% FICHIER OLR … … 60 58 ind_olr2=std(ind_olr1); 61 59 clear ind_olr1; 62 63 60 64 61 % intialisation des tableaux … … 160 157 end; 161 158 end; 162 163 159 164 160 % on fait la somme des valeurs de uwnd925_filtreentiel sur le nbre d occurences … … 322 318 end; 323 319 324 325 320 % on fait la somme des valeurs de vwnd925_filtreentiel sur le nbre d occurences 326 321 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) -
trunk/src/mode_sahelien/composite_wind925_EOF1_30100.m
r17 r18 41 41 ncload('uwnd925.2006.nc'); uwnd925(:,:,:,28)=uwnd(siz:sizbis,1,:,:);clear uwnd; 42 42 43 44 43 clear siz1 siz1bis; 45 44 46 45 uwnd=permute(uwnd925,[1 4 2 3]); 47 46 clear uwnd925; 48 49 47 50 48 %%%%%%%%%%%%%%% FICHIER OLR … … 60 58 ind_olr2=std(ind_olr1); 61 59 clear ind_olr1; 62 63 60 64 61 % intialisation des tableaux … … 160 157 end; 161 158 end; 162 163 159 164 160 % on fait la somme des valeurs de uwnd925_filtreentiel sur le nbre d occurences … … 322 318 end; 323 319 324 325 320 % on fait la somme des valeurs de vwnd925_filtreentiel sur le nbre d occurences 326 321 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 384 379 end; 385 380 386 387 381 print -depsc2 composite_wind925_EOF1_ER30100_part1.eps; 388 382 389 383 figure(2); 390 384 orient('landscape'); 391 392 385 393 386 for e=1:4; … … 440 433 end; 441 434 442 443 435 print -depsc2 composite_wind925_EOF1_ER30100_part2.eps; 444 436 … … 495 487 end; 496 488 497 498 489 print -depsc2 composite_wind925_EOF1_ER30100_part3.eps; 499 490 -
trunk/src/mode_sahelien/composite_wind925_EOF234_1030.m
r17 r18 41 41 ncload('uwnd925.2006.nc'); uwnd925(:,:,:,28)=uwnd(siz:sizbis,1,:,:);clear uwnd; 42 42 43 44 43 clear siz1 siz1bis; 45 44 46 45 uwnd=permute(uwnd925,[1 4 2 3]); 47 46 clear uwnd925; 48 49 47 50 48 %%%%%%%%%%%%%%% FICHIER OLR … … 60 58 ind_olr2=std(ind_olr1); 61 59 clear ind_olr1; 62 63 60 64 61 % intialisation des tableaux … … 160 157 end; 161 158 end; 162 163 159 164 160 % on fait la somme des valeurs de uwnd925_filtreentiel sur le nbre d occurences … … 222 218 clear ind_olr1; 223 219 224 225 220 % intialisation des tableaux 226 221 vwnd925_filtre_compomax=NaN*ones(60,31,49,144); … … 322 317 end; 323 318 324 325 319 % on fait la somme des valeurs de vwnd925_filtreentiel sur le nbre d occurences 326 320 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 384 378 end; 385 379 386 387 380 print -depsc2 composite_wind925_EOF234_1030_part1.eps; 388 381 389 382 figure(2); 390 383 orient('landscape'); 391 392 384 393 385 for e=1:4; … … 440 432 end; 441 433 442 443 434 print -depsc2 composite_wind925_EOF234_1030_part2.eps; 444 435 -
trunk/src/mode_sahelien/composite_wind925_EOF23_30100.m
r17 r18 41 41 ncload('uwnd925.2006.nc'); uwnd925(:,:,:,28)=uwnd(siz:sizbis,1,:,:);clear uwnd; 42 42 43 44 43 clear siz1 siz1bis; 45 44 46 45 uwnd=permute(uwnd925,[1 4 2 3]); 47 46 clear uwnd925; 48 49 47 50 48 %%%%%%%%%%%%%%% OLR … … 60 58 ind_olr2=std(ind_olr1); 61 59 clear ind_olr1; 62 63 60 64 61 % intialisation des tableaux … … 160 157 end; 161 158 end; 162 163 159 164 160 % on fait la somme des valeurs de uwnd925_filtreentiel sur le nbre d occurences … … 222 218 clear ind_olr1; 223 219 224 225 220 % intialisation des tableaux 226 221 vwnd925_filtre_compomax=NaN*ones(60,31,49,144); … … 322 317 end; 323 318 324 325 319 % on fait la somme des valeurs de vwnd925_filtreentiel sur le nbre d occurences 326 320 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 384 378 end; 385 379 386 387 380 print -depsc2 composite_wind925_EOF23_30100_part1.eps; 388 381 389 382 figure(2); 390 383 orient('landscape'); 391 392 384 393 385 for e=1:4; … … 440 432 end; 441 433 442 443 434 print -depsc2 composite_wind925_EOF23_30100_part2.eps; 444 435 … … 495 486 end; 496 487 497 498 488 print -depsc2 composite_wind925_EOF23_30100_part3.eps; 499 489 -
trunk/src/mode_sahelien/composite_wind925_eof23_2590.m
r17 r18 41 41 ncload('uwnd925.2006.nc'); uwnd925(:,:,:,28)=uwnd(siz:sizbis,1,:,:);clear uwnd; 42 42 43 44 43 clear siz1 siz1bis; 45 44 46 45 uwnd=permute(uwnd925,[1 4 2 3]); 47 46 clear uwnd925; 48 49 47 50 48 %%%%%%%%%%%%%%% FICHIER OLR … … 59 57 ind_olr2=std(ind_olr1); 60 58 clear ind_olr1; 61 62 59 63 60 % intialisation des tableaux … … 159 156 end; 160 157 end; 161 162 158 163 159 % on fait la somme des valeurs de uwnd925_filtreentiel sur le nbre d occurences … … 220 216 clear ind_olr1; 221 217 222 223 218 % intialisation des tableaux 224 219 vwnd925_filtre_compomax=NaN*ones(60,31,49,144); … … 320 315 end; 321 316 322 323 317 % on fait la somme des valeurs de vwnd925_filtreentiel sur le nbre d occurences 324 318 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 327 321 vwnd925_filtre_compomax1(:,:,:)=nanmean(vwnd925_filtre_compomax); 328 322 clear vwnd925_filtre_compomax; 329 330 323 331 324 figure(1); … … 393 386 end; 394 387 395 396 388 print -depsc2 composite_eof23_2090_u925_hgt_part1.eps; 397 389 … … 460 452 end; 461 453 462 463 454 print -depsc2 composite_eof23_2090_u925_hgt_part2.eps; 464 455 … … 527 518 end; 528 519 529 530 520 print -depsc2 composite_eof23_2090_u925_hgt_part3.eps; 531 521 -
trunk/src/mode_sahelien/composite_wind925_rossby.m
r17 r18 148 148 end; 149 149 150 151 150 % on fait la somme des valeurs de uwnd925_filtreentiel sur le nbre d occurences 152 151 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 298 297 end; 299 298 end; 300 301 299 302 300 % on fait la somme des valeurs de vwnd925_filtreentiel sur le nbre d occurences … … 370 368 end; 371 369 372 373 370 print -depsc2 composite_rossby_eof2a_u925_hgt.eps; 374 371 -
trunk/src/mode_sahelien/composite_wind925_sahel.m
r17 r18 39 39 clear uwnd700; 40 40 41 42 41 %%%%%%%%%%%%%%% FICHIER OLR 43 42 % ouverture / lecture du fichier … … 61 60 ind_olr2=std(ind_olr1); 62 61 clear ind_olr1; 63 64 62 65 63 % intialisation des tableaux … … 162 160 end; 163 161 164 165 162 % on fait la somme des valeurs de uwnd700_filtreentiel sur le nbre d occurences 166 163 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) … … 324 321 end; 325 322 end; 326 327 323 328 324 % on fait la somme des valeurs de vwnd700_filtreentiel sur le nbre d occurences -
trunk/src/mode_sahelien/lagged_correlations_eof2_eof3_1030.m
r16 r18 93 93 set (gca,'tick','out'); 94 94 95 96 95 figure(1); 97 96 print -depsc2 lag_correlations_eof2_eof3_1030.eps; -
trunk/src/mode_sahelien/olr_eof_jjas2006_rossby.m
r17 r18 48 48 siz27=siz26+365; siz27bis=siz27+a; % 2005 49 49 siz28=siz27+365; siz28bis=siz28+a; % 2006 50 51 50 52 51 olr(:,1,:,:)=ind_olr(:,:,siz1:siz1bis); % 1979 -
trunk/src/mode_sahelien/olr_eof_jjas2006_rossby_1025.m
r17 r18 15 15 % on selectionne les mois d ete (attention aux annees bixestiles!) 16 16 17 17 siz1=31+28+31+30+31+1; 18 18 siz1bis=31+28+31+30+31+30+31+31+30; 19 19 a=siz1bis-siz1; … … 48 48 siz27=siz26+365; siz27bis=siz27+a; % 2005 49 49 siz28=siz27+365; siz28bis=siz28+a; % 2006 50 51 50 52 51 olr(:,1,:,:)=ind_olr(:,:,siz1:siz1bis); % 1979 -
trunk/src/mode_sahelien/olr_eof_jjas2006_rossby_1030.m
r17 r18 15 15 % on selectionne les mois d ete (attention aux annees bixestiles!) 16 16 17 17 siz1=31+28+31+30+31+1; 18 18 siz1bis=31+28+31+30+31+30+31+31+30; 19 19 a=siz1bis-siz1; … … 48 48 siz27=siz26+365; siz27bis=siz27+a; % 2005 49 49 siz28=siz27+365; siz28bis=siz28+a; % 2006 50 51 50 52 51 olr(:,1,:,:)=ind_olr(:,:,siz1:siz1bis); % 1979 -
trunk/src/mode_sahelien/olr_eof_jjas2006_rossby_2590.m
r17 r18 15 15 % on selectionne les mois d ete (attention aux annees bixestiles!) 16 16 17 17 siz1=31+28+31+30+31+1; 18 18 siz1bis=31+28+31+30+31+30+31+31+30; 19 19 a=siz1bis-siz1; … … 48 48 siz27=siz26+365; siz27bis=siz27+a; % 2005 49 49 siz28=siz27+365; siz28bis=siz28+a; % 2006 50 51 50 52 51 olr(:,1,:,:)=ind_olr(:,:,siz1:siz1bis); % 1979 -
trunk/src/mode_sahelien/olr_eof_jjas2006_rossby_30100.m
r17 r18 15 15 % on selectionne les mois d ete (attention aux annees bixestiles!) 16 16 17 17 siz1=31+28+31+30+31+1; 18 18 siz1bis=31+28+31+30+31+30+31+31+30; 19 19 a=siz1bis-siz1; … … 48 48 siz27=siz26+365; siz27bis=siz27+a; % 2005 49 49 siz28=siz27+365; siz28bis=siz28+a; % 2006 50 51 50 52 51 olr(:,1,:,:)=ind_olr(:,:,siz1:siz1bis); % 1979 -
trunk/src/mode_sahelien/olr_eof_jjas_2006.m
r17 r18 14 14 % on selectionne les mois d ete (attention aux annees bixestiles!) 15 15 16 16 siz1=31+28+31+30+31+1; 17 17 siz1bis=31+28+31+30+31+30+31+31+30; 18 18 a=siz1bis-siz1; … … 47 47 siz27=siz26+365; siz27bis=siz27+a; % 2005 48 48 siz28=siz27+365; siz28bis=siz28+a; % 2006 49 50 49 51 50 olr(:,1,:,:)=ind_olr(:,:,siz1:siz1bis); % 1979 -
trunk/src/mode_sahelien/palette.m
r17 r18 25 25 palette2(1,:) = [0.9 0 0.9]; % Mauve 26 26 27 28 27 colormap(palette2); -
trunk/src/mode_sahelien/prepare_olr_total.m
r17 r18 13 13 olr(a,:,:)=coco; 14 14 end 15 16 17 15 18 16 cd('mode_sahelien'); … … 60 58 61 59 % on met sous la forme longitude latitude jour 62 63 60 % on selectionne les mois d ete (attention aux annees bixestiles!) 64 61 … … 94 91 siz28=siz27+365; siz28bis=siz28+a; % 2000 95 92 96 97 93 olr(:,:,:,1)=olr1(:,:,siz1:siz1bis); % 1979 98 94 olr(:,:,:,2)=olr1(:,:,siz2:siz2bis); % 1980 … … 123 119 olr(:,:,:,27)=olr1(:,:,siz27:siz27bis); % 2005 124 120 olr(:,:,:,28)=olr1(:,:,siz28:siz28bis); % 2006 125 126 121 127 122 clear siz1 siz1bis; -
trunk/src/mode_sahelien/read_pluie_ird.m
r17 r18 96 96 end; 97 97 end; 98 99 98 100 99 pluie(:,:,:,5)=v(:,:,siz1:siz1bis); % 1983 -
trunk/src/mode_sahelien/read_pluie_ird_25.m
r17 r18 96 96 end; 97 97 end; 98 99 98 100 99 pluie(:,:,:,5)=v(:,:,siz1:siz1bis); % 1983 -
trunk/src/mode_sahelien/reconstitution.m
r17 r18 67 67 cosa3=reshape(cosa3,122,22); 68 68 69 70 71 69 % On selectionne la demaine geographique 72 70 olr_centre=reshape(olr_centre,2684,17,25); … … 129 127 end; 130 128 131 132 133 129 figure(1); 134 130 print -depsc2 reconstitution_OLR_JJAS_part1.eps; -
trunk/src/mode_sahelien/reconstitution_1025.m
r17 r18 71 71 cosa24=reshape(cosa24,122,28); 72 72 73 74 73 fid=fopen('sahel_scav.dat','r'); 75 74 v=fread(fid,3538,'float'); -
trunk/src/mode_sahelien/reconstitution_1030.m
r17 r18 71 71 cosa24=reshape(cosa24,122,28); 72 72 73 74 75 76 73 % On selectionne la demaine geographique 77 74 olr=reshape(olr,3416,25,17); -
trunk/src/mode_sahelien/reconstitution_1030_std.m
r17 r18 65 65 cosa24=reshape(cosa24,122,28); 66 66 67 68 67 for a=1:28; 69 68 if a<12; -
trunk/src/mode_sahelien/reconstitution_30100_eof1.m
r17 r18 22 22 cosa1(:,:)=cosa(1,:,:); 23 23 cosa1=reshape(cosa1,3416,25,17); 24 25 24 26 25 cosa2(:,:)=cosa(2,:,:); … … 48 47 lon2=17; 49 48 taille=(lat2-lat1+1)*(lon2-lon1+1); 50 51 49 52 50 cosa1=cosa1(:,lon1:lon2,lat1:lat2); -
trunk/src/mode_sahelien/reconstitution_30100_eof23.m
r17 r18 44 44 lon2=17; 45 45 taille=(lat2-lat1+1)*(lon2-lon1+1); 46 47 46 48 47 cosa23=cosa23(:,lon1:lon2,lat1:lat2); -
trunk/src/mode_sahelien/reconstitution_ER_rossby_30100.m
r17 r18 47 47 lon2=17; 48 48 taille=(lat2-lat1+1)*(lon2-lon1+1); 49 50 49 51 50 cosa12=cosa12(:,lon1:lon2,lat1:lat2); -
trunk/src/mode_sahelien/reconstitution_ER_rossby_30100_bis.m
r17 r18 34 34 35 35 cosa_all=reshape(olr,3416,25,17); 36 37 36 38 37 % On selectionne la demaine geographique … … 75 74 cosa_all=reshape(cosa_all,122,28); 76 75 77 78 76 olr=reshape(olr,3416,25,17); 79 77 cosa4=olr(:,lon1:lon2,lat1:lat2); -
trunk/src/mode_sahelien/time_serie_ER_rossby_1030.m
r17 r18 56 56 axis([1 122 -100 100]); 57 57 58 59 58 set (gca,'YTick',[-100:20:100],'YTickLabel',[-100:20:100],'fontname','Arial','fontsize',6); 60 59 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); -
trunk/src/mode_sahelien/time_serie_ER_rossby_1030_complet.m
r17 r18 10 10 figure(3); 11 11 orient('landscape') 12 13 12 14 13 load eof_rossby_jjas_2006_1030.mat … … 55 54 axis([1 122 -100 100]); 56 55 57 58 56 set (gca,'YTick',[-100:20:100],'YTickLabel',[-100:20:100],'fontname','Arial','fontsize',6); 59 57 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); -
trunk/src/mode_sahelien/time_serie_ER_rossby_30100.m
r17 r18 56 56 axis([1 122 -100 100]); 57 57 58 59 58 set (gca,'YTick',[-100:20:100],'YTickLabel',[-100:20:100],'fontname','Arial','fontsize',6); 60 59 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); 61 -
trunk/src/mode_sahelien/time_serie_ER_rossby_30100_complet.m
r17 r18 10 10 figure(3); 11 11 orient('landscape') 12 13 12 14 13 load eof_rossby_jjas_2006_30100.mat … … 55 54 axis([1 122 -100 100]); 56 55 57 58 56 set (gca,'YTick',[-100:20:100],'YTickLabel',[-100:20:100],'fontname','Arial','fontsize',6); 59 57 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); -
trunk/src/mode_sahelien/time_serie_histog_phase1.m
r17 r18 73 73 sum_kara_max=sum((type_temps_compomax1)'); 74 74 75 76 75 % on traduit le nbre d occurence en frequences relatives 77 76 kara_max=zeros(21,4); -
trunk/src/mode_sahelien/time_serie_histog_phase2.m
r17 r18 73 73 sum_kara_max=sum((type_temps_compomax1)'); 74 74 75 76 75 % on traduit le nbre d occurence en frequences relatives 77 76 kara_max=zeros(21,4); -
trunk/src/mode_sahelien/time_serie_histog_phase4.m
r17 r18 91 91 end 92 92 93 94 93 for a=1:4; 95 94 subplot(4,4,a+12); -
trunk/src/mode_sahelien/time_serie_histog_phase5.m
r17 r18 91 91 end 92 92 93 94 93 for a=1:4; 95 94 subplot(4,4,a); -
trunk/src/mode_sahelien/time_serie_histog_phase6.m
r17 r18 91 91 end 92 92 93 94 93 for a=1:4; 95 94 subplot(4,4,a+4); -
trunk/src/mode_sahelien/time_serie_histog_phase8.m
r17 r18 73 73 sum_kara_max=sum((type_temps_compomax1)'); 74 74 75 76 75 % on traduit le nbre d occurence en frequences relatives 77 76 kara_max=zeros(21,4); … … 91 90 end 92 91 end 93 94 92 95 93 for a=1:4; -
trunk/src/mode_sahelien/time_serie_jjas.m
r17 r18 6 6 olr=olr'; 7 7 PC1=olr*C; 8 9 8 10 9 %%%%%%% cronique de la VP1 … … 61 60 axis([1 122 -3 3]); 62 61 63 64 62 set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); 65 63 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); -
trunk/src/mode_sahelien/time_serie_jjas_10_25.m
r17 r18 6 6 olr=olr'; 7 7 PC1=olr*C; 8 9 8 10 9 %%%%%%% cronique de la VP1 … … 60 59 axis([1 122 -3 3]); 61 60 62 63 61 set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); 64 62 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); -
trunk/src/mode_sahelien/time_serie_jjas_25_90.m
r17 r18 6 6 olr=olr'; 7 7 PC1=olr*C; 8 9 8 10 9 %%%%%%% cronique de la VP1 … … 61 60 axis([1 122 -3 3]); 62 61 63 64 62 set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); 65 63 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); -
trunk/src/mode_sahelien/time_serie_jjas_30100.m
r17 r18 6 6 olr=olr'; 7 7 PC1=olr*C; 8 9 8 10 9 %%%%%%% cronique de la VP1 … … 61 60 axis([1 122 -3 3]); 62 61 63 64 62 set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); 65 63 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); -
trunk/src/sauvegrads.m
r16 r18 89 89 fclose(fid); 90 90 91 92 91 % init des variables 93 92 [sxp syp]=size(posit); … … 105 104 error('La latitude doit ??tre comprise entre -90 et 90') 106 105 end 107 108 106 109 107 % Init date du jour pour enregistrement
Note: See TracChangeset
for help on using the changeset viewer.