57 |
|
painCave.isFatal = 1; |
58 |
|
simError(); |
59 |
|
} else { |
60 |
< |
surfaceTension= simParams->getSurfaceTension(); |
60 |
> |
surfaceTension= simParams->getSurfaceTension()* OOPSEConstant::surfaceTensorConvert * OOPSEConstant::energyConvert; |
61 |
|
} |
62 |
|
|
63 |
|
} |
175 |
|
double area = hmat(0,0) * hmat(1, 1); |
176 |
|
|
177 |
|
double conservedQuantity = totalEnergy + thermostat_kinetic + thermostat_potential + |
178 |
< |
barostat_kinetic + barostat_potential - surfaceTension * area; |
178 |
> |
barostat_kinetic + barostat_potential - surfaceTension * area/ OOPSEConstant::energyConvert; |
179 |
|
|
180 |
|
return conservedQuantity; |
181 |
|
|