[gmx-developers] introduce effective potential

Berk Hess gmx3 at hotmail.com
Tue Apr 20 14:35:35 CEST 2004

>Could you tell me the subroutine which is needed to be done to get
>correct the pressure?

The shift forces fr->fshift are updated in do_dih_fup,
so if you only apply a scale factor to ddphi before this call
everything is consistent.

>At the moment, whenever I try with each value of q and epsilon, I have to
>compile the gromacs. It is very time consumming and of course it is a very
>bad way. How can I introduce these parameters into the *mdp file?

You could do it with an enviroment variable.
For that you need to include the following code:

static double q=-1;
char *ptr;

if (q == -1)
  ptr = getenv("GMX_Q");
  if (ptr) {
  } else {
    q = 1;


MSN Search, for accurate results! http://search.msn.nl

More information about the gromacs.org_gmx-developers mailing list