diff --git a/src/ASPHERE/compute_temp_asphere.cpp b/src/ASPHERE/compute_temp_asphere.cpp index e39b82652d..14a737742d 100755 --- a/src/ASPHERE/compute_temp_asphere.cpp +++ b/src/ASPHERE/compute_temp_asphere.cpp @@ -70,12 +70,14 @@ void ComputeTempAsphere::init() fix_dof += modify->fix[i]->dof(igroup); recount(); + tempbias = 0; + tbias = NULL; if (id_bias) { tempbias = 1; int i = modify->find_compute(id_bias); if (i < 0) error->all("Could not find compute ID for temperature bias"); tbias = modify->compute[i]; - } else tempbias = 0; + } calculate_inertia(); } diff --git a/src/compute_temp.cpp b/src/compute_temp.cpp index 8b2a8931c2..451247f910 100644 --- a/src/compute_temp.cpp +++ b/src/compute_temp.cpp @@ -58,6 +58,8 @@ void ComputeTemp::init() fix_dof += modify->fix[i]->dof(igroup); recount(); + tempbias = 0; + tbias = NULL; if (id_bias) { tempbias = 1; int i = modify->find_compute(id_bias); diff --git a/src/compute_temp_com.cpp b/src/compute_temp_com.cpp index 3b808c7cf2..2cfbf72164 100644 --- a/src/compute_temp_com.cpp +++ b/src/compute_temp_com.cpp @@ -66,6 +66,7 @@ void ComputeTempCOM::init() recount(); masstotal = group->mass(igroup); + tbias = NULL; if (id_bias) { int i = modify->find_compute(id_bias); if (i < 0) error->all("Could not find compute ID for temperature bias"); diff --git a/src/compute_temp_deform.cpp b/src/compute_temp_deform.cpp index cca5e0e5b5..40ba16136e 100644 --- a/src/compute_temp_deform.cpp +++ b/src/compute_temp_deform.cpp @@ -73,6 +73,7 @@ void ComputeTempDeform::init() fix_dof += modify->fix[i]->dof(igroup); recount(); + tbias = NULL; if (id_bias) { int i = modify->find_compute(id_bias); if (i < 0) error->all("Could not find compute ID for temperature bias"); diff --git a/src/compute_temp_partial.cpp b/src/compute_temp_partial.cpp index eaacee47e4..e9b17ed930 100644 --- a/src/compute_temp_partial.cpp +++ b/src/compute_temp_partial.cpp @@ -67,6 +67,7 @@ void ComputeTempPartial::init() fix_dof += modify->fix[i]->dof(igroup); recount(); + tbias = NULL; if (id_bias) { int i = modify->find_compute(id_bias); if (i < 0) error->all("Could not find compute ID for temperature bias"); diff --git a/src/compute_temp_ramp.cpp b/src/compute_temp_ramp.cpp index 91f6409816..f068f233bb 100644 --- a/src/compute_temp_ramp.cpp +++ b/src/compute_temp_ramp.cpp @@ -132,6 +132,7 @@ void ComputeTempRamp::init() fix_dof += modify->fix[i]->dof(igroup); recount(); + tbias = NULL; if (id_bias) { int i = modify->find_compute(id_bias); if (i < 0) error->all("Could not find compute ID for temperature bias"); diff --git a/src/compute_temp_region.cpp b/src/compute_temp_region.cpp index efab6219de..873d1ae2eb 100644 --- a/src/compute_temp_region.cpp +++ b/src/compute_temp_region.cpp @@ -63,6 +63,7 @@ void ComputeTempRegion::init() { dof = 0; + tbias = NULL; if (id_bias) { int i = modify->find_compute(id_bias); if (i < 0) error->all("Could not find compute ID for temperature bias"); diff --git a/src/compute_temp_sphere.cpp b/src/compute_temp_sphere.cpp index 4145db51c7..4e381b11e4 100644 --- a/src/compute_temp_sphere.cpp +++ b/src/compute_temp_sphere.cpp @@ -62,12 +62,14 @@ void ComputeTempSphere::init() fix_dof += modify->fix[i]->dof(igroup); recount(); + tempbias = 0; + tbias = NULL; if (id_bias) { tempbias = 1; int i = modify->find_compute(id_bias); if (i < 0) error->all("Could not find compute ID for temperature bias"); tbias = modify->compute[i]; - } else tempbias = 0; + } if (atom->mass) { double *mass = atom->mass;