Changeset 91 for trunk/src/SIMULS_IRCAAM


Ignore:
Timestamp:
03/16/14 20:25:53 (10 years ago)
Author:
pinsard
Message:

fix thanks to coding rules

Location:
trunk/src/SIMULS_IRCAAM
Files:
34 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/SIMULS_IRCAAM/anomalies_olr_AFR_1979_2000.m

    r26 r91  
    107107olr(:,30,:,:)=rlut; clear rlut; 
    108108 
    109 % selction du domaine geographique 
     109% sélection du domaine géographique 
    110110olr1=zeros(122,30,64,128); 
    111111olr1(:,:,:,1:64)=olr(:,:,:,65:128); 

  • trunk/src/SIMULS_IRCAAM/carte_eof_df.m

    r85 r91  
    123123 
    124124fullfilename=[IRCAAM_ID,'eof_',ircaam_dataset,'_', simulation, '_', df, '.mat']; 
    125 try  
     125try 
    126126 status=load(fullfilename,'C','L'); 
    127127catch 

  • trunk/src/SIMULS_IRCAAM/climato_olr_AFR.m

    r24 r91  
    9393olr(:,29,:,:)=rlut; clear rlut; 
    9494 
    95 % selction du domaine geographique 
     95% sélection du domaine géographique 
    9696olr1=zeros(122,29,64,128); 
    9797olr1(:,:,:,1:64)=olr(:,:,:,65:128); 

  • trunk/src/SIMULS_IRCAAM/climato_olr_CTL.m

    r22 r91  
    9292olr(:,29,:,:)=rlut; clear rlut; 
    9393 
    94 % selction du domaine geographique 
     94% sélection du domaine géographique 
    9595olr1=zeros(122,29,64,128); 
    9696olr1(:,:,:,1:64)=olr(:,:,:,65:128); 

  • trunk/src/SIMULS_IRCAAM/climato_olr_noaa.m

    r22 r91  
    77 
    88% taille 12327x73x144 1er juin 1974-31 decembre2007 
    9 % on enleve les annees 1974 a 1978 nombre de jours a enlever au debut 214+365+366+365+365 
    10 % on enleve les annees 2007 a 2008 nombre de jours a enlever au debut 12327-(365+60) 
     9% on enleve les années 1974 a 1978 nombre de jours a enlever au debut 214+365+366+365+365 
     10% on enleve les années 2007 a 2008 nombre de jours a enlever au debut 12327-(365+60) 
    1111 
    1212debut=214+365+366+365+365; 
    … …  
    2929clear olr; 
    3030 
    31 % on permture le temps et l'espce 
    32 % nouvelles dimentions 
     31% on permute le temps et l'espace 
     32% nouvelles dimensions 
    3333% 21 144 siz 
    3434olr1=permute(olr1, [2 3 1]); 
    … …  
    3636% on met sous la forme longitude latitude jour 
    3737 
    38 % on selectionne les mois d ete (attention aux annees bixestiles!) 
     38% on selectionne les mois d'été (attention aux années bissextiles!) 
    3939 
    4040siz1=31+28+31+30+31+1; 

  • trunk/src/SIMULS_IRCAAM/composite_hgt500_AFR_eof1_1030.m

    r26 r91  
    1111ncload('zg_d.AfNQIVIV.197106-09.nc'); 
    1212 
    13 % selction du domaine geographique 
     13% sélection du domaine géographique 
    1414hgt500b=zeros(122,29,64,128); 
    1515hgt500b(:,:,:,1:64)=hgt500(:,:,:,65:128); 
    … …  
    7272      % 1ere condition (laisser une marge de 10) 
    7373      if (b>20 && b<112); 
    74 % cas 1 : max deltas apres negetifs et deltas avant postifs 
     74% cas 1 : max deltas aprÚs negetifs et deltas avant postifs 
    7575         if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && hgt500_eof(b,a)>(hgt500bis)); 
    7676            compteurmax=compteurmax+1; 
    … …  
    108108            hgt500_filtre_compomax(compteurmax,30,:,:)=hgt500(b+9,a,:,:); 
    109109            hgt500_filtre_compomax(compteurmax,31,:,:)=hgt500(b+10,a,:,:); 
    110 % cas 2 : min deltas apres positifs et deltas avant negatifs 
     110% cas 2 : min deltas aprÚs positifs et deltas avant negatifs 
    111111         elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && hgt500_eof(b,a)<(hgt500bis)); 
    112112            compteurmin=compteurmin+1; 
    … …  
    153153 
    154154% on fait la somme des valeurs de hgt500_filtreentiel sur le nbre d occurences 
    155 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 
     155% on obtent un tableau de la forme (jours-avant-aprÚs,lon,lat) 
    156156hgt500_filtre_compomin1(:,:,:)=mean(hgt500_filtre_compomin); 
    157157clear hgt500_filtre_compomin; 

  • trunk/src/SIMULS_IRCAAM/composite_hgt500_AFR_eof1_30100.m

    r26 r91  
    1111ncload('zg_d.AfNQIVIV.197106-09.nc'); 
    1212 
    13 % selction du domaine geographique 
     13% sélection du domaine géographique 
    1414hgt500b=zeros(122,29,64,128); 
    1515hgt500b(:,:,:,1:64)=hgt500(:,:,:,65:128); 
    … …  
    7272      % 1ere condition (laisser une marge de 10) 
    7373      if (b>20 && b<112); 
    74 % cas 1 : max deltas apres negetifs et deltas avant postifs 
     74% cas 1 : max deltas aprÚs negatifs et deltas avant postifs 
    7575         if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && hgt500_eof(b,a)>(hgt500bis)); 
    7676            compteurmax=compteurmax+1; 
    … …  
    108108            hgt500_filtre_compomax(compteurmax,30,:,:)=hgt500(b+9,a,:,:); 
    109109            hgt500_filtre_compomax(compteurmax,31,:,:)=hgt500(b+10,a,:,:); 
    110 % cas 2 : min deltas apres positifs et deltas avant negatifs 
     110% cas 2 : min deltas aprÚs positifs et deltas avant negatifs 
    111111         elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && hgt500_eof(b,a)<(hgt500bis)); 
    112112            compteurmin=compteurmin+1; 
    … …  
    153153 
    154154% on fait la somme des valeurs de hgt500_filtreentiel sur le nbre d occurences 
    155 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 
     155% on obtent un tableau de la forme (jours-avant-aprÚs,lon,lat) 
    156156hgt500_filtre_compomin1(:,:,:)=mean(hgt500_filtre_compomin); 
    157157clear hgt500_filtre_compomin; 

  • trunk/src/SIMULS_IRCAAM/composite_hgt500_AFR_eof2_1030.m

    r26 r91  
    1111ncload('zg_d.AfNQIVIV.197106-09.nc'); 
    1212 
    13 % selction du domaine geographique 
     13% sélection du domaine géographique 
    1414hgt500b=zeros(122,29,64,128); 
    1515hgt500b(:,:,:,1:64)=hgt500(:,:,:,65:128); 
    … …  
    7272      % 1ere condition (laisser une marge de 10) 
    7373      if (b>20 && b<112); 
    74 % cas 1 : max deltas apres negetifs et deltas avant postifs 
     74% cas 1 : max deltas aprÚs negetifs et deltas avant postifs 
    7575         if (delta(b+1)<0 && delta(b)>0 && hgt500_eof(b,a)>(hgt500bis)); 
    7676            compteurmax=compteurmax+1; 
    … …  
    108108            hgt500_filtre_compomax(compteurmax,30,:,:)=hgt500(b+9,a,:,:); 
    109109            hgt500_filtre_compomax(compteurmax,31,:,:)=hgt500(b+10,a,:,:); 
    110 % cas 2 : min deltas apres positifs et deltas avant negatifs 
     110% cas 2 : min deltas aprÚs positifs et deltas avant negatifs 
    111111         elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && hgt500_eof(b,a)<(-hgt500bis)); 
    112112            compteurmin=compteurmin+1; 
    … …  
    168168 
    169169% on fait la somme des valeurs de hgt500_filtreentiel sur le nbre d occurences 
    170 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 
     170% on obtent un tableau de la forme (jours-avant-aprÚs,lon,lat) 
    171171hgt500_filtre_compomin1(:,:,:)=mean(hgt500_filtre_compomin); 
    172172clear hgt500_filtre_compomin; 

  • trunk/src/SIMULS_IRCAAM/composite_hgt500_AFR_eof2_1030_min.m

    r26 r91  
    1111ncload('zg_d.AfNQIVIV.197106-09.nc'); 
    1212 
    13 % selction du domaine geographique 
     13% sélection du domaine géographique 
    1414hgt500b=zeros(122,29,64,128); 
    1515hgt500b(:,:,:,1:64)=hgt500(:,:,:,65:128); 

  • trunk/src/SIMULS_IRCAAM/composite_hgt500_AFR_eof2_30100.m

    r26 r91  
    1111ncload('zg_d.AfNQIVIV.197106-09.nc'); 
    1212 
    13 % selction du domaine geographique 
     13% sélection du domaine géographique 
    1414hgt500b=zeros(122,29,64,128); 
    1515hgt500b(:,:,:,1:64)=hgt500(:,:,:,65:128); 
    … …  
    7272      % 1ere condition (laisser une marge de 10) 
    7373      if (b>20 && b<112); 
    74 % cas 1 : max deltas apres negetifs et deltas avant postifs 
     74% cas 1 : max deltas aprÚs negetifs et deltas avant postifs 
    7575         if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && hgt500_eof(b,a)>(hgt500bis)); 
    7676            compteurmax=compteurmax+1; 
    … …  
    108108            hgt500_filtre_compomax(compteurmax,30,:,:)=hgt500(b+9,a,:,:); 
    109109            hgt500_filtre_compomax(compteurmax,31,:,:)=hgt500(b+10,a,:,:); 
    110 % cas 2 : min deltas apres positifs et deltas avant negatifs 
     110% cas 2 : min deltas aprÚs positifs et deltas avant negatifs 
    111111         elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && hgt500_eof(b,a)<(hgt500bis)); 
    112112            compteurmin=compteurmin+1; 
    … …  
    153153 
    154154% on fait la somme des valeurs de hgt500_filtreentiel sur le nbre d occurences 
    155 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 
     155% on obtent un tableau de la forme (jours-avant-aprÚs,lon,lat) 
    156156hgt500_filtre_compomin1(:,:,:)=mean(hgt500_filtre_compomin); 
    157157clear hgt500_filtre_compomin; 

  • trunk/src/SIMULS_IRCAAM/composite_hgt500_TROP_eof1_1030.m

    r26 r91  
    1111ncload('zg_d.TrNQIVIV.197106-09.nc'); 
    1212 
    13 % selction du domaine geographique 
     13% sélection du domaine géographique 
    1414hgt500b=zeros(122,29,64,128); 
    1515hgt500b(:,:,:,1:64)=hgt500(:,:,:,65:128); 
    … …  
    7272      % 1ere condition (laisser une marge de 10) 
    7373      if (b>20 && b<112); 
    74 % cas 1 : max deltas apres negetifs et deltas avant postifs 
     74% cas 1 : max deltas aprÚs negetifs et deltas avant postifs 
    7575         if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && hgt500_eof(b,a)>(hgt500bis)); 
    7676            compteurmax=compteurmax+1; 
    … …  
    108108            hgt500_filtre_compomax(compteurmax,30,:,:)=hgt500(b+9,a,:,:); 
    109109            hgt500_filtre_compomax(compteurmax,31,:,:)=hgt500(b+10,a,:,:); 
    110 % cas 2 : min deltas apres positifs et deltas avant negatifs 
     110% cas 2 : min deltas aprÚs positifs et deltas avant negatifs 
    111111         elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && hgt500_eof(b,a)<(hgt500bis)); 
    112112            compteurmin=compteurmin+1; 
    … …  
    153153 
    154154% on fait la somme des valeurs de hgt500_filtreentiel sur le nbre d occurences 
    155 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 
     155% on obtent un tableau de la forme (jours-avant-aprÚs,lon,lat) 
    156156hgt500_filtre_compomin1(:,:,:)=mean(hgt500_filtre_compomin); 
    157157clear hgt500_filtre_compomin; 

  • trunk/src/SIMULS_IRCAAM/composite_hgt500_TROP_eof1_30100.m

    r26 r91  
    1111ncload('zg_d.TrNQIVIV.197106-09.nc','lon','lat'); 
    1212 
    13 % selction du domaine geographique 
     13% sélection du domaine géographique 
    1414hgt500b=zeros(122,29,64,128); 
    1515hgt500b(:,:,:,1:64)=hgt500(:,:,:,65:128); 
    … …  
    7272      % 1ere condition (laisser une marge de 10) 
    7373      if (b>20 && b<112); 
    74 % cas 1 : max deltas apres negetifs et deltas avant postifs 
     74% cas 1 : max deltas aprÚs negetifs et deltas avant postifs 
    7575         if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && hgt500_eof(b,a)>(hgt500bis)); 
    7676            compteurmax=compteurmax+1; 
    … …  
    108108            hgt500_filtre_compomax(compteurmax,30,:,:)=hgt500(b+9,a,:,:); 
    109109            hgt500_filtre_compomax(compteurmax,31,:,:)=hgt500(b+10,a,:,:); 
    110 % cas 2 : min deltas apres positifs et deltas avant negatifs 
     110% cas 2 : min deltas aprÚs positifs et deltas avant negatifs 
    111111         elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && hgt500_eof(b,a)<(hgt500bis)); 
    112112            compteurmin=compteurmin+1; 
    … …  
    153153 
    154154% on fait la somme des valeurs de hgt500_filtreentiel sur le nbre d occurences 
    155 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 
     155% on obtent un tableau de la forme (jours-avant-aprÚs,lon,lat) 
    156156hgt500_filtre_compomin1(:,:,:)=mean(hgt500_filtre_compomin); 
    157157clear hgt500_filtre_compomin; 

  • trunk/src/SIMULS_IRCAAM/composite_hgt500_TROP_eof2_1030.m

    r26 r91  
    1111ncload('zg_d.TrNQIVIV.197106-09.nc'); 
    1212 
    13 % selction du domaine geographique 
     13% sélection du domaine géographique 
    1414hgt500b=zeros(122,29,64,128); 
    1515hgt500b(:,:,:,1:64)=hgt500(:,:,:,65:128); 
    … …  
    7272      % 1ere condition (laisser une marge de 10) 
    7373      if (b>20 && b<112); 
    74 % cas 1 : max deltas apres negetifs et deltas avant postifs 
     74% cas 1 : max deltas aprÚs negetifs et deltas avant postifs 
    7575         if (delta(b+1)<0 && delta(b)>0 && hgt500_eof(b,a)>(hgt500bis)); 
    7676            compteurmax=compteurmax+1; 
    … …  
    108108            hgt500_filtre_compomax(compteurmax,30,:,:)=hgt500(b+9,a,:,:); 
    109109            hgt500_filtre_compomax(compteurmax,31,:,:)=hgt500(b+10,a,:,:); 
    110 % cas 2 : min deltas apres positifs et deltas avant negatifs 
     110% cas 2 : min deltas aprÚs positifs et deltas avant negatifs 
    111111         elseif (delta(b+1)>0 && delta(b)<0 && hgt500_eof(b,a)<(hgt500bis)); 
    112112            compteurmin=compteurmin+1; 
    … …  
    153153 
    154154% on fait la somme des valeurs de hgt500_filtreentiel sur le nbre d occurences 
    155 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 
     155% on obtent un tableau de la forme (jours-avant-aprÚs,lon,lat) 
    156156hgt500_filtre_compomin1(:,:,:)=mean(hgt500_filtre_compomin); 
    157157clear hgt500_filtre_compomin; 

  • trunk/src/SIMULS_IRCAAM/composite_hgt500_TROP_eof2_30100.m

    r26 r91  
    1111ncload('zg_d.TrNQIVIV.197106-09.nc','lon','lat'); 
    1212 
    13 % selction du domaine geographique 
     13% sélection du domaine géographique 
    1414hgt500b=zeros(122,29,64,128); 
    1515hgt500b(:,:,:,1:64)=hgt500(:,:,:,65:128); 
    … …  
    7272      % 1ere condition (laisser une marge de 10) 
    7373      if (b>20 && b<112); 
    74 % cas 1 : max deltas apres negetifs et deltas avant postifs 
     74% cas 1 : max deltas aprÚs negetifs et deltas avant postifs 
    7575         if (delta(b+1)<0 && delta(b)>0 && hgt500_eof(b,a)>(hgt500bis)); 
    7676            compteurmax=compteurmax+1; 
    … …  
    108108            hgt500_filtre_compomax(compteurmax,30,:,:)=hgt500(b+9,a,:,:); 
    109109            hgt500_filtre_compomax(compteurmax,31,:,:)=hgt500(b+10,a,:,:); 
    110 % cas 2 : min deltas apres positifs et deltas avant negatifs 
     110% cas 2 : min deltas aprÚs positifs et deltas avant negatifs 
    111111         elseif (delta(b+1)>0 && delta(b)<0 && hgt500_eof(b,a)<(hgt500bis)); 
    112112            compteurmin=compteurmin+1; 
    … …  
    153153 
    154154% on fait la somme des valeurs de hgt500_filtreentiel sur le nbre d occurences 
    155 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 
     155% on obtent un tableau de la forme (jours-avant-aprÚs,lon,lat) 
    156156hgt500_filtre_compomin1(:,:,:)=mean(hgt500_filtre_compomin); 
    157157clear hgt500_filtre_compomin; 

  • trunk/src/SIMULS_IRCAAM/composite_olr_AFR_eof1_1030.m

    r26 r91  
    9393olr(:,29,:,:)=rlut; clear rlut; 
    9494 
    95 % selction du domaine geographique 
     95% sélection du domaine géographique 
    9696olr1=zeros(122,29,64,128); 
    9797olr1(:,:,:,1:64)=olr(:,:,:,65:128); 
    … …  
    133133      % 1ere condition (laisser une marge de 10) 
    134134      if (b>20 && b<112); 
    135 % cas 1 : max deltas apres negetifs et deltas avant postifs 
     135% cas 1 : max deltas aprÚs negetifs et deltas avant postifs 
    136136         if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 
    137137            compteurmax=compteurmax+1; 
    … …  
    169169            olr_filtre_compomax(compteurmax,30,:,:)=olr1(b+9,a,:,:); 
    170170            olr_filtre_compomax(compteurmax,31,:,:)=olr1(b+10,a,:,:); 
    171 % cas 2 : min deltas apres positifs et deltas avant negatifs 
     171% cas 2 : min deltas aprÚs positifs et deltas avant negatifs 
    172172         elseif (delta(b+1)>0 && delta(b)<0 && olr_eof(b,a)<(-olrbis)); 
    173173            compteurmin=compteurmin+1; 
    … …  
    224224 
    225225% on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 
    226 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 
     226% on obtent un tableau de la forme (jours-avant-aprÚs,lon,lat) 
    227227olr_filtre_compomin1(:,:,:)=nanmean(olr_filtre_compomin); 
    228228clear olr_filtre_compomin; 

  • trunk/src/SIMULS_IRCAAM/composite_olr_AFR_eof1_30100.m

    r26 r91  
    9393olr(:,29,:,:)=rlut; clear rlut; 
    9494 
    95 % selction du domaine geographique 
     95% sélection du domaine géographique 
    9696olr1=zeros(122,29,64,128); 
    9797olr1(:,:,:,1:64)=olr(:,:,:,65:128); 
    … …  
    133133      % 1ere condition (laisser une marge de 10) 
    134134      if (b>20 && b<112); 
    135 % cas 1 : max deltas apres negetifs et deltas avant postifs 
     135% cas 1 : max deltas aprÚs negetifs et deltas avant postifs 
    136136         if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 
    137137            compteurmax=compteurmax+1; 
    … …  
    169169            olr_filtre_compomax(compteurmax,30,:,:)=olr1(b+9,a,:,:); 
    170170            olr_filtre_compomax(compteurmax,31,:,:)=olr1(b+10,a,:,:); 
    171 % cas 2 : min deltas apres positifs et deltas avant negatifs 
     171% cas 2 : min deltas aprÚs positifs et deltas avant negatifs 
    172172         elseif (delta(b+1)>0 && delta(b)<0 && olr_eof(b,a)<(-olrbis)); 
    173173            compteurmin=compteurmin+1; 
    … …  
    224224 
    225225% on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 
    226 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 
     226% on obtent un tableau de la forme (jours-avant-aprÚs,lon,lat) 
    227227olr_filtre_compomin1(:,:,:)=nanmean(olr_filtre_compomin); 
    228228clear olr_filtre_compomin; 

  • trunk/src/SIMULS_IRCAAM/composite_olr_AFR_eof2_1030.m

    r26 r91  
    9393olr(:,29,:,:)=rlut; clear rlut; 
    9494 
    95 % selction du domaine geographique 
     95% sélection du domaine géographique 
    9696olr1=zeros(122,29,64,128); 
    9797olr1(:,:,:,1:64)=olr(:,:,:,65:128); 
    … …  
    133133      % 1ere condition (laisser une marge de 10) 
    134134      if (b>20 && b<112); 
    135 % cas 1 : max deltas apres negetifs et deltas avant postifs 
     135% cas 1 : max deltas aprÚs negetifs et deltas avant postifs 
    136136         if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 
    137137            compteurmax=compteurmax+1; 
    … …  
    169169            olr_filtre_compomax(compteurmax,30,:,:)=olr1(b+9,a,:,:); 
    170170            olr_filtre_compomax(compteurmax,31,:,:)=olr1(b+10,a,:,:); 
    171 % cas 2 : min deltas apres positifs et deltas avant negatifs 
     171% cas 2 : min deltas aprÚs positifs et deltas avant negatifs 
    172172         elseif (delta(b+1)>0 && delta(b)<0 && olr_eof(b,a)<(-olrbis)); 
    173173            compteurmin=compteurmin+1; 
    … …  
    228228 
    229229% on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 
    230 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 
     230% on obtent un tableau de la forme (jours-avant-aprÚs,lon,lat) 
    231231olr_filtre_compomin1(:,:,:)=nanmean(olr_filtre_compomin); 
    232232clear olr_filtre_compomin; 

  • trunk/src/SIMULS_IRCAAM/composite_olr_AFR_eof2_30100.m

    r26 r91  
    9393olr(:,29,:,:)=rlut; clear rlut; 
    9494 
    95 % selction du domaine geographique 
     95% sélection du domaine géographique 
    9696olr1=zeros(122,29,64,128); 
    9797olr1(:,:,:,1:64)=olr(:,:,:,65:128); 
    … …  
    133133      % 1ere condition (laisser une marge de 10) 
    134134      if (b>20 && b<112); 
    135 % cas 1 : max deltas apres negetifs et deltas avant postifs 
     135% cas 1 : max deltas aprÚs negetifs et deltas avant postifs 
    136136         if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 
    137137            compteurmax=compteurmax+1; 
    … …  
    169169            olr_filtre_compomax(compteurmax,30,:,:)=olr1(b+9,a,:,:); 
    170170            olr_filtre_compomax(compteurmax,31,:,:)=olr1(b+10,a,:,:); 
    171 % cas 2 : min deltas apres positifs et deltas avant negatifs 
     171% cas 2 : min deltas aprÚs positifs et deltas avant negatifs 
    172172         elseif (delta(b+1)>0 && delta(b)<0 && olr_eof(b,a)<(-olrbis)); 
    173173            compteurmin=compteurmin+1; 
    … …  
    224224 
    225225% on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 
    226 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 
     226% on obtent un tableau de la forme (jours-avant-aprÚs,lon,lat) 
    227227olr_filtre_compomin1(:,:,:)=nanmean(olr_filtre_compomin); 
    228228clear olr_filtre_compomin; 

  • trunk/src/SIMULS_IRCAAM/composite_olr_AS_eof1_1030.m

    r26 r91  
    9393olr(:,29,:,:)=rlut; clear rlut; 
    9494 
    95 % selction du domaine geographique 
     95% sélection du domaine géographique 
    9696olr1=zeros(122,29,64,128); 
    9797olr1(:,:,:,1:64)=olr(:,:,:,65:128); 
    … …  
    133133      % 1ere condition (laisser une marge de 10) 
    134134      if (b>20 && b<112); 
    135 % cas 1 : max deltas apres negetifs et deltas avant postifs 
     135% cas 1 : max deltas aprÚs negetifs et deltas avant postifs 
    136136         if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 
    137137            compteurmax=compteurmax+1; 
    … …  
    169169            olr_filtre_compomax(compteurmax,30,:,:)=olr1(b+9,a,:,:); 
    170170            olr_filtre_compomax(compteurmax,31,:,:)=olr1(b+10,a,:,:); 
    171 % cas 2 : min deltas apres positifs et deltas avant negatifs 
     171% cas 2 : min deltas aprÚs positifs et deltas avant negatifs 
    172172         elseif (delta(b+1)>0 && delta(b)<0  && olr_eof(b,a)<(-olrbis)); 
    173173            compteurmin=compteurmin+1; 
    … …  
    224224 
    225225% on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 
    226 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 
     226% on obtent un tableau de la forme (jours-avant-aprÚs,lon,lat) 
    227227olr_filtre_compomin1(:,:,:)=nanmean(olr_filtre_compomin); 
    228228clear olr_filtre_compomin; 

  • trunk/src/SIMULS_IRCAAM/composite_olr_AS_eof1_30100.m

    r26 r91  
    9393olr(:,29,:,:)=rlut; clear rlut; 
    9494 
    95 % selction du domaine geographique 
     95% sélection du domaine géographique 
    9696olr1=zeros(122,29,64,128); 
    9797olr1(:,:,:,1:64)=olr(:,:,:,65:128); 
    … …  
    133133      % 1ere condition (laisser une marge de 10) 
    134134      if (b>20 && b<112); 
    135 % cas 1 : max deltas apres negetifs et deltas avant postifs 
     135% cas 1 : max deltas aprÚs negetifs et deltas avant postifs 
    136136         if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 
    137137            compteurmax=compteurmax+1; 
    … …  
    169169            olr_filtre_compomax(compteurmax,30,:,:)=olr1(b+9,a,:,:); 
    170170            olr_filtre_compomax(compteurmax,31,:,:)=olr1(b+10,a,:,:); 
    171 % cas 2 : min deltas apres positifs et deltas avant negatifs 
     171% cas 2 : min deltas aprÚs positifs et deltas avant negatifs 
    172172         elseif (delta(b+1)>0 && delta(b) &&  olr_eof(b,a)<(-olrbis)); 
    173173            compteurmin=compteurmin+1; 
    … …  
    224224 
    225225% on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 
    226 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 
     226% on obtent un tableau de la forme (jours-avant-aprÚs,lon,lat) 
    227227olr_filtre_compomin1(:,:,:)=nanmean(olr_filtre_compomin); 
    228228clear olr_filtre_compomin; 

  • trunk/src/SIMULS_IRCAAM/composite_olr_AS_eof2_1030.m

    r26 r91  
    9393olr(:,29,:,:)=rlut; clear rlut; 
    9494 
    95 % selction du domaine geographique 
     95% sélection du domaine géographique 
    9696olr1=zeros(122,29,64,128); 
    9797olr1(:,:,:,1:64)=olr(:,:,:,65:128); 
    … …  
    133133      % 1ere condition (laisser une marge de 10) 
    134134      if (b>20 && b<112); 
    135 % cas 1 : max deltas apres negetifs et deltas avant postifs 
     135% cas 1 : max deltas aprÚs negetifs et deltas avant postifs 
    136136         if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 
    137137            compteurmax=compteurmax+1; 
    … …  
    169169            olr_filtre_compomax(compteurmax,30,:,:)=olr1(b+9,a,:,:); 
    170170            olr_filtre_compomax(compteurmax,31,:,:)=olr1(b+10,a,:,:); 
    171 % cas 2 : min deltas apres positifs et deltas avant negatifs 
     171% cas 2 : min deltas aprÚs positifs et deltas avant negatifs 
    172172         elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && olr_eof(b,a)<(-olrbis)); 
    173173            compteurmin=compteurmin+1; 
    … …  
    224224 
    225225% on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 
    226 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 
     226% on obtent un tableau de la forme (jours-avant-aprÚs,lon,lat) 
    227227olr_filtre_compomin1(:,:,:)=nanmean(olr_filtre_compomin); 
    228228clear olr_filtre_compomin; 

  • trunk/src/SIMULS_IRCAAM/composite_olr_AS_eof2_30100.m

    r26 r91  
    9393olr(:,29,:,:)=rlut; clear rlut; 
    9494 
    95 % selction du domaine geographique 
     95% sélection du domaine géographique 
    9696olr1=zeros(122,29,64,128); 
    9797olr1(:,:,:,1:64)=olr(:,:,:,65:128); 
    … …  
    133133      % 1ere condition (laisser une marge de 10) 
    134134      if (b>20 && b<112); 
    135 % cas 1 : max deltas apres negetifs et deltas avant postifs 
     135% cas 1 : max deltas aprÚs negetifs et deltas avant postifs 
    136136         if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 
    137137            compteurmax=compteurmax+1; 
    … …  
    169169            olr_filtre_compomax(compteurmax,30,:,:)=olr1(b+9,a,:,:); 
    170170            olr_filtre_compomax(compteurmax,31,:,:)=olr1(b+10,a,:,:); 
    171 % cas 2 : min deltas apres positifs et deltas avant negatifs 
     171% cas 2 : min deltas aprÚs positifs et deltas avant negatifs 
    172172         elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && olr_eof(b,a)<(-olrbis)); 
    173173            compteurmin=compteurmin+1; 
    … …  
    224224 
    225225% on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 
    226 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 
     226% on obtent un tableau de la forme (jours-avant-aprÚs,lon,lat) 
    227227olr_filtre_compomin1(:,:,:)=nanmean(olr_filtre_compomin); 
    228228clear olr_filtre_compomin; 

  • trunk/src/SIMULS_IRCAAM/composite_olr_CTL_eof1_1030.m

    r26 r91  
    9393olr(:,29,:,:)=rlut; clear rlut; 
    9494 
    95 % selction du domaine geographique 
     95% sélection du domaine géographique 
    9696olr1=zeros(122,29,64,128); 
    9797olr1(:,:,:,1:64)=olr(:,:,:,65:128); 

  • trunk/src/SIMULS_IRCAAM/composite_olr_CTL_eof1_30100.m

    r26 r91  
    9393olr(:,29,:,:)=rlut; clear rlut; 
    9494 
    95 % selction du domaine geographique 
     95% sélection du domaine géographique 
    9696olr1=zeros(122,29,64,128); 
    9797olr1(:,:,:,1:64)=olr(:,:,:,65:128); 

  • trunk/src/SIMULS_IRCAAM/composite_olr_CTL_eof2_1030.m

    r26 r91  
    9393olr(:,29,:,:)=rlut; clear rlut; 
    9494 
    95 % selction du domaine geographique 
     95% sélection du domaine géographique 
    9696olr1=zeros(122,29,64,128); 
    9797olr1(:,:,:,1:64)=olr(:,:,:,65:128); 

  • trunk/src/SIMULS_IRCAAM/composite_olr_CTL_eof2_30100.m

    r26 r91  
    9393olr(:,29,:,:)=rlut; clear rlut; 
    9494 
    95 % selction du domaine geographique 
     95% sélection du domaine géographique 
    9696olr1=zeros(122,29,64,128); 
    9797olr1(:,:,:,1:64)=olr(:,:,:,65:128); 

  • trunk/src/SIMULS_IRCAAM/composite_olr_TR_eof1_1030.m

    r26 r91  
    9393olr(:,29,:,:)=rlut; clear rlut; 
    9494 
    95 % selction du domaine geographique 
     95% sélection du domaine géographique 
    9696olr1=zeros(122,29,64,128); 
    9797olr1(:,:,:,1:64)=olr(:,:,:,65:128); 

  • trunk/src/SIMULS_IRCAAM/composite_olr_TR_eof1_30100.m

    r26 r91  
    9696olr(:,30,:,:)=rlut; clear rlut; 
    9797 
    98 % selction du domaine geographique 
     98% sélection du domaine géographique 
    9999olr1=zeros(122,30,64,128); 
    100100olr1(:,:,:,1:64)=olr(:,:,:,65:128); 

  • trunk/src/SIMULS_IRCAAM/composite_olr_TR_eof2_1030.m

    r26 r91  
    9696olr(:,30,:,:)=rlut; clear rlut; 
    9797 
    98 % selction du domaine geographique 
     98% sélection du domaine géographique 
    9999olr1=zeros(122,30,64,128); 
    100100olr1(:,:,:,1:64)=olr(:,:,:,65:128); 

  • trunk/src/SIMULS_IRCAAM/composite_olr_TR_eof2_30100.m

    r26 r91  
    9696olr(:,30,:,:)=rlut; clear rlut; 
    9797 
    98 % selction du domaine geographique 
     98% sélection du domaine géographique 
    9999olr1=zeros(122,30,64,128); 
    100100olr1(:,:,:,1:64)=olr(:,:,:,65:128); 

  • trunk/src/SIMULS_IRCAAM/correlations_eofs_olr_jas.m

    r26 r91  
    77 
    88% taille 12327x73x144 1er juin 1974-31 decembre2007 
    9 % on enleve les annees 1974 a 1978 nombre de jours a enlever au debut 214+365+366+365+365 
    10 % on enleve les annees 2007 a 2008 nombre de jours a enlever au debut 12327-(365+60) 
     9% on enleve les années 1974 a 1978 nombre de jours a enlever au debut 214+365+366+365+365 
     10% on enleve les années 2007 a 2008 nombre de jours a enlever au debut 12327-(365+60) 
    1111 
    1212debut=214+365+366+365+365; 
    … …  
    3030 
    3131% on permture le temps et l'espce 
    32 % nouvelles dimentions 
     32% nouvelles dimensions 
    3333% 21 144 siz 
    3434olr1=permute(olr1, [2 3 1]); 
    … …  
    5050% on met sous la forme longitude latitude jour 
    5151 
    52 % on selectionne les mois d ete (attention aux annees bixestiles!) 
     52% on sélectionne les mois d'été (attention aux années bixestiles!) 
    5353 
    5454siz1=31+28+31+30+31+30+1; 

  • trunk/src/SIMULS_IRCAAM/correlations_eofs_olr_jas_20.m

    r26 r91  
    88 
    99% taille 12327x73x144 1er juin 1974-31 decembre2007 
    10 % on enleve les annees 1974 a 1978 nombre de jours a enlever au debut 214+365+366+365+365 
    11 % on enleve les annees 2007 a 2008 nombre de jours a enlever au debut 12327-(365+60) 
     10% on enleve les années 1974 a 1978 nombre de jours a enlever au debut 214+365+366+365+365 
     11% on enleve les années 2007 a 2008 nombre de jours a enlever au debut 12327-(365+60) 
    1212 
    1313debut=214+365+366+365+365; 
    … …  
    3030clear olr; 
    3131 
    32 % on permture le temps et l'espce 
    33 % nouvelles dimentions 
     32% on permture le temps et l'espace 
     33% nouvelles dimensions 
    3434% 21 144 siz 
    3535olr1=permute(olr1, [2 3 1]); 
    … …  
    5151% on met sous la forme longitude latitude jour 
    5252 
    53 % on selectionne les mois d ete (attention aux annees bixestiles!) 
     53% on selectionne les mois d ete (attention aux années bixestiles!) 
    5454 
    5555siz1=31+28+31+30+31+30+1; 

  • trunk/src/SIMULS_IRCAAM/prepare_olr_AFR_annuel.m

    r29 r91  
    9797olr(:,30,:,:)=rlut; clear rlut; 
    9898 
    99 % selction du domaine geographique 
     99% sélection du domaine géographique 
    100100olr1=zeros(122,30,64,128); 
    101101olr1(:,:,:,1:64)=olr(:,:,:,65:128); 

  • trunk/src/SIMULS_IRCAAM/prepare_pluie_AFR_annuel.m

    r29 r91  
    9797pluie(:,30,:,:)=pr; clear pr; 
    9898 
    99 % selction du domaine geographique 
     99% sélection du domaine géographique 
    100100pluie1=zeros(122,30,64,128); 
    101101pluie1(:,:,:,1:64)=pluie(:,:,:,65:128); 
Note: See TracChangeset for help on using the changeset viewer.