[gmx-users] rerun of trajectory!

[Paul Barrett]

It looks right...

I use 

  tpbconv -s post_pr.tpr -f main.trr -o main_restart

  nohup mdrun -v  -s main_restart -o main_restart -x main_restart -c
main_restart -g main_restartlog >& main_restart.job

which works fine

Im sure you have already checked the following:
1) are you executing the command from the same working directory as
traj.trr?
2) you did call the previous run traj.trr and not traj.xtc or traj.trj
or somethingelse.trr?
3) if traj.trr is in the directory where you are running it try

 ls -l traj.trr 

and see if it has a non-zero file size.

If all those are fine I cant think what the problem is.


sorry

Paul Barrett

>dear gromacs users,
>
>i'm trying to simulate a protein solvated in water.I ran the simulation
>for 1ns.This was successfully done.Now i want to continue the MD for
>another 9ns with the trr file of the 1ns simulation.
>So i tried making a continuation input file using tpbconv.i.e
>
>>tpbconv -s b4md_emsd_saksol.tpr -f traj.trr -n md_emsd_saksol.ndx -o
>traj_out.tpr
>
>But the command does not seem to work. My .trr file does not seem to be
>getting recognized. i'm getting the following error after running the tpbconv
>command:
>
>Reading toplogy and shit from b4md_emsd_saksol.tpr
>Reading file b4md_emsd_saksol.tpr, VERSION 3.1.4 (single precision)
>Reading file b4md_emsd_saksol.tpr, VERSION 3.1.4 (single precision)
>Modifying ir->bUncStart to TRUE
>
>READING COORDS, VELS AND BOX FROM TRAJECTORY traj.trr...
>
>Fatal error: File traj.trr not found
>
>Do i have to give any other option in the tpbconv command.
>Can anyone please tell me how i can make a continutaion input file for
>further rerun of the MD.
>Or if i take the output gro file created from the previous 1ns simulation
>and use it as input for the grompp command for the 9ns simulation, will it be
>equivalent to rerunning from where i left off. Or will the velocities and
>energies calculated in the trr file of the 1ns simulation change now as
>fresh initail velocities will be assigned for the new run.
>
>will really appreciate if someone can tell me how to go about it and what
>exactly can be done to solve the problem.
>
>Thanks & Regards ,
>
>Swetha Vijayakrishnan
>
>
>
>
>
>
>
>
>*******************************************************************************
>
>
>Swetha Vijayakrishnan
>
>Project Assistant
>C/O Dr.Shekhar.C.Mande
>Centre for DNA Fingerprinting and Diagnostics(CDFD),
>ECIL Road,Nacharam,
>Hyderabad - 500076
>
>Phone:: Direct: 91-040-7171442
>        Indirect: 91-040-7151344/46/47/56 (lab extn 1400)
>
>Fax  :: 91-040-(7155479-7155610>
>
>Email:: swetha at cdfd.org.in
>        swethasenorita at yahoo.com
>
>
>_______________________________________________
>gmx-users mailing list
>gmx-users at gromacs.org
>http://www.gromacs.org/mailman/listinfo/gmx-users
>Please don't post (un)subscribe requests to the list. Use the 
>www interface or send it to gmx-users-request at gromacs.org.
>--------------------------------------------------------------------------------
>Return-Path: <gmx-users-admin at gromacs.org>
>Received: from oxmail.ox.ac.uk by nutmeg.biop.ox.ac.uk (MX V4.2 AXP) with SMTP;
>          Wed, 14 May 2003 11:29:59 +0100
>Received: from hawk.theophys.kth.se ([130.237.25.19]) by oxmail.ox.ac.uk with
>          smtp (Exim 4.14) id 19FtWJ-0002jm-Dw for barrett at biop.ox.ac.uk; Wed,
>          14 May 2003 11:29:59 +0100
>Received: (qmail 20169 invoked from network); 14 May 2003 10:30:02 -0000
>Received: from localhost (HELO hawk.theophys.kth.se) (127.0.0.1) by 0 with
>          SMTP; 14 May 2003 10:30:02 -0000
>Received: (qmail 20124 invoked from network); 14 May 2003 10:29:15 -0000
>Received: from unknown (HELO www.cdfd.org.in) (210.212.212.2) by 0 with SMTP;
>          14 May 2003 10:29:15 -0000
>Received: from localhost (swetha at localhost) by www.cdfd.org.in
>          (SGI-8.9.3/8.9.3) with ESMTP id PAA84395 for <gmx-users at gromacs.org>;
>          Wed, 14 May 2003 15:59:57 -0700 (PDT)
>Message-ID: <Pine.SGI.4.44.0305141523160.13923218-100000 at www.cdfd.org.in>