Changeset 21
- Timestamp:
- 01/06/09 16:22:54 (15 years ago)
- Location:
- trunk/src
- Files:
-
- 62 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/SIMULS_IRCAAM/anomalies_olr_AFR_1979_2000.m
r18 r21 111 111 olr1(:,:,:,1:64)=olr(:,:,:,65:128); 112 112 olr1(:,:,:,65:128)=olr(:,:,:,1:64); 113 clear olr 113 clear olr; 114 114 olr=olr1; 115 115 clear olr1; … … 120 120 % 122 30 15 24 121 121 olr1=olr(31:122,:,:,:); 122 clear olr 122 clear olr; 123 123 olr(:,:,:)=mean(olr1); 124 124 clear olr1; -
trunk/src/SIMULS_IRCAAM/climato_olr_AFR.m
r16 r21 97 97 olr1(:,:,:,1:64)=olr(:,:,:,65:128); 98 98 olr1(:,:,:,65:128)=olr(:,:,:,1:64); 99 clear olr 99 clear olr; 100 100 olr=olr1; 101 101 clear olr1; -
trunk/src/SIMULS_IRCAAM/climato_olr_CTL.m
r16 r21 98 98 olr1(:,:,:,1:64)=olr(:,:,:,65:128); 99 99 olr1(:,:,:,65:128)=olr(:,:,:,1:64); 100 clear olr 100 clear olr; 101 101 olr=olr1; 102 102 clear olr1; -
trunk/src/SIMULS_IRCAAM/composite_hgt500_AFR_eof1_1030.m
r17 r21 73 73 if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && hgt500_eof(b,a)>(hgt500bis)); 74 74 compteurmax=compteurmax+1; 75 hgt500_filtre_compomax(compteurmax,1,:,:)=hgt500(b-20,a,:,:);75 hgt500_filtre_compomax(compteurmax,1,:,:)=hgt500(b-20,a,:,:); 76 76 hgt500_filtre_compomax(compteurmax,2,:,:)=hgt500(b-19,a,:,:); 77 77 hgt500_filtre_compomax(compteurmax,3,:,:)=hgt500(b-18,a,:,:); … … 200 200 if e==4; 201 201 co=colorbar; 202 set (co,'ytick',[-80:40:80],'yticklabel',[-80:40:80],'fontname','Arial','fontsize',5);202 set (co,'ytick',[-80:40:80],'yticklabel',[-80:40:80],'fontname','Arial','fontsize',5); 203 203 end; 204 204 hold on; … … 208 208 axis([-180 180 20 80]); 209 209 if e==7; 210 title('HGT500 composite EOF1 1030j simulation AFR 210 title('HGT500 composite EOF1 1030j simulation AFR','fontsize',8) ; 211 211 end 212 212 text(-270,40,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 255 255 if e==4; 256 256 co=colorbar; 257 set (co,'ytick',[-80:40:80],'yticklabel',[-80:40:80],'fontname','Arial','fontsize',5);257 set (co,'ytick',[-80:40:80],'yticklabel',[-80:40:80],'fontname','Arial','fontsize',5); 258 258 end; 259 259 hold on; … … 263 263 axis([-180 180 20 80]); 264 264 if e==7; 265 265 title('HGT500 composite EOF1 1030j simulation AFR','fontsize',8) ; 266 266 end 267 267 text(-270,40,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_hgt500_AFR_eof1_30100.m
r17 r21 73 73 if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && hgt500_eof(b,a)>(hgt500bis)); 74 74 compteurmax=compteurmax+1; 75 hgt500_filtre_compomax(compteurmax,1,:,:)=hgt500(b-20,a,:,:);75 hgt500_filtre_compomax(compteurmax,1,:,:)=hgt500(b-20,a,:,:); 76 76 hgt500_filtre_compomax(compteurmax,2,:,:)=hgt500(b-19,a,:,:); 77 77 hgt500_filtre_compomax(compteurmax,3,:,:)=hgt500(b-18,a,:,:); … … 200 200 if e==4; 201 201 co=colorbar; 202 set (co,'ytick',[-80:40:80],'yticklabel',[-80:40:80],'fontname','Arial','fontsize',5);202 set (co,'ytick',[-80:40:80],'yticklabel',[-80:40:80],'fontname','Arial','fontsize',5); 203 203 end; 204 204 hold on; … … 208 208 axis([-180 180 20 80]); 209 209 if e==7; 210 title('HGT500 composite EOF1 30100j simulation AFR 210 title('HGT500 composite EOF1 30100j simulation AFR','fontsize',8) ; 211 211 end 212 212 text(-270,40,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 255 255 if e==4; 256 256 co=colorbar; 257 set (co,'ytick',[-80:40:80],'yticklabel',[-80:40:80],'fontname','Arial','fontsize',5);257 set (co,'ytick',[-80:40:80],'yticklabel',[-80:40:80],'fontname','Arial','fontsize',5); 258 258 end; 259 259 hold on; … … 263 263 axis([-180 180 20 80]); 264 264 if e==7; 265 265 title('HGT500 composite EOF1 30100j simulation AFR','fontsize',8) ; 266 266 end 267 267 text(-270,40,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_hgt500_AFR_eof2_1030.m
r17 r21 73 73 if (delta(b+1)<0 && delta(b)>0 && hgt500_eof(b,a)>(hgt500bis)); 74 74 compteurmax=compteurmax+1; 75 hgt500_filtre_compomax(compteurmax,1,:,:)=hgt500(b-20,a,:,:);75 hgt500_filtre_compomax(compteurmax,1,:,:)=hgt500(b-20,a,:,:); 76 76 hgt500_filtre_compomax(compteurmax,2,:,:)=hgt500(b-19,a,:,:); 77 77 hgt500_filtre_compomax(compteurmax,3,:,:)=hgt500(b-18,a,:,:); … … 220 220 if e==4; 221 221 co=colorbar; 222 set (co,'ytick',[-50:5:50],'yticklabel',[-50:5:50],'fontname','Arial','fontsize',5);222 set (co,'ytick',[-50:5:50],'yticklabel',[-50:5:50],'fontname','Arial','fontsize',5); 223 223 end; 224 224 hold on; … … 228 228 axis([-180 180 20 80]); 229 229 if e==7; 230 title('HGT500 composite EOF2 1030j simulation AFR 230 title('HGT500 composite EOF2 1030j simulation AFR','fontsize',8) ; 231 231 end 232 232 text(-270,40,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 280 280 if e==4; 281 281 co=colorbar; 282 set (co,'ytick',[-40:20:40],'yticklabel',[-40:20:40],'fontname','Arial','fontsize',5);282 set (co,'ytick',[-40:20:40],'yticklabel',[-40:20:40],'fontname','Arial','fontsize',5); 283 283 end; 284 284 hold on; … … 288 288 axis([-180 180 20 80]); 289 289 if e==7; 290 290 title('HGT500 composite EOF2 1030j simulation AFR','fontsize',8) ; 291 291 end 292 292 text(-270,40,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_hgt500_AFR_eof2_1030_min.m
r17 r21 226 226 if e==4; 227 227 co=colorbar; 228 set (co,'ytick',[-50:5:50],'yticklabel',[-50:5:50],'fontname','Arial','fontsize',5);228 set (co,'ytick',[-50:5:50],'yticklabel',[-50:5:50],'fontname','Arial','fontsize',5); 229 229 end; 230 230 hold on; … … 234 234 axis([-180 180 20 80]); 235 235 if e==7; 236 title('HGT500 composite EOF2 1030j simulation AFR 236 title('HGT500 composite EOF2 1030j simulation AFR','fontsize',8) ; 237 237 end 238 238 text(-270,40,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 288 288 if e==4; 289 289 co=colorbar; 290 set (co,'ytick',[-40:20:40],'yticklabel',[-40:20:40],'fontname','Arial','fontsize',5);290 set (co,'ytick',[-40:20:40],'yticklabel',[-40:20:40],'fontname','Arial','fontsize',5); 291 291 end; 292 292 hold on; … … 296 296 axis([-180 180 20 80]); 297 297 if e==7; 298 298 title('HGT500 composite EOF2 1030j simulation AFR','fontsize',8) ; 299 299 end 300 300 text(-270,40,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_hgt500_AFR_eof2_30100.m
r17 r21 73 73 if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && hgt500_eof(b,a)>(hgt500bis)); 74 74 compteurmax=compteurmax+1; 75 hgt500_filtre_compomax(compteurmax,1,:,:)=hgt500(b-20,a,:,:);75 hgt500_filtre_compomax(compteurmax,1,:,:)=hgt500(b-20,a,:,:); 76 76 hgt500_filtre_compomax(compteurmax,2,:,:)=hgt500(b-19,a,:,:); 77 77 hgt500_filtre_compomax(compteurmax,3,:,:)=hgt500(b-18,a,:,:); … … 200 200 if e==4; 201 201 co=colorbar; 202 set (co,'ytick',[-80:40:80],'yticklabel',[-80:40:80],'fontname','Arial','fontsize',5);202 set (co,'ytick',[-80:40:80],'yticklabel',[-80:40:80],'fontname','Arial','fontsize',5); 203 203 end; 204 204 hold on; … … 208 208 axis([-180 180 20 80]); 209 209 if e==7; 210 title('HGT500 composite EOF2 30100j simulation AFR','fontsize',8) ;210 title('HGT500 composite EOF2 30100j simulation AFR','fontsize',8) ; 211 211 end 212 212 text(-270,40,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 255 255 if e==4; 256 256 co=colorbar; 257 set (co,'ytick',[-80:40:80],'yticklabel',[-80:40:80],'fontname','Arial','fontsize',5);257 set (co,'ytick',[-80:40:80],'yticklabel',[-80:40:80],'fontname','Arial','fontsize',5); 258 258 end; 259 259 hold on; … … 263 263 axis([-180 180 20 80]); 264 264 if e==7; 265 265 title('HGT500 composite EOF2 30100j simulation AFR','fontsize',8) ; 266 266 end 267 267 text(-270,40,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_hgt500_TROP_eof1_1030.m
r17 r21 73 73 if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && hgt500_eof(b,a)>(hgt500bis)); 74 74 compteurmax=compteurmax+1; 75 hgt500_filtre_compomax(compteurmax,1,:,:)=hgt500(b-20,a,:,:);75 hgt500_filtre_compomax(compteurmax,1,:,:)=hgt500(b-20,a,:,:); 76 76 hgt500_filtre_compomax(compteurmax,2,:,:)=hgt500(b-19,a,:,:); 77 77 hgt500_filtre_compomax(compteurmax,3,:,:)=hgt500(b-18,a,:,:); … … 200 200 if e==4; 201 201 co=colorbar; 202 set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6);202 set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6); 203 203 end; 204 204 hold on; … … 208 208 axis([-180 180 20 80]); 209 209 if e==7; 210 title('HGT500 composite EOF1 1030j simulation TROP','fontsize',8) ;210 title('HGT500 composite EOF1 1030j simulation TROP','fontsize',8) ; 211 211 end 212 212 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 255 255 if e==4; 256 256 co=colorbar; 257 set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6);257 set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6); 258 258 end; 259 259 hold on; … … 263 263 axis([-180 180 20 80]); 264 264 if e==7; 265 265 title('HGT500 composite EOF1 1030j simulation TROP','fontsize',8) ; 266 266 end 267 267 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_hgt500_TROP_eof1_30100.m
r17 r21 73 73 if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && hgt500_eof(b,a)>(hgt500bis)); 74 74 compteurmax=compteurmax+1; 75 hgt500_filtre_compomax(compteurmax,1,:,:)=hgt500(b-20,a,:,:);75 hgt500_filtre_compomax(compteurmax,1,:,:)=hgt500(b-20,a,:,:); 76 76 hgt500_filtre_compomax(compteurmax,2,:,:)=hgt500(b-19,a,:,:); 77 77 hgt500_filtre_compomax(compteurmax,3,:,:)=hgt500(b-18,a,:,:); … … 200 200 if e==4; 201 201 co=colorbar; 202 set (co,'ytick',[-80:40:80],'yticklabel',[-80:40:80],'fontname','Arial','fontsize',5);202 set (co,'ytick',[-80:40:80],'yticklabel',[-80:40:80],'fontname','Arial','fontsize',5); 203 203 end; 204 204 hold on; … … 208 208 axis([-180 180 20 80]); 209 209 if e==7; 210 title('HGT500 composite EOF1 30100j simulation TROP','fontsize',8) ;210 title('HGT500 composite EOF1 30100j simulation TROP','fontsize',8) ; 211 211 end 212 212 text(-270,40,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 255 255 if e==4; 256 256 co=colorbar; 257 set (co,'ytick',[-80:40:80],'yticklabel',[-80:40:80],'fontname','Arial','fontsize',5);257 set (co,'ytick',[-80:40:80],'yticklabel',[-80:40:80],'fontname','Arial','fontsize',5); 258 258 end; 259 259 hold on; … … 263 263 axis([-180 180 20 80]); 264 264 if e==7; 265 265 title('HGT500 composite EOF1 30100j simulation TROP','fontsize',8) ; 266 266 end 267 267 text(-270,40,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_hgt500_TROP_eof2_1030.m
r17 r21 73 73 if (delta(b+1)<0 && delta(b)>0 && hgt500_eof(b,a)>(hgt500bis)); 74 74 compteurmax=compteurmax+1; 75 hgt500_filtre_compomax(compteurmax,1,:,:)=hgt500(b-20,a,:,:);75 hgt500_filtre_compomax(compteurmax,1,:,:)=hgt500(b-20,a,:,:); 76 76 hgt500_filtre_compomax(compteurmax,2,:,:)=hgt500(b-19,a,:,:); 77 77 hgt500_filtre_compomax(compteurmax,3,:,:)=hgt500(b-18,a,:,:); … … 200 200 if e==4; 201 201 co=colorbar; 202 set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6);202 set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6); 203 203 end; 204 204 hold on; … … 208 208 axis([-180 180 20 80]); 209 209 if e==7; 210 title('HGT500 composite EOF2 1030j simulation TROP','fontsize',8) ;210 title('HGT500 composite EOF2 1030j simulation TROP','fontsize',8) ; 211 211 end 212 212 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 255 255 if e==4; 256 256 co=colorbar; 257 set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6);257 set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6); 258 258 end; 259 259 hold on; … … 263 263 axis([-180 180 20 80]); 264 264 if e==7; 265 265 title('HGT500 composite EOF2 1030j simulation TROP','fontsize',8) ; 266 266 end 267 267 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_hgt500_TROP_eof2_30100.m
r17 r21 73 73 if (delta(b+1)<0 && delta(b)>0 && hgt500_eof(b,a)>(hgt500bis)); 74 74 compteurmax=compteurmax+1; 75 hgt500_filtre_compomax(compteurmax,1,:,:)=hgt500(b-20,a,:,:);75 hgt500_filtre_compomax(compteurmax,1,:,:)=hgt500(b-20,a,:,:); 76 76 hgt500_filtre_compomax(compteurmax,2,:,:)=hgt500(b-19,a,:,:); 77 77 hgt500_filtre_compomax(compteurmax,3,:,:)=hgt500(b-18,a,:,:); … … 200 200 if e==4; 201 201 co=colorbar; 202 set (co,'ytick',[-80:40:80],'yticklabel',[-80:40:80],'fontname','Arial','fontsize',5);202 set (co,'ytick',[-80:40:80],'yticklabel',[-80:40:80],'fontname','Arial','fontsize',5); 203 203 end; 204 204 hold on; … … 208 208 axis([-180 180 20 80]); 209 209 if e==7; 210 title('HGT500 composite EOF2 30100j simulation TROP','fontsize',8) ;210 title('HGT500 composite EOF2 30100j simulation TROP','fontsize',8) ; 211 211 end 212 212 text(-270,40,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 255 255 if e==4; 256 256 co=colorbar; 257 set (co,'ytick',[-80:40:80],'yticklabel',[-80:40:80],'fontname','Arial','fontsize',5);257 set (co,'ytick',[-80:40:80],'yticklabel',[-80:40:80],'fontname','Arial','fontsize',5); 258 258 end; 259 259 hold on; … … 263 263 axis([-180 180 20 80]); 264 264 if e==7; 265 265 title('HGT500 composite EOF2 30100j simulation TROP','fontsize',8) ; 266 266 end 267 267 text(-270,40,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_olr_AFR_eof1_1030.m
r18 r21 91 91 92 92 ncload('rlut_d.AfNQIVIV.200006-09.nc','rlut','lat','lon'); 93 olr(:,29,:,:)=rlut; 94 clear rlut; 93 olr(:,29,:,:)=rlut; clear rlut; 95 94 96 95 % selction du domaine geographique … … 135 134 if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 136 135 compteurmax=compteurmax+1; 137 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:);136 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:); 138 137 olr_filtre_compomax(compteurmax,2,:,:)=olr1(b-19,a,:,:); 139 138 olr_filtre_compomax(compteurmax,3,:,:)=olr1(b-18,a,:,:); … … 277 276 if e==4; 278 277 co=colorbar; 279 set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6);278 set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6); 280 279 end; 281 280 hold on; … … 285 284 axis([-180 180 -20 40]); 286 285 if e==7; 287 title('OLR composite EOF1 1030j simulation AFR 286 title('OLR composite EOF1 1030j simulation AFR','fontsize',8) ; 288 287 end 289 288 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 339 338 if e==4; 340 339 co=colorbar; 341 set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6);340 set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6); 342 341 end; 343 342 hold on; … … 347 346 axis([-180 180 -20 40]); 348 347 if e==7; 349 348 title('OLR composite EOF1 1030j simulation AFR','fontsize',8) ; 350 349 end 351 350 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_olr_AFR_eof1_30100.m
r18 r21 134 134 if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 135 135 compteurmax=compteurmax+1; 136 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:);136 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:); 137 137 olr_filtre_compomax(compteurmax,2,:,:)=olr1(b-19,a,:,:); 138 138 olr_filtre_compomax(compteurmax,3,:,:)=olr1(b-18,a,:,:); … … 284 284 axis([-180 180 -20 40]); 285 285 if e==7; 286 title('OLR composite EOF1 30100j simulation AFR','fontsize',8) ;286 title('OLR composite EOF1 30100j simulation AFR','fontsize',8) ; 287 287 end 288 288 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 346 346 axis([-180 180 -20 40]); 347 347 if e==7; 348 348 title('OLR composite EOF1 30100j simulation AFR','fontsize',8) ; 349 349 end 350 350 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_olr_AFR_eof2_1030.m
r18 r21 134 134 if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 135 135 compteurmax=compteurmax+1; 136 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:);136 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:); 137 137 olr_filtre_compomax(compteurmax,2,:,:)=olr1(b-19,a,:,:); 138 138 olr_filtre_compomax(compteurmax,3,:,:)=olr1(b-18,a,:,:); … … 288 288 axis([-180 180 -20 40]); 289 289 if e==7; 290 title('OLR composite EOF2 1030j simulation AFR','fontsize',8) ;290 title('OLR composite EOF2 1030j simulation AFR','fontsize',8) ; 291 291 end 292 292 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 351 351 axis([-180 180 -20 40]); 352 352 if e==7; 353 353 title('OLR composite EOF2 1030j simulation AFR','fontsize',8) ; 354 354 end 355 355 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_olr_AFR_eof2_30100.m
r18 r21 134 134 if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 135 135 compteurmax=compteurmax+1; 136 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:);136 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:); 137 137 olr_filtre_compomax(compteurmax,2,:,:)=olr1(b-19,a,:,:); 138 138 olr_filtre_compomax(compteurmax,3,:,:)=olr1(b-18,a,:,:); … … 284 284 axis([-180 180 -20 40]); 285 285 if e==7; 286 title('OLR composite EOF2 30100j simulation AFR','fontsize',8) ;286 title('OLR composite EOF2 30100j simulation AFR','fontsize',8) ; 287 287 end 288 288 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 346 346 axis([-180 180 -20 40]); 347 347 if e==7; 348 348 title('OLR composite EOF2 30100j simulation AFR','fontsize',8) ; 349 349 end 350 350 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_olr_AS_eof1_1030.m
r18 r21 285 285 axis([-180 180 -20 40]); 286 286 if e==7; 287 title('OLR composite EOF1 1030j simulation AS','fontsize',8) ;287 title('OLR composite EOF1 1030j simulation AS','fontsize',8) ; 288 288 end 289 289 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 347 347 axis([-180 180 -20 40]); 348 348 if e==7; 349 349 title('OLR composite EOF1 1030j simulation AS','fontsize',8) ; 350 350 end 351 351 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_olr_AS_eof1_30100.m
r18 r21 135 135 if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 136 136 compteurmax=compteurmax+1; 137 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:);137 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:); 138 138 olr_filtre_compomax(compteurmax,2,:,:)=olr1(b-19,a,:,:); 139 139 olr_filtre_compomax(compteurmax,3,:,:)=olr1(b-18,a,:,:); … … 285 285 axis([-180 180 -20 40]); 286 286 if e==7; 287 title('OLR composite EOF1 30100j simulation AS','fontsize',8) ;287 title('OLR composite EOF1 30100j simulation AS','fontsize',8) ; 288 288 end 289 289 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 347 347 axis([-180 180 -20 40]); 348 348 if e==7; 349 349 title('OLR composite EOF1 30100j simulation AS','fontsize',8) ; 350 350 end 351 351 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_olr_AS_eof2_1030.m
r18 r21 135 135 if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 136 136 compteurmax=compteurmax+1; 137 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:);137 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:); 138 138 olr_filtre_compomax(compteurmax,2,:,:)=olr1(b-19,a,:,:); 139 139 olr_filtre_compomax(compteurmax,3,:,:)=olr1(b-18,a,:,:); … … 285 285 axis([-180 180 -20 40]); 286 286 if e==7; 287 title('OLR composite EOF2 1030j simulation AS','fontsize',8) ;287 title('OLR composite EOF2 1030j simulation AS','fontsize',8) ; 288 288 end 289 289 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 347 347 axis([-180 180 -20 40]); 348 348 if e==7; 349 349 title('OLR composite EOF2 1030j simulation AS','fontsize',8) ; 350 350 end 351 351 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_olr_AS_eof2_30100.m
r18 r21 135 135 if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 136 136 compteurmax=compteurmax+1; 137 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:);137 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:); 138 138 olr_filtre_compomax(compteurmax,2,:,:)=olr1(b-19,a,:,:); 139 139 olr_filtre_compomax(compteurmax,3,:,:)=olr1(b-18,a,:,:); … … 285 285 axis([-180 180 -20 40]); 286 286 if e==7; 287 title('OLR composite EOF2 30100j simulation AS 287 title('OLR composite EOF2 30100j simulation AS','fontsize',8) ; 288 288 end 289 289 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 347 347 axis([-180 180 -20 40]); 348 348 if e==7; 349 349 title('OLR composite EOF2 30100j simulation AS','fontsize',8) ; 350 350 end 351 351 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_olr_CTL_eof1_1030.m
r18 r21 137 137 if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 138 138 compteurmax=compteurmax+1; 139 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:);139 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:); 140 140 olr_filtre_compomax(compteurmax,2,:,:)=olr1(b-19,a,:,:); 141 141 olr_filtre_compomax(compteurmax,3,:,:)=olr1(b-18,a,:,:); … … 287 287 axis([-180 180 -20 40]); 288 288 if e==7; 289 title('OLR composite EOF1 1030j simulation CTL','fontsize',8) ;289 title('OLR composite EOF1 1030j simulation CTL','fontsize',8) ; 290 290 end 291 291 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 349 349 axis([-180 180 -20 40]); 350 350 if e==7; 351 351 title('OLR composite EOF1 1030j simulation CTL','fontsize',8) ; 352 352 end 353 353 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_olr_CTL_eof1_30100.m
r18 r21 136 136 if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 137 137 compteurmax=compteurmax+1; 138 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:);138 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:); 139 139 olr_filtre_compomax(compteurmax,2,:,:)=olr1(b-19,a,:,:); 140 140 olr_filtre_compomax(compteurmax,3,:,:)=olr1(b-18,a,:,:); … … 286 286 axis([-180 180 -20 40]); 287 287 if e==7; 288 title('OLR composite EOF1 30100j simulation CTL','fontsize',8) ;288 title('OLR composite EOF1 30100j simulation CTL','fontsize',8) ; 289 289 end 290 290 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 348 348 axis([-180 180 -20 40]); 349 349 if e==7; 350 350 title('OLR composite EOF1 30100j simulation CTL','fontsize',8) ; 351 351 end 352 352 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_olr_CTL_eof2_1030.m
r18 r21 136 136 if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 137 137 compteurmax=compteurmax+1; 138 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:);138 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:); 139 139 olr_filtre_compomax(compteurmax,2,:,:)=olr1(b-19,a,:,:); 140 140 olr_filtre_compomax(compteurmax,3,:,:)=olr1(b-18,a,:,:); … … 286 286 axis([-180 180 -20 40]); 287 287 if e==7; 288 title('OLR composite EOF2 1030j simulation CTL','fontsize',8) ;288 title('OLR composite EOF2 1030j simulation CTL','fontsize',8) ; 289 289 end 290 290 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 348 348 axis([-180 180 -20 40]); 349 349 if e==7; 350 350 title('OLR composite EOF2 1030j simulation CTL','fontsize',8) ; 351 351 end 352 352 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_olr_CTL_eof2_30100.m
r18 r21 136 136 if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 137 137 compteurmax=compteurmax+1; 138 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:);138 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:); 139 139 olr_filtre_compomax(compteurmax,2,:,:)=olr1(b-19,a,:,:); 140 140 olr_filtre_compomax(compteurmax,3,:,:)=olr1(b-18,a,:,:); … … 286 286 axis([-180 180 -20 40]); 287 287 if e==7; 288 title('OLR composite EOF2 30100j simulation CTL 288 title('OLR composite EOF2 30100j simulation CTL','fontsize',8) ; 289 289 end 290 290 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 348 348 axis([-180 180 -20 40]); 349 349 if e==7; 350 350 title('OLR composite EOF2 30100j simulation CTL','fontsize',8) ; 351 351 end 352 352 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_olr_TR_eof1_1030.m
r18 r21 135 135 if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 136 136 compteurmax=compteurmax+1; 137 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:);137 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:); 138 138 olr_filtre_compomax(compteurmax,2,:,:)=olr1(b-19,a,:,:); 139 139 olr_filtre_compomax(compteurmax,3,:,:)=olr1(b-18,a,:,:); … … 283 283 axis([-180 180 -20 40]); 284 284 if e==7; 285 title('OLR composite EOF1 1030j simulation TR','fontsize',8) ;285 title('OLR composite EOF1 1030j simulation TR','fontsize',8) ; 286 286 end 287 287 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 345 345 axis([-180 180 -20 40]); 346 346 if e==7; 347 347 title('OLR composite EOF1 1030j simulation TR','fontsize',8) ; 348 348 end 349 349 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_olr_TR_eof1_30100.m
r17 r21 143 143 if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>1); 144 144 compteurmax=compteurmax+1; 145 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:);145 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:); 146 146 olr_filtre_compomax(compteurmax,2,:,:)=olr1(b-19,a,:,:); 147 147 olr_filtre_compomax(compteurmax,3,:,:)=olr1(b-18,a,:,:); … … 291 291 axis([-180 180 -20 40]); 292 292 if e==7; 293 title('OLR composite EOF1 30100j simulation TR','fontsize',8) ;293 title('OLR composite EOF1 30100j simulation TR','fontsize',8) ; 294 294 end 295 295 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 353 353 axis([-180 180 -20 40]); 354 354 if e==7; 355 355 title('OLR composite EOF1 30100j simulation TR','fontsize',8) ; 356 356 end 357 357 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_olr_TR_eof2_1030.m
r18 r21 143 143 if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>1); 144 144 compteurmax=compteurmax+1; 145 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:);145 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:); 146 146 olr_filtre_compomax(compteurmax,2,:,:)=olr1(b-19,a,:,:); 147 147 olr_filtre_compomax(compteurmax,3,:,:)=olr1(b-18,a,:,:); … … 291 291 axis([-180 180 -20 40]); 292 292 if e==7; 293 title('OLR composite EOF2 30100j simulation TR','fontsize',8) ;293 title('OLR composite EOF2 30100j simulation TR','fontsize',8) ; 294 294 end 295 295 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 353 353 axis([-180 180 -20 40]); 354 354 if e==7; 355 355 title('OLR composite EOF2 30100j simulation TR','fontsize',8) ; 356 356 end 357 357 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/composite_olr_TR_eof2_30100.m
r18 r21 143 143 if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>1); 144 144 compteurmax=compteurmax+1; 145 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:);145 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:); 146 146 olr_filtre_compomax(compteurmax,2,:,:)=olr1(b-19,a,:,:); 147 147 olr_filtre_compomax(compteurmax,3,:,:)=olr1(b-18,a,:,:); … … 291 291 axis([-180 180 -20 40]); 292 292 if e==7; 293 title('OLR composite EOF2 30100j simulation TR','fontsize',8) ;293 title('OLR composite EOF2 30100j simulation TR','fontsize',8) ; 294 294 end 295 295 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); … … 353 353 axis([-180 180 -20 40]); 354 354 if e==7; 355 355 title('OLR composite EOF2 30100j simulation TR','fontsize',8) ; 356 356 end 357 357 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); -
trunk/src/SIMULS_IRCAAM/prepare_olr_filtre_AFR.m
r16 r21 110 110 olr1(:,:,:,1:64)=olr(:,:,:,65:128); 111 111 olr1(:,:,:,65:128)=olr(:,:,:,1:64); 112 clear olr 112 clear olr; 113 113 olr=olr1; 114 114 clear olr1; -
trunk/src/SIMULS_IRCAAM/prepare_olr_filtre_CTL.m
r16 r21 112 112 olr1(:,:,:,1:64)=olr(:,:,:,65:128); 113 113 olr1(:,:,:,65:128)=olr(:,:,:,1:64); 114 clear olr 114 clear olr; 115 115 olr=olr1; 116 116 clear olr1; -
trunk/src/SIMULS_IRCAAM/prepare_olr_filtre_CTL_clim.m
r16 r21 112 112 olr1(:,:,:,1:64)=olr(:,:,:,65:128); 113 113 olr1(:,:,:,65:128)=olr(:,:,:,1:64); 114 clear olr 114 clear olr; 115 115 olr=olr1; 116 116 clear olr1; -
trunk/src/SIMULS_IRCAAM/time_serie_afr_1030.m
r18 r21 35 35 subplot(4,3,a-22); 36 36 end 37 clear koro1 37 clear koro1; 38 38 debut=(a-1)*122+1; 39 39 fin=a*122; -
trunk/src/SIMULS_IRCAAM/time_serie_afr_30100.m
r18 r21 34 34 subplot(4,3,a-22); 35 35 end 36 clear koro1 36 clear koro1; 37 37 debut=(a-1)*122+1; 38 38 fin=a*122; -
trunk/src/SIMULS_IRCAAM/time_serie_as_30100.m
r18 r21 36 36 subplot(4,3,a-22); 37 37 end 38 clear koro1 38 clear koro1; 39 39 debut=(a-1)*122+1; 40 40 fin=a*122; -
trunk/src/SIMULS_IRCAAM/time_serie_ctl_1030.m
r18 r21 35 35 subplot(4,3,a-22); 36 36 end 37 clear koro1 37 clear koro1; 38 38 debut=(a-1)*122+1; 39 39 fin=a*122; -
trunk/src/SIMULS_IRCAAM/time_serie_ctl_30100.m
r18 r21 35 35 subplot(4,3,a-22); 36 36 end 37 clear koro1 37 clear koro1; 38 38 debut=(a-1)*122+1; 39 39 fin=a*122; -
trunk/src/SIMULS_IRCAAM/time_serie_noaa_annuel_20.m
r18 r21 1 1 % 2 2 3 clear; 3 4 close all; -
trunk/src/SIMULS_IRCAAM/time_serie_tr_1030.m
r18 r21 34 34 subplot(4,3,a-22); 35 35 end 36 clear koro1 36 clear koro1; 37 37 debut=(a-1)*122+1; 38 38 fin=a*122; -
trunk/src/SIMULS_IRCAAM/time_serie_tr_30100.m
r18 r21 35 35 subplot(4,3,a-22); 36 36 end 37 clear koro1 37 clear koro1; 38 38 debut=(a-1)*122+1; 39 39 fin=a*122; -
trunk/src/mode_sahelien/composite_olr_phase2.m
r18 r21 149 149 clear olr; 150 150 olr2=reshape(olr1,3416,21,144); 151 olr3(:,:)=mean(olr2); clear olr2; 151 olr3(:,:)=mean(olr2); 152 clear olr2; 152 153 153 154 load phase2.txt; 154 155 phase2=reshape(phase2,122,28,3); 155 156 olr_filtre_compomax=zeros(90,21,144); 156 157 compteurmax=0; 157 158 -
trunk/src/mode_sahelien/composite_olr_phase3.m
r18 r21 154 154 load phase3.txt; 155 155 phase3=reshape(phase3,122,28,3); 156 156 olr_filtre_compomax=zeros(90,21,144); 157 157 compteurmax=0; 158 158 -
trunk/src/mode_sahelien/composite_olr_phase4.m
r18 r21 149 149 clear olr; 150 150 olr2=reshape(olr1,3416,21,144); 151 olr3(:,:)=mean(olr2); clear olr2; 151 olr3(:,:)=mean(olr2); 152 clear olr2; 152 153 153 154 load phase4.txt; 154 155 phase4=reshape(phase4,122,28,3); 155 156 olr_filtre_compomax=zeros(90,21,144); 156 157 compteurmax=0; 157 158 -
trunk/src/mode_sahelien/composite_olr_phase5.m
r18 r21 154 154 load phase5.txt; 155 155 phase5=reshape(phase5,122,28,3); 156 156 olr_filtre_compomax=zeros(90,21,144); 157 157 compteurmax=0; 158 158 -
trunk/src/mode_sahelien/composite_olr_phase6.m
r18 r21 154 154 load phase6.txt; 155 155 phase6=reshape(phase6,122,28,3); 156 156 olr_filtre_compomax=zeros(90,21,144); 157 157 compteurmax=0; 158 158 -
trunk/src/mode_sahelien/composite_olr_phase7.m
r18 r21 149 149 clear olr; 150 150 olr2=reshape(olr1,3416,21,144); 151 olr3(:,:)=mean(olr2); clear olr2; 151 olr3(:,:)=mean(olr2); 152 clear olr2; 152 153 153 154 load phase7.txt; 154 155 phase7=reshape(phase7,122,28,3); 155 156 olr_filtre_compomax=zeros(90,21,144); 156 157 compteurmax=0; 157 158 -
trunk/src/mode_sahelien/composite_olr_phase8.m
r18 r21 154 154 load phase8.txt; 155 155 phase8=reshape(phase8,122,28,3); 156 156 olr_filtre_compomax=zeros(90,21,144); 157 157 compteurmax=0; 158 158 -
trunk/src/mode_sahelien/time_serie_ER_rossby_1030.m
r18 r21 13 13 for b=1:28; 14 14 year(a,b)=1979+(b-1); 15 end 16 end 15 end; 16 end; 17 17 year=reshape(year,3416,1); 18 18 … … 21 21 for b=1:28; 22 22 day(a,b)=31+28+31+30+31+a; 23 end 24 end 23 end; 24 end; 25 25 day=reshape(day,3416,1); 26 26 … … 58 58 set (gca,'YTick',[-100:20:100],'YTickLabel',[-100:20:100],'fontname','Arial','fontsize',6); 59 59 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); 60 -
trunk/src/mode_sahelien/time_serie_ER_rossby_1030_complet.m
r18 r21 28 28 figure(3); 29 29 subplot(4,3,a-22); 30 end 30 end; 31 31 clear koro1 32 32 debut=(a-1)*122+1; … … 49 49 elseif a==24; 50 50 text(0,183,'time series EOF1 and EOF2 OLR ER 10-30d'); 51 end 51 end; 52 52 title(a+1978,'FontSize',8); 53 53 set (gca,'tickDir','out'); … … 56 56 set (gca,'YTick',[-100:20:100],'YTickLabel',[-100:20:100],'fontname','Arial','fontsize',6); 57 57 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); 58 end 58 end; 59 59 60 60 figure(1); -
trunk/src/mode_sahelien/time_serie_ER_rossby_30100.m
r18 r21 13 13 for b=1:28; 14 14 year(a,b)=1979+(b-1); 15 end 16 end 15 end; 16 end; 17 17 year=reshape(year,3416,1); 18 18 … … 21 21 for b=1:28; 22 22 day(a,b)=31+28+31+30+31+a; 23 end 24 end 23 end; 24 end; 25 25 day=reshape(day,3416,1); 26 26 -
trunk/src/mode_sahelien/time_serie_ER_rossby_30100_complet.m
r18 r21 28 28 figure(3); 29 29 subplot(4,3,a-22); 30 end 30 end; 31 31 clear koro1 32 32 debut=(a-1)*122+1; … … 49 49 elseif a==24; 50 50 text(0,183,'time series EOF1 and EOF2 OLR ER 30-100d'); 51 end 51 end; 52 52 title(a+1978,'FontSize',8); 53 53 set (gca,'tickDir','out'); … … 56 56 set (gca,'YTick',[-100:20:100],'YTickLabel',[-100:20:100],'fontname','Arial','fontsize',6); 57 57 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); 58 end 58 end; 59 59 60 60 figure(1); -
trunk/src/mode_sahelien/time_serie_histog_phase1.m
r18 r21 88 88 for b=1:4 89 89 kara_diff(a,b)=kara_max(a,b)-kara_moy(b); 90 end 91 end 90 end; 91 end; 92 92 93 93 for a=1:4; … … 105 105 else 106 106 bar(jour,diff,'r'); 107 end 107 end; 108 108 set (gca,'XTick',[-10:2:10],'XTickLabel',[-10:2:10],'fontname','Arial','fontsize',6); 109 109 set (gca,'YTick',[-20:10:20],'YTickLabel',[-20:5:20],'fontname','Arial','fontsize',6); … … 120 120 title('Blocking','Fontsize',8,'Fontname','Arial'); 121 121 end; 122 end 122 end; -
trunk/src/mode_sahelien/time_serie_histog_phase2.m
r18 r21 88 88 for b=1:4 89 89 kara_diff(a,b)=kara_max(a,b)-kara_moy(b); 90 end 91 end 90 end; 91 end; 92 92 93 93 for a=1:4; … … 105 105 else 106 106 bar(jour,diff,'r'); 107 end 107 end; 108 108 set (gca,'XTick',[-10:2:10],'XTickLabel',[-10:2:10],'fontname','Arial','fontsize',6); 109 109 set (gca,'YTick',[-20:10:20],'YTickLabel',[-20:5:20],'fontname','Arial','fontsize',6); … … 120 120 title('Blocking','Fontsize',8,'Fontname','Arial'); 121 121 end; 122 end 122 end; -
trunk/src/mode_sahelien/time_serie_histog_phase3.m
r17 r21 23 23 % boucle sur les annees 24 24 for a=1:28; 25 % 1ere condition (laisser une marge de 10)26 25 for b=1:122; 26 % 1ere condition (laisser une marge de 10) 27 27 if (b>11 && b<112); 28 28 % condition 1 on est a un maximum deltas apres sont negetifs et … … 88 88 for b=1:4 89 89 kara_diff(a,b)=kara_max(a,b)-kara_moy(b); 90 end 91 end 90 end; 91 end; 92 92 93 93 for a=1:4; … … 105 105 else 106 106 bar(jour,diff,'r'); 107 end 107 end; 108 108 set (gca,'XTick',[-10:2:10],'XTickLabel',[-10:2:10],'fontname','Arial','fontsize',6); 109 109 set (gca,'YTick',[-20:10:20],'YTickLabel',[-20:5:20],'fontname','Arial','fontsize',6); -
trunk/src/mode_sahelien/time_serie_histog_phase4.m
r18 r21 88 88 for b=1:4 89 89 kara_diff(a,b)=kara_max(a,b)-kara_moy(b); 90 end 91 end 90 end; 91 end; 92 92 93 93 for a=1:4; … … 105 105 else 106 106 bar(jour,diff,'r'); 107 end 107 end; 108 108 set (gca,'XTick',[-10:2:10],'XTickLabel',[-10:2:10],'fontname','Arial','fontsize',6); 109 109 set (gca,'YTick',[-20:10:20],'YTickLabel',[-20:5:20],'fontname','Arial','fontsize',6); … … 120 120 title('Blocking','Fontsize',8,'Fontname','Arial'); 121 121 end; 122 end 122 end; -
trunk/src/mode_sahelien/time_serie_histog_phase5.m
r18 r21 88 88 for b=1:4 89 89 kara_diff(a,b)=kara_max(a,b)-kara_moy(b); 90 end 91 end 90 end; 91 end; 92 92 93 93 for a=1:4; … … 105 105 else 106 106 bar(jour,diff,'r'); 107 end 107 end; 108 108 set (gca,'XTick',[-10:2:10],'XTickLabel',[-10:2:10],'fontname','Arial','fontsize',6); 109 109 set (gca,'YTick',[-20:10:20],'YTickLabel',[-20:5:20],'fontname','Arial','fontsize',6); … … 120 120 title('Blocking','Fontsize',8,'Fontname','Arial'); 121 121 end; 122 end 122 end; -
trunk/src/mode_sahelien/time_serie_histog_phase6.m
r18 r21 88 88 for b=1:4 89 89 kara_diff(a,b)=kara_max(a,b)-kara_moy(b); 90 end 91 end 90 end; 91 end; 92 92 93 93 for a=1:4; … … 105 105 else 106 106 bar(jour,diff,'r'); 107 end 107 end; 108 108 set (gca,'XTick',[-10:2:10],'XTickLabel',[-10:2:10],'fontname','Arial','fontsize',6); 109 109 set (gca,'YTick',[-20:10:20],'YTickLabel',[-20:5:20],'fontname','Arial','fontsize',6); … … 120 120 title('Blocking','Fontsize',8,'Fontname','Arial'); 121 121 end; 122 end 122 end; -
trunk/src/mode_sahelien/time_serie_histog_phase7.m
r17 r21 88 88 for b=1:4 89 89 kara_diff(a,b)=kara_max(a,b)-kara_moy(b); 90 end 91 end 90 end; 91 end; 92 92 93 93 for a=1:4; … … 105 105 else 106 106 bar(jour,diff,'r'); 107 end 107 end; 108 108 set (gca,'XTick',[-10:2:10],'XTickLabel',[-10:2:10],'fontname','Arial','fontsize',6); 109 109 set (gca,'YTick',[-20:10:20],'YTickLabel',[-20:5:20],'fontname','Arial','fontsize',6); … … 120 120 title('Blocking','Fontsize',8,'Fontname','Arial'); 121 121 end; 122 end 122 end; -
trunk/src/mode_sahelien/time_serie_histog_phase8.m
r18 r21 88 88 for b=1:4 89 89 kara_diff(a,b)=kara_max(a,b)-kara_moy(b); 90 end 91 end 90 end; 91 end; 92 92 93 93 for a=1:4; … … 105 105 else 106 106 bar(jour,diff,'r'); 107 end 107 end; 108 108 set (gca,'XTick',[-10:2:10],'XTickLabel',[-10:2:10],'fontname','Arial','fontsize',6); 109 109 set (gca,'YTick',[-20:10:20],'YTickLabel',[-20:5:20],'fontname','Arial','fontsize',6); … … 117 117 elseif a==3; 118 118 title('NAO-','Fontsize',8,'Fontname','Arial'); 119 119 else 120 120 title('Blocking','Fontsize',8,'Fontname','Arial'); 121 121 end; 122 end 122 end; -
trunk/src/mode_sahelien/time_serie_jjas.m
r18 r21 20 20 time_serie_olr2(a)=time_serie_olr2(a)/dt_olr2; 21 21 time_serie_olr3(a)=time_serie_olr3(a)/dt_olr3; 22 end 22 end; 23 23 24 24 save eof1_1030.txt time_serie_olr1 -ascii; … … 37 37 figure(3); 38 38 subplot(4,3,a-22); 39 end 40 clear koro1 39 end; 40 clear koro1; 41 41 debut=(a-1)*122+1; 42 42 fin=a*122; … … 49 49 y=zeros(122,1); 50 50 plot(x,y,'k:'); 51 title(a+1978,'FontSize',8) 51 title(a+1978,'FontSize',8); 52 52 if a==2; 53 53 text(0,5.5,'serie temporelle VP1'); 54 54 elseif a==13; 55 55 text(0,5.5,'serie temporelle VP1'); 56 56 elseif a==24; 57 57 text(0,5.5,'serie temporelle VP1'); 58 end 58 end; 59 59 set (gca,'tickDir','out'); 60 60 axis([1 122 -3 3]); 61 61 62 set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6);63 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6);62 set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); 63 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); 64 64 end; -
trunk/src/mode_sahelien/time_serie_jjas_10_25.m
r18 r21 20 20 time_serie_olr2(a)=time_serie_olr2(a)/dt_olr2; 21 21 time_serie_olr3(a)=time_serie_olr3(a)/dt_olr3; 22 end 22 end; 23 23 24 24 save eof1_olr_rossby_jjas_1025.txt time_serie_olr1 -ascii; … … 37 37 figure(3); 38 38 subplot(4,3,a-22); 39 end 40 clear koro1 39 end; 40 clear koro1; 41 41 debut=(a-1)*122+1; 42 42 fin=a*122; … … 51 51 if a==2; 52 52 text(0,5.5,'serie temporelle VP1 filtrage rossby 10-25j'); 53 53 elseif a==13; 54 54 text(0,5.5,'serie temporelle VP1 filtrage rossby 10-25j'); 55 55 elseif a==24; … … 60 60 61 61 set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); 62 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6);62 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); 63 63 end; 64 64 -
trunk/src/mode_sahelien/time_serie_jjas_25_90.m
r18 r21 20 20 time_serie_olr2(a)=time_serie_olr2(a)/dt_olr2; 21 21 time_serie_olr3(a)=time_serie_olr3(a)/dt_olr3; 22 end 22 end; 23 23 24 24 save eof1_olr_rossby_jjas_2590.txt time_serie_olr1 -ascii; … … 37 37 figure(3); 38 38 subplot(4,3,a-22); 39 end 40 clear koro1 39 end; 40 clear koro1; 41 41 debut=(a-1)*122+1; 42 42 fin=a*122; … … 52 52 if a==2; 53 53 text(0,5.5,'serie temporelle VP1 filtrage rossby 25-90j'); 54 54 elseif a==13; 55 55 text(0,5.5,'serie temporelle VP1 filtrage rossby 25-90j'); 56 56 elseif a==24; … … 61 61 62 62 set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); 63 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6);63 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); 64 64 end; 65 65 -
trunk/src/mode_sahelien/time_serie_jjas_30100.m
r18 r21 20 20 time_serie_olr2(a)=time_serie_olr2(a)/dt_olr2; 21 21 time_serie_olr3(a)=time_serie_olr3(a)/dt_olr3; 22 end 22 end; 23 23 24 24 save eof1_30100.txt time_serie_olr1 -ascii; … … 37 37 figure(3); 38 38 subplot(4,3,a-22); 39 end 40 clear koro1 39 end; 40 clear koro1; 41 41 debut=(a-1)*122+1; 42 42 fin=a*122; … … 49 49 y=zeros(122,1); 50 50 plot(x,y,'k:'); 51 title(a+1978,'FontSize',8) 51 title(a+1978,'FontSize',8); 52 52 if a==2; 53 53 text(0,5.5,'serie temporelle VP1 filtrage 30100'); 54 54 elseif a==13; 55 55 text(0,5.5,'serie temporelle VP1 filtrage 30100'); 56 56 elseif a==24; 57 57 text(0,5.5,'serie temporelle VP1 filtrage 30100'); 58 end 58 end; 59 59 set (gca,'tickDir','out'); 60 60 axis([1 122 -3 3]); 61 61 62 set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6);63 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6);62 set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); 63 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); 64 64 end; 65 65
Note: See TracChangeset
for help on using the changeset viewer.