[AMBER] Benign typo for a constant in barostats.F90

From: Neale, Christopher Andrew <cneale.lanl.gov>
Date: Mon, 19 Dec 2016 21:52:00 +0000

I think there's a minor typo in Amber16 src/pmemd/src/barostats.F90:

  ! p*dV (6.02204d-2 / 4184.0d0 converts from bar*A^3/particle to kcal/mol)
  pv_work = pres0 * (uc_volume - orig_vol) * 6.02204d-2 / 4184.0d0

Where I presume 6.02204d-2 was intended to be 6.02214d-2 ?

There's in some cases a truncated version of Avogadro's constant, e.g.:
barostats.F90: ! converts dyne.A^2/cm to kcal/mol (2.390057e-27 * 6.022e23)

but I can't see any reason to have a value of 6.02204d-2 coming from 1e-25 N/m * Avogadro's constant.

Not complaining. Just posting so it can get fixed.
AMBER mailing list
Received on Mon Dec 19 2016 - 14:00:03 PST
Custom Search