1 #ifndef REACTION_H
2 #define REACTION_H
3 
4 #include "chemical.hpp"
5 using namespace std;
6 
7 class reaction {
8 private:
9 
10   int m;
11   double *n, k0, E, Hr, *safe_n, *safe_a;
12   // double tmp;
13 //   char file[41], line[31];
14   chemical ** list;
15 //   ifstream in;
16   // ofstream logf;
17   // terminator *end;
18 
19   int find_chemical ( const string & chem_name ) const;
20 
21 public:
22   // reaction(){};
23   reaction ( const string & , int ,  chemical ** );
24   ~reaction();
25   double *a;									//contains the molar coefficients
26   double dHr(double);					//returns heat of rection at T, in kJ/mol
27   double rate(double, double*);		//returns rate of reaction aT and C[], in mol/s.m3
28 
29 //   void show_name(){cout<<name;}
30 };
31 
32 #endif
33