Re: RAMD Setting parameter TclForces from script failed

From: Vlad Cojocaru (vlad.cojocaru_at_mpi-muenster.mpg.de)
Date: Tue Mar 29 2011 - 04:59:25 CDT

Francesco,

Does the colvar module work with othe tcl forces scripts ?

Vlad

On 03/29/2011 09:51 AM, Francesco Pietra wrote:
> Hello Jim:
> It starts to become clear this morning (hopefully). It seems to me
> that the problem is the simultaneous presence of colvars and ramd (a
> situation that I have emphasized on previous posts, without having any
> reason for that). Both involve individual atom requests, which is not
> allowed in parallel simulations.
>
> Thus, sourcing ramd script from only the ramd portion of namd conf
> file, and having "numsteps" (probably the same having "run") at the
> end of namd conf file (after ramd), the error is:
>
> FATAL ERROR: Due to a design error, GlobalMasterServer does not
> support individual atom requests from multiple global force clients on
> parallel runs.
>
> Here the relevant portions of the log file, where I have shortened the
> portion related to colvars:
>
> arm++> cpu topology info is gathered in 0.001 seconds.
> Info: NAMD CVS-2011-02-14 for Linux-x86_64
> Info:
> Info: Please visit http://www.ks.uiuc.edu/Research/namd/
> Info: and send feedback or bug reports to namd_at_ks.uiuc.edu
> Info:
> Info: Please cite Phillips et al., J. Comp. Chem. 26:1781-1802 (2005)
> Info: in all publications reporting results obtained with NAMD.
> Info:
> Info: Based on Charm++/Converse 60202 for net-linux-x86_64-iccstatic
> Info: Built Mon Feb 14 02:16:25 CST 2011 by jim on lisboa.ks.uiuc.edu
> Info: 1 NAMD CVS-2011-02-14 Linux-x86_64 4 tya64 francesco
> Info: Running on 4 processors.
> Info: CPU topology information available.
> Info: Charm++/Converse parallel runtime startup completed at 0.078851 s
> Info: 1.63774 MB of memory in use based on CmiMemoryUsage
> Info: Configuration file is ramd_acc-03.conf
> Info: Working in the current directory
> /home/francesco/work_hasic1a_GLUP_ASPP_allHSP_S-S/9.ramd_04-1.1
> TCL: Suspending until startup complete.
> Info: EXTENDED SYSTEM FILE ./press-02.restart.xsc
> Info: SIMULATION PARAMETERS:
> Info: TIMESTEP 1
> Info: NUMBER OF STEPS 100
> Info: STEPS PER CYCLE 20
> Info: PERIODIC CELL BASIS 1 122.1 0 0
> Info: PERIODIC CELL BASIS 2 0 119.2 0
> Info: PERIODIC CELL BASIS 3 0 0 148.634
> Info: PERIODIC CELL CENTER -0.0430445 -1.34946 -7.98849
> Info: WRAPPING ALL CLUSTERS AROUND PERIODIC BOUNDARIES ON OUTPUT.
> Info: WRAPPING TO IMAGE NEAREST TO PERIODIC CELL CENTER.
> Info: LOAD BALANCER Centralized
> Info: LOAD BALANCING STRATEGY New Load Balancers -- DEFAULT
> Info: LDB PERIOD 4000 steps
> Info: FIRST LDB TIMESTEP 100
> Info: LAST LDB TIMESTEP -1
> Info: LDB BACKGROUND SCALING 1
> Info: HOM BACKGROUND SCALING 1
> Info: PME BACKGROUND SCALING 1
> Info: MAX SELF PARTITIONS 20
> Info: MAX PAIR PARTITIONS 8
> Info: SELF PARTITION ATOMS 154
> Info: SELF2 PARTITION ATOMS 154
> Info: PAIR PARTITION ATOMS 318
> Info: PAIR2 PARTITION ATOMS 637
> Info: MIN ATOMS PER PATCH 100
> Info: VELOCITY FILE ./press-02.restart.vel
> Info: CENTER OF MASS MOVING INITIALLY? NO
> Info: DIELECTRIC 1
> Info: EXCLUDE SCALED ONE-FOUR
> Info: 1-4 ELECTROSTATICS SCALED BY 1
> Info: MODIFIED 1-4 VDW PARAMETERS WILL BE USED
> Info: DCD FILENAME ./ramd_acc-03.dcd
> Info: DCD FREQUENCY 100
> Info: DCD FIRST STEP 100
> Info: DCD FILE WILL CONTAIN UNIT CELL DATA
> Info: NO EXTENDED SYSTEM TRAJECTORY OUTPUT
> Info: NO VELOCITY DCD OUTPUT
> Info: OUTPUT FILENAME ./ramd_acc-03
> Info: RESTART FILENAME ./ramd_acc-03.restart
> Info: RESTART FREQUENCY 10
> Info: BINARY RESTART FILES WILL BE USED
> Info: SWITCHING ACTIVE
> Info: SWITCHING ON 10
> Info: SWITCHING OFF 12
> Info: PAIRLIST DISTANCE 13.5
> Info: PAIRLIST SHRINK RATE 0.01
> Info: PAIRLIST GROW RATE 0.01
> Info: PAIRLIST TRIGGER 0.3
> Info: PAIRLISTS PER CYCLE 2
> Info: PAIRLISTS ENABLED
> Info: MARGIN 3
> Info: HYDROGEN GROUP CUTOFF 2.5
> Info: PATCH DIMENSION 19
> Info: ENERGY OUTPUT STEPS 10
> Info: CROSSTERM ENERGY INCLUDED IN DIHEDRAL
> Info: TIMING OUTPUT STEPS 100
> Info: TCL GLOBAL FORCES ACTIVE
> Info: TCL GLOBAL FORCES SCRIPT
> /usr/local/namd_2.8nb/lib/ramd/scripts/ramd-4.1_script.tcl
> Info: COLLECTIVE VARIABLES CALCULATION REQUESTED
> Info: COLLECTIVE VARIABLES CONFIGURATION ./CLA-ligands_colvars_press01.in
> Info: COLLECTIVE VARIABLES RESTART INFORMATION
> ./press-02.restart.colvars.state
> Info: LANGEVIN DYNAMICS ACTIVE
> Info: LANGEVIN TEMPERATURE 310
> Info: LANGEVIN DAMPING COEFFICIENT IS 5 INVERSE PS
> Info: LANGEVIN DYNAMICS NOT APPLIED TO HYDROGENS
> Info: LANGEVIN PISTON PRESSURE CONTROL ACTIVE
> Info: TARGET PRESSURE IS 1.01325 BAR
> Info: OSCILLATION PERIOD IS 200 FS
> Info: DECAY TIME IS 100 FS
> Info: PISTON TEMPERATURE IS 310 K
> Info: PRESSURE CONTROL IS GROUP-BASED
> Info: INITIAL STRAIN RATE IS 1.30719e-05 3.88275e-05 2.39687e-06
> Info: CELL FLUCTUATION IS ANISOTROPIC
> Info: CONSTANT AREA PRESSURE CONTROL ACTIVE
> Info: PARTICLE MESH EWALD (PME) ACTIVE
> Info: PME TOLERANCE 1e-06
> Info: PME EWALD COEFFICIENT 0.257952
> Info: PME INTERPOLATION ORDER 4
> Info: PME GRID DIMENSIONS 128 120 160
> Info: PME MAXIMUM GRID SPACING 1
> Info: Attempting to read FFTW data from
> FFTW_NAMD_CVS-2011-02-14_Linux-x86_64.txt
> Info: Optimizing 6 FFT steps. 1... 2... 3... 4... 5... 6... Done.
> Info: Writing FFTW data to FFTW_NAMD_CVS-2011-02-14_Linux-x86_64.txt
> Info: FULL ELECTROSTATIC EVALUATION FREQUENCY 5
> Info: USING VERLET I (r-RESPA) MTS SCHEME.
> Info: C1 SPLITTING OF LONG RANGE ELECTROSTATICS
> Info: PLACING ATOMS IN PATCHES BY HYDROGEN GROUPS
> Info: RIGID BONDS TO HYDROGEN : WATER
> Info: ERROR TOLERANCE : 1e-06
> Info: MAX ITERATIONS : 100
> Info: RIGID WATER USING SETTLE ALGORITHM
> Info: RANDOM NUMBER SEED 15343
> Info: USE HYDROGEN BONDS? NO
> Info: COORDINATE PDB ./complex.r.pdb
> Info: STRUCTURE FILE ./complex.r.psf
> Info: PARAMETER file: CHARMM format!
> Info: PARAMETERS ./par_all27_prot_lipid.prm
> Info: PARAMETERS ./AP1.r.prm
> Info: USING ARITHMETIC MEAN TO COMBINE L-J SIGMA PARAMETERS
> Info: BINARY COORDINATES ./press-02.restart.coor
> Info: SUMMARY OF PARAMETERS:
> Info: 195 BONDS
> Info: 485 ANGLES
> Info: 608 DIHEDRAL
> Info: 54 IMPROPER
> Info: 6 CROSSTERM
> Info: 130 VDW
> Info: 0 VDW_PAIRS
> Info: TIME FOR READING PSF FILE: 2.7509
> Info: TIME FOR READING PDB FILE: 0.902384
> Info:
> Info: Reading from binary file ./press-02.restart.coor
> Info: ****************************
> Info: STRUCTURE SUMMARY:
> Info: 219393 ATOMS
> Info: 156441 BONDS
> Info: 119832 ANGLES
> Info: 82263 DIHEDRALS
> Info: 3446 IMPROPERS
> Info: 1264 CROSSTERMS
> Info: 0 EXCLUSIONS
> Info: 188277 RIGID BONDS
> Info: 469902 DEGREES OF FREEDOM
> Info: 77450 HYDROGEN GROUPS
> Info: 4 ATOMS IN LARGEST HYDROGEN GROUP
> Info: 77450 MIGRATION GROUPS
> Info: 4 ATOMS IN LARGEST MIGRATION GROUP
> Info: TOTAL MASS = 1.34691e+06 amu
> Info: TOTAL CHARGE = -0.000173653 e
> Info: MASS DENSITY = 1.03392 g/cm^3
> Info: ATOM DENSITY = 0.101417 atoms/A^3
> Info: *****************************
> Info:
> Info: Entering startup at 3.82823 s, 57.3364 MB of memory in use
> Info: Startup phase 0 took 0.192449 s, 57.336 MB of memory in use
> Info: Startup phase 1 took 1.98068 s, 95.3553 MB of memory in use
> Info: Startup phase 2 took 0.056417 s, 97.0344 MB of memory in use
> Info: Startup phase 3 took 0.0479629 s, 97.0343 MB of memory in use
> Info: PATCH GRID IS 6 (PERIODIC) BY 6 (PERIODIC) BY 7 (PERIODIC)
> Info: PATCH GRID IS 1-AWAY BY 1-AWAY BY 1-AWAY
> Info: Reading from binary file ./press-02.restart.vel
> Info: REMOVING COM VELOCITY 0.00292499 0.0305455 0.0250787
> Info: LARGEST PATCH (51) HAS 943 ATOMS
> Info: Startup phase 4 took 0.242286 s, 124.273 MB of memory in use
> Info: PME using 4 and 4 processors for FFT and reciprocal sum.
> Info: PME GRID LOCATIONS: 0 1 2 3
> Info: PME TRANS LOCATIONS: 0 1 2 3
> Info: Optimizing 4 FFT steps. 1... 2... 3... 4... Done.
> Info: Startup phase 5 took 0.0777831 s, 129.158 MB of memory in use
> Info: Startup phase 6 took 0.124777 s, 109.334 MB of memory in use
> LDB: Central LB being created...
> Info: Startup phase 7 took 1.33913 s, 109.481 MB of memory in use
> Info: CREATING 16022 COMPUTE OBJECTS
> Info: useSync: 1 useProxySync: 0
> Info: NONBONDED TABLE R-SQUARED SPACING: 0.0625
> Info: NONBONDED TABLE SIZE: 769 POINTS
> TCL: RAMD:
> TCL: RAMD: -------------------------------------------------------------------
> TCL: RAMD: Random Acceleration Molecular Dynamics Simulation version 4.1
> TCL: RAMD: -------------------------------------------------------------------
> TCL: RAMD:
> TCL: RAMD: mdSteps = 0
> TCL: RAMD: accel = 0.3
> TCL: RAMD: ramdSeed = 14257
> TCL: RAMD: forceOutFreq = 50
> TCL: RAMD: maxDist = 50
> TCL: RAMD: firstProtAtom = 1
> TCL: RAMD: rMinRamd = 0.02
> TCL: RAMD: debugLevel = 0
> TCL: RAMD: mdStart = no
> TCL: RAMD: ramdSteps = 50
> TCL: RAMD: firstRamdAtom = 219371
> TCL: RAMD: lastRamdAtom = 219393
> TCL: RAMD: lastProtAtom = 20137
> TCL: RAMD: Pure RAMD simulation is performed
> TCL: RAMD:
> TCL: RAMD: Atoms subject to the random acceleration are: 219371 219372
> 219373 219374 219375 219376 219377 219378 219379 219380 219381 219382
> 219383 219384 219385 219386 219387 219388 219389 219390 219391 219392
> 219393
> colvars: ----------------------------------------------------------------------
> colvars: Initializing the collective variables module, version 2011-01-03.
> colvars: # analysis = off [default]
> colvars: # colvarsTrajFrequency = 100
> colvars: # colvarsRestartFrequency = 200
> colvars: # colvarsTrajAppend = off [default]
> colvars: The restart output state file will be
> "./ramd_acc-03.restart.colvars.state".
> colvars: The final output state file will be "./ramd_acc-03.colvars.state".
> colvars: The trajectory file will be "./ramd_acc-03.colvars.traj".
> colvars: ----------------------------------------------------------------------
> colvars: Initializing a new collective variable.
> colvars: # name = colvar1
> colvars: Initializing a new "distance" component.
> colvars: # componentCoeff = 1 [default]
> colvars: # componentExp = 1 [default]
> colvars: # period = 0 [default]
> colvars: # wrapAround = 0 [default]
> colvars: # forceNoPBC = off [default]
> colvars: # oneSiteSystemForce = off [default]
> colvars: Initializing atom group "group1".
> colvars: Atom group "group1" defined, 1 initialized: total mass = 35.45.
> colvars: Initializing atom group "group2".
> colvars: Atom group "group2" defined, 2 initialized: total mass = 28.014.
> colvars: All components initialized.
> colvars: # width = 1
> colvars: # lowerBoundary = -10
> colvars: # lowerWallConstant = 100
> colvars: # lowerWall = -5
> colvars: # upperBoundary = 10
> colvars: # upperWallConstant = 100
> colvars: # upperWall = 5
> colvars: # expandBoundaries = off [default]
> colvars: # extendedLagrangian = off [default]
> colvars: # outputValue = on [default]
> colvars: # outputVelocity = off [default]
> colvars: # outputSystemForce = on
> colvars: # outputAppliedForce = off [default]
> colvars: ----------------------------------------------------------------------
> ..................................................
> ..................................................
> colvars: # name = colvar9
> colvars: Initializing a new "distance" component.
> colvars: # componentCoeff = 1 [default]
> colvars: # componentExp = 1 [default]
> colvars: # period = 0 [default]
> colvars: # wrapAround = 0 [default]
> colvars: # forceNoPBC = off [default]
> colvars: # oneSiteSystemForce = off [default]
> colvars: Initializing atom group "group1".
> colvars: Atom group "group1" defined, 1 initialized: total mass = 35.45.
> colvars: Initializing atom group "group2".
> colvars: Atom group "group2" defined, 1 initialized: total mass = 15.999.
> colvars: All components initialized.
> colvars: # width = 1
> colvars: # lowerBoundary = -10
> colvars: # lowerWallConstant = 100
> colvars: # lowerWall = -5
> colvars: # upperBoundary = 10
> colvars: # upperWallConstant = 100
> colvars: # upperWall = 5
> colvars: # expandBoundaries = off [default]
> colvars: # extendedLagrangian = off [default]
> colvars: # outputValue = on [default]
> colvars: # outputVelocity = off [default]
> colvars: # outputSystemForce = on
> colvars: # outputAppliedForce = off [default]
> colvars: ----------------------------------------------------------------------
> colvars: Collective variables initialized, 9 in total.
> colvars: ----------------------------------------------------------------------
> colvars: Initializing a new "harmonic" instance.
> colvars: # name = CLA_harmonic
> colvars: # colvars = { colvar1, colvar2, colvar3, colvar4, colvar5,
> colvar6, colvar7, colvar8, colvar9 }
> colvars: # forceConstant = 5
> colvars: # centers = { 3.6, 4.9, 3, 3.2, 3, 4, 3.5, 3.1, 5 }
> colvars: # targetCenters = { 3.6, 4.9, 3, 3.2, 3, 4, 3.5, 3.1, 5 } [default]
> colvars: # targetForceConstant = 0 [default]
> colvars: ----------------------------------------------------------------------
> colvars: Collective variables biases initialized, 1 in total.
> colvars: ----------------------------------------------------------------------
> colvars: Restarting from file "./press-02.restart.colvars.state".
> colvars: Restarting collective variable "colvar1" from value: 3.48985
> colvars: Restarting collective variable "colvar2" from value: 4.50437
> colvars: Restarting collective variable "colvar3" from value: 3.32463
> colvars: Restarting collective variable "colvar4" from value: 3.04155
> colvars: Restarting collective variable "colvar5" from value: 2.99285
> colvars: Restarting collective variable "colvar6" from value: 3.1526
> colvars: Restarting collective variable "colvar7" from value: 3.20845
> colvars: Restarting collective variable "colvar8" from value: 3.11906
> colvars: Restarting collective variable "colvar9" from value: 5.10255
> colvars: Restarting harmonic bias "CLA_harmonic".
> colvars: ----------------------------------------------------------------------
> colvars: Collective variables module initialized.
> colvars: ----------------------------------------------------------------------
> Info: Startup phase 8 took 0.171692 s, 125.178 MB of memory in use
> Info: Startup phase 9 took 0.000221014 s, 133.126 MB of memory in use
> Info: Finished startup at 8.06162 s, 133.126 MB of memory in use
>
> FATAL ERROR: Due to a design error, GlobalMasterServer does not
> support individual atom requests from multiple global force clients on
> parallel runs.
>
>
>
> If the colvar module in the conf file is commented out, ramd runs fine
> and the final pdb (coor) file is also obtained. Conversely, without
> ramd, colvars run fine. It seems to be a matter for developers unless
> (hopefully) it is simply mishandling from my side. Attached hehe
> please find the whole conf file (with colvars commented out), which is
> where errors from my side could be detected.
>
> # Job
> # System mod21-popc-AP1 TIP3-water-solvated, NaCl isotonic
> # ramd with CLA-ligands colvars, PME following press-02.conf
>
> # forcefield
> paratypecharmm on
> parameters ./par_all27_prot_lipid.prm
> parameters ./AP1.r.prm
>
>
> # molecules
> structure ./complex.r.psf
> coordinates ./complex.r.pdb
> bincoordinates ./press-02.restart.coor
> binvelocities ./press-02.restart.vel
> extendedSystem ./press-02.restart.xsc
>
> # constraints
> # colvars on # provide a colvarsConfig (and colvarsINput if restart)
> # colvarsConfig ./CLA-ligands_colvars_press01.in
> # colvarsInput ./press-02.restart.colvars.state
>
>
> # constant temp control
> langevin on
> langevinTemp 310
> langevinDamping 5
> langevinHydrogen off
>
> # constant pressure control
> useGroupPressure yes
> useFlexibleCell yes
> useConstantRatio no
> LangevinPiston on
> LangevinPistonTarget 1.01325
> LangevinPistonPeriod 200.
> LangevinPistonDecay 100.
> LangevinPistonTemp 310
> UseConstantArea yes
>
> # integrator
> timestep 1.0 ;# 1fs/step
> nonbondedFreq 1 ;# nonbonded forces every step
> fullElectFrequency 5 ;# PME only every five step
> stepspercycle 20 ;# redo pairlist every 20 steps
>
> # Approximations_1
> # rigidBonds all ;# needed for 2fs/step
> rigidBonds water
> rigidTolerance 0.000001
> exclude scaled1-4
> cutoff 12.0
> switching on
> switchdist 10.0
> pairlistdist 13.5 ;# cutoff +3.5
> margin 3
> 1-4scaling 1.0
> PME yes
> # cellBasisVector1 122.1 0. 0.
> # cellBasisVector2 0. 119.2 0.
> # cellBasisVector3 0. 0. 153.99
> # Don't set the periodic cell basis if you have also specified an .xsc
> # restart file
> cellOrigin -0.043044526129961014 -1.3494617938995361 -7.98848533630371
> PMEGridSpacing 1.0
>
>
> # output
> outputName ./ramd_acc-03
> outputEnergies 10 # multiple of fullElectFrequency or viceversa
> restartfreq 10
> DCDfreq 100
> binaryrestart yes
> binaryoutput no
> wrapNearest on
> wrapAll on
>
>
> #############################################################
> #*** Random Acceleration Molecular Dynamics
> *************************************
>
> source /usr/local/namd_2.8nb/lib/ramd/scripts/ramd-4.1.tcl
> #*** sources the wrapper script ramd-4.1.tcl;
> #*** please change the directory '../scripts/' to '' ( the correct path );
> #*** directory '' should contain the scripts: ramd-4.1.tcl,
> ramd-4.1_script.tcl, and vectors.tcl
>
> ramd debugLevel 0
> #*** activates verbose output if set to something else than 0
>
> ramd ramdSteps 50
> #*** specifies the number of steps in 1 ramd stint;
> #*** defaults to 50
>
> ramd accel 0.3
> #*** specifies the acceleration to be applied;
> #*** defaults to 0.25 kcal/mol*A*amu
>
> ramd rMinRamd 0.02
> #*** specifies the minimum distance to be travelled by the ligand in 1
> ramd stint;
> #*** defaults to 0.01 Angstr
>
> ramd forceOutFreq 50
> #*** every 'forceOutFreq' steps detailed output of forces will be written;
> #*** defaults to 0 (no detailed output)
>
> ramd maxDist 50
> #*** specifies the distance between the COMs of the ligand and the
> protein when the simulation is stopped
> #*** defaults to 50 Angstr (myoglobin max with 43 A)
>
> ramd firstProtAtom 1
> #*** specifies the index of the first protein atom
> #*** defaults to 1 (assumes first atom in the system corresponds to
> first protein atom
>
> ramd lastProtAtom 20137
> #*** specifies the index of the last protein atom
> #*** required; simulation exits if this parameter is not set
>
> ramd firstRamdAtom 219371
> #*** specifies the index of the first ligand atom
> #*** required; simulation exits if this parameter is not set
>
> ramd lastRamdAtom 219393
> #*** specifies the index of the last ligand atom
> #*** required; simulation exits if this parameter is not set
>
> ramd ramdSeed 14257
> #*** specifies the seed for the random number generator (for the
> generation of acceleration directions)
> #*** defaults to 14253
> #*** please change if you wish to run different trajectories
> ######################
> # Run protocol (steps multiple of stepspercycle)
> seed 15343
> numsteps 100
>
>
> thanks
> francesco
>
>
> On Mon, Mar 28, 2011 at 11:52 PM, Jim Phillips<jim_at_ks.uiuc.edu> wrote:
>> You only want the tclForce on that's in the ramd.tcl file. You need to look
>> at the RAMD example since I'm not familiar with how it works.
>>
>> -Jim
>>
>>
>> On Mon, 28 Mar 2011, Francesco Pietra wrote:
>>
>>> On Mon, Mar 28, 2011 at 9:33 PM, Jim Phillips<jim_at_ks.uiuc.edu> wrote:
>>>> It looks like you have "tclForces on" after run or minimize in your
>>>> config
>>>> file, so possibly your tclForces script was not active during the RAMD
>>>> run
>>>> at all.
>>> Hi Jim:
>>> Yes, "tclforces on" is in the sourced ramd script, which comes after
>>> the "run" command in namd conf file. I had checked that against namd
>>> 2.7 with external namd: there was the same order of statements.
>>>
>>> Now, if "tclforces on" is located in namd conf file, at the top:
>>>
>>>
>>> # forcefield
>>> paratypecharmm on
>>> parameters ./par_all27_prot_lipid.prm
>>> parameters ./AP1.r.prm
>>> tclForces on
>>>
>>> the error is:
>>>
>>> Info: Configuration file is ramd_acc-03.conf
>>> Info: Working in the current directory
>>> /home/francesco/work_hasic1a_GLUP_ASPP_allHSP_S-S/9.ramd_04-1.1
>>> TCL: Suspending until startup complete.
>>> ERROR: 'tclForcesScript' is a required configuration option
>>> ERROR: when 'tclForces' is set
>>> ERROR: tclForcesScript defines: Tcl script for global forces
>>> FATAL ERROR: ERROR(S) IN THE CONFIGURATION FILE
>>>
>>>
>>>
>>> Where should "tclforces on" be now located?
>>>
>>> thanks
>>> francesco
>>>
>>>
>>>
>>>> -Jim
>>>>
>>>> On Mon, 28 Mar 2011, Francesco Pietra wrote:
>>>>
>>>>> Hi:
>>>>> I carried out a trial RAMD simulation with namd.2.8 (from binary
>>>>> Linux64) charmmff for a protein-organic-ligand in a lipidic membrane.
>>>>>
>>>>> The simulation ended correctly at the bare 100 chosen steps, the
>>>>> ligand did not move - as expected from 1000 steps acc 0.3 ts 1fs.
>>>>> Everything seems OK except the error below. I vaguely remember that
>>>>> Vlad warned me in the past to compile namd in order that ramd exists
>>>>> correctly. Was that the same errior message or anything new about
>>>>> which I should now care?
>>>>>
>>>>> thanks
>>>>> francesco pietra
>>>>>
>>>>> **************************
>>>>>
>>>>> ENERGY: 1000 10573.6218 17958.2104 9301.2411
>>>>> 779.3272 -763083.0117 55270.6864 0.0000
>>>>> 1.9652 144208.1552 -524989.8044 308.8680
>>>>> -669197.9596 -524727.0908 309.0692 -173.1971
>>>>> -32.1881 2166358.3052 -35.4825 -31.7813
>>>>>
>>>>> WRITING EXTENDED SYSTEM TO RESTART FILE AT STEP 1000
>>>>> WRITING COORDINATES TO DCD FILE AT STEP 1000
>>>>> WRITING COORDINATES TO RESTART FILE AT STEP 1000
>>>>> FINISHED WRITING RESTART COORDINATES
>>>>> The last position output (seq=1000) takes 0.134 seconds, 744.610 MB of
>>>>> memory in use
>>>>> WRITING VELOCITIES TO RESTART FILE AT STEP 1000
>>>>> FINISHED WRITING RESTART VELOCITIES
>>>>> The last velocity output (seq=1000) takes 0.021 seconds, 744.609 MB of
>>>>> memory in use
>>>>> TCL: Setting parameter TclForces to on
>>>>> FATAL ERROR: Setting parameter TclForces from script failed!
>>>>>
>>>>> FATAL ERROR: Setting parameter TclForces from script failed!
>>>>>
>>>>> ------------- Processor 1 Exiting: Called CmiAbort ------------
>>>>> Reason: FATAL ERROR: Setting parameter TclForces from script failed!
>>>>>
>>>>>
>>>>> FATAL ERROR: Setting parameter TclForces from script failed!
>>>>>
>>>>> ------------- Processor 3 Exiting: Called CmiAbort ------------
>>>>> Reason: FATAL ERROR: Setting parameter TclForces from script failed!
>>>>>
>>>>>
>>>>> ------------- Processor 0 Exiting: Called CmiAbort ------------
>>>>> Reason: FATAL ERROR: Setting parameter TclForces from script failed!
>>>>>
>>>>>
>>>>> FATAL ERROR: Setting parameter TclForces from script failed!
>>>>>
>>>>> ------------- Processor 2 Exiting: Called CmiAbort ------------
>>>>> Reason: FATAL ERROR: Setting parameter TclForces from script failed!
>>>>> **************************
>>>>>

-- 
Dr. Vlad Cojocaru
Max Planck Institute for Molecular Biomedicine
Department of Cellular and Developmental Biology
Roentgenstrasse 20
48149 Muenster, Germany
tel: +49-251-70365-324
fax: +49-251-70365-399
email: vlad.cojocaru[at]mpi-muenster.mpg.de

This archive was generated by hypermail 2.1.6 : Mon Dec 31 2012 - 23:20:02 CST