| 83 |
|
|
| 84 |
|
double Thermo::getTemperature(){ |
| 85 |
|
|
| 86 |
< |
const double kb = 1.88E-3; // boltzman's constant in kcal/(mol K) |
| 86 |
> |
const double kb = 1.9872179E-3; // boltzman's constant in kcal/(mol K) |
| 87 |
|
double temperature; |
| 88 |
|
|
| 89 |
|
int ndf = 3 * entry_plug->n_atoms + 3 * entry_plug->n_oriented |
| 95 |
|
|
| 96 |
|
double Thermo::getPressure(){ |
| 97 |
|
|
| 98 |
< |
const double conv_Pa_atm = 9.901E-6; // convert Pa -> atm |
| 99 |
< |
const double conv_internal_Pa = 1.661E-7; //convert amu/(fs^2 A) -> Pa |
| 100 |
< |
const double conv_A_m = 1.0E-10; //convert A -> m |
| 98 |
> |
// const double conv_Pa_atm = 9.901E-6; // convert Pa -> atm |
| 99 |
> |
// const double conv_internal_Pa = 1.661E-7; //convert amu/(fs^2 A) -> Pa |
| 100 |
> |
// const double conv_A_m = 1.0E-10; //convert A -> m |
| 101 |
|
|
| 102 |
|
return 0.0; |
| 103 |
|
} |