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