Hi Una,

Do the xaver.pdb and the md.trr correspond in terms of atom numbers?
The JACOBI error has to do with the fitting. The 'broken molecule'
thing is a warning, issued because you give a .pdb file which does not
contain information regarding connectivity, which coulod otherwise be
used to make broken molecules whole. But you shouldn't have broken
molecules if you didn't run trjconv on the trajectory with certain the
-pbc options.

Cheers,

Tsjerk

On 5/29/07, Una Bjarnadottir <[EMAIL PROTECTED]> wrote:
Dear users,

I'm using g_covar to analyse my trajectories and saw here on the list
that the one should use the average structure as the reference structre
in the -s option in g_covar.

I made my average structure with g_rmsf but when I ran g_covar

g_covar -f md.trr -s xaver.pdb -ascii covar.dat -ref

I got these comments and error:

Note: the fit and analysis group are identical,
      while the fit is mass weighted and the analysis is not.
      Making the fit non mass weighted.


Warning: can not make broken molecules whole without a run input file,
         don't worry, mdrun doesn't write broken molecules

Calculating the average structure ...
trn version: GMX_trn_file (single precision)
Reading frame       0 time    0.000   Fatal error: Error: Too many
iterations in routine JACOBI

Why does it say "broken molecule" when I have full trajectory and an
average structure what am I doing wrong to cause the error?

Many thanks to all commets.

All my best, Una

_______________________________________________
gmx-users mailing list    [email protected]
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php



--
Tsjerk A. Wassenaar, Ph.D.
Junior UD (post-doc)
Biomolecular NMR, Bijvoet Center
Utrecht University
Padualaan 8
3584 CH Utrecht
The Netherlands
P: +31-30-2539931
F: +31-30-2537623
_______________________________________________
gmx-users mailing list    [email protected]
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Reply via email to