| 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 | 
  | 
 |