- Timestamp:
- 01/06/09 10:53:18 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/mode_sahelien/composite_vwnd925_EOF12_ER30100_dat.m
r15 r16 1 % 1 2 % initialisation 2 clear; close all;3 3 clear; 4 close all; 4 5 5 6 % etape 1 :on selectionne les mois d ete (attention aux annees bixestiles!) … … 8 9 vwnd925=zeros(122,49,144,28); 9 10 10 siz=31+28+31+30+31+1; sizbis=31+28+31+30+31+30+31+31+30; 11 siz=31+28+31+30+31+1; 12 sizbis=31+28+31+30+31+30+31+31+30; 11 13 ncload('vwnd925.1979.nc'); vwnd925(:,:,:,1)=vwnd(siz:sizbis,1,:,:); clear vwnd; 12 14 ncload('vwnd925.1980.nc'); vwnd925(:,:,:,2)=vwnd(siz+1:sizbis+1,1,:,:); clear vwnd; … … 48 50 49 51 %%%%%%%%%%%%%%% FICHIER OLR 50 % ouverture / lecture du fichier 52 % ouverture / lecture du fichier 51 53 load eof12_rossby_30100.txt 52 54 ind_olr=reshape(eof12_rossby_30100,122,28); … … 55 57 % on calcule l ecart type sur la serie entiere de l OLR 56 58 ind_olr1=reshape(ind_olr, k,1); 57 ind_olr2=std(ind_olr1); clear ind_olr1;58 59 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 59 ind_olr2=std(ind_olr1); 60 clear ind_olr1; 61 60 62 % GEOPOTENTIEL 61 63 % intialisation des tableaux 62 vwnd925_filtre_compomax=NaN*ones(60,31,49,144); 64 vwnd925_filtre_compomax=NaN*ones(60,31,49,144); 63 65 vwnd925_filtre_compomin=NaN*ones(60,31,49,144); 64 66 65 compteurmax=0; 67 compteurmax=0; 66 68 compteurmin=0; 67 69 68 for a=1:28; % boucle sur les annees 69 clear delta; 70 % boucle sur les annees 71 for a=1:28; 72 clear delta; 70 73 delta=zeros(122,1); 71 for b=1:122; % boucle sur les jours : on prepare avant 74 % boucle sur les jours : on prepare avant 75 for b=1:122; 72 76 if b>1; 73 77 delta(b)=ind_olr(b,a)-ind_olr(b-1,a); 74 78 end; 75 79 end; 76 for b=1:122; % boucle sur les jours 77 if (b>22 && b<112); % 1ere condition (laisser une marge de 10) 78 % cas 1 : max deltas apres negetifs et deltas avant postifs 79 if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && ind_olr(b,a)>ind_olr2); 80 % boucle sur les jours 81 for b=1:122; 82 % 1ere condition (laisser une marge de 10) 83 if (b>22 && b<112); 84 % cas 1 : max deltas apres negetifs et deltas avant postifs 85 if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && ind_olr(b,a)>ind_olr2); 80 86 compteurmax=compteurmax+1; 81 87 vwnd925_filtre_compomax(compteurmax,1,:,:)=vwnd(b-21,a,:,:); … … 113 119 vwnd925_filtre_compomax(compteurmax,31,:,:)=vwnd(b+9,a,:,:); 114 120 115 % cas 2 : min deltas apres positifs et deltas avant negatifs 121 % cas 2 : min deltas apres positifs et deltas avant negatifs 116 122 elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && ind_olr(b,a)<(-ind_olr2)); 117 123 compteurmin=compteurmin+1; … … 149 155 vwnd925_filtre_compomin(compteurmin,30,:,:)=vwnd(b+8,a,:,:); 150 156 vwnd925_filtre_compomin(compteurmin,31,:,:)=vwnd(b+9,a,:,:); 151 157 152 158 end; 153 159 end; 154 end; 160 end; 155 161 end; 156 162 … … 158 164 % on fait la somme des valeurs de vwnd925_filtreentiel sur le nbre d occurences 159 165 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 160 vwnd925_filtre_compomin1(:,:,:)=nanmean(vwnd925_filtre_compomin); clear vwnd925_filtre_compomin; 161 vwnd925_filtre_compomax1(:,:,:)=nanmean(vwnd925_filtre_compomax); clear vwnd925_filtre_compomax; 162 163 x=[-10:1:10]'; %' 164 165 figure(1); 166 vwnd925_filtre_compomin1(:,:,:)=nanmean(vwnd925_filtre_compomin); 167 clear vwnd925_filtre_compomin; 168 vwnd925_filtre_compomax1(:,:,:)=nanmean(vwnd925_filtre_compomax); 169 clear vwnd925_filtre_compomax; 170 171 x=[-10:1:10]'; 172 173 figure(1); 166 174 orient('landscape'); 167 175 lat=lat';
Note: See TracChangeset
for help on using the changeset viewer.