git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@13402 f3b2605a-c512-4ea7-a41b-209d697bcdaa

This commit is contained in:
sjplimp
2015-04-11 23:17:54 +00:00
parent 29e574b720
commit 4ae6abbea5

View File

@ -197,7 +197,10 @@ void ComputeTempProfile::dof_compute()
adjust_dof_fix();
natoms_temp = group->count(igroup);
dof = domain->dimension * natoms_temp;
dof -= extra_dof + fix_dof;
// subtract additional d*Nbins DOF, as in Evans and Morriss paper
dof -= extra_dof + fix_dof + domain->dimension*nbins;
if (dof > 0) tfactor = force->mvv2e / (dof * force->boltz);
else tfactor = 0.0;
}