Changed the check on initial and final temperature to <= 0 for both.

This commit is contained in:
Stefan Paquay
2017-03-28 11:22:10 -04:00
parent 3dfe4505dd
commit ee55a98103

View File

@ -133,7 +133,7 @@ FixNH::FixNH(LAMMPS *lmp, int narg, char **arg) :
t_target = t_start;
t_stop = force->numeric(FLERR,arg[iarg+2]);
t_period = force->numeric(FLERR,arg[iarg+3]);
if (t_start < 0.0 || t_stop <= 0.0)
if (t_start <= 0.0 || t_stop <= 0.0)
error->all(FLERR,
"Target temperature for fix nvt/npt/nph cannot be 0.0");
iarg += 4;