Changeset 16 for trunk/src/mode_sahelien/olr_eof_1025.m
- Timestamp:
- 01/06/09 10:53:18 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/mode_sahelien/olr_eof_1025.m
r10 r16 1 % 2 1 3 %%%%%%%%%%%%%%% FICHIER OLR 2 4 % ouverture / lecture du fichier … … 8 10 v=fread(fid,siz,'float'); 9 11 fclose(fid); 10 olr1=reshape(v,25,17,10500); clear v;11 12 olr1=reshape(v,25,17,10500); 13 clear v; 12 14 13 15 % on selectionne les mois d ete (attention aux annees bixestiles!) … … 101 103 clear siz28 siz28bis; 102 104 103 104 105 105 clear olr1; 106 106 olr=reshape(olr,425,3416); 107 107 108 108 % calcul de la matrice de covariance 109 109 vaudou=zeros(425,425); … … 111 111 for b=1:425; 112 112 clear e d covaria; 113 e=olr(b,:); 114 d=olr(a,:); 113 e=olr(b,:); 114 d=olr(a,:); 115 115 covaria=cov(e,d); 116 116 vaudou(a,b)=covaria(2,1); … … 121 121 [C,L]=eig(vaudou); 122 122 123 % spatialisation du premier mode124 125 126 % serie temporelle de la premiere vP127 128 123 delete('eof_1025_jjas.mat'); 129 124 save('eof_1025_jjas.mat','C','L','olr')
Note: See TracChangeset
for help on using the changeset viewer.