[gmx-users] Parallel pulling with Gromacs 4.0.7: COMM mode problem

[Aykut Erbas]

Hi everyone

I have a pulling code running in parallel (32 processors) with Gromacs 4.0.7
The systems composes of a surface and a protein

with Gromacs 3.3.3 (on single processors ) when I set the COMM mode for 
the surface, the surface stands still,
and of course for the pulling code (reference groups is also surface) , 
the position of the protein is given /wrt the surface.
However for G4, no matter which group you set for COMM mode, the surface 
moves in the opposite direction of motion and you have the pulling 
coordinates in absolute coord. which is not what I want.

I want a pilling output /wrt the surface as in G3

Has anyone experienced smtg like this?
suggestions will be appreciated

Thanks
Aykut


*.mdp file
DIAM=surface
//////////////////////////////////////////

integrator               = md
nsteps                    = 30000000
dt                           = 0.002
nstlist                     = 50
nstxout                  = 20000
nstvout                  = 20000
nstfout                   = 20000
nstxtcout               = 4000
nstlog                    = 20000
nstenergy              = 20000
constraints             = hbonds
ns_type                 = grid
coulombtype         = pme
pme_order             = 4
fourierspacing           = 0.12
rlist                        = 0.8
rvdw                     = 0.8
rcoulomb               = 0.8
energygrps            = DIAM Protein SOL
tcoupl                   = berendsen
tc_grps                  = DIAM Protein SOL
tau_t                    = 0.1 0.1 0.1
ref_t                    = 300 300 300
gen_vel                  = yes
Pcoupl                   = berendsen
Pcoupltype               = semiisotropic
compressibility          = 2.5E-8 4.5E-5
tau_p                    = 1.0 1.0
ref_p                    = 1.0 1.0
comm_mode                = angular
nstcomm                  = 1
comm_grps                = DIAM

;PULLING
pull                     = umbrella
pull_start               = yes
pull_geometry            = position
pull_dim                 = Y Y N
pull_nstxout             = 200
pull_nstfout             = 200
pull_ngroups             = 1
pull_group0              = DIAM
pull_group1              = AA1
pull_vec1                = 1.0 0.0 0.0
pull_init1               = 2.0 0.0 1.0
pull_rate1               = 0.001
pull_k1                  = 50
///////////////////////////////////////////