Hi all.<br><br>First of all, I would really want want to apologize if I somehow missunderstood any previous message from anyone, and made myself a bit annoying because of this fact.<br><br>Second, I would like to say thanks again to eferybody for all help in these matters.
<br><br><div><span class="gmail_quote">On 10/1/07, <b class="gmail_sendername">Tsjerk Wassenaar</b> <<a href="mailto:tsjerkw@gmail.com" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">tsjerkw@gmail.com
</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">I'd say the easiest is to get into gmx_editconf and look at the option
<br>"-grasp" or "-mead" which put the atomic charge and vdw radius in the<br>b-factor and occupancy field of a .pdb file when given a .tpr file.<br>It's trivial to make it write charges and masses instead. Would also
<br>be a good first exercise for C/Gromacs :p</blockquote><div><br>Good idea! Even I don't choose to sitck with many file type conversions for the sake of the disk space, that specific program now seems to be a good start to learn explicitly how all the topology data is passed from one place to the other. And that would help me to learn how to read the topology files (since it seems to be only a matter now of reading atomic charges and masses).
<br></div><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">I'm not sure what you mean with "shape" and "that integer". There's
<br>nothing much redundant in the .xtc file and there are no integers<br>having to do with box shape. In fact, there is no "shape"...</blockquote><div><br>Well, by "shape" I meant "cubic" and "truncated octahedral", for example. But the integer was my mistake (old dlpoly adaptions that just came to my mind, together with the first integer in the header of all xdrf files and that I could not find a proper meaning at time), thought for a moment that the first integer in the .xtc header mean the "shape of the box". Just played a bit later with the manual and recorded that it's only the 9 box parameters that are important to define that within gromacs. Sorry for that mistake.
<br></div></div><br><div>> > So, I would have atom names,<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Yes</blockquote><div><br>
>From .gro or .pdb file in the standard way. <br></div><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">> number of molecules types,<br>
No<br>> number of molecules of each type,
<br>No<br>> number of atoms of each molecule type,<br>No</blockquote><div><br>I think those three can become yes in a really easy way: some kind of "state machine" would easily do the trick and "count" them properly.
<br></div><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">> coordinates,<br>Yes</blockquote><div><br>From .xtc or .pdb file. <br></div>
<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">charges<br>Yes, if you use the right editconf option<br>and masses on a .pdb file, correct?
<br>No.</blockquote><div><br>Only modifying editconf to print those in a adapted .pdb file, throw a file conversion and much disk space "wasted". Or I can learn exactly from editconf how to read it from .tpr or .top (I guess it takes from these files the atomic charges when needed) and adapt into my codes.
<br></div><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">> But, still, where could I then easilly get the box type and sizes, as<br>
> well as simulation times?<br>The .xtc file, the .edr file...</blockquote><div><br>So .xtc, since its already under control. ;) <br></div><br></div><div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Fine :-) That's the sort of reason that you'd want to have to keep<br>Fortran for the present problem.</blockquote><div><br>So we agree here. Three ways so: converting the trajectories to another format (ugly for my taste, and too much disk space needed), linking to a library (like done in reading .xtc files) or linking to an object. ;)
<br></div><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Like I said, the masses are easy to infer from the atom names. You don't<br>
need to read them... you know anything starting with N has mass 14.011<br>
(or whatever it is). Or you can read share/gromacs/top/atommass.dat to<br>do the lookups.</blockquote><div><br>That is the kind of stuff I really don't like. For instance, all codes I written for all this time long aimed at being abolutelly general. So, if anything starting by N means nitrogen, how to deal when someone using the same program (even in the group) chooses to say that Nb is a "different type of nitrogen" while another says it is "niobium"? The program must be able to deal with such problems that arise (and believe, I've already got locked at this one once, and that's why I don't go such a way anymore), and the only way I see it doing so is directly reading the topologies, instead of "knowing" it in advance. ;)
<br></div><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">> Sorry for not noticing the box shape in the .xtc file. I passed over<br>> that integer in the reading, dismissing it as only some kind of old
<br>> feature useless now. Really sorry for that point.<br>><br>> Which files should I look for in order to get the proper way of reading<br>> masses and atomic charges from a .tpr file? And which objects will need
<br>> to be included in the linking stage? Is there any "outsider" program<br>> (fortran would be perfect, but I don't think it will be too common too)<br>> that works directly with the .tpr files?
<br>
<br>tpxio.c will have relevant functions. Or, follow the execution of mdrun<br>until it reads the .tpr file and see how it handles these things.</blockquote><div><br>Thanks a lot! That final information will be really usefull for me to sort out how to properly link those routines in my codes. :)
<br></div><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">If you use the editconf suggestion to get charges, you don't need to
<br>read a .tpr file.</blockquote><div><br>But then or I get locked into the problem of generality of having too much disk usage. But that idea will be also usefull to learn how to read the masses and charges from .tpr files.
<br></div><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">> Thanks a lot everybody for all the help, and sorry for taking so much of<br>
> your time,<br><br>Help is fine, but repeating myself is tiresome :-)<br></blockquote><div> <br><span class="gmail_quote"></span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
> Wait a second, you mean that editconf would allow me to have both
<br>> charges and masses on a .pdb file directly from command line (having<br>> other matters asking for my presence absolutelly now, so please forgive<br>> me if this message becomes somehow strange)<br>No, nobody said that. Please read free advice carefully, since there's
<br>no better way to annoy the giver than to not be understood when they've<br>expressed an idea clearly.</blockquote><div><br>I'm absolutelly really sorry for that. I can only apologize myself that I commited such a reading mistake on your email message. Please forgive me for being so boring at that moment.
<br><br>Thanks a lot to everybody for all help with this matter, and really sorry for taking so much of your time with these matters.<br><br>If and when everything goes smoothly here, I'll probably post something related, in order to help others that come with the same question. Possibly a small test program to read those from .tpr files and print them on screen.
<br></div></div>
<br></div>
Again, thanks a lot for all help.<br><br>Sincerally yours,<br><br>Jones<br>