git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@591 f3b2605a-c512-4ea7-a41b-209d697bcdaa
This commit is contained in:
@ -71,7 +71,7 @@ void FixHeat::end_of_step()
|
|||||||
|
|
||||||
double vsub[3],vcm[3];
|
double vsub[3],vcm[3];
|
||||||
|
|
||||||
double ke = group->ke(igroup);
|
double ke = group->ke(igroup)*force->ftm2v;
|
||||||
group->vcm(igroup,masstotal,vcm);
|
group->vcm(igroup,masstotal,vcm);
|
||||||
double vcmsq = vcm[0]*vcm[0] + vcm[1]*vcm[1] + vcm[2]*vcm[2];
|
double vcmsq = vcm[0]*vcm[0] + vcm[1]*vcm[1] + vcm[2]*vcm[2];
|
||||||
double escale = (ke + heat - 0.5*vcmsq*masstotal)/(ke - 0.5*vcmsq*masstotal);
|
double escale = (ke + heat - 0.5*vcmsq*masstotal)/(ke - 0.5*vcmsq*masstotal);
|
||||||
|
|||||||
Reference in New Issue
Block a user