From: Jose Borreguero (borreguero_at_gmail.com)
Date: Tue Apr 27 2010 - 14:04:11 CDT
Dear NAMD users,
I have a trajectory (dcd format) resulting of merging 10 diffferent
trajectories, each with each own box size. Each snapshot contains its own
box size, though.
Now I want to evaluate the energy of each snapshot. I used PME and I am
passing an extended file system (box.xsc) defining a box size corresponding
to the first of the ten trajectories.
My question is whether namd will read the box size of the extended system
file and assume it is the same for every single snapshot, or if it will read
the box size of each snapshot.
The last behavior is the one I need so that namd correctly applies the
periodic boundary conditions and the PME calculations.
Here's my configuration file:
amber on
parmfile ./top
numsteps 1
switching off
scnb 2
exclude scaled1-4
1-4scaling 0.833333
outputname junk
temperature 0
COMmotion yes
cutoff 10
dielectric 1.0
extendedSystem ./box.xsc
PME on
PMEGridSizeX 60
PMEGridSizeY 54
PMEGridSizeZ 60
switchdist 10
pairInteraction on
pairInteractionGroup1 1
pairInteractionFile ./all.pdb
pairInteractionSelf on
coordinates ./all.pdb
set ts 0
coorfile open dcd dcd
while { ![coorfile read] } {
firstTimestep $ts
run 0
incr ts 1
}
coorfile close
Best regards,
-Jose Borreguero
This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:54:03 CST