chem/KPP/util/write_decomp/integr_edit.c

   1 #include <stdio.h>
   2 #include <stdlib.h>
   3
   4 #define NAMELEN 4096
   5
   6
   7 /* replace decomp routine in KPP Integr file */
   8
   9
  10 main( int argc, char *argv[] )
  11 {
  12
  13   char mechname[NAMELEN];
  14   char intfname[NAMELEN], incfname[NAMELEN], tfname[NAMELEN];
  15   char inln[NAMELEN], incln[NAMELEN];
  16   char callln[NAMELEN], endln[NAMELEN] ;
  17   char cp_command[NAMELEN];
  18
  19   int copyit, add_sub;
  20
  21   FILE * intf;
  22   FILE * incf;
  23   FILE * ofile;
  24
  25
  26   argv++ ;
  27
  28   if ( *argv )
  29     {
  30   strcpy( mechname, *argv );
  31     }
  32   else
  33     {
  34      printf("ERROR: USAGE: integr_edit mech_name\n");
  35      exit(11);
  36     }
  37
  38   sprintf ( intfname, "module_kpp_%s_Integr.F", mechname);
  39   sprintf ( incfname, "decomp_%s.inc", mechname);
  40   sprintf ( tfname, "%s_new", intfname );
  41
  42
  43
  44   sprintf(  cp_command,"cp %s %s",tfname, intfname);
  45
  46   intf = fopen( intfname , "r" );
  47   incf = fopen( incfname , "r" );
  48   ofile = fopen( tfname , "w" );
  49
  50
  51
  52   sprintf ( callln , "   CALL %s_KppDecomp\0", mechname );
  53   sprintf ( endln , "END MODULE");
  54
  55      /* loop over lines in Integr file */
  56        while ( fgets ( inln , 4096 , intf ) != NULL ){
  57
  58            copyit=1;
  59
  60
  61            /* replace call to decomp routine */
  62
  63            if ( !strncmp (inln, callln, strlen(callln)-1) ) {
  64
  65           printf("   integr_edit: replacing  %s \n", inln);
  66
  67           fprintf(ofile, "!!!  use direct adressing in decomp \n");
  68           fprintf(ofile, "!!! %s", inln);
  69           fprintf(ofile, "CALL decomp_%s ( A, ising )\n", mechname );
  70
  71
  72           add_sub=1;
  73           copyit=0;
  74            }
  75
  76
  77            /* add decomp routine w. direct referncing */
  78
  79           if ( !strncmp (inln, endln, strlen(endln)-1) ) {
  80
  81             if ( add_sub ){
  82             printf("  %s ", inln );
  83                while ( fgets ( incln , 4096 , incf ) != NULL ){
  84
  85                 fprintf(ofile, "%s", incln);
  86
  87                }
  88
  89                 fprintf(ofile, " \n\n\n");
  90             }
  91           }
  92
  93
  94
  95
  96        /* copy line from original file */
  97          if ( copyit ) {
  98
  99            fprintf(ofile, "%s", inln);
 100
 101          }
 102
 103
 104      }
 105
 106
 107        if ( ! add_sub  ) {
 108          printf(" integr_edit: Kept previous version. \n ");
 109        }
 110
 111
 112
 113   fclose( intf );
 114   fclose( incf );
 115   fclose( ofile );
 116
 117
 118    system(cp_command);
 119
 120   exit (0);
 121 }