--- trunk/OOPSE/libmdtools/Restraints.cpp 2004/05/22 18:16:18 1187 +++ trunk/OOPSE/libmdtools/Restraints.cpp 2004/05/22 20:55:04 1189 @@ -140,6 +140,7 @@ double Restraints::Calc_Restraint_Forces(vectorgetA(A); Calc_rVal( pos, i ); Calc_body_thetaVal( A, i ); - vecParticles[i]->getZangle(omega); - Calc_body_omegaVal( A, omega ); + omegaPass = vecParticles[i]->getZangle(); + Calc_body_omegaVal( A, omegaPass ); if (omega > PI || omega < -PI) cout << "oops... " << omega << "\n";