Changeset 16 for trunk/src/mode_sahelien/reconstitution_ER_rossby_1030.m
- Timestamp:
- 01/06/09 10:53:18 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/mode_sahelien/reconstitution_ER_rossby_1030.m
r15 r16 1 % 1 2 2 3 clear; … … 10 11 11 12 cosa=zeros(3,3416,425); 12 for a=1:3; % boucle sur les VP 13 % boucle sur les VP 14 for a=1:3; 13 15 for b=1:3416 % boucle sur le temps 14 16 for c=1:425 % boucle sur l espace … … 18 20 end; 19 21 20 cosa1(:,:)=cosa(1,:,:); 21 cosa1=reshape(cosa1,3416,25,17); 22 cosa1(:,:)=cosa(1,:,:); 23 cosa1=reshape(cosa1,3416,25,17); 22 24 23 cosa2(:,:)=cosa(2,:,:); 24 cosa2=reshape(cosa2,3416,25,17); 25 cosa2(:,:)=cosa(2,:,:); 26 cosa2=reshape(cosa2,3416,25,17); 25 27 26 cosabis=cosa(1:2,:,:); 27 cosa12(:,:)=sum(cosabis); clear cosabis; 28 cosa12=reshape(cosa12,3416,25,17); 28 cosabis=cosa(1:2,:,:); 29 cosa12(:,:)=sum(cosabis); 30 clear cosabis; 31 cosa12=reshape(cosa12,3416,25,17); 29 32 30 31 32 figure(1); orient('landscape') 33 figure(2); orient('landscape') 34 figure(3); orient('landscape') 33 figure(1); 34 orient('landscape') 35 figure(2); 36 orient('landscape') 37 figure(3); 38 orient('landscape') 35 39 36 40 % On selectionne la demaine geographique 37 41 38 42 % la latitude : 7.5N (8) 12.5N (10) 39 % la longitude : 10W (9) 10E(17) 40 lat1=8; lat2=10; 41 lon1=9; lon2=17; 43 % la longitude : 10W (9) 10E(17) 44 lat1=8; 45 lat2=10; 46 lon1=9; 47 lon2=17; 42 48 taille=(lat2-lat1+1)*(lon2-lon1+1); 43 49 44 45 c osa12=cosa12(:,lon1:lon2,lat1:lat2); clear cosi2;46 cosa12=reshape(cosa12,3416,taille); 50 cosa12=cosa12(:,lon1:lon2,lat1:lat2); 51 clear cosi2; 52 cosa12=reshape(cosa12,3416,taille); 47 53 cosa12=cosa12'; 48 cosa12=mean(cosa12); 54 cosa12=mean(cosa12); 49 55 cosa12=cosa12'; 50 56 save eof12_rossby_1030.txt cosa12 -ascii; … … 52 58 53 59 % On selectionne la demaine geographique 54 olr=reshape(olr,3416,25,17); 55 cosa4=olr(:,lon1:lon2,lat1:lat2); clear olr; 56 cosa4=reshape(cosa4,3416,taille); clear taille; 60 olr=reshape(olr,3416,25,17); 61 cosa4=olr(:,lon1:lon2,lat1:lat2); 62 clear olr; 63 cosa4=reshape(cosa4,3416,taille); 64 clear taille; 57 65 cosa4=cosa4'; 58 66 cosa4=mean(cosa4); … … 62 70 for a=1:28; 63 71 if a<12; 64 figure(1); subplot(4,3,a); 72 figure(1); 73 subplot(4,3,a); 65 74 elseif a<24 66 figure(2); subplot(4,3,a-11); 75 figure(2); 76 subplot(4,3,a-11); 67 77 else 68 figure(3); subplot(4,3,a-23); 78 figure(3); 79 subplot(4,3,a-23); 69 80 end 70 81 y2=cosa12(:,a); 71 82 y4=cosa4(:,a); 72 x=[1:122]'; 73 plot(x,y2); hold on; 83 x=[1:122]'; 84 plot(x,y2); 85 hold on; 74 86 plot(x,y4,'k'); 75 87 y5=zeros(122,1); 76 88 plot(x,y5,'k:'); 77 89 78 90 title(a+1978,'Fontsize',8); 79 91 if a==2; … … 97 109 set (gca,'XTick',[10:10:122],'XTickLabel',[10:10:122],'fontname','Arial','fontsize',6); 98 110 99 end; 111 end; 100 112 101 figure(1); print -depsc2 reconstitution1_eof12_filtreER_1030.eps; 102 figure(2); print -depsc2 reconstitution2_eof12_filtreER_1030.eps; 103 figure(3); print -depsc2 reconstitution3_eof12_filtreER_1030.eps; 104 113 figure(1); 114 print -depsc2 reconstitution1_eof12_filtreER_1030.eps; 115 figure(2); 116 print -depsc2 reconstitution2_eof12_filtreER_1030.eps; 117 figure(3); 118 print -depsc2 reconstitution3_eof12_filtreER_1030.eps;
Note: See TracChangeset
for help on using the changeset viewer.