Changeset 16 for trunk/src/mode_sahelien/reconstitution.m
- Timestamp:
- 01/06/09 10:53:18 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/mode_sahelien/reconstitution.m
r15 r16 1 clear; close all; 1 % 2 3 clear; 4 close all; 2 5 3 6 load eof_rossby_jjas.mat … … 6 9 7 10 cosa=zeros(3,2684,425); 8 for a=1:3; % boucle sur les VP 11 % boucle sur les VP 12 for a=1:3; 9 13 for b=1:2684 % boucle sur le temps 10 14 for c=1:425 % boucle sur l espace … … 14 18 end; 15 19 16 cosi3(:,:)=sum(cosa); cosi3=reshape(cosi3,2684,17,25); 20 cosi3(:,:)=sum(cosa); 21 cosi3=reshape(cosi3,2684,17,25); 17 22 18 cosabis=cosa(1:2,:,:); 19 cosi2(:,:)=sum(cosabis); clear cosabis; 20 cosi2=reshape(cosi2,2684,17,25); 23 cosabis=cosa(1:2,:,:); 24 cosi2(:,:)=sum(cosabis); 25 clear cosabis; 26 cosi2=reshape(cosi2,2684,17,25); 21 27 22 cosa1(:,:)=cosa(1,:,:); clear cosa; 23 cosi1=reshape(cosa1,2684,17,25); clear cosa1; 28 cosa1(:,:)=cosa(1,:,:); 29 clear cosa; 30 cosi1=reshape(cosa1,2684,17,25); 31 clear cosa1; 24 32 25 26 figure(1); orient('landscape') 27 figure(2); orient('landscape') 28 figure(3); orient('landscape') 33 figure(1); 34 orient('landscape') 35 figure(2); 36 orient('landscape') 37 figure(3); 38 orient('landscape') 29 39 30 40 % On selectionne la demaine geographique 31 41 32 42 taille=(12-7+1)*(17-9+1); 33 cosa1=cosi1(:,7:12,9:17); clear cosi1; 34 cosa1=reshape(cosa1,2684,taille); 43 cosa1=cosi1(:,7:12,9:17); 44 clear cosi1; 45 cosa1=reshape(cosa1,2684,taille); 35 46 cosa1=cosa1'; 36 cosa1=mean(cosa1); 47 cosa1=mean(cosa1); 37 48 cosa1=cosa1'; 38 49 cosa1=reshape(cosa1,122,22); 39 50 40 cosa2=cosi2(:,7:12,9:17); clear cosi2; 41 cosa2=reshape(cosa2,2684,taille); 51 cosa2=cosi2(:,7:12,9:17); 52 clear cosi2; 53 cosa2=reshape(cosa2,2684,taille); 42 54 cosa2=cosa2'; 43 cosa2=mean(cosa2); 55 cosa2=mean(cosa2); 44 56 cosa2=cosa2'; 45 57 save mode_rossby_jjas.txt cosa2 -ascii; … … 47 59 cosa2=reshape(cosa2,122,22); 48 60 49 cosa3=cosi3(:,7:12,9:17); clear cosi3; 50 cosa3=reshape(cosa3,2684,taille); 61 cosa3=cosi3(:,7:12,9:17); 62 clear cosi3; 63 cosa3=reshape(cosa3,2684,taille); 51 64 cosa3=cosa3'; 52 cosa3=mean(cosa3); 65 cosa3=mean(cosa3); 53 66 cosa3=cosa3'; 54 67 cosa3=reshape(cosa3,122,22); … … 57 70 58 71 % On selectionne la demaine geographique 59 olr_centre=reshape(olr_centre,2684,17,25); 60 cosa4=olr_centre(:,7:12,9:17); clear olr_centre; 61 cosa4=reshape(cosa4,2684,taille); clear taille; 72 olr_centre=reshape(olr_centre,2684,17,25); 73 cosa4=olr_centre(:,7:12,9:17); 74 clear olr_centre; 75 cosa4=reshape(cosa4,2684,taille); 76 clear taille; 62 77 cosa4=cosa4'; 63 78 cosa4=mean(cosa4); … … 67 82 for a=1:22; 68 83 if a<12; 69 figure(1); subplot(4,3,a); 84 figure(1); 85 subplot(4,3,a); 70 86 else 71 figure(2); subplot(4,3,a-11); 87 figure(2); 88 subplot(4,3,a-11); 72 89 end 73 90 y1=cosa1(:,a); … … 75 92 y3=cosa3(:,a); 76 93 y4=cosa4(:,a); 77 x=[1:122]'; 94 x=[1:122]'; 78 95 plot(x,y1,'r'); 79 96 hold on; … … 83 100 y5=zeros(122,1); 84 101 plot(x,y5,'k:'); 85 102 86 103 title(a+1978,'Fontsize',8); 87 104 if a==2; … … 110 127 set (gca,'XTick',[10:10:90],'XTickLabel',[10:10:90],'fontname','Arial','fontsize',6); 111 128 112 end; 129 end; 113 130 114 131
Note: See TracChangeset
for help on using the changeset viewer.