adapted for volatile abinit format

ab2tds/TDS_Reading.py

@@ -758,7 +758,12 @@ def ReadCastep(filename, filename_extra=None,
                     items=string.split(sl[linecount+(imode/5)  ]  )
                     frequencies[iwave,imode] = float(string.atof(items[imode%5]))*energy_scaling*hartree2cm
 
-                while( "; " not in sl[linecount] ):
+
+                    
+                while( "Eigendisplacements" not in sl[linecount]  ):
+                    linecount+=1
+                    
+                while( ("; " not in sl[linecount][:2])      and ("- " not in sl[linecount][:2])    ):
                     linecount+=1
                     
                 for imode in range(NofBranches):
@@ -778,9 +783,11 @@ def ReadCastep(filename, filename_extra=None,
                     linecount+=1
                     
                     print " LETTO UN VETTORE ", (eigenvectors[iwave,imode].imag*eigenvectors[iwave,imode].imag + eigenvectors[iwave,imode].real*eigenvectors[iwave,imode].real).sum()
-                    
-                    while( "; " not in sl[linecount] ):
-                        linecount+=1
+
+                    if imode != NofBranches-1:
+                        while( ("; " not in sl[linecount][:2])      and ("- " not in sl[linecount][:2])    ):
+
+                            linecount+=1