Colvars Error?

From: David Minh (daveminh_at_gmail.com)
Date: Fri Oct 26 2012 - 11:13:21 CDT

Hello NAMD users and developers,

I am having a strange issue related to colvars. I am trying to do a
coorfile loop in my NAMD configuration file:

coorfile open dcd {begin1.dcd}
while { ![coorfile read] } {
  reinitvels $temperature
  run 100
}
coorfile close

It works fine if I remove the lines:

colvars on
colvarsConfig gorgeCylinder.colvars

Here are contents of gorgeCylinder.colvars:

colvar {
  name DistXY
  upperWall 6.0
  upperWallConstant 10.0
  distanceXY {
  main {
    atomsFile {../../../ligands/AAA/ligand.grid.pdb}
    atomsCol B
    atomsColValue 1.0
  }
  ref {
    dummyAtom (52.548195, 57.459400, 49.528568)
  }
  axis (-0.294390, -0.561701, 0.773192)

  }
}

colvar {
  name DistZ
  lowerWall 3
  upperWall 25
  lowerWallConstant 10.0
  upperWallConstant 10.0
  distanceZ {
  main {
    atomsFile {../../../ligands/AAA/ligand.grid.pdb}
    atomsCol B
    atomsColValue 1.0
  }
  ref {
    dummyAtom (52.548195, 57.459400, 49.528568)
  }
  axis (-0.294390, -0.561701, 0.773192)

  }
}

Basically it keeps the molecule within a cylinder.

If I include the colvars lines in the NAMD configuration file, the code
runs for 100 steps and then says:

FATAL ERROR: Got the wrong number of atoms
------------- Processor 0 Exiting: Called CmiAbort ------------
Reason: FATAL ERROR: Got the wrong number of atoms

Charm++ fatal error:
FATAL ERROR: Got the wrong number of atoms

Any ideas on what could cause this? The error doesn't really make sense
because it works fine without the colvars lines. Is the error really
related to colvarsInput? I can't seem to put a colvarsInput into the tcl
loop.

Thanks,
David Minh

This archive was generated by hypermail 2.1.6 : Mon Dec 31 2012 - 23:22:12 CST