chem/KPP/util/wkc/get_kpp_chem_specs.c

   1 #include <stdio.h>
   2 #include <stdlib.h>
   3 #include <string.h>
   4 #include <strings.h>
   5 #include <dirent.h>
   6 #include <unistd.h>
   7 #include <sys/stat.h>
   8
   9
  10 #include "protos.h"
  11 #include "protos_kpp.h"
  12 #include "registry.h"
  13 #include "data.h"
  14 #include "kpp_data.h"
  15
  16 #define DEBUGR 0
  17
  18 int
  19 get_kpp_chem_specs ( char* kpp_dirname )
  20 {
  21
  22
  23 knode_t * q , * member  ;
  24 DIR * dir;
  25 struct dirent * entry;
  26 struct stat dir_stat;
  27 char fulldirname[NAMELEN], spcfilename[NAMELEN];
  28 char inln[NAMELEN], kpp_spec[NAMELEN];
  29 FILE * spcFile;
  30 int   in_comment, got_it;
  31
  32
  33
  34   /* http://users.actcom.co.il/~choo/lupg/tutorials/handling-files/handling-files.html#directory_struct */
  35
  36
  37
  38     dir = opendir(kpp_dirname);
  39     if (!dir) {
  40         fprintf(stderr, "WARNING from gen_kpp: Cannot read directory: %s \n", kpp_dirname);
  41         perror("");
  42         return;
  43     }
  44
  45
  46     /* loop through sub directories in KPP directory */
  47
  48  while ((entry = readdir(dir))) {
  49    if (entry->d_name ) {
  50
  51          if ( strcmp(entry->d_name, ".") == 0)
  52           continue;
  53          if ( strcmp(entry->d_name, "..") == 0)
  54            continue;
  55
  56
  57
  58       sprintf(  fulldirname, "%s/%s", kpp_dirname, entry->d_name);
  59
  60       printf("%s \n", fulldirname );
  61
  62        /* check if the given entry is a directory. */
  63         if (stat(fulldirname, &dir_stat) == -1) {
  64           fprintf(stderr, "WAA\n\n");
  65            perror("WARNING from gen_kpp: ");
  66             continue;
  67         }
  68
  69
  70         /* check if KPP species file is present. */
  71
  72          sprintf(  spcfilename, "%s/%s/%s.spc", kpp_dirname, entry->d_name, entry->d_name);
  73
  74
  75          ;fprintf(stderr, " spcfilename: %s \n",spcfilename);
  76
  77       spcFile = fopen (spcfilename, "r" );
  78
  79         if ( spcFile == NULL ) {
  80           fprintf(stderr,"WARNING from gen_kpp: File %s not found. Skipping. \n", spcfilename);
  81            continue;
  82         }
  83
  84    printf(" Found %s \n", spcfilename );
  85
  86
  87
  88
  89    /*----------------------------------------------------*/
  90
  91
  92         /* put KPP packagename into linked list */
  93
  94          q = new_knode( );
  95          q->next = NULL ;
  96          strcpy( q->name, entry->d_name  );
  97          add_knode_to_end( q , &(KPP_packs) ) ;
  98
  99          /* loop over lines in KPP species file */
 100          while ( fgets ( inln , NAMELEN , spcFile ) != NULL ){
 101               if ( DEBUGR == 1 ){    printf("%s  ", inln); }
 102          /* strip from comments (loop through letters) */
 103          int n=0;
 104          int nn = 0;
 105          int j;
 106          in_comment = 0;
 107          got_it = 0;
 108
 109             for(j = 0; j < NAMELEN ; j++) kpp_spec[j]='\0';
 110           while ( inln[n] !=  '\0' ){
 111            if ( inln[n] == '{') in_comment=1;
 112             if ( in_comment == 0 ) {
 113                     if (inln[n] == '=' || inln[n] == '#') {
 114                         got_it=1;
 115                     }
 116                       if ( got_it == 0  && inln[n] !=  ' '){
 117                     /* printf("%c %i \n ", inln[n], in_comment ); */
 118
 119                        kpp_spec[nn]=inln[n];
 120                        nn++;
 121
 122                       }
 123           }
 124
 125           if (inln[n] == '}') in_comment=0;
 126           n++;
 127
 128           }
 129
 130           /* printf("spec: %s \n ", kpp_spec);  */
 131
 132           if (kpp_spec[0] != '\0'  && got_it == 1 )  {
 133
 134             if ( DEBUGR == 1 ){
 135                   printf("spec: %s \n ", kpp_spec);
 136                   fprintf(stderr," p, name   %s %s \n", q->name, kpp_spec );
 137                }
 138
 139            member = new_knode( ) ;
 140            strcpy( member->name , kpp_spec )   ;
 141            member->next = NULL ;
 142            add_knode_to_end( member , &(q->members) ) ;
 143
 144           }
 145          }
 146
 147
 148
 149          fclose(spcFile);
 150
 151     }
 152
 153  }
 154
 155   return(0) ;
 156 }
 157