Setting molecule COM to 0 after moving atoms

This commit is contained in:
Anders Hafreager
2017-04-26 20:10:18 +02:00
parent 2f32fb7f8b
commit 67fced37c8

View File

@ -595,6 +595,9 @@ void FixGCMC::init()
onemols[imol]->x[i][1] -= onemols[imol]->com[1];
onemols[imol]->x[i][2] -= onemols[imol]->com[2];
}
onemols[imol]->com[0] = 0;
onemols[imol]->com[1] = 0;
onemols[imol]->com[2] = 0;
} else gas_mass = atom->mass[ngcmc_type];