diff --git a/src/fix_gravity.cpp b/src/fix_gravity.cpp index 37e605f74d..24bc7b4b3c 100644 --- a/src/fix_gravity.cpp +++ b/src/fix_gravity.cpp @@ -248,12 +248,12 @@ void FixGravity::post_force(int vflag) if (varflag != CONSTANT) { modify->clearstep_compute(); if (mstyle == EQUAL) magnitude = input->variable->compute_equal(mvar); - if (vstyle == EQUAL) magnitude = input->variable->compute_equal(vvar); - if (pstyle == EQUAL) magnitude = input->variable->compute_equal(pvar); - if (tstyle == EQUAL) magnitude = input->variable->compute_equal(tvar); - if (xstyle == EQUAL) magnitude = input->variable->compute_equal(xvar); - if (ystyle == EQUAL) magnitude = input->variable->compute_equal(yvar); - if (zstyle == EQUAL) magnitude = input->variable->compute_equal(zvar); + if (vstyle == EQUAL) vert = input->variable->compute_equal(vvar); + if (pstyle == EQUAL) phi = input->variable->compute_equal(pvar); + if (tstyle == EQUAL) theta = input->variable->compute_equal(tvar); + if (xstyle == EQUAL) xdir = input->variable->compute_equal(xvar); + if (ystyle == EQUAL) ydir = input->variable->compute_equal(yvar); + if (zstyle == EQUAL) zdir = input->variable->compute_equal(zvar); modify->addstep_compute(update->ntimestep + 1); set_acceleration();