<br><br><div class="gmail_quote">On Sun, Nov 21, 2010 at 2:05 PM, Spitaleri Andrea <span dir="ltr"><<a href="mailto:spitaleri.andrea@hsr.it">spitaleri.andrea@hsr.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hi,<br>
yes sure. Basically I do:<br>
<br>
1. mdrun -s runA_ -deffnm runA_ -replex 5000 -multi 5 -> i get the first 5x25ns of remd simulation (five xtc every 5ps and five trr every 20ps, for each replica). I check those files by gmxcheck and they are fine. no errors.<br>
<br>
2. for i in 'seq 1 5'; tpbconv -s runA_$i -nsteps 25000000 -o runB_$i -> extension the simulation to 50ns total<br>
<br>
3. mdrun -s runB_ -replex 5000 -multi 5 -deffnm runB_ -cpi runA_ -> at end some (2 or 1 on the 5 xtc file) of the xtc are corrupted (from gmxcheck) whereas the trr are fine. These are from 25ns to 50ns.<br>
<br>
>From the log file I do not see any errors. Everything seems fine. I have free room space in the hd too :)<br>
<br>
I am just wondering whether the problem is in the xtc options (precision and writing step)<br></blockquote><div><br></div><div>I doubt it that it has anything to do with your xtc options.</div><div><br></div><div>Are all you xtc corrupted or only some? Are those which are corrupted all corrupted on the same frame or different ones?</div>
<div><br></div><div>Roland</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div class="im">
<br>
------------------------------------------------------------<br>
Andrea Spitaleri PhD<br>
Dulbecco Telethon Institute<br>
Center of Genomics, BioInformatics and BioStatistics<br>
c/o DIBIT Scientific Institute<br>
Biomolecular NMR, 1B4<br>
Via Olgettina 58<br>
20132 Milano (Italy)<br>
<a href="http://sites.google.com/site/andreaspitaleri/" target="_blank">http://sites.google.com/site/andreaspitaleri/</a><br>
Tel: 0039-0226434348/5622/3497/4922<br>
Fax: 0039-0226434153<br>
------------------------------------------------------------<br>
</div>________________________________________<br>
Da: <a href="mailto:gmx-users-bounces@gromacs.org">gmx-users-bounces@gromacs.org</a> [<a href="mailto:gmx-users-bounces@gromacs.org">gmx-users-bounces@gromacs.org</a>] per conto di Mark Abraham [<a href="mailto:Mark.Abraham@anu.edu.au">Mark.Abraham@anu.edu.au</a>]<br>
Inviato: domenica 21 novembre 2010 17.03<br>
A: Discussion list for GROMACS users<br>
Oggetto: Re: [gmx-users] xtc corrupted during REMD<br>
<div><div></div><div class="h5"><br>
On 21/11/2010 2:34 AM, Spitaleri Andrea wrote:<br>
> Hi there,<br>
> I am encountering a weird problem with a REMD simulation using 4.5.3. The total simulation is 50ns with 5 replica, and I do in two runs: 25ns and then continuing to 50ns (walltime queue). The first run is okay, the continue run (the last 25ns) randomly make some xtc files corrupted (from gmxcheck I get the Magic Number Error).<br>
<br>
I don't understand how the simulation can continue writing the .xtc<br>
files when you are getting magic number errors from gmxcheck. We need to<br>
see command lines for your workflow, please :-)<br>
<br>
Mark<br>
<br>
> It is strange since the respective trr files are okay and the simulation is still going (it is not blowing up from the log, not step.pdb files, not crash). The only difference is that I am writing the xtc often respect to the trr file and just the complex not the solvent:<br>
><br>
> nstxout = 10000 ; coordinates every 20ps<br>
> nstvout = 0 ; velocity every 0ps<br>
> nstfout = 0 ; forces every 0 ps<br>
> nstlog = 2500 ; energies log every 5ps<br>
> nstenergy = 2500 ; energies every 5ps<br>
> nstxtcout = 2500 ; coordinates every 5ps to xtc<br>
> xtc-precision = 2500 ;<br>
> xtc-grps = complex;<br>
><br>
><br>
> Since the error is happening only for the continuing run, I am just wondering if there is any reason for this.<br>
><br>
> thanks for any help<br>
><br>
><br>
> and<br>
><br>
> ------------------------------------------------------------<br>
> Andrea Spitaleri PhD<br>
> Dulbecco Telethon Institute<br>
> Center of Genomics, BioInformatics and BioStatistics<br>
> c/o DIBIT Scientific Institute<br>
> Biomolecular NMR, 1B4<br>
> Via Olgettina 58<br>
> 20132 Milano (Italy)<br>
> <a href="http://sites.google.com/site/andreaspitaleri/" target="_blank">http://sites.google.com/site/andreaspitaleri/</a><br>
> Tel: 0039-0226434348/5622/3497/4922<br>
> Fax: 0039-0226434153<br>
> ------------------------------------------------------------<br>
><br>
><br>
> -----------------------------------------------------------------------------------<br>
> SOSTIENI ANCHE TU LA RICERCA DEL SAN RAFFAELE.<br>
> NON C'E' CURA SENZA RICERCA.<br>
> Per donazioni: ccp 42437681 intestato a Fondazione Arete' Onlus del San Raffaele.<br>
> Per informazioni: tel. 02.2643.4461 - <a href="http://www.sanraffaele.org" target="_blank">www.sanraffaele.org</a><br>
<br>
--<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the<br>
www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
<br>
-----------------------------------------------------------------------------------<br>
SOSTIENI ANCHE TU LA RICERCA DEL SAN RAFFAELE.<br>
NON C'E' CURA SENZA RICERCA.<br>
Per donazioni: ccp 42437681 intestato a Fondazione Arete' Onlus del San Raffaele.<br>
Per informazioni: tel. 02.2643.4461 - <a href="http://www.sanraffaele.org" target="_blank">www.sanraffaele.org</a><br>
--<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the<br>
www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
<br>
<br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>ORNL/UT Center for Molecular Biophysics <a href="http://cmb.ornl.gov">cmb.ornl.gov</a><br>865-241-1537, ORNL PO BOX 2008 MS6309<br>