source: trunk/LMDZ.TITAN/libf/chimtitan/comp.c @ 1058

/* comp: Compounds characteristics. */
/* GCCM */

#include "titan.h"

void comp_(char CORPS[][10], double *MASS)
{
   int   i;
   char  corps[100][10];

   for( i = 0; i <= NC; i++)
   {
     strcpy( corps[i], CORPS[i] );
     corps[i][strcspn(CORPS[i], " ")] = '\0';
   }

   for( i = 0; i <= NC-1; i++ )
   {
      if( strcmp(corps[i], "CH4") == 0 )
      {
         MASS[i] = 16.04e0;
      }
      if( strcmp(corps[i], "H") == 0 )
      {
         MASS[i] = 1.01e0;
      }
      if( strcmp(corps[i], "H2") == 0 )
      {
         MASS[i] = 2.0158e0;
      }
      if( strcmp(corps[i], "CH") == 0 )
      {
         MASS[i] = 13.02e0;
      }
      if( ( strcmp( corps[i], "CH2" ) == 0 ) || ( strcmp( corps[i], "CH2s" ) == 0 ) )
      {
         MASS[i] = 14.03e0;
      }
      if( strcmp(corps[i], "CH3") == 0 )
      {
         MASS[i] = 15.03e0;
      }
      if( strcmp(corps[i], "C") == 0 )
      {
         MASS[i] = 12.01e0;
      }
      if( strcmp(corps[i], "C2") == 0 )
      {
         MASS[i] = 24.02e0;
      }
      if( strcmp(corps[i], "C2H") == 0 )
      {
         MASS[i] = 25.03e0;
      }
      if( strcmp(corps[i], "C2H3") == 0 )
      {
         MASS[i] = 27.05e0;
      }
      if( strcmp(corps[i], "C2H4") == 0 )
      {
         MASS[i] = 28.05e0;
      }
      if( strcmp(corps[i], "C2H2") == 0 )
      {
         MASS[i] = 26.04e0;
      }
      if( strcmp(corps[i], "C2H5") == 0 )
      {
         MASS[i] = 29.06e0;
      }
      if( strcmp(corps[i], "C2H6") == 0 )
      {
         MASS[i] = 30.07e0;
      }
      if( strcmp(corps[i], "C3H2") == 0 )
      {
         MASS[i] = 38.05e0;
      }
      if( strcmp(corps[i], "C3H3") == 0 )
      {
         MASS[i] = 39.06e0;
      }
      if( ( strcmp(corps[i], "CH2CCH2") == 0 ) || ( strcmp(corps[i], "CH3CCH") == 0 ) )
      {
         MASS[i] = 40.07e0;
      }
      if( strcmp(corps[i], "C3H5") == 0 )
      {
         MASS[i] = 41.07e0;
      }
      if( strcmp(corps[i], "C3H6") == 0 )
      {
         MASS[i] = 42.08e0;
      }
      if( strcmp(corps[i], "C3H7") == 0 )
      {
         MASS[i] = 43.09e0;
       }
      if( strcmp(corps[i], "C3H8") == 0 )
      {
         MASS[i] = 44.11e0;
      }
      if( strcmp(corps[i], "C4H") == 0 )
      {
         MASS[i] = 49.05e0;
      }
      if( ( strcmp(corps[i], "C4H2") == 0 )||( strcmp(corps[i], "C4H2s") == 0 ) )
      {
         MASS[i] = 50.06e0;
      }
      if( strcmp(corps[i], "C4H3") == 0 )
      {
         MASS[i] = 51.07e0;
      }
      if( strcmp(corps[i], "C4H4") == 0 )
      {
         MASS[i] = 52.08e0;
      }
      if( strcmp(corps[i], "C4H5") == 0 )
      {
         MASS[i] = 53.07e0;
      }
      if( strcmp(corps[i], "C4H6") == 0 )
      {
         MASS[i] = 54.09e0;
      }
      if( strcmp(corps[i], "C4H10") == 0 )
      {
         MASS[i] = 58.13e0;
      }
      if( strcmp(corps[i], "C6H") == 0 )
      {
         MASS[i] = 73.07e0;
      }
      if( strcmp(corps[i], "C6H2") == 0 )
      {
         MASS[i] = 74.08e0;
      }
      if( strcmp(corps[i], "C8H2") == 0 )
      {
         MASS[i] = 98.10e0;
      }
      if( strcmp( corps[i], "AC6H6" ) == 0 )
      {
         MASS[i] = 78.1136e0;
      }
      if( ( strcmp( corps[i], "C6H5" ) == 0 ) || ( strcmp( corps[i], "AC6H5" ) == 0 ) )
      {
         MASS[i] = 77.1136e0;
      }
      if( strcmp( corps[i], "C6H6" ) == 0 )
      {
         MASS[i] = 78.1136e0;
      }
      if( strcmp(corps[i], "N2") == 0 )
      {
         MASS[i] = 28.0134e0;
      }
      if( strcmp(corps[i], "N4S") == 0 )
      {
         MASS[i] = 14.01e0;
      }
      if( strcmp(corps[i], "NH") == 0 )
      {
         MASS[i] = 15.01e0;
      }
      if( strcmp(corps[i], "CN") == 0 )
      {
         MASS[i] = 26.02e0;
      }
      if( strcmp(corps[i], "HCN") == 0 )
      {
         MASS[i] = 27.04e0;
      }
      if( strcmp(corps[i], "H2CN") == 0 )
      {
         MASS[i] = 28.05e0;
      }
      if( strcmp(corps[i], "C2N") == 0 )         /* C2N */
      {
         MASS[i] = 39.05e0;
      }
      if( strcmp( corps[i], "CHCN" ) == 0 )
      {
         MASS[i]   = 39.05e0;
      }
      if( strcmp( corps[i], "CH2CN" ) == 0 )
      {
         MASS[i]   = 40.04e0;
      }
      if( strcmp( corps[i], "CH3CN" ) == 0 )
      {
         MASS[i]   = 41.05e0;
      }
      if( strcmp( corps[i], "C2H3CN" ) == 0 )
      {
         MASS[i]   = 53.06e0;
      }
      if( strcmp(corps[i], "NCCN") == 0 )        /* NCCN */
      {
         MASS[i] = 52.04e0;
      }
      if( strcmp(corps[i], "C3N") == 0 )         /* C3N */
      {
         MASS[i] = 50.04e0;
      }
      if( strcmp(corps[i], "HC3N") == 0 )        /* HC3N */
      {
         MASS[i] = 51.05e0;
      }
      if( strcmp( corps[i], "C4N2" ) == 0 )
      {
         MASS[i]   = 76.1e0;
      }
      if( strcmp(corps[i], "H2O") == 0 )       
      {
         MASS[i] = 18.02e0;
      }
      if( ( strcmp(corps[i], "O3P") == 0 ) || ( strcmp(corps[i], "O1D") == 0 ) )      
      {
         MASS[i] = 16.0e0;
      }
      if( strcmp(corps[i], "OH") == 0 )        
      {
         MASS[i] = 17.01e0;
      }
      if( strcmp(corps[i], "CO") == 0 )          
      {
         MASS[i] = 28.01e0;
      }
      if( strcmp(corps[i], "HCO") == 0 )        
      {
         MASS[i] = 29.02e0;
      }
      if( strcmp(corps[i], "CO2") == 0 )       
      {
         MASS[i] = 44.01e0;
      }
      if( strcmp(corps[i], "CH2CO") == 0 )    
      {
         MASS[i] = 42.04e0;
      }
      if( strcmp(corps[i], "CH2O") == 0 )    
      {
         MASS[i] = 30.03e0;
      }
      if( ( strcmp(corps[i], "CH2OH") == 0 ) || ( strcmp(corps[i], "CH3O") == 0 ) )  
      {
         MASS[i] = 31.04e0;
      }
      if( strcmp(corps[i], "CH3OH") == 0 )       
      {
         MASS[i] = 32.042e0;
      }
   }
}