Ignore:
Timestamp:
01/06/09 10:53:18 (16 years ago)
Author:
pinsard
Message:

remove trailing blanks, split lines with multiple statements

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 
     3clear; 
     4close all; 
    85 
    96cd('..'); 
     
    3027olr1=zeros(siz,37,144); 
    3128olr1(:,:,1:72)=olr(:,:,73:144); 
    32 olr1(:,:,73:144)=olr(:,:,1:72); clear olr; 
     29olr1(:,:,73:144)=olr(:,:,1:72); 
     30clear olr; 
    3331 
    3432% on permture le temps et l'espce 
     
    4240      vecteur=zeros(siz,1); 
    4341      for c=1:siz; 
    44            vecteur(c)=olr1(a,b,c);     
     42           vecteur(c)=olr1(a,b,c); 
    4543      end 
    4644      std_c=nanmean(vecteur); 
    4745      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 
    5048   end 
    51 end;    
    52  
    53 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
     49end; 
     50 
    5451% on met sous la forme longitude latitude jour 
    55  
    5652 
    5753% on selectionne les mois d ete (attention aux annees bixestiles!) 
     
    153149%olr=zeros(37,144,92,28); 
    154150 
    155  
    156 olr1=permute(olr,[3 1 2 4]); clear olr; 
    157 olr(:,:,:)=nanmean(olr1);  
    158  
     151olr1=permute(olr,[3 1 2 4]); 
     152clear olr; 
     153olr(:,:,:)=nanmean(olr1); 
    159154 
    160155load eof3_noaa_annuel.txt 
     
    167162correlations_eof3=zeros(37,144); 
    168163 
    169 lat=lat(5:41)';  
     164lat=lat(5:41)'; 
    170165lon=lon-180; 
    171166 
     
    177172           y(c)=olr(a,b,c); 
    178173        end 
    179         clear coco; coco=corrcoef(y,eof1_noaa_annuel); 
     174        clear coco; 
     175        coco=corrcoef(y,eof1_noaa_annuel); 
    180176        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); 
    182179        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); 
    184182        correlations_eof3(a,b)=coco(2,1); 
    185183 
    186     end     
     184    end 
    187185end 
    188186 
    189187subplot(3,1,1); 
    190 correlations_eof1(1,1)=1.0001; correlations_eof1(1,2)=-1.0001; 
     188correlations_eof1(1,1)=1.0001; 
     189correlations_eof1(1,2)=-1.0001; 
    191190 
    192191contourf(lon,lat,correlations_eof1,[-1:0.1:1],'LineStyle','none'); 
     
    206205subplot(3,1,2); 
    207206palette 
    208 correlations_eof2(1,1)=1.0001; correlations_eof2(1,2)=-1.0001; 
     207correlations_eof2(1,1)=1.0001; 
     208correlations_eof2(1,2)=-1.0001; 
    209209contourf(lon,lat,correlations_eof2,[-1:0.1:1],'LineStyle','none'); 
    210210colorbar; 
     
    224224subplot(3,1,3); 
    225225palette; 
    226 correlations_eof3(1,1)=1.0001; correlations_eof3(1,2)=-1.0001; 
     226correlations_eof3(1,1)=1.0001; 
     227correlations_eof3(1,2)=-1.0001; 
    227228contourf(lon,lat,correlations_eof3,[-1:0.1:1],'LineStyle','none'); 
    228229set (gca,'tickDir','out'); 
Note: See TracChangeset for help on using the changeset viewer.