Thank you again for your time and help. Performing "rerun" on the original system passes without warnings with the following output:
---------------------------------------------------------------------------------------------------------------------------- Option Filename Type Description ------------------------------------------------------------ -s inp/dyn/dppc_dupc_chol_034x051x015_295K.tpr Input Run input file: tpr tpb tpa -o /scratch/tmp/dhako_01/projects/_rafts_cg/crd/dynmc_restart/dppc_dupc_chol_034x051x015_295K_testrerun.trr Output Full precision trajectory: trr trj cpt -x traj.xtc Output, Opt. Compressed trajectory (portable xdr format) -cpi state.cpt Input, Opt. Checkpoint file -cpo state.cpt Output, Opt. Checkpoint file -c crd/dynmc_end/dppc_dupc_chol_034x051x015_295K_testrerun.gro Output Structure file: gro g96 pdb etc. -e dat/enr/dppc_dupc_chol_034x051x015_295K_testrerun.edr Output Energy file -g out/dppc_dupc_chol_034x051x015_295K_testrerun.log Output Log file -dhdl dhdl.xvg Output, Opt. xvgr/xmgr file -field field.xvg Output, Opt. xvgr/xmgr file -table table.xvg Input, Opt. xvgr/xmgr file -tablep tablep.xvg Input, Opt. xvgr/xmgr file -tableb table.xvg Input, Opt. xvgr/xmgr file -rerun /scratch/tmp/dhako_01/projects/_rafts_cg/crd/dynmc_restart/dppc_dupc_chol_034x051x015_295K.trr Input, Opt! Trajectory: xtc trr trj gro g96 pdb cpt -tpi tpi.xvg Output, Opt. xvgr/xmgr file -tpid tpidist.xvg Output, Opt. xvgr/xmgr file -ei sam.edi Input, Opt. ED sampling input -eo sam.edo Output, Opt. ED sampling output -j wham.gct Input, Opt. General coupling stuff -jo bam.gct Output, Opt. General coupling stuff -ffout gct.xvg Output, Opt. xvgr/xmgr file -devout deviatie.xvg Output, Opt. xvgr/xmgr file -runav runaver.xvg Output, Opt. xvgr/xmgr file -px pullx.xvg Output, Opt. xvgr/xmgr file -pf pullf.xvg Output, Opt. xvgr/xmgr file -mtx nm.mtx Output, Opt. Hessian matrix -dn dipole.ndx Output, Opt. Index file Option Type Value Description ------------------------------------------------------ -[no]h bool no Print help info and quit -[no]version bool no Print version info and quit -nice int 0 Set the nicelevel -deffnm string Set the default filename for all file options -xvg enum xmgrace xvg plot formatting: xmgrace, xmgr or none -[no]pd bool no Use particle decompostion -dd vector 0 0 0 Domain decomposition grid, 0 is optimize -nt int 0 Number of threads to start (0 is guess) -npme int -1 Number of separate nodes to be used for PME, -1 is guess -ddorder enum interleave DD node order: interleave, pp_pme or cartesian -[no]ddcheck bool yes Check for all bonded interactions with DD -rdd real 0 The maximum distance for bonded interactions with DD (nm), 0 is determine from initial coordinates -rcon real 0 Maximum distance for P-LINCS (nm), 0 is estimate -dlb enum auto Dynamic load balancing (with DD): auto, no or yes -dds real 0.8 Minimum allowed dlb scaling of the DD cell size -gcom int -1 Global communication frequency -[no]v bool no Be loud and noisy -[no]compact bool yes Write a compact log file -[no]seppot bool no Write separate V and dVdl terms for each interaction type and node to the log file(s) -pforce real -1 Print all forces larger than this (kJ/mol nm) -[no]reprod bool no Try to avoid optimizations that affect binary reproducibility -cpt real 15 Checkpoint interval (minutes) -[no]cpnum bool no Keep and number checkpoint files -[no]append bool yes Append to previous output files when continuing from checkpoint instead of adding the simulation part number to all file names -maxh real -1 Terminate after 0.99 times this time (hours) -multi int 0 Do multiple simulations in parallel -replex int 0 Attempt replica exchange every # steps -reseed int -1 Seed for replica exchange, -1 is generate a seed -[no]ionize bool no Do a simulation including the effect of an X-Ray bombardment on your system Reading file inp/dyn/dppc_dupc_chol_034x051x015_295K.tpr, VERSION 4.5.1 (single precision) Starting 8 threads Making 2D domain decomposition 4 x 2 x 1 WARNING: This run will generate roughly 8779 Mb of data starting md rerun 'DPPC/DUPC/CHOLESTEROL BILAYER', reading coordinates from input trajectory '/scratch/tmp/dhako_01/projects/_rafts_cg/crd/dynmc_restart/dppc_dupc_chol_034x051x015_295K.trr' trn version: GMX_trn_file (single precision) Reading frame 0 time 0.000 Reading frame 1 time 1000.000 Reading frame 2 time 2000.000 Reading frame 3 time 3000.000 Reading frame 4 time 4000.000 Reading frame 5 time 5000.000 Reading frame 6 time 6000.000 Reading frame 7 time 7000.000 Reading frame 8 time 8000.000 Reading frame 9 time 9000.000 Reading frame 10 time 10000.000 Reading frame 11 time 11000.000 Reading frame 12 time 12000.000 Reading frame 13 time 13000.000 Reading frame 14 time 14000.000 Reading frame 15 time 15000.000 Reading frame 16 time 16000.000 Reading frame 17 time 17000.000 Reading frame 18 time 18000.000 Reading frame 19 time 19000.000 Reading frame 20 time 20000.000 Reading frame 30 time 30000.000 Reading frame 40 time 40000.000 Reading frame 50 time 50000.000 Reading frame 60 time 60000.000 Reading frame 70 time 70000.000 Reading frame 80 time 80000.000 Reading frame 90 time 90000.000 Reading frame 100 time 100000.000 Reading frame 110 time 110000.000 Reading frame 120 time 120000.000 Reading frame 130 time 130000.000 Reading frame 140 time 140000.000 Reading frame 150 time 150000.000 Reading frame 160 time 160000.000 Reading frame 170 time 170000.000 Reading frame 180 time 180000.000 Reading frame 190 time 190000.000 Reading frame 200 time 200000.000 Reading frame 300 time 300000.000 Reading frame 400 time 400000.000 Reading frame 500 time 500000.000 Reading frame 600 time 600000.000 Reading frame 700 time 700000.000 Reading frame 800 time 800000.000 Reading frame 900 time 900000.000 Reading frame 1000 time 1000000.000 Reading frame 1100 time 1100000.000 Reading frame 1200 time 1200000.000 Reading frame 1300 time 1300000.000 Reading frame 1400 time 1400000.000 Reading frame 1500 time 1500000.000 Reading frame 1600 time 1600000.000 Reading frame 1700 time 1700000.000 Reading frame 1800 time 1800000.000 Reading frame 1900 time 1900000.000 Reading frame 2000 time 2000000.000 Reading frame 3000 time 3000000.000 Reading frame 4000 time 4000000.000 Reading frame 5000 time 5000000.000 Reading frame 6000 time 6000000.000 Reading frame 7000 time 7000000.000 Reading frame 8000 time 8000000.000 Reading frame 9000 time 9000000.000 Reading frame 10000 time 10000000.000 Reading frame 11000 time 11000000.000 Reading frame 12000 time 12000000.000 Last frame 12000 time 12000000.000 Average load imbalance: 30.8 % Part of the total run time spent waiting due to load imbalance: 1.7 % NOTE: 6 % of the run time was spent communicating energies, you might want to use the -gcom option of mdrun Parallel run - timing based on wallclock. NODE (s) Real (s) (%) Time: 1157.505 1157.505 100.0 19:17 (Mnbf/s) (GFlops) (ns/day) (hour/ns) Performance: 12.993 1.855 895719.837 0.000 ---------------------------------------------------------------------------------------------------------------------------- The subsequent g_energy also runs smoothly. One gets no warning also if passing the ordered *ordered.tpr and the new index (*ordered.ndx) files to the "mdrun -rerun" of the original system. This indeed leaves with the ordered trajectory file as the only problematic point which might cause the warnings. Interestingly, the ordered *.trr file (product of trjorder command) is nearly twice as small as the original *.trr file. Since the coordinates in the ordered *.trr file look fine then probably it is the velocities which might be missing in it. Can the warnings come from this issue ? In this case the g_energy should not be able to calculate the kinetic energy. If this is the case may one still rely on the numbers obtained by g_energy for the coulomb and vdw terms (despite the warning messages) ? Thanks a lot for all the suggestions. > Date: Wed, 24 Jul 2013 07:12:54 -0400 > From: jalem...@vt.edu > To: gmx-users@gromacs.org > Subject: Re: [gmx-users] Calculate interaction energy dynamically > > On 7/24/13 5:49 AM, Davit Hakobyan wrote: > > This is smelling a bit buggy to me, but then too, I've never tried a rerun on > a > reordered trajectory. As a test, can you produce an .edr file from a rerun > of > the original trajectory (i.e. not manipulated by trjorder)? > > -Justin > > -- > ================================================== > > Justin A. Lemkul, Ph.D. > Postdoctoral Fellow > > Department of Pharmaceutical Sciences > School of Pharmacy > Health Sciences Facility II, Room 601 > University of Maryland, Baltimore > 20 Penn St. > Baltimore, MD 21201 > > jalem...@outerbanks.umaryland.edu | (410) 706-7441 > > ================================================== > -- > gmx-users mailing list gmx-users@gromacs.org > http://lists.gromacs.org/mailman/listinfo/gmx-users > * Please search the archive at > http://www.gromacs.org/Support/Mailing_Lists/Search before posting! > * Please don't post (un)subscribe requests to the list. Use the > www interface or send it to gmx-users-requ...@gromacs.org. > * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists -- gmx-users mailing list gmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting! * Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists