chem/KPP/kpp/kpp-2.1/drv/main.c

   1 void INTEGRATE( double TIN, double TOUT );
   2
   3 void main()
   4 {
   5 KPP_REAL rtols;
   6 KPP_REAL Tstart, Tend, DT;
   7 FILE * fp;
   8 int i;
   9
  10 /* ---- TIME VARIABLES ------------------ */
  11
  12   rtols = 1e-3;
  13   for( i = 0; i < NVAR; i++ ) {
  14     RTOL[i] = rtols;
  15     ATOL[i] = 1E-18;
  16   }
  17
  18   Initialize();
  19
  20   fp = fopen("Extrapd.m", "w");
  21   fprintf( fp, "ed=[\n");
  22
  23   Tstart = 3600*12;
  24   Tend = Tstart + 3600*24*5;
  25   STEPMIN = 0.001;
  26   STEPMAX = 900;
  27   DT = 3600.;
  28   TEMP = 236.21;
  29
  30 /* -- BELOW THIS LIMIT USE TAYLOR INSTEAD OF EXP --- */
  31 /* ********** TIME LOOP **************************** */
  32
  33   TIME = Tstart;
  34   while (TIME <= Tend) {
  35
  36     printf("\nMonitor: ");
  37     for( i = 0; i < NMONITOR; i++ )
  38       printf( "%12.8g  ", C[ MONITOR[i] ]/CFACTOR );
  39
  40     fprintf( fp, "\n%6.1f ", (TIME-Tstart)/3600.0 );
  41     for( i = 0; i < NLOOKAT; i++ )
  42       fprintf( fp, "%24.16e ", C[ LOOKAT[i] ]/CFACTOR );
  43
  44     INTEGRATE( DT );
  45   }
  46
  47 /* *********** END TIME LOOP *********************** */
  48
  49   fprintf(fp, "\n];");
  50   fclose( fp );
  51   printf("\n");
  52 }