Changeset 91 for trunk


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

fix thanks to coding rules

Location:
trunk
Files:
70 edited

Legend:

Unmodified
Added
Removed

  • trunk/adm/sphinx/conf.py

    r83 r91  
    2020# 
    2121#   * find a way to produce olr_split.sh.1 even if graphviz inside source 
    22 #   * find a way to avoird hard coded man pages name here (should be in 
     22#   * find a way to avoid hard coded man pages name here (should be in 
    2323#     adm/makefile) 
    2424#   * use a specific index.rst for man troff 

  • 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); 

  • trunk/src/cartemonde.m

    r86 r91  
    2020% Data of continents and borders are in :file:`monde.mat`. 
    2121% 
    22 % A appeler à la fin de votre programme 
     22% A appeler à la fin de votre programme 
    2323% 
    2424% Ajuster la fenetre spatiale au besoin avec la 
    … …  
    3030% :: 
    3131% 
    32 % >> ............ (début de votre programme) 
     32% >> ............ (début de votre programme) 
    3333% >> cartemonde; 
    3434% >> axis ([40 100 0 50]); 

  • trunk/src/cartemonde1.m

    r86 r91  
    1 % CARTEMONDE1 trace les contours des continents et les frontieres 
     1% CARTEMONDE1 trace les contours des continents et les frontiÚres 
    22 
    33%+ 
    … …  
    1010% 
    1111% --------------------------------------------------- 
    12 % trace les contours des continents et les frontieres 
     12% trace les contours des continents et les frontiÚres 
    1313% --------------------------------------------------- 
    1414% 
    … …  
    1616% =========== 
    1717% 
    18 % trace les contours des continents et les frontieres. 
     18% trace les contours des continents et les frontiÚres. 
    1919% 
    2020% Data of continents and borders are in :file:`monde.mat`. 
    2121% 
    22 % A appeler à la fin de votre programme 
     22% À appeler à la fin de votre programme 
    2323% 
    24 % Ajuster la fenetre spatiale au besoin avec la 
    25 %  fonction AXIS([longmin longmax latmin latmax]) 
     24% Ajuster la fenêtre spatiale au besoin avec la 
     25% fonction AXIS([longmin longmax latmin latmax]) 
    2626% 
    2727% EXAMPLES 
    2828% ======== 
    2929% 
    30 % :: 
     30% .. code-block:: matlab 
    3131% 
    32 % >> ............ (début de votre programme) 
    33 % >> cartemonde; 
    34 % >> axis ([40 100 0 50]); 
     32%    cartemonde; 
     33%    axis ([40 100 0 50]); 
    3534% 
    3635% SEE ALSO 
    … …  
    7473end 
    7574 
    76 %  tracage des frontieres 
     75%  traçage des frontiÚres 
    7776 [nnl,nnc]=size(npfronti); 
    7877 aa=1; 
    … …  
    8483%        gg=plot (fronti(aa:cc,1),fronti(aa:cc,2)); 
    8584%        set (gg,'Color',[0.5 0.5 0.5],'LineWidth',0.4) 
    86    bb=npfronti(iii); 
    87  end 
     85    bb=npfronti(iii); 
     86end 

  • trunk/src/ircaam_profile.sh

    r90 r91  
    1515# ======== 
    1616# 
    17 # Online usage:: 
    18 # 
    19 #  $ . ./ircaam_profile.sh -d directory -i indir -o outdir -t tempdir 
    20 # 
    21 # In :file:`${HOME}/.profile`, add the following line:: 
    22 # 
    23 #  . ./ircaam_profile.sh -d directory -i indir -o outdir -t tempdir 
     17# Online usage: 
     18# 
     19# .. code-block:: bash 
     20# 
     21#    . ./ircaam_profile.sh -d directory -i indir -o outdir -t tempdir 
     22# 
     23# In :file:`${HOME}/.profile`, add the following line: 
     24# 
     25# .. code-block:: bash 
     26# 
     27#    . ./ircaam_profile.sh -d directory -i indir -o outdir -t tempdir 
    2428# 
    2529# DESCRIPTION 
    … …  
    4852# ======== 
    4953# 
    50 # For fplod, on aedon.locean-ipsl.upmc.fr:: 
    51 # 
    52 #  $ cd /usr/home/fplod/incas/ircaam/ircaam_ws/ 
    53 #  $ . ./src/ircaam_profile.sh \ 
     54# For fplod, on aedon.locean-ipsl.upmc.fr: 
     55# 
     56# .. code-block:: bash 
     57# 
     58#    cd /usr/home/fplod/incas/ircaam/ircaam_ws/ 
     59#    . ./src/ircaam_profile.sh \ 
    5460#    -d $(pwd) \ 
    5561#    -i /usr/temp/${LOGNAME}/ircaam_d/ \ 
    … …  
    5763#    -t /usr/temp/${LOGNAME}/log/ 
    5864# 
    59 # For fplod on zeus.locean-ipsl.umpc.fr:: 
    60 # 
    61 #  $ cd ${HOME}/incas/ircaam/ircaam_ws/ 
    62 #  $ . ./src/ircaam_profile.sh \ 
     65# For fplod on zeus.locean-ipsl.umpc.fr: 
     66# 
     67# .. code-block:: bash 
     68# 
     69#    cd ${HOME}/incas/ircaam/ircaam_ws/ 
     70#    . ./src/ircaam_profile.sh \ 
    6371#    -d $(pwd) \ 
    6472#    -i /usr/temp/${LOGNAME}/ircaam_d/ \ 
    … …  
    7583# ++ besoin de posix 
    7684# 
    77 # ++ pas de MANPATH defini par défaut sur zeus 
     85# ++ pas de MANPATH defini par défaut sur zeus 
    7886# 
    7987# ++ climserv 
    … …  
    129137usage=" Usage : ${command} -d directory -i indir -o outdir -t tempdir" 
    130138# 
    131 while [ ${#} -gt 0 ] 
    132 do 
    133     case ${1} in 
    134         -d) 
    135             # directory for application choosen by user (see svn checkout command used) 
    136             directory=${2} 
    137             shift 
    138         ;; 
    139         -i) 
    140             # directory for inputs choosen by user 
    141             indir=${2} 
    142             shift 
    143         ;; 
    144         -o) 
    145             # directory for outputs choosen by user 
    146             outdir=${2} 
    147             shift 
    148         ;; 
    149         -t) 
    150             # directory for temporary outputs choosen by user 
    151             tempdir=${2} 
    152             shift 
    153         ;; 
    154         *) 
    155             # other choice 
    156             echo "eee : unknown option ${1}" 
    157             echo "${usage}" 
     139pb=0 
     140# 
     141minargcount=8 
     142maxargcount=8 
     143narg=${#} 
     144# 
     145if [ ${narg} -lt ${minargcount} ] 
     146then 
     147    echo "eee : not enought arguments (${narg} vs ${minargcount})" 
     148    echo "${usage}" 
     149    # nb : no exit because this file should be launched by login process 
     150    pb=1 
     151fi 
     152unset minargcount 
     153# 
     154if [ ${narg} -gt ${maxargcount} ] 
     155then 
     156    echo "eee : too many arguments (${narg} vs ${maxargcount})" 
     157    echo "${usage}" 
     158    # nb : no exit because this file should be launched by login process 
     159    pb=1 
     160fi 
     161unset maxargcount 
     162unset narg 
     163# 
     164if [ ${pb} -eq 0 ] 
     165then 
     166    while [ ${#} -gt 0 ] 
     167    do 
     168        case ${1} in 
     169            -d) 
     170                # directory for application choosen by user (see svn checkout command used) 
     171                directory=${2} 
     172                shift 
     173            ;; 
     174            -i) 
     175                # directory for inputs choosen by user 
     176                indir=${2} 
     177                shift 
     178            ;; 
     179            -o) 
     180                # directory for outputs choosen by user 
     181                outdir=${2} 
     182                shift 
     183            ;; 
     184            -t) 
     185                # directory for temporary outputs choosen by user 
     186                tempdir=${2} 
     187                shift 
     188            ;; 
     189            *) 
     190                # other choice 
     191                echo "eee : unknown option ${1}" 
     192                echo "${usage}" 
     193                # nb : no exit because this file should be launched by login process 
     194            ;; 
     195        esac 
     196        # next flag 
     197        shift 
     198    done 
     199    unset usage 
     200fi 
     201if [ ${pb} -eq 0 ] 
     202then 
     203    # 
     204    # check for ${directory} 
     205    if [ ! -d ${directory} ] 
     206    then 
     207        echo " eee : ${directory} not found" 
     208        # nb : no exit because this file should be launched by login process 
     209    fi 
     210    # 
     211    # check for permission on directory 
     212    if [ ! -x ${directory} ] 
     213    then 
     214        echo " eee : ${directory} not reachable" 
     215        # nb : no exit because this file should be launched by login process 
     216    fi 
     217    # 
     218    IRCAAM=${directory} 
     219    export IRCAAM 
     220    unset directory 
     221    # 
     222    # add IRCAAM tools to PATH 
     223    # if not already done 
     224    suppath=$(echo ${IRCAAM} | tr -s "/") 
     225    echo ${PATH} | grep -q "${suppath}:" 
     226    test_path=${?} 
     227    if [ ${test_path} -ne 0 ] 
     228    then 
     229        PATH=${suppath}:${PATH} 
     230        export PATH 
     231    else 
     232        # option bavarde oui/non pas encore implantée ++ 
     233        echo "${command} : iii : ${suppath} already in \${PATH}" 
     234    fi 
     235    unset test_path 
     236    unset suppath 
     237    # 
     238    # add ${IRCAAM}/../bin tools to PATH 
     239    # if not already done 
     240    suppath=$(echo ${IRCAAM}/../bin | tr -s "/") 
     241    echo ${PATH} | grep -q "${suppath}:" 
     242    test_path=${?} 
     243    if [ ${test_path} -ne 0 ] 
     244    then 
     245        PATH=${suppath}:${PATH} 
     246        export PATH 
     247    else 
     248        # option bavarde oui/non pas encore implantée ++ 
     249        echo "${command} : iii : ${suppath} already in \${PATH}" 
     250    fi 
     251    unset test_path 
     252    unset suppath 
     253    # 
     254    # add IRCAAM manuals to MANPATH 
     255    # if not already done 
     256    suppath=$(echo ${IRCAAM}/doc/manuals/man | tr -s "/") 
     257    echo ${MANPATH} | grep -q ${suppath} 
     258    test_manpath=${?} 
     259    if [ ${test_manpath} -ne 0 ] 
     260    then 
     261        MANPATH=${suppath}:${MANPATH} 
     262        export MANPATH 
     263    else 
     264        # option bavarde oui/non pas encore implantée ++ 
     265        echo "${command} : iii : ${suppath} already in \${MANPATH}" 
     266    fi 
     267    unset test_manpath 
     268    unset suppath 
     269    # 
     270    IRCAAM_LOG=${tempdir} 
     271    export IRCAAM_LOG 
     272    unset tempdir 
     273    if [ ! -d ${IRCAAM_LOG} ] 
     274    then 
     275        mkdir -p ${IRCAAM_LOG} 
     276        status=${?} 
     277        if [ ${status} -ne 0 ] 
     278        then 
     279            echo "${command} : eee : can not create \${IRCAAM_LOG}" 
    158280            # nb : no exit because this file should be launched by login process 
    159         ;; 
    160     esac 
    161     # next flag 
    162     shift 
    163 done 
    164 unset usage 
    165 # 
    166 # check for ${directory} 
    167 if [ ! -d ${directory} ] 
    168 then 
    169     echo " eee : ${directory} not found" 
    170     # nb : no exit because this file should be launched by login process 
     281        else 
     282            echo "${command} : iii : creation of \${IRCAAM_LOG}" 
     283        fi 
     284        unset status 
     285    fi 
     286    # check for permission on IRCAAM_LOG 
     287    if [ ! -x ${IRCAAM_LOG} ] 
     288    then 
     289        echo " eee : ${IRCAAM_LOG} not reachable" 
     290        # nb : no exit because this file should be launched by login process 
     291    fi 
     292    # 
     293    # check for permission on IRCAAM_LOG 
     294    if [ ! -w ${IRCAAM_LOG} ] 
     295    then 
     296        echo " eee : ${IRCAAM_LOG} not writable" 
     297        # nb : no exit because this file shouldreachable be launched 
     298        # by login process 
     299    fi 
     300    # 
     301    EDITOR=vi 
     302    export EDITOR 
     303    # 
     304    # io directories 
     305    IRCAAM_ID=${indir} 
     306    export IRCAAM_ID 
     307    unset indir 
     308    if [ ! -d ${IRCAAM_ID} ] 
     309    then 
     310        mkdir -p ${IRCAAM_ID} 
     311        echo "${command} : iii : creation of \${IRCAAM_ID}" 
     312    fi 
     313    # check for permission on IRCAAM_ID 
     314    if [ ! -x ${IRCAAM_ID} ] 
     315    then 
     316        echo " eee : ${IRCAAM_ID} not reachable" 
     317        # nb : no exit because this file should be launched by login process 
     318    fi 
     319    # 
     320    IRCAAM_OD=${outdir} 
     321    export IRCAAM_OD 
     322    unset outdir 
     323    if [ ! -d ${IRCAAM_OD} ] 
     324    then 
     325        mkdir -p ${IRCAAM_OD} 
     326        echo "${command} : iii : creation of \${IRCAAM_OD}" 
     327    fi 
     328    # check for permission on IRCAAM_OD 
     329    if [ ! -x ${IRCAAM_OD} ] 
     330    then 
     331        echo " eee : ${IRCAAM_OD} not reachable" 
     332        # nb : no exit because this file should be launched by login process 
     333    fi 
     334    if [ ! -w ${IRCAAM_OD} ] 
     335    then 
     336        echo " eee : ${IRCAAM_OD} not writable" 
     337        # nb : no exit because this file should be launched by login process 
     338    fi 
    171339fi 
    172 # 
    173 # check for permission on directory 
    174 if [ ! -x ${directory} ] 
    175 then 
    176     echo " eee : ${directory} not reachable" 
    177     # nb : no exit because this file should be launched by login process 
    178 fi 
    179 # 
    180 IRCAAM=${directory} 
    181 export IRCAAM 
    182 unset directory 
    183 # 
    184 # add IRCAAM tools to PATH 
    185 # if not already done 
    186 suppath=$(echo ${IRCAAM} | tr -s "/") 
    187 echo ${PATH} | grep -q "${suppath}:" 
    188 test_path=${?} 
    189 if [ ${test_path} -ne 0 ] 
    190 then 
    191     PATH=${suppath}:${PATH} 
    192     export PATH 
    193 else 
    194     # option bavarde oui/non pas encore implantée ++ 
    195     echo "${command} : iii : ${suppath} already in \${PATH}" 
    196 fi 
    197 unset test_path 
    198 unset suppath 
    199 # 
    200 # add ${IRCAAM}/../bin tools to PATH 
    201 # if not already done 
    202 suppath=$(echo ${IRCAAM}/../bin | tr -s "/") 
    203 echo ${PATH} | grep -q "${suppath}:" 
    204 test_path=${?} 
    205 if [ ${test_path} -ne 0 ] 
    206 then 
    207     PATH=${suppath}:${PATH} 
    208     export PATH 
    209 else 
    210     # option bavarde oui/non pas encore implantée ++ 
    211     echo "${command} : iii : ${suppath} already in \${PATH}" 
    212 fi 
    213 unset test_path 
    214 unset suppath 
    215 # 
    216 # add IRCAAM manuals to MANPATH 
    217 # if not already done 
    218 suppath=$(echo ${IRCAAM}/doc/manuals/man | tr -s "/") 
    219 echo ${MANPATH} | grep -q ${suppath} 
    220 test_manpath=${?} 
    221 if [ ${test_manpath} -ne 0 ] 
    222 then 
    223     MANPATH=${suppath}:${MANPATH} 
    224     export MANPATH 
    225 else 
    226     # option bavarde oui/non pas encore implantée ++ 
    227     echo "${command} : iii : ${suppath} already in \${MANPATH}" 
    228 fi 
    229 unset test_manpath 
    230 unset suppath 
    231 # 
    232 IRCAAM_LOG=${tempdir} 
    233 export IRCAAM_LOG 
    234 unset tempdir 
    235 if [ ! -d ${IRCAAM_LOG} ] 
    236 then 
    237     mkdir -p ${IRCAAM_LOG} 
    238     status=${?} 
    239     if [ ${status} -ne 0 ] 
    240     then 
    241         echo "${command} : eee : can not create \${IRCAAM_LOG}" 
    242         # nb : no exit because this file should be launched by login process 
    243     else 
    244         echo "${command} : iii : creation of \${IRCAAM_LOG}" 
    245     fi 
    246     unset status 
    247 fi 
    248 # check for permission on IRCAAM_LOG 
    249 if [ ! -x ${IRCAAM_LOG} ] 
    250 then 
    251     echo " eee : ${IRCAAM_LOG} not reachable" 
    252     # nb : no exit because this file should be launched by login process 
    253 fi 
    254 # 
    255 # check for permission on IRCAAM_LOG 
    256 if [ ! -w ${IRCAAM_LOG} ] 
    257 then 
    258     echo " eee : ${IRCAAM_LOG} not writable" 
    259     # nb : no exit because this file shouldreachable be launched by login process 
    260 fi 
    261 # 
    262 EDITOR=vi 
    263 export EDITOR 
    264 # 
    265 # io directories 
    266 IRCAAM_ID=${indir} 
    267 export IRCAAM_ID 
    268 unset indir 
    269 if [ ! -d ${IRCAAM_ID} ] 
    270 then 
    271     mkdir -p ${IRCAAM_ID} 
    272     echo "${command} : iii : creation of \${IRCAAM_ID}" 
    273 fi 
    274 # check for permission on IRCAAM_ID 
    275 if [ ! -x ${IRCAAM_ID} ] 
    276 then 
    277     echo " eee : ${IRCAAM_ID} not reachable" 
    278     # nb : no exit because this file should be launched by login process 
    279 fi 
    280 # 
    281 IRCAAM_OD=${outdir} 
    282 export IRCAAM_OD 
    283 unset outdir 
    284 if [ ! -d ${IRCAAM_OD} ] 
    285 then 
    286     mkdir -p ${IRCAAM_OD} 
    287     echo "${command} : iii : creation of \${IRCAAM_OD}" 
    288 fi 
    289 # check for permission on IRCAAM_OD 
    290 if [ ! -x ${IRCAAM_OD} ] 
    291 then 
    292     echo " eee : ${IRCAAM_OD} not reachable" 
    293     # nb : no exit because this file should be launched by login process 
    294 fi 
    295 if [ ! -w ${IRCAAM_OD} ] 
    296 then 
    297     echo " eee : ${IRCAAM_OD} not writable" 
    298     # nb : no exit because this file should be launched by login process 
    299 fi 
     340unset pb 
    300341# 
    301342# end 

  • trunk/src/mode_sahelien/composite_olr_EOF12_ER1030.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; 
    … …  
    3131 
    3232% on permture le temps et l'espce 
    33 % nouvelles dimentions 
     33% nouvelles dimensions 
    3434% 21 144 siz 
    3535olr1=permute(olr1, [2 3 1]); 
    … …  
    5050 
    5151% on met sous la forme longitude latitude jour 
    52 % on selectionne les mois d ete (attention aux annees bixestiles!) 
     52% on selectionne les mois d ete (attention aux années bixestiles!) 
    5353 
    5454siz1=31+28+31+30+31+1; 
    … …  
    168168compteurmin=0; 
    169169 
    170 % boucle sur les annees 
     170% boucle sur les années 
    171171for a=1:28; 
    172172   clear delta; 

  • trunk/src/mode_sahelien/composite_olr_EOF12_ER30100.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; 
    … …  
    3131 
    3232% on permture le temps et l'espce 
    33 % nouvelles dimentions 
     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+1; 
    … …  
    169169compteurmin=0; 
    170170 
    171 % boucle sur les annees 
     171% boucle sur les années 
    172172for a=1:28; 
    173173   clear delta; 

  • trunk/src/mode_sahelien/composite_olr_EOF12_ER30100_dat.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; 
    … …  
    3131 
    3232% on permture le temps et l'espce 
    33 % nouvelles dimentions 
     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+1; 
    … …  
    164164compteurmin=0; 
    165165 
    166 % boucle sur les annees 
     166% boucle sur les années 
    167167for a=1:28; 
    168168   clear delta; 

  • trunk/src/mode_sahelien/composite_olr_EOF23_1030.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; 
    … …  
    3131 
    3232% on permture le temps et l'espce 
    33 % nouvelles dimentions 
     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+1; 
    … …  
    172172compteurmin=0; 
    173173 
    174 % boucle sur les annees 
     174% boucle sur les années 
    175175for a=1:28; 
    176176   clear delta; 

  • trunk/src/mode_sahelien/composite_olr_EOF23_30100.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; 
    … …  
    3131 
    3232% on permture le temps et l'espce 
    33 % nouvelles dimentions 
     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+1; 
    … …  
    164164compteurmin=0; 
    165165 
    166 % boucle sur les annees 
     166% boucle sur les années 
    167167for a=1:28; 
    168168   clear delta; 

  • trunk/src/mode_sahelien/composite_olr_EOF24_1030.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; 
    … …  
    3131 
    3232% on permture le temps et l'espce 
    33 % nouvelles dimentions 
     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+1; 
    … …  
    171171compteurmin=0; 
    172172 
    173 % boucle sur les annees 
     173% boucle sur les années 
    174174for a=1:28; 
    175175   clear delta; 

  • trunk/src/mode_sahelien/composite_olr_EOF24_1030_passbande.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; 
    … …  
    3131 
    3232% on permture le temps et l'espce 
    33 % nouvelles dimentions 
     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+1; 
    … …  
    154154clear v; 
    155155 
    156 % on selectionne les annees de 1979 a 2006 
     156% on selectionne les années de 1979 a 2006 
    157157ind_olr1=ind_olr(:,1:28); 
    158158clear ind_olr; 
    … …  
    179179compteurmin=0; 
    180180 
    181 % boucle sur les annees 
     181% boucle sur les années 
    182182for a=1:28; 
    183183   clear delta; 

  • trunk/src/mode_sahelien/composite_olr_ER1030.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 2001 a 2008 nombre de jours a enlever au debut 12327-(365*3+366+365*3+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 2001 a 2008 nombre de jours a enlever au debut 12327-(365*3+366+365*3+60) 
    1212 
    1313debut=214+365+366+365+365; 
    … …  
    3131 
    3232% on permture le temps et l'espce 
    33 % nouvelles dimentions 
     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+1; 
    … …  
    165165compteurmin=0; 
    166166 
    167 % boucle sur les annees 
     167% boucle sur les années 
    168168for a=1:28; 
    169169   clear delta; 

  • trunk/src/mode_sahelien/composite_olr_eof_df.m

    r85 r91  
    121121 
    122122% taille 12327x73x144 1er juin 1974-31 decembre2007 
    123 % on enleve les annees 1974 a 1978 nombre de jours a enlever au debut 214+365+366+365+365 
    124 % on enleve les annees 2007 a 2008 nombre de jours a enlever au debut 12327-(365+60) 
     123% on enleve les années 1974 a 1978 nombre de jours a enlever au debut 214+365+366+365+365 
     124% on enleve les années 2007 a 2008 nombre de jours a enlever au debut 12327-(365+60) 
    125125 
    126126debut=214+365+366+365+365; 
    … …  
    164164% on met sous la forme longitude latitude jour 
    165165 
    166 % on selectionne les mois d ete (attention aux annees bixestiles!) 
     166% on selectionne les mois d ete (attention aux années bixestiles!) 
    167167 
    168168siz1=31+28+31+30+31+1; 
    … …  
    285285compteurmin=0; 
    286286 
    287 % boucle sur les annees 
     287% boucle sur les années 
    288288for a=1:28; 
    289289   clear delta; 

  • trunk/src/mode_sahelien/composite_olr_phase1_old.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 2001 a 2008 nombre de jours a enlever au debut 12327-(365*3+366+365*3+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 2001 a 2008 nombre de jours a enlever au debut 12327-(365*3+366+365*3+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 selectionne les mois d ete (attention aux années bixestiles!) 
    5353 
    5454siz1=31+28+31+30+31+1; 
    … …  
    155155compteurmax=0; 
    156156 
    157 % boucle sur les annees 
     157% boucle sur les années 
    158158for a=1:28; 
    159159   % boucle sur les jours 

  • trunk/src/mode_sahelien/composite_olr_phase_num.m

    r85 r91  
    101101 
    102102% taille 12327x73x144 1er juin 1974-31 decembre2007 
    103 % on enleve les annees 1974 a 1978 nombre de jours a enlever au debut 214+365+366+365+365 
    104 % on enleve les annees 2001 a 2008 nombre de jours a enlever au debut 12327-(365*3+366+365*3+60) 
     103% on enleve les années 1974 a 1978 nombre de jours a enlever au debut 214+365+366+365+365 
     104% on enleve les années 2001 a 2008 nombre de jours a enlever au debut 12327-(365*3+366+365*3+60) 
    105105 
    106106debut=214+365+366+365+365; 
    … …  
    124124 
    125125% on permture le temps et l'espce 
    126 % nouvelles dimentions 
     126% nouvelles dimensions 
    127127% 21 144 siz 
    128128olr1=permute(olr1, [2 3 1]); 
    … …  
    144144% on met sous la forme longitude latitude jour 
    145145 
    146 % on selectionne les mois d ete (attention aux annees bixestiles!) 
     146% on selectionne les mois d ete (attention aux années bixestiles!) 
    147147 
    148148siz1=31+28+31+30+31+1; 
    … …  
    260260compteurmax=0; 
    261261 
    262 % boucle sur les annees 
     262% boucle sur les années 
    263263for a=1:28; 
    264264   % boucle sur les jours 

  • trunk/src/mode_sahelien/composite_olr_sahel_1025.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; 
    … …  
    3131 
    3232% on permture le temps et l'espce 
    33 % nouvelles dimentions 
     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+1; 
    … …  
    164164compteurmin=0; 
    165165 
    166 % boucle sur les annees 
     166% boucle sur les années 
    167167for a=1:28; 
    168168   clear delta; 

  • trunk/src/mode_sahelien/composite_olr_sahel_pb_3j.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; 
    … …  
    3131 
    3232% on permture le temps et l'espce 
    33 % nouvelles dimentions 
     33% nouvelles dimensions 
    3434% 41 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+1; 
    … …  
    156156clear v; 
    157157 
    158 % on selectionne les annees de 1979 a 2006 
     158% on selectionne les années de 1979 a 2006 
    159159ind_olr=ind_olr(:,1:28); 
    160160k=122*28; 
    … …  
    176176compteurmax=0; 
    177177compteurmin=0; 
    178 % boucle sur les annees 
     178% boucle sur les années 
    179179for a=1:28; 
    180180   clear delta; 

  • trunk/src/mode_sahelien/composite_pluie_ird_sahel.m

    r26 r91  
    583583       axis([-180 180 -10 40]); 
    584584       if e==7; 
    585            title('Pluie IRD : diffrences composites min -max (mode Sahelien)','fontsize',8); 
     585           title('Pluie IRD : differences composites min -max (mode Sahelien)','fontsize',8); 
    586586       end; 
    587587       text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 

  • trunk/src/mode_sahelien/prepare_olr_filtre.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; 
    … …  
    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+1; 

  • trunk/src/mode_sahelien/prepare_olr_total.m

    r30 r91  
    1515 
    1616% taille 12327x73x144 1er juin 1974-31 decembre2007 
    17 % on enleve les annees 1974 a 1978 nombre de jours a enlever au debut 214+365+366+365+365 
    18 % on enleve les annees 2007 a 2008 nombre de jours a enlever au debut 12327-(365+60) 
     17% on enleve les années 1974 a 1978 nombre de jours a enlever au debut 214+365+366+365+365 
     18% on enleve les années 2007 a 2008 nombre de jours a enlever au debut 12327-(365+60) 
    1919 
    2020debut=214+365+366+365+365; 
    … …  
    5757 
    5858% on met sous la forme longitude latitude jour 
    59 % on selectionne les mois d ete (attention aux annees bixestiles!) 
     59% on selectionne les mois d ete (attention aux années bixestiles!) 
    6060 
    6161siz1=31+28+31+30+31+1; 

  • trunk/src/mode_sahelien/reconstitution.m

    r30 r91  
    4242orient('landscape') 
    4343 
    44 % On selectionne la demaine geographique 
     44% On sélectionne le domaine géographique 
    4545 
    4646taille=(12-7+1)*(17-9+1); 
    … …  
    7171cosa3=reshape(cosa3,122,22); 
    7272 
    73 % On selectionne la demaine geographique 
     73% On sélectionne le domaine géographique 
    7474olr_centre=reshape(olr_centre,2684,17,25); 
    7575cosa4=olr_centre(:,7:12,9:17); 

  • trunk/src/mode_sahelien/reconstitution_1025.m

    r30 r91  
    4747orient('landscape') 
    4848 
    49 % On selectionne la demaine geographique 
     49% On sélectionne le domaine géographique 
    5050 
    5151% la latitude : 5N (7) 12.5N (10) 

  • trunk/src/mode_sahelien/reconstitution_1030.m

    r30 r91  
    4747orient('landscape') 
    4848 
    49 % On selectionne la demaine geographique 
     49% On sélectionne le domaine géographique 
    5050 
    5151% la latitude : 7.5N (8) 12.5N (10) 
    … …  
    7575cosa24=reshape(cosa24,122,28); 
    7676 
    77 % On selectionne la demaine geographique 
     77% On sélectionne la domaine géographique 
    7878olr=reshape(olr,3416,25,17); 
    7979cosa4=olr(:,lon1:lon2,lat1:lat2); 

  • trunk/src/mode_sahelien/reconstitution_1030_std.m

    r25 r91  
    3939orient('landscape') 
    4040 
    41 % On selectionne la demaine geographique 
     41% On sélectionne la domaine géographique 
    4242 
    4343% la latitude : 7.5N (8) 12.5N (10) 

  • trunk/src/mode_sahelien/reconstitution_30100_eof1.m

    r30 r91  
    4242orient('landscape') 
    4343 
    44 % On selectionne la demaine geographique 
     44% On sélectionne le domaine géographique 
    4545 
    4646% la latitude : 7.5N (8) 12.5N (10) 
    … …  
    6161cosa1=reshape(cosa1,122,28); 
    6262 
    63 % On selectionne la demaine geographique 
     63% On sélectionne le domaine géographique 
    6464olr=reshape(olr,3416,25,17); 
    6565cosa4=olr(:,lon1:lon2,lat1:lat2); 

  • trunk/src/mode_sahelien/reconstitution_30100_eof23.m

    r30 r91  
    3939orient('landscape') 
    4040 
    41 % On selectionne la demaine geographique 
     41% On sélectionne le domaine géographique 
    4242 
    4343% la latitude : 7.5N (8) 12.5N (10) 
    … …  
    6060cosa23=reshape(cosa23,122,28); 
    6161 
    62 % On selectionne la demaine geographique 
     62% On sélectionne le domaine géographique 
    6363 
    6464cosa1=cosa1(:,lon1:lon2,lat1:lat2); 

  • trunk/src/mode_sahelien/reconstitution_ER_rossby_30100_bis.m

    r30 r91  
    3939cosa_all=reshape(olr,3416,25,17); 
    4040 
    41 % On selectionne la demaine geographique 
     41% On sélectionne le domaine géographique 
    4242 
    4343% la latitude : 7.5N (8) 12.5N (10) 

  • trunk/src/mode_sahelien/reconstitution_ER_rossby_df.m

    r86 r91  
    138138orient('landscape') 
    139139 
    140 % On selectionne la demaine geographique 
     140% On sélectionne le domaine géographique 
    141141 
    142142% la latitude : 7.5N (8) 12.5N (10) 
    … …  
    161161cosa12=reshape(cosa12,122,28); 
    162162 
    163 % On selectionne la demaine geographique 
     163% On sélectionne le domaine géographique 
    164164olr=reshape(olr,3416,25,17); 
    165165cosa4=olr(:,lon1:lon2,lat1:lat2); 

  • trunk/src/mode_sahelien/time_serie_jjas_df.m

    r88 r91  
    5353% nb of figure depends on number of years to be plot : so figure(3) contains  
    5454% something only if nb years > 23 (or so), so Postscript file of figure(3)  
    55 % shoul not be saved if  nb years is less than 23 ! 
     55% should not be saved if  nb years is less than 23 ! 
    5656% 
    5757% EVOLUTIONS 

  • trunk/src/palette.m

    r86 r91  
    3838palette2=zeros(20,3); 
    3939palette2(20,:)= [0.2 0 0];          % noir 
    40 palette2(19,:)= [0.4 0 0];          % Marron foncé 
     40palette2(19,:)= [0.4 0 0];          % Marron foncé 
    4141palette2(18,:)= [0.6 0 0];          % Marron rouge 
    4242palette2(17,:)= [0.8 0 0];          % Rouge 

  • trunk/src/palette1.m

    r86 r91  
    3838palette2=zeros(20,3); 
    3939palette2(20,:)= [0.2 0 0];          % noir 
    40 palette2(19,:)= [0.4 0 0];          % Marron foncé 
     40palette2(19,:)= [0.4 0 0];          % Marron foncé 
    4141palette2(18,:)= [0.6 0 0];          % Marron rouge 
    4242palette2(17,:)= [0.8 0 0];          % Rouge 

  • trunk/src/palette3.m

    r86 r91  
    3939palette2=zeros(16,3); 
    4040%palette2(20,:)= [0.2 0 0];          % noir 
    41 %palette2(19,:)= [0.4 0 0];          % Marron foncé 
     41%palette2(19,:)= [0.4 0 0];          % Marron foncé 
    4242%palette2(18,:)= [0.6 0 0];          % Marron rouge 
    4343%palette2(17,:)= [0.8 0 0];          % Rouge 

  • trunk/src/palettecomplet.m

    r86 r91  
    3939palette2=zeros(20,3); 
    4040palette2(20,:)= [0.2 0 0];          % noir 
    41 palette2(19,:)= [0.4 0 0];          % Marron foncé 
     41palette2(19,:)= [0.4 0 0];          % Marron foncé 
    4242palette2(18,:)= [0.6 0 0];          % Marron rouge 
    4343palette2(17,:)= [0.8 0 0];          % Rouge 

  • trunk/src/palettejerome30.m

    r86 r91  
    3838palette=zeros(20,3); 
    3939palette(20,:)= [0.3 0 0];  % noir 
    40 palette(19,:)= [0.5 0 0];  % Marron foncé 
     40palette(19,:)= [0.5 0 0];  % Marron foncé 
    4141palette(18,:)= [0.7 0 0];  % Marron rouge 
    4242palette(17,:)= [0.9 0 0];  % Rouge 

  • trunk/src/refdataget.sh

    r89 r91  
    9393system=$(uname) 
    9494case "${system}" in 
    95    AIX|IRIX64) 
    96       echo " www : no specific posix checking" 
    97    ;; 
    98    *) 
    99       set -o posix 
    100    ;; 
     95    AIX|IRIX64) 
     96        echo " www : no specific posix checking" 
     97    ;; 
     98    *) 
     99        set -o posix 
     100    ;; 
    101101esac 
    102102unset system 
    … …  
    116116if [ "${IRCAAM_LOG}" = "" ] 
    117117then 
    118    echo " eee : \${IRCAAM_LOG} not defined" 
    119    exit 1 
     118    echo " eee : \${IRCAAM_LOG} not defined" 
     119    exit 1 
    120120fi 
    121121# 
    … …  
    123123if [ ! -d ${IRCAAM_LOG} ] 
    124124then 
    125    echo " eee : ${IRCAAM_LOG} not found" 
    126    exit 1 
     125    echo " eee : ${IRCAAM_LOG} not found" 
     126    exit 1 
    127127fi 
    128128# 
    … …  
    130130if [ ! -x ${IRCAAM_LOG} ] 
    131131then 
    132    echo " eee : ${IRCAAM_LOG} not reachable" 
    133    exit 1 
     132    echo " eee : ${IRCAAM_LOG} not reachable" 
     133    exit 1 
    134134fi 
    135135# 
    … …  
    137137if [ ! -w ${IRCAAM_LOG} ] 
    138138then 
    139    echo " eee : ${IRCAAM_LOG} not writable" 
    140    exit 1 
     139    echo " eee : ${IRCAAM_LOG} not writable" 
     140    exit 1 
    141141fi 
    142142# 
    … …  
    150150if [ ! -w ${IRCAAM_ID} ] 
    151151then 
    152    echo " eee : ${IRCAAM_ID} not writable" 
    153    exit 1 
     152    echo " eee : ${IRCAAM_ID} not writable" 
     153    exit 1 
    154154fi 
    155155 
    … …  
    165165for simu in ${list_simu} 
    166166do 
    167    case ${simu} in 
    168       CtIV|CtCl|AfNQIVIV) 
    169          dirref="${dirrefb}/Sim/" 
    170       ;; 
    171       AsNQIVIV) 
    172          dirref="${dirrefb}/Sim_As/" 
    173       ;; 
    174       TrNQIVIV) 
    175          dirref="${dirrefb}/Sim_Tr/" 
    176       ;; 
    177       *) 
    178          echo "eee : wrong \${simu}" 
    179          exit 1 
    180       ;; 
    181    esac 
    182    dirref="${dirref}/${simu}/" 
    183    for param in ${list_param} 
    184    do 
     167    case ${simu} in 
     168        CtIV|CtCl|AfNQIVIV) 
     169            dirref="${dirrefb}/Sim/" 
     170        ;; 
     171        AsNQIVIV) 
     172            dirref="${dirrefb}/Sim_As/" 
     173        ;; 
     174        TrNQIVIV) 
     175            dirref="${dirrefb}/Sim_Tr/" 
     176        ;; 
     177        *) 
     178            echo "eee : wrong \${simu}" 
     179            exit 1 
     180        ;; 
     181    esac 
     182    dirref="${dirref}/${simu}/" 
     183    for param in ${list_param} 
     184    do 
    185185      year_i=${year_min} 
    186186      while [ ${year_i} -le ${year_max} ] 
    187187      do 
    188          year_c=$(printf "%4.4d" ${year_i}') 
    189          fileref=${param}.${simu}.${year_c}06-09.nc 
    190          if  [ -f ${IRCAAM_ID}/${fileref} ] 
    191          then 
    192             echo "iii : ${IRCAAM_ID}/${fileref} exist" ${log} 2>&1 
    193             echo "iii : nothing done" >> ${log} 2>&1 
    194          else 
    195             wget --no-verbose --user=ircaam-dods --password=${dods_passwd} \ 
    196                -P ${IRCAAM_ID} \ 
    197                ${dirref}/${fileref} >> ${log} 2>&1 
    198             # ++ ok ? 
    199             # ++ si climserv faire ln -s ${dirref}/${fileref} ${IRCAAM_ID}/ 
    200          fi 
    201          year_i=$(( ${year_i} + 1 )) 
     188          year_c=$(printf "%4.4d" ${year_i}') 
     189          fileref=${param}.${simu}.${year_c}06-09.nc 
     190          if  [ -f ${IRCAAM_ID}/${fileref} ] 
     191          then 
     192              echo "iii : ${IRCAAM_ID}/${fileref} exist" ${log} 2>&1 
     193              echo "iii : nothing done" >> ${log} 2>&1 
     194          else 
     195              wget --no-verbose --user=ircaam-dods --password=${dods_passwd} \ 
     196              -P ${IRCAAM_ID} \ 
     197              ${dirref}/${fileref} >> ${log} 2>&1 
     198              # ++ ok ? 
     199              # ++ si climserv faire ln -s ${dirref}/${fileref} ${IRCAAM_ID}/ 
     200          fi 
     201          year_i=$(( ${year_i} + 1 )) 
    202202      done 
    203203      unset year_i 
    204    done 
    205    unset param 
     204    done 
     205    unset param 
    206206done 
    207207# OLR 
    … …  
    210210if  [ -f ${IRCAAM_ID}/${fileref} ] 
    211211then 
    212    echo "iii : ${IRCAAM_ID}/${fileref} exist"  >> ${log} 2>&1 
    213    echo "iii : nothing done" >> ${log} 2>&1 
     212    echo "iii : ${IRCAAM_ID}/${fileref} exist"  >> ${log} 2>&1 
     213    echo "iii : nothing done" >> ${log} 2>&1 
    214214else 
    215    wget --no-verbose -P ${IRCAAM_ID} ${dirref}/${fileref} >> ${log} 
     215    wget --no-verbose -P ${IRCAAM_ID} ${dirref}/${fileref} >> ${log} 
    216216   # ++ ok ? 
    217217fi 

  • trunk/src/sauvegrads.m

    r86 r91  
    7676h=[]; 
    7777if nargin < 3 
    78         error('Il n''y a pas assez d''arguments.'); 
     78    error('Il n''y a pas assez d''arguments.'); 
    7979elseif nargin > 3 
    80         error('Il y a trop d''arguments.') 
     80    error('Il y a trop d''arguments.') 
    8181end 
    8282 
Note: See TracChangeset for help on using the changeset viewer.