This email reminded me about a converter that I made for Materials Studio xtd and trj files to MMTK/nMoldyn&#39;s nc format. I have attached it for anyone interested. <br><br>It shows how to use the &quot;snapshot&quot; feature to build up the nc file and how to insert time information into that snapshot. <br>

<br>If anyone is actually going to use it for MS-&gt;MMTK stuff, you&#39;ll need the lxml package for python and you will need to convert your files to ASCII using trj2ascii.exe. (trj2ascii.exe is a small program found in the MS root directory on Windows.) You need to provide both the ascii trj file and the xtd file. It worked fine for MS4.1, but I haven&#39;t tested it with other versions.<br>
<br>Jay<br><br><div class="gmail_quote">On Fri, Jan 23, 2009 at 9:30 AM, Konrad Hinsen <span dir="ltr">&lt;<a href="mailto:hinsen@cnrs-orleans.fr" target="_blank">hinsen@cnrs-orleans.fr</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div>On Jan 23, 2009, at 2:50, Yi Liu wrote:<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
But I still do not understand why G3 works but G4 not? I essentially<br>
created them in the same way. I considered the other differences<br>
between these two systems:<br>
(1) the number of atoms should not matter since the good G3 systems<br>
have more atoms;<br>
(2) the number of residuals (G3 only has one residual but G4 have<br>
many). How this could cause the problem?<br>
<br>
I can not get more hints by comparing the description strings in these<br>
two cdl files. Do you have any idea about it?<br>
</blockquote>
<br></div>
Since it is ncgen that complains about the description string, the problem cannot be related to its content. For ncgen, it is just some string. So there is either a syntactical problem with representing the string in cdl, or a bug in ncgen. I haven&#39;t use ncgen or cdl very much, so I cannot say more about this.<div>

<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
I attempted to extract the new description line generated from<br>
SnapshotGenerator (convert nc to cdl and cut it out), and replaced the<br>
description line in my cdl files and used the correct string length.<br>
Nevertheless, it did not work.<br>
</blockquote>
<br></div>
Something that you could try is running ncdump on a valid trajectory and then ncgen on the result. If that produces an error or an incorrect file, something is necessarily wrong with ncdump and/or ncgen.<div>
<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Side notes: The way I got the current description line is quite<br>
&quot;dirty&quot;. I actually did a short run using NAMD, and convert the output<br>
dcd to nc (then cdl) using dcd_to_nc converter. Then I am able to<br>
extract the description line from cdl file and plug into my cdl files.<br>
</blockquote>
<br></div>
Then why not take the dcd_to_nc converter and modify it to read your LAMMPS file instead of the DCD file? All you need to know to do that is the LAMMPS format, but apparently you have all that information. If there is something you don&#39;t understand in the dcd_to_nc converter, just contact me.<br>


<br>
This seems the best solution to me because in the end you&#39;d have a general one-step LAMMPS-to-nc converter that is much easier to use.<div><br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
I thought about the version of netCDF. I recently updated to 4.0. Even<br>
though there are two new binary formats are introduced from 4.0, it<br>
should be &quot;classic format&quot; by default. And I can get G3 working using<br>
netCDF 4.0.<br>
</blockquote>
<br></div>
There shouldn&#39;t be any problem with netCDF 4, as long as you use the classic format. The 64-bit format should work as well without changes, but you need ScientificPython 2.9 to be able to create such files. I haven&#39;t tested HDF5 format support yet, so I cannot say if it works.<div>

<br>
<br>
Konrad.<br>
--<br>
---------------------------------------------------------------------<br>
Konrad Hinsen<br>
Centre de Biophysique Moléculaire, CNRS Orléans<br>
Synchrotron Soleil - Division Expériences<br>
Saint Aubin - BP 48<br>
91192 Gif sur Yvette Cedex, France<br>
Tel. +33-1 69 35 97 15<br>
E-Mail: <a href="mailto:hinsen@cnrs-orleans.fr" target="_blank">hinsen@cnrs-orleans.fr</a><br></div>
Web: <a href="http://dirac.cnrs-orleans.fr/%7Ehinsen/" target="_blank">http://dirac.cnrs-orleans.fr/~hinsen/</a><br>
---------------------------------------------------------------------<div><div></div><div><br>
<br>
<br>
<br>
<br>
_______________________________________________<br>
mmtk maillist &nbsp;- &nbsp;<a href="mailto:mmtk@starship.python.net" target="_blank">mmtk@starship.python.net</a><br>
<a href="http://starship.python.net/mailman/listinfo/mmtk" target="_blank">http://starship.python.net/mailman/listinfo/mmtk</a><br>
</div></div></blockquote></div><br>