Hi David,<br><br>Though I see some sense in having the proposed hack, it should be used with care since in most cases the additional groups will be ions or small subsets of atoms, which are better not to be coupled to a heat bath independently. I would advocate not to rely on such an option, but rather check the pdb for non-standard groups or use the protein/non-protein groups.
<br><br>JM2C,<br><br>Tsjerk<br><br><div><span class="gmail_quote">On 2/14/06, <b class="gmail_sendername">David van der Spoel</b> <<a href="mailto:spoel@xray.bmc.uu.se">spoel@xray.bmc.uu.se</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;">
David Mathog wrote:<br>> Tamas Horvath wrote:<br>><br>>> Check out what groups you have. The values at tau_t and ref_t refer to<br>>><br>> those<br>><br>>> groups.<br>>> In the default case, in the demo, those refer to "protein" and "SOL"
<br>>><br>> groups.<br>><br>> I went another way on this one. Since the values in each column<br>> are the same (at least in the demos) I modified readir.c so that it<br>> just filled in the missing tau_t and ref_t values with those from
<br>> the last column. After emitting a warning of course. Then it ran<br>> to completion.<br>><br>Sounds like a useful hack, since it is quite rare that one used<br>different values. Can you send me the hacked
readir.c code?<br>A general problem with automating stuff is that people may stop<br>thinking... Let's be honest if you want to get something done and get a<br>messagebox with the options Cancel OK, then you click OK (I do...).
<br>However for scientific programs the user should know what (s)he's doing...<br>> Which brings up the next mystery. When ngmx runs on the output<br>> the filter screen has<br>><br>> Group-1: system,protein,protein-h,c-alpha,backbone,
<br>> mainchain,mainchain_Cb, and MainChain+H<br>> Group-2: sidechain,sidechain-H,prot-masses,nonprotein,CA,SOL,Other<br>><br>> In the gmxdemo with cpeptide there was only one Group-1, which<br>> contained all of the entries from Groups 1 and 2 except CA.
<br>><br>> What significance does one group vs. two have? Any? Is it just that<br>> the program thinks the menu is too long and splits it at 15 entries?<br>><br>It's just that, they won't fit on the screen otherwise.
<br><br><br>If you have more hacks that you think are of general interest please<br>send them to me.<br>Thanks,<br><br>--<br>David.<br>________________________________________________________________________<br>David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,
<br>Dept. of Cell and Molecular Biology, Uppsala University.<br>Husargatan 3, Box 596, 75124 Uppsala, Sweden<br>phone: 46 18 471 4205 fax: 46 18 511 755<br><a href="mailto:spoel@xray.bmc.uu.se">spoel@xray.bmc.uu.se
</a> <a href="mailto:spoel@gromacs.org">spoel@gromacs.org</a> <a href="http://folding.bmc.uu.se">http://folding.bmc.uu.se</a><br>++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++<br><br>_______________________________________________
<br>gmx-users mailing list <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br><a href="http://www.gromacs.org/mailman/listinfo/gmx-users">http://www.gromacs.org/mailman/listinfo/gmx-users</a><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/mailing_lists/users.php">http://www.gromacs.org/mailing_lists/users.php
</a><br></blockquote></div><br><br clear="all"><br>-- <br><br>Tsjerk A. Wassenaar, M.Sc.<br>Groningen Biomolecular Sciences and Biotechnology Institute (GBB)<br>Dept. of Biophysical Chemistry<br>University of Groningen<br>
Nijenborgh 4<br>9747AG Groningen, The Netherlands<br>+31 50 363 4336<br>