Changeset 16 for trunk/src/SIMULS_IRCAAM/correlations_eofs_olr_jas.m
- Timestamp:
- 01/06/09 10:53:18 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/SIMULS_IRCAAM/correlations_eofs_olr_jas.m
r10 r16 1 %%%%%%%%%%%%%%% FICHIER OLR 2 % ouverture / lecture du fichier 3 clear 4 close all 5 6 % initialisation 7 clear; close all; 1 % 2 3 clear; 4 close all; 8 5 9 6 cd('..'); … … 30 27 olr1=zeros(siz,37,144); 31 28 olr1(:,:,1:72)=olr(:,:,73:144); 32 olr1(:,:,73:144)=olr(:,:,1:72); clear olr; 29 olr1(:,:,73:144)=olr(:,:,1:72); 30 clear olr; 33 31 34 32 % on permture le temps et l'espce … … 42 40 vecteur=zeros(siz,1); 43 41 for c=1:siz; 44 vecteur(c)=olr1(a,b,c); 42 vecteur(c)=olr1(a,b,c); 45 43 end 46 44 std_c=nanmean(vecteur); 47 45 for c=1:siz; 48 olr1(a,b,c)=olr1(a,b,c)-std_c; 49 end 46 olr1(a,b,c)=olr1(a,b,c)-std_c; 47 end 50 48 end 51 end; 52 53 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 49 end; 50 54 51 % on met sous la forme longitude latitude jour 55 56 52 57 53 % on selectionne les mois d ete (attention aux annees bixestiles!) … … 153 149 %olr=zeros(37,144,92,28); 154 150 155 156 olr1=permute(olr,[3 1 2 4]); clear olr; 157 olr(:,:,:)=nanmean(olr1); 158 151 olr1=permute(olr,[3 1 2 4]); 152 clear olr; 153 olr(:,:,:)=nanmean(olr1); 159 154 160 155 load eof3_noaa_annuel.txt … … 167 162 correlations_eof3=zeros(37,144); 168 163 169 lat=lat(5:41)'; 164 lat=lat(5:41)'; 170 165 lon=lon-180; 171 166 … … 177 172 y(c)=olr(a,b,c); 178 173 end 179 clear coco; coco=corrcoef(y,eof1_noaa_annuel); 174 clear coco; 175 coco=corrcoef(y,eof1_noaa_annuel); 180 176 correlations_eof1(a,b)=coco(2,1); 181 clear coco; coco=corrcoef(y,eof2_noaa_annuel); 177 clear coco; 178 coco=corrcoef(y,eof2_noaa_annuel); 182 179 correlations_eof2(a,b)=coco(2,1); 183 clear coco; coco=corrcoef(y,eof3_noaa_annuel); 180 clear coco; 181 coco=corrcoef(y,eof3_noaa_annuel); 184 182 correlations_eof3(a,b)=coco(2,1); 185 183 186 end 184 end 187 185 end 188 186 189 187 subplot(3,1,1); 190 correlations_eof1(1,1)=1.0001; correlations_eof1(1,2)=-1.0001; 188 correlations_eof1(1,1)=1.0001; 189 correlations_eof1(1,2)=-1.0001; 191 190 192 191 contourf(lon,lat,correlations_eof1,[-1:0.1:1],'LineStyle','none'); … … 206 205 subplot(3,1,2); 207 206 palette 208 correlations_eof2(1,1)=1.0001; correlations_eof2(1,2)=-1.0001; 207 correlations_eof2(1,1)=1.0001; 208 correlations_eof2(1,2)=-1.0001; 209 209 contourf(lon,lat,correlations_eof2,[-1:0.1:1],'LineStyle','none'); 210 210 colorbar; … … 224 224 subplot(3,1,3); 225 225 palette; 226 correlations_eof3(1,1)=1.0001; correlations_eof3(1,2)=-1.0001; 226 correlations_eof3(1,1)=1.0001; 227 correlations_eof3(1,2)=-1.0001; 227 228 contourf(lon,lat,correlations_eof3,[-1:0.1:1],'LineStyle','none'); 228 229 set (gca,'tickDir','out');
Note: See TracChangeset
for help on using the changeset viewer.