Index: trunk/LMDZ.TITAN/libf/chimtitan/chimie.c
===================================================================
--- trunk/LMDZ.TITAN/libf/chimtitan/chimie.c	(revision 1968)
+++ trunk/LMDZ.TITAN/libf/chimtitan/chimie.c	(revision 2099)
@@ -13,14 +13,7 @@
    double  ai,ao,ei,eo,ki,ko,m,ti,to;
    FILE   *out;
-   char   corps[100][10];
    static char   reaction[NREAC+1][12][10]={
 #include VERCHIM
    "",       "",       "",       "",     "",  "","","","","","",""};
-
-   for( i = 0; i <= NC; i++)
-   {
-     strcpy( corps[i], CORPS[i] );
-     corps[i][strcspn(CORPS[i], " ")] = '\0';
-   }
 
    for( i = 0; i <= NC-1; i++ )
@@ -40,5 +33,5 @@
           && (strcmp(reaction[i][j],""))     ) /* Which compound ? */
          {
-            while( strcmp(reaction[i][j],corps[k]) )     /* Compound in reaction j, column i */
+            while( strcmp(reaction[i][j],CORPS[k]) )     /* Compound in reaction j, column i */
             {
                if( k == NC+1 )
Index: trunk/LMDZ.TITAN/libf/chimtitan/comp.c
===================================================================
--- trunk/LMDZ.TITAN/libf/chimtitan/comp.c	(revision 1968)
+++ trunk/LMDZ.TITAN/libf/chimtitan/comp.c	(revision 2099)
@@ -8,5 +8,4 @@
 {
    int   i,j;
-   char  corps[NC+1][10];
 
    double m,ma,epsa,sig,siga,p;
@@ -19,10 +18,4 @@
    siga = 3.798e0;         /* Lennard-Jones length of background gas 1/10 nm */
    epsa = 71.4e0;       /* Lennard-Jones energy of background gas */
-
-   for( i = 0; i <= NC; i++)
-   {
-     strcpy( corps[i], CORPS[i] );
-     corps[i][strcspn(CORPS[i], " ")] = '\0';
-   }
 
    for( i = 0; i <= NC-1; i++)
@@ -33,5 +26,5 @@
    for( i = 0; i <= NC-1; i++ )
    {
-      if( strcmp(corps[i], "CH4") == 0 )
+      if( strcmp(CORPS[i], "CH4") == 0 )
       {
          MASS[i] = 16.04e0;
@@ -42,9 +35,9 @@
            / ( CT[j] * sig * omega(TEMP[j],epsa,148.6e0) );
       }
-      if( strcmp(corps[i], "H") == 0 )
+      if( strcmp(CORPS[i], "H") == 0 )
       {
          MASS[i] = 1.01e0;
       }
-      if( strcmp(corps[i], "H2") == 0 )
+      if( strcmp(CORPS[i], "H2") == 0 )
       {
          MASS[i] = 2.0158e0;
@@ -55,5 +48,5 @@
            / ( CT[j] * sig * omega(TEMP[j],epsa,59.7e0) );
       }
-      if( strcmp(corps[i], "CH") == 0 )
+      if( strcmp(CORPS[i], "CH") == 0 )
       {
          MASS[i] = 13.02e0;
@@ -63,5 +56,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( ( strcmp( corps[i], "CH2" ) == 0 ) || ( strcmp( corps[i], "CH2s" ) == 0 ) )
+      if( ( strcmp( CORPS[i], "CH2" ) == 0 ) || ( strcmp( CORPS[i], "CH2s" ) == 0 ) )
       {
          MASS[i] = 14.03e0;
@@ -71,5 +64,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "CH3") == 0 )
+      if( strcmp(CORPS[i], "CH3") == 0 )
       {
          MASS[i] = 15.03e0;
@@ -79,5 +72,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "C") == 0 )
+      if( strcmp(CORPS[i], "C") == 0 )
       {
          MASS[i] = 12.01e0;
@@ -87,5 +80,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "C2") == 0 )
+      if( strcmp(CORPS[i], "C2") == 0 )
       {
          MASS[i] = 24.02e0;
@@ -95,5 +88,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "C2H") == 0 )
+      if( strcmp(CORPS[i], "C2H") == 0 )
       {
          MASS[i] = 25.03e0;
@@ -103,5 +96,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "C2H3") == 0 )
+      if( strcmp(CORPS[i], "C2H3") == 0 )
       {
          MASS[i] = 27.05e0;
@@ -111,5 +104,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "C2H4") == 0 )
+      if( strcmp(CORPS[i], "C2H4") == 0 )
       {
          MASS[i] = 28.05e0;
@@ -120,5 +113,5 @@
            / ( CT[j] * sig * omega(TEMP[j],epsa,224.7e0) );
       }
-      if( strcmp(corps[i], "C2H2") == 0 )
+      if( strcmp(CORPS[i], "C2H2") == 0 )
       {
          MASS[i] = 26.04e0;
@@ -129,5 +122,5 @@
            / ( CT[j] * sig * omega(TEMP[j],epsa,231.8e0) );
       }
-      if( strcmp(corps[i], "C2H5") == 0 )
+      if( strcmp(CORPS[i], "C2H5") == 0 )
       {
          MASS[i] = 29.06e0;
@@ -137,5 +130,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "C2H6") == 0 )
+      if( strcmp(CORPS[i], "C2H6") == 0 )
       {
          MASS[i] = 30.07e0;
@@ -146,5 +139,5 @@
            / ( CT[j] * sig * omega(TEMP[j],epsa,215.7e0) );
       }
-      if( strcmp(corps[i], "C3H2") == 0 )
+      if( strcmp(CORPS[i], "C3H2") == 0 )
       {
          MASS[i] = 38.05e0;
@@ -154,5 +147,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "C3H3") == 0 )
+      if( strcmp(CORPS[i], "C3H3") == 0 )
       {
          MASS[i] = 39.06e0;
@@ -162,5 +155,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( ( strcmp(corps[i], "CH2CCH2") == 0 ) || ( strcmp(corps[i], "CH3CCH") == 0 ) )
+      if( ( strcmp(CORPS[i], "CH2CCH2") == 0 ) || ( strcmp(CORPS[i], "CH3CCH") == 0 ) )
       {
          MASS[i] = 40.07e0;
@@ -171,5 +164,5 @@
            / ( CT[j] * sig * omega(TEMP[j],epsa,251.8e0) );
       }
-      if( strcmp(corps[i], "C3H5") == 0 )
+      if( strcmp(CORPS[i], "C3H5") == 0 )
       {
          MASS[i] = 41.07e0;
@@ -179,5 +172,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "C3H6") == 0 )
+      if( strcmp(CORPS[i], "C3H6") == 0 )
       {
          MASS[i] = 42.08e0;
@@ -188,5 +181,5 @@
            / ( CT[j] * sig * omega(TEMP[j],epsa,248.9e0) );
       }
-      if( strcmp(corps[i], "C3H7") == 0 )
+      if( strcmp(CORPS[i], "C3H7") == 0 )
       {
          MASS[i] = 43.09e0;
@@ -196,5 +189,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
        }
-      if( strcmp(corps[i], "C3H8") == 0 )
+      if( strcmp(CORPS[i], "C3H8") == 0 )
       {
          MASS[i] = 44.11e0;
@@ -205,5 +198,5 @@
            / ( CT[j] * sig * omega(TEMP[j],epsa,237.1e0) );
       }
-      if( strcmp(corps[i], "C4H") == 0 )
+      if( strcmp(CORPS[i], "C4H") == 0 )
       {
          MASS[i] = 49.05e0;
@@ -213,5 +206,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( ( strcmp(corps[i], "C4H2") == 0 )||( strcmp(corps[i], "C4H2s") == 0 ) )
+      if( ( strcmp(CORPS[i], "C4H2") == 0 )||( strcmp(CORPS[i], "C4H2s") == 0 ) )
       {
          MASS[i] = 50.06e0;
@@ -221,5 +214,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "C4H3") == 0 )
+      if( strcmp(CORPS[i], "C4H3") == 0 )
       {
          MASS[i] = 51.07e0;
@@ -229,5 +222,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "C4H4") == 0 )
+      if( strcmp(CORPS[i], "C4H4") == 0 )
       {
          MASS[i] = 52.08e0;
@@ -237,5 +230,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "C4H5") == 0 )
+      if( strcmp(CORPS[i], "C4H5") == 0 )
       {
          MASS[i] = 53.07e0;
@@ -245,5 +238,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "C4H6") == 0 )
+      if( strcmp(CORPS[i], "C4H6") == 0 )
       {
          MASS[i] = 54.09e0;
@@ -253,5 +246,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "C4H10") == 0 )
+      if( strcmp(CORPS[i], "C4H10") == 0 )
       {
          MASS[i] = 58.13e0;
@@ -262,5 +255,5 @@
            / ( CT[j] * sig * omega(TEMP[j],epsa,531.4e0) );
       }
-      if( strcmp(corps[i], "C6H") == 0 )
+      if( strcmp(CORPS[i], "C6H") == 0 )
       {
          MASS[i] = 73.07e0;
@@ -270,5 +263,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "C6H2") == 0 )
+      if( strcmp(CORPS[i], "C6H2") == 0 )
       {
          MASS[i] = 74.08e0;
@@ -278,5 +271,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "C8H2") == 0 )
+      if( strcmp(CORPS[i], "C8H2") == 0 )
       {
          MASS[i] = 98.10e0;
@@ -286,5 +279,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp( corps[i], "AC6H6" ) == 0 )
+      if( strcmp( CORPS[i], "AC6H6" ) == 0 )
       {
          MASS[i] = 78.1136e0;
@@ -294,5 +287,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( ( strcmp( corps[i], "C6H5" ) == 0 ) || ( strcmp( corps[i], "AC6H5" ) == 0 ) )
+      if( ( strcmp( CORPS[i], "C6H5" ) == 0 ) || ( strcmp( CORPS[i], "AC6H5" ) == 0 ) )
       {
          MASS[i] = 77.1136e0;
@@ -302,5 +295,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp( corps[i], "C6H6" ) == 0 )
+      if( strcmp( CORPS[i], "C6H6" ) == 0 )
       {
          MASS[i] = 78.1136e0;
@@ -310,9 +303,9 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "N2") == 0 )
+      if( strcmp(CORPS[i], "N2") == 0 )
       {
          MASS[i] = 28.0134e0;
       }
-      if( strcmp(corps[i], "N4S") == 0 )
+      if( strcmp(CORPS[i], "N4S") == 0 )
       {
          MASS[i] = 14.01e0;
@@ -322,5 +315,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "NH") == 0 )
+      if( strcmp(CORPS[i], "NH") == 0 )
       {
          MASS[i] = 15.01e0;
@@ -330,5 +323,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "CN") == 0 )
+      if( strcmp(CORPS[i], "CN") == 0 )
       {
          MASS[i] = 26.02e0;
@@ -338,5 +331,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "HCN") == 0 )
+      if( strcmp(CORPS[i], "HCN") == 0 )
       {
          MASS[i] = 27.04e0;
@@ -347,5 +340,5 @@
            / ( CT[j] * sig * omega(TEMP[j],epsa,569.1e0) );
       }
-      if( strcmp(corps[i], "H2CN") == 0 )
+      if( strcmp(CORPS[i], "H2CN") == 0 )
       {
          MASS[i] = 28.05e0;
@@ -355,5 +348,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "C2N") == 0 )         /* C2N */
+      if( strcmp(CORPS[i], "C2N") == 0 )         /* C2N */
       {
          MASS[i] = 39.05e0;
@@ -363,5 +356,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp( corps[i], "CHCN" ) == 0 )
+      if( strcmp( CORPS[i], "CHCN" ) == 0 )
       {
          MASS[i]   = 39.05e0;
@@ -371,5 +364,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp( corps[i], "CH2CN" ) == 0 )
+      if( strcmp( CORPS[i], "CH2CN" ) == 0 )
       {
          MASS[i]   = 40.04e0;
@@ -379,5 +372,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp( corps[i], "CH3CN" ) == 0 )
+      if( strcmp( CORPS[i], "CH3CN" ) == 0 )
       {
          MASS[i]   = 41.05e0;
@@ -387,5 +380,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp( corps[i], "C2H3CN" ) == 0 )
+      if( strcmp( CORPS[i], "C2H3CN" ) == 0 )
       {
          MASS[i]   = 53.06e0;
@@ -395,5 +388,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "NCCN") == 0 )        /* NCCN */
+      if( strcmp(CORPS[i], "NCCN") == 0 )        /* NCCN */
       {
          MASS[i] = 52.04e0;
@@ -404,5 +397,5 @@
            / ( CT[j] * sig * omega(TEMP[j],epsa,348.6e0) );
       }
-      if( strcmp(corps[i], "C3N") == 0 )         /* C3N */
+      if( strcmp(CORPS[i], "C3N") == 0 )         /* C3N */
       {
          MASS[i] = 50.04e0;
@@ -412,5 +405,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "HC3N") == 0 )        /* HC3N */
+      if( strcmp(CORPS[i], "HC3N") == 0 )        /* HC3N */
       {
          MASS[i] = 51.05e0;
@@ -420,5 +413,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp( corps[i], "C4N2" ) == 0 )
+      if( strcmp( CORPS[i], "C4N2" ) == 0 )
       {
          MASS[i]   = 76.1e0;
@@ -428,5 +421,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "H2O") == 0 )       
+      if( strcmp(CORPS[i], "H2O") == 0 )       
       {
          MASS[i] = 18.02e0;
@@ -437,5 +430,5 @@
            / ( CT[j] * sig * omega(TEMP[j],epsa,809.1e0) );
       }
-      if( ( strcmp(corps[i], "O3P") == 0 ) || ( strcmp(corps[i], "O1D") == 0 ) )      
+      if( ( strcmp(CORPS[i], "O3P") == 0 ) || ( strcmp(CORPS[i], "O1D") == 0 ) )      
       {
          MASS[i] = 16.0e0;
@@ -445,5 +438,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "OH") == 0 )        
+      if( strcmp(CORPS[i], "OH") == 0 )        
       {
          MASS[i] = 17.01e0;
@@ -453,5 +446,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "HO2") == 0 )
+      if( strcmp(CORPS[i], "HO2") == 0 )
       {
          MASS[i] = 33.01e0;
@@ -461,5 +454,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "H2O2") == 0 )
+      if( strcmp(CORPS[i], "H2O2") == 0 )
       {
          MASS[i] = 33.01e0;
@@ -469,5 +462,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "O2") == 0 )
+      if( strcmp(CORPS[i], "O2") == 0 )
       {
          MASS[i] = 32.0e0;
@@ -477,5 +470,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "O3") == 0 )
+      if( strcmp(CORPS[i], "O3") == 0 )
       {
          MASS[i] = 32.0e0;
@@ -485,5 +478,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "CO") == 0 )          
+      if( strcmp(CORPS[i], "CO") == 0 )          
       {
          MASS[i] = 28.01e0;
@@ -494,5 +487,5 @@
            / ( CT[j] * sig * omega(TEMP[j],epsa,91.7e0) );
       }
-      if( strcmp(corps[i], "HCO") == 0 )        
+      if( strcmp(CORPS[i], "HCO") == 0 )        
       {
          MASS[i] = 29.02e0;
@@ -502,5 +495,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "CO2") == 0 )       
+      if( strcmp(CORPS[i], "CO2") == 0 )       
       {
          MASS[i] = 44.01e0;
@@ -511,5 +504,5 @@
            / ( CT[j] * sig * omega(TEMP[j],epsa,195.2e0) );
       }
-      if( strcmp(corps[i], "CH2CO") == 0 )    
+      if( strcmp(CORPS[i], "CH2CO") == 0 )    
       {
          MASS[i] = 42.04e0;
@@ -519,5 +512,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "CH2O") == 0 )    
+      if( strcmp(CORPS[i], "CH2O") == 0 )    
       {
          MASS[i] = 30.03e0;
@@ -527,5 +520,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( ( strcmp(corps[i], "CH2OH") == 0 ) || ( strcmp(corps[i], "CH3O") == 0 ) )  
+      if( ( strcmp(CORPS[i], "CH2OH") == 0 ) || ( strcmp(CORPS[i], "CH3O") == 0 ) )  
       {
          MASS[i] = 31.04e0;
@@ -535,5 +528,5 @@
             MD[i][j] = sqrt( p * TEMP[j] * m ) / ( CT[j] * sig );
       }
-      if( strcmp(corps[i], "CH3OH") == 0 )       
+      if( strcmp(CORPS[i], "CH3OH") == 0 )       
       {
          MASS[i] = 32.042e0;
Index: trunk/LMDZ.TITAN/libf/chimtitan/gptitan.c
===================================================================
--- trunk/LMDZ.TITAN/libf/chimtitan/gptitan.c	(revision 1968)
+++ trunk/LMDZ.TITAN/libf/chimtitan/gptitan.c	(revision 2099)
@@ -19,5 +19,5 @@
 	      int *htoh2, double *surfhaze)
 {
-  char   outlog[100],corps[100][10];
+  char   outlog[100];
   int    i,j,k,l;
   int    ireac,ncom1,ncom2;
@@ -58,10 +58,4 @@
 #endif
   /**/
-
-  for( i = 0; i <= NC; i++)
-    {
-      strcpy( corps[i], CORPS[i] );
-      corps[i][strcspn(CORPS[i], " ")] = '\0';
-    }
    
   strcpy( outlog, "chimietitan" );
@@ -105,8 +99,8 @@
   for( i = 0; i <= ST-1; i++ ) 
     {
-      if( strcmp(corps[i], "H"  ) == 0 ) ih   = i;
-      if( strcmp(corps[i], "H2" ) == 0 ) ih2  = i;
-      if( strcmp(corps[i], "N2" ) == 0 ) in2  = i;
-      if( strcmp(corps[i], "N4S") == 0 ) in4s = i;
+      if( strcmp(CORPS[i], "H"  ) == 0 ) ih   = i;
+      if( strcmp(CORPS[i], "H2" ) == 0 ) ih2  = i;
+      if( strcmp(CORPS[i], "N2" ) == 0 ) in2  = i;
+      if( strcmp(CORPS[i], "N4S") == 0 ) in4s = i;
     }
 
@@ -160,5 +154,5 @@
 	 fprintf(out,"%d %e %e %e\n",j,(RA[j]-R0),NB[j],TEMP[j]);
 	 for (i=0;i<=NREAC-1;i++) fprintf(out,"%d %e\n",i,KRATE[i][j]);
-	 for (i=0;i<=ST-1;i++) fprintf(out,"%10s %e\n",corps[i],Y[i][j]);
+	 for (i=0;i<=ST-1;i++) fprintf(out,"%10s %e\n",CORPS[i],Y[i][j]);
 	 fclose( out );
 	 }
@@ -231,5 +225,5 @@
 	  if( (*aerprod) == 1 )
 	    {
-	      aer(corps,TEMP,NB,Y,&j,k_dep,faer,
+	      aer(CORPS,TEMP,NB,Y,&j,k_dep,faer,
 		  &dyc2h2,&dyhc3n,&dyhcn,&dynccn,&dych3cn,&dyc2h3cn,utilaer,
 		  mmolaer,productaer,csurn,csurh);
@@ -245,20 +239,20 @@
 		 printf("AERPROD : LAT = %g - J = %d\n",(*LAT),j);
 		 if(fabs(dyc2h2*NB[j])>fabs(fp[utilaer[2]]/10.))
-		 printf("fp(%s) =%e; dyc2h2 =%e\n",corps[utilaer[2]],
+		 printf("fp(%s) =%e; dyc2h2 =%e\n",CORPS[utilaer[2]],
 		 fp[utilaer[2]],dyc2h2*NB[j]);
 		 if(fabs(dyhcn*NB[j])>fabs(fp[utilaer[5]]/10.))
-		 printf("fp(%s) =%e; dyhcn  =%e\n",corps[utilaer[5]],
+		 printf("fp(%s) =%e; dyhcn  =%e\n",CORPS[utilaer[5]],
 		 fp[utilaer[5]],dyhcn*NB[j]);
 		 if(fabs(dyhc3n*NB[j])>fabs(fp[utilaer[6]]/10.))
-		 printf("fp(%s) =%e; dyhc3n =%e\n",corps[utilaer[6]],
+		 printf("fp(%s) =%e; dyhc3n =%e\n",CORPS[utilaer[6]],
 		 fp[utilaer[6]],dyhc3n*NB[j]);
 		 if(fabs(dynccn*NB[j])>fabs(fp[utilaer[13]]/10.))
-		 printf("fp(%s) =%e; dynccn =%e\n",corps[utilaer[13]],
+		 printf("fp(%s) =%e; dynccn =%e\n",CORPS[utilaer[13]],
 		 fp[utilaer[13]],dynccn*NB[j]);
 		 if(fabs(dych3cn*NB[j])>fabs(fp[utilaer[14]]/10.))
-		 printf("fp(%s) =%e; dych3cn=%e\n",corps[utilaer[14]],
+		 printf("fp(%s) =%e; dych3cn=%e\n",CORPS[utilaer[14]],
 		 fp[utilaer[14]],dych3cn*NB[j]);
 		 if(fabs(dyc2h3cn*NB[j])>fabs(fp[utilaer[15]]/10.))
-		 printf("fp(%s) =%e; dyc2h3cn=%e\n",corps[utilaer[15]],
+		 printf("fp(%s) =%e; dyc2h3cn=%e\n",CORPS[utilaer[15]],
 		 fp[utilaer[15]],dyc2h3cn*NB[j]);
 	      */
@@ -289,5 +283,5 @@
 	  if( (*htoh2) == 1 )
 	    {
-              heterohtoh2(corps,TEMP,NB,Y,surfhaze,&j,&dyh,&dyh2,utilaer);
+              heterohtoh2(CORPS,TEMP,NB,Y,surfhaze,&j,&dyh,&dyh2,utilaer);
 	      /* dyh <= 0 / 1.0 en adsor., 1 en reac. */
 
@@ -295,7 +289,7 @@
 		 printf("HTOH2 : LAT = %g - J = %d\n",(*LAT),j);
 		 if(fabs(dyh*NB[j])>fabs(fp[utilaer[0]]/10.))
-		 printf("fp(%s) = %e; dyh  = %e\n",corps[utilaer[0]],fp[utilaer[0]],dyh*NB[j]);
+		 printf("fp(%s) = %e; dyh  = %e\n",CORPS[utilaer[0]],fp[utilaer[0]],dyh*NB[j]);
 		 if(fabs(dyh2*NB[j])>fabs(fp[utilaer[1]]/10.))
-		 printf("fp(%s) = %e; dyh2 = %e\n",corps[utilaer[1]],fp[utilaer[1]],dyh2*NB[j]);
+		 printf("fp(%s) = %e; dyh2 = %e\n",CORPS[utilaer[1]],fp[utilaer[1]],dyh2*NB[j]);
 	      */
 
@@ -527,5 +521,5 @@
 		    out = fopen( outlog, "a" );
 		    fprintf( out, "Latitude %g;", (*LAT));
-		    fprintf(out, " alt:%e; %s %e %e ; %e %e\n",(RA[j]-R0),corps[i],ym1[i],Y[i][j],time,delta);
+		    fprintf(out, " alt:%e; %s %e %e ; %e %e\n",(RA[j]-R0),CORPS[i],ym1[i],Y[i][j],time,delta);
 		    fclose( out );
 		    }
@@ -546,5 +540,5 @@
                   out = fopen( outlog, "a" );
                   fprintf( out, "WARNING %s mixing ratio is %e %e at %d\n",
-                           corps[i], ym1[i], Y[i][j], j );
+                           CORPS[i], ym1[i], Y[i][j], j );
                   for( k = 0; k <= NLEV-1; k++ ) fprintf( out, "%d %e %e\n",k,ym1[i],Y[i][k] );
                   fclose( out );
@@ -669,5 +663,5 @@
 	    if( (*aerprod) == 1 )
 	      {
-		aer(corps,TEMP,NB,Y,&j,k_dep,faer,
+		aer(CORPS,TEMP,NB,Y,&j,k_dep,faer,
 		    &dyc2h2,&dyhc3n,&dyhcn,&dynccn,&dych3cn,&dyc2h3cn,utilaer,
 		    mmolaer,productaer,csurn,csurh);
@@ -683,20 +677,20 @@
 		   printf("AERPROD : LAT = %g - J = %d\n",(*LAT),j);
 		   if(fabs(dyc2h2*NB[j])>fabs(fp[utilaer[2]]/10.))
-		   printf("fp(%s) =%e; dyc2h2 =%e\n",corps[utilaer[2]],
+		   printf("fp(%s) =%e; dyc2h2 =%e\n",CORPS[utilaer[2]],
 		   fp[utilaer[2]],dyc2h2*NB[j]);
 		   if(fabs(dyhcn*NB[j])>fabs(fp[utilaer[5]]/10.))
-		   printf("fp(%s) =%e; dyhcn  =%e\n",corps[utilaer[5]],
+		   printf("fp(%s) =%e; dyhcn  =%e\n",CORPS[utilaer[5]],
 		   fp[utilaer[5]],dyhcn*NB[j]);
 		   if(fabs(dyhc3n*NB[j])>fabs(fp[utilaer[6]]/10.))
-		   printf("fp(%s) =%e; dyhc3n =%e\n",corps[utilaer[6]],
+		   printf("fp(%s) =%e; dyhc3n =%e\n",CORPS[utilaer[6]],
 		   fp[utilaer[6]],dyhc3n*NB[j]);
 		   if(fabs(dynccn*NB[j])>fabs(fp[utilaer[13]]/10.))
-		   printf("fp(%s) =%e; dynccn =%e\n",corps[utilaer[13]],
+		   printf("fp(%s) =%e; dynccn =%e\n",CORPS[utilaer[13]],
 		   fp[utilaer[13]],dynccn*NB[j]);
 		   if(fabs(dych3cn*NB[j])>fabs(fp[utilaer[14]]/10.))
-		   printf("fp(%s) =%e; dych3cn=%e\n",corps[utilaer[14]],
+		   printf("fp(%s) =%e; dych3cn=%e\n",CORPS[utilaer[14]],
 		   fp[utilaer[14]],dych3cn*NB[j]);
 		   if(fabs(dyc2h3cn*NB[j])>fabs(fp[utilaer[15]]/10.))
-		   printf("fp(%s) =%e; dyc2h3cn=%e\n",corps[utilaer[15]],
+		   printf("fp(%s) =%e; dyc2h3cn=%e\n",CORPS[utilaer[15]],
 		   fp[utilaer[15]],dyc2h3cn*NB[j]);
 		*/
@@ -727,5 +721,5 @@
 	    if( (*htoh2) == 1 )
 	      {
-		heterohtoh2(corps,TEMP,NB,Y,surfhaze,&j,&dyh,&dyh2,utilaer);
+		heterohtoh2(CORPS,TEMP,NB,Y,surfhaze,&j,&dyh,&dyh2,utilaer);
 		/* dyh <= 0 / 1.0 en adsor., 1 en reac. */
 
@@ -733,7 +727,7 @@
 		   printf("HTOH2 : LAT = %g - J = %d\n",(*LAT),j);
 		   if(fabs(dyh*NB[j])>fabs(fp[utilaer[0]]/10.))
-		   printf("fp(%s) = %e; dyh  = %e\n",corps[utilaer[0]],fp[utilaer[0]],dyh*NB[j]);
+		   printf("fp(%s) = %e; dyh  = %e\n",CORPS[utilaer[0]],fp[utilaer[0]],dyh*NB[j]);
 		   if(fabs(dyh2*NB[j])>fabs(fp[utilaer[1]]/10.))
-		   printf("fp(%s) = %e; dyh2 = %e\n",corps[utilaer[1]],fp[utilaer[1]],dyh2*NB[j]);
+		   printf("fp(%s) = %e; dyh2 = %e\n",CORPS[utilaer[1]],fp[utilaer[1]],dyh2*NB[j]);
 		*/
 
@@ -817,5 +811,5 @@
 			  out = fopen( outlog, "a" );
 			  fprintf( out, "Latitude: %g; declin:%e;", (*LAT), (*DECLIN) );
-			  fprintf(out, " alt:%e; %s %e %e ; %e %e\n",(RA[j]-R0),corps[i],ym1[i],Y[i][j],time,delta);
+			  fprintf(out, " alt:%e; %s %e %e ; %e %e\n",(RA[j]-R0),CORPS[i],ym1[i],Y[i][j],time,delta);
 			  fclose( out );
 			  }
@@ -836,5 +830,5 @@
 		      out = fopen( outlog, "a" );
 		      fprintf( out, "WARNING %s mixing ratio is %e %e at %d\n",
-			       corps[i], ym1[i][j], Y[i][j], j );
+			       CORPS[i], ym1[i][j], Y[i][j], j );
 		      for( k = 0; k <= NLEV-1; k++ ) fprintf( out, "%d %e %e\n",k,ym1[i][j],Y[i][k] );
 		      fclose( out );
Index: trunk/LMDZ.TITAN/libf/chimtitan/tractitan.c
===================================================================
--- trunk/LMDZ.TITAN/libf/chimtitan/tractitan.c	(revision 1968)
+++ trunk/LMDZ.TITAN/libf/chimtitan/tractitan.c	(revision 2099)
@@ -7,5 +7,5 @@
      double Y0[][NLEV], double *FIN )
 {
-   char   outlog[100],corps[100][10];
+   char   outlog[100];
    int    i,j,k,l;
    double  annee,**tau,**ym1;
@@ -15,10 +15,4 @@
    FILE   *out;
 
-   for( i = 0; i <= NC; i++)
-   {
-     strcpy( corps[i], CORPS[i] );
-     corps[i][strcspn(CORPS[i], " ")] = '\0';
-   }
-   
    time     = ts = 0.0e0;
    annee    = 9.46728e8;
@@ -54,9 +48,9 @@
         tau[i][j] = 1.e6;   /* autres corps = 1.e6 s, donc rappel tres fort */
         
-        if( strcmp(corps[i],"C2H2") == 0 ) 
+        if( strcmp(CORPS[i],"C2H2") == 0 ) 
            tau[i][j] = annee*pow( 10., 2.+1.*(100.-(RB[j]-R0))/200. );
-        if( strcmp(corps[i],"C2H6") == 0 ) 
+        if( strcmp(CORPS[i],"C2H6") == 0 ) 
            tau[i][j] = annee*pow( 10., 1.+1.*(200.-(RB[j]-R0))/300. );  
-        if( strcmp(corps[i],"HCN") == 0 )
+        if( strcmp(CORPS[i],"HCN") == 0 )
         { 
           if( (RB[j]-R0) >= 350. ) 
@@ -65,5 +59,5 @@
              tau[i][j] = annee*10.;  
         }
-        if( strcmp(corps[i],"C4H2") == 0 )
+        if( strcmp(CORPS[i],"C4H2") == 0 )
         {
           if( (RB[j]-R0) >= 300. ) 
@@ -113,5 +107,5 @@
                   if( conv > 0.1 )
                   {
-                     fprintf(out, "%d %s %e %e\n",j,corps[i],ym1[i][j],Y[i][j]); 
+                     fprintf(out, "%d %s %e %e\n",j,CORPS[i],ym1[i][j],Y[i][j]); 
                   }
 */
@@ -131,5 +125,5 @@
 /*
                   fprintf( out, "WARNING %s mixing ratio is %e %e at %d",
-                           corps[i], ym1[i][j], Y[i][j], j );
+                           CORPS[i], ym1[i][j], Y[i][j], j );
                   fclose( out );  
 */
