[gmx-users] Can you please tell me the exact command? Re: Restart a stopped simulation after a power failure? : Not working in Gromacs 5.1.4
Mark Abraham
mark.j.abraham at gmail.com
Thu May 18 17:38:38 CEST 2017
Hi,
Use
gmx mdrun
gmx mdrun -cpi
so that by default you get appending and the filenames haven't changed.
Later, when you're comfortable, try
gmx mdrun -deffnm whatever
gmx mdrun -deffnm whatever -cpi
but really there's very few cases where deffnm is actually useful. You have
a directory name for organizing your work - do only one simulation in it
and avoid problems from trying to name each file descriptively.
Your problems likely stem from trying to do
gmx mdrun -deffnm whatever
gmx mdrun -s whatever -cpi
but it's hard to guess what you've done from partial reports. If so, you
can do it better by doing
gmx mdrun -deffnm whatever -cpi
Now, you probably need to do
gmx trjcat -f whatever traj -o combined
to get all your output into one file
Hope this helps!
Mark
On Thu, May 18, 2017 at 5:31 PM Adarsh V. K. <adarsh_p130085bt at nitc.ac.in>
wrote:
> Dear Dr. Mark,
>
> Thank you for the prompt reply. I have been struggling with the issue for
> couple month time. Would you please advise me how to use " -deffnm " flag
> in commands?
>
> # gmx mdrun -s md_0_1.tpr -cpi md_0_1_prev.cpt -v -deffnm
>
> Will the above command will successful restart a stopped simulation and
> append?, If no.... could you please type the exact command I should use to
> solve the issue.
>
> Sincerely,
>
> Adarsh V. K.
> Research Scholar (PhD)
> National Institute of Technology, Calicut ( NIT C )
> Kerala state, India
>
> --------------------------------------------------------------------------
> On Thu, May 18, 2017 at 7:13 PM, Mark Abraham <mark.j.abraham at gmail.com>
> wrote:
>
> > Hi,
> >
> > Some of your commands and outputs are inconsistent. You can't get
> > md_0_1.log without e.g. using -deffnm, and you'll only get md.log if you
> > don't use deffnm. But if you change that when trying to do your
> > continuation, mdrun won't know how to append because you haven't told it
> > that the old name. Be consistent. Your whole trajectory is present in all
> > your output files, but you're going to have to concatenate them yourself,
> > e.g. with gmx trjcat since you didn't tell mdrun a consistent name for
> the
> > files.
> >
> > GROMACS 2016 likely refuses to do your continuation, precisely so that
> you
> > are prompted to say what you actually want, rather than mdrun trying to
> be
> > clever and getting it wrong.
> >
> > Mark
> >
> > On Thu, May 18, 2017 at 11:03 AM Adarsh V. K. <
> adarsh_p130085bt at nitc.ac.in>
> > wrote:
> >
> >> Dear gmx users,
> >>
> >> I use Gromacs 5.1.4 for protein ligand simulation. How can I restart a
> >> stopped simulation after a power failure? (used 8 processor cores +
> >> Graphics card GTX 780Ti).
> >>
> >> I have attached the log file along with this mail.
> >>
> >> It appeared that no details appended after restarting simulation (after
> >> power failure at 3.5 ns), Even Though the 'terminal' Ubuntu 16.04 shows
> >> that the simulation successfully restarted from 3.5 ns (used 8 processor
> >> cores + Graphics card GTX 780Ti) and completed the entire 8 ns
> >> simulation...!!. Frames...350 -> 3500 ps (not appended the log file
> after
> >> the restart from 3500ps to 8000ps).
> >>
> >> I used the command (gromacs 5.1.4),
> >> # gmx mdrun -s md_0_1.tpr -cpi md_0_1_prev.cpt -v
> >>
> >> I also tried, # gmx mdrun -s md_0_1.tpr -cpi md_0_1_prev.cpt -append
> >>
> >> ------------------------------------------------------------
> >> ------------------------------------
> >> md_0_1.log : (final)
> >> -----------------
> >> Step Time Lambda
> >> 1765000 3530.00000 0.00000
> >>
> >> Energies (kJ/mol)
> >> G96Angle Proper Dih. Improper Dih. LJ-14
> >> Coulomb-14
> >> 6.73096e+03 3.13198e+03 2.15931e+03 1.31375e+03
> >> 4.51549e+04
> >> LJ (SR) Disper. corr. Coulomb (SR) Coul. recip.
> >> Potential
> >> 3.66659e+05 -7.80170e+03 -2.71529e+06 4.12416e+03
> >> -2.29382e+06
> >> Kinetic En. Total Energy Temperature Pres. DC (bar) Pressure
> >> (bar)
> >> 4.12474e+05 -1.88135e+06 3.00209e+02 -7.62778e+01
> >> 2.10104e+01
> >> Constr. rmsd
> >> 2.37431e-05
> >> ------------------------------------------------------------
> >> -------------------------------------
> >> md.log (final)
> >> ----------
> >> Step Time Lambda
> >> 4000000 8000.00000 0.00000
> >>
> >> Writing checkpoint, step 4000000 at Thu May 18 13:55:01 2017
> >>
> >>
> >> Energies (kJ/mol)
> >> G96Angle Proper Dih. Improper Dih. LJ-14
> >> Coulomb-14
> >> 6.44959e+03 3.18212e+03 2.06343e+03 1.53953e+03
> >> 4.49878e+04
> >> LJ (SR) Disper. corr. Coulomb (SR) Coul. recip.
> >> Potential
> >> 3.66131e+05 -7.78489e+03 -2.71259e+06 4.78708e+03
> >> -2.29124e+06
> >> Kinetic En. Total Energy Temperature Pres. DC (bar) Pressure
> >> (bar)
> >> 4.13852e+05 -1.87738e+06 3.01212e+02 -7.59495e+01
> >> 5.34060e+01
> >> Constr. rmsd
> >> 2.33151e-05
> >> ---------
> >> Force evaluation time GPU/CPU: 12.005 ms/12.263 ms = 0.979
> >> For optimal performance this ratio should be close to 1!
> >>
> >> Core t (s) Wall t (s) (%)
> >> Time: 312346.797 39374.380 793.3
> >> 10h56:14
> >> (ns/day) (hour/ns)
> >> Performance: 9.967 2.408
> >> Finished mdrun on rank 0 Thu May 18 13:55:02 2017
> >>
> >> ------------------------------------------------------------
> >> -----------------------------------
> >> But while using the command it shows only 3.5 ns of simulation. How to
> >> solve this problem
> >>
> >> gmx trjconv -s md_0_1.tpr -f md_0_1.xtc -o md_0_1_noPBC.xtc -pbc mol -ur
> >> compact
> >>
> >> Select 0 ("System")
> >> ------------------------------------------------------------
> >> -------------------------------------
> >>
> >> Later I repeated the whole 8 ns simulation from beginning.... It
> >> completed with out any issues.
> >>
> >> Can any body tell me what happens during the power failure and why
> >> re-starting a stopped simulation is not working.
> >>
> >> Is it any modification in 're-start' command is required?
> >>
> >> Regards,
> >> Adarsh V. K.
> >>
> >>
> --
> Gromacs Users mailing list
>
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> send a mail to gmx-users-request at gromacs.org.
>
More information about the gromacs.org_gmx-users
mailing list