Dear All
I performed the tutorial lesson on Molecular Dynamics, tmoldyn_01, with Abinit-8.0.7. The calculation was successfully completed but post processing of the MD results with diag_moldyn.py got the following error:
...
[thanusit@physics-clusters tmoldyn_01]$ python diag_moldyn.py
Vous n'avez pas de fichier netcdf pour diagnostics !
...
I wonder if this is due to the script cannot cope with the change of file extension introduced in abinit-8.0.7, i.e. from _HIST to _HIST.nc, or other issues?
(Please note that my Abinit-8.0.7 was built with links against python-2.7.5, netcdf-4.3.1, hdf-5.18.12, openmpi-1.6.5.)
Regards,
Thanusit
diag_moldyn.py gets error for MD ouput files of Abinit-8.0.7 [SOLVED]
Moderators: MMNSchmitt, gonze
Re: diag_moldyn.py gets error for MD ouput files of Abinit-8
I wonder if this is due to the script cannot cope with the change of file extension introduced in abinit-8.0.7, i.e. from _HIST to _HIST.nc, or other issues?
Yes, could you try this patch?
=== modified file 'doc/tutorial/lesson_paral_moldyn/diag_moldyn.py'
--- doc/tutorial/lesson_paral_moldyn/diag_moldyn.py 2015-01-24 08:15:18 +0000
+++ doc/tutorial/lesson_paral_moldyn/diag_moldyn.py 2016-06-29 11:22:09 +0000
@@ -43,7 +43,7 @@
ok_OUT=-1;ok_MOLDYN=-1
if len(OUT_list) > 0 :
OUT_list.sort() ; OUT_list.reverse()
- fic_HIST=OUT_list[0][1].replace('_OUT.nc','_HIST')
+ fic_HIST=OUT_list[0][1].replace('_OUT.nc','_HIST.nc')
if os.path.exists(fic_HIST) :
fiche = fic_HIST
ok_OUT=0
@@ -74,7 +74,7 @@
TypeFichier=1
fichier1=file
fichier2=''
- elif file.find('_HIST') == len(file)-5 :
+ elif file.find('_HIST.nc') == len(file)-8 :
TypeFichier=2
fichier1=file
fichier2=file.replace('_HIST','_OUT.nc')
Re: diag_moldyn.py gets error for MD ouput files of Abinit-8
Dear Gmatteo
Thank you for your reply.
I tried the patch you suggested. This fixes the above error but a new one came up, as shown below.
It looks like the '.nc' extension is appended to '_OUT.nc' so that the script is looking for '_OUT.nc.nc' which doesn't exist.
Regards,
Thanusit
Thank you for your reply.
I tried the patch you suggested. This fixes the above error but a new one came up, as shown below.
Code: Select all
[thanusit@physics-clusters tmoldyn_01]$ python diag_moldyn_patched.py
Traceback (most recent call last):
File "diag_moldyn_patched.py", line 213, in <module>
FDR = lireNetcdf()
File "diag_moldyn_patched.py", line 60, in lireNetcdf
ConvNetCdf_Ascii(fiche)
File "diag_moldyn_patched.py", line 122, in ConvNetCdf_Ascii
ncfile2 = netcdf.netcdf_file(fichier2, 'r')
File "/share/apps/python2.7/gcc-4.7.2/lib/python2.7/site-packages/scipy/io/netcdf.py", line 198, in __init__
self.fp = open(self.filename, '%sb' % mode)
IOError: [Errno 2] No such file or directory: 'tmoldyn_out_OUT.nc.nc'
Exception AttributeError: "'netcdf_file' object has no attribute 'fp'" in <bound method netcdf_file.close of <scipy.io.netcdf.netcdf_file object at 0xfb0610>> ignored
It looks like the '.nc' extension is appended to '_OUT.nc' so that the script is looking for '_OUT.nc.nc' which doesn't exist.
Regards,
Thanusit
Re: diag_moldyn.py gets error for MD ouput files of Abinit-8
I'm attaching another version of the script that seems to work.
I had to change the file extension by adding a `.log` at the end in order to upload the file.
HTH
M
I had to change the file extension by adding a `.log` at the end in order to upload the file.
HTH
M
- Attachments
-
diag_moldyn.py.log- (7.75 KiB) Downloaded 439 times
Re: diag_moldyn.py gets error for MD ouput files of Abinit-8
Dear Gmatteo
That works.
Thank you very much.
Thanusit
That works.
Thank you very much.
Thanusit