SAMBA-ilum 1.0.0.541__tar.gz → 1.0.0.543__tar.gz

{samba_ilum-1.0.0.541 → samba_ilum-1.0.0.543}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: SAMBA_ilum
-Version: 1.0.0.541
+Version: 1.0.0.543
 Summary: ...
 Author: Augusto de Lelis Araujo
 Author-email: augusto-lelis@outlook.com

{samba_ilum-1.0.0.541 → samba_ilum-1.0.0.543}/SAMBA_ilum.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: SAMBA_ilum
-Version: 1.0.0.541
+Version: 1.0.0.543
 Summary: ...
 Author: Augusto de Lelis Araujo
 Author-email: augusto-lelis@outlook.com

{samba_ilum-1.0.0.541 → samba_ilum-1.0.0.543}/SAMBA_ilum.egg-info/SOURCES.txt

@@ -31,6 +31,7 @@ samba_ilum/src/kpoints.py
 samba_ilum/src/lattice_plot3d.py
 samba_ilum/src/make_files.py
 samba_ilum/src/output.py
+samba_ilum/src/poscar_fix_direct_coord.py
 samba_ilum/src/potcar.py
 samba_ilum/src/vdW_DF.py
 samba_ilum/src/xy-scan.py

{samba_ilum-1.0.0.541 → samba_ilum-1.0.0.543}/samba_ilum/__main__.py

@@ -16,7 +16,7 @@ dir_samba = os.path.dirname(os.path.realpath(__file__))
 print(f'{dir_samba}')
 #--------------------
 
-version = '1.0.0.541'
+version = '1.0.0.543'
 
 print(" ")
 print("=============================================================")
@@ -59,6 +59,13 @@ except Exception as e:
     print(" ")
 
 
+#=============================================================
+# Fixing the coordinates of POSCAR files in direct form ======
+#=============================================================
+# exec(open(dir_codes + '/poscar_fix_direct_coord.py').read())
+#=============================================================
+
+
 # ------------------------------------------------------------------------------
 # Checking if the "run.input" file exists --------------------------------------
 # ------------------------------------------------------------------------------

{samba_ilum-1.0.0.541 → samba_ilum-1.0.0.543}/samba_ilum/src/HeteroStructure_Generator.py

@@ -1707,6 +1707,7 @@ print(f'{n_structures} remaining cells')
 
 #=====================================================================
 # Filtering/Excluding Similar POSCAR Files ===========================
+# within the temporary subdirectories of the POSCAR files ============
 #=====================================================================
 
 n += 1
@@ -1825,7 +1826,7 @@ print(f'{n_structures} remaining cells')
 
 
 #=====================================================================
-# Deleting temporary directory of POSCAR files =======================
+# Deleting temporary subdirectories of POSCAR files ==================
 #=====================================================================
 number = 0
 #--------------------------
@@ -1845,114 +1846,125 @@ for i in range(len(files0)):
 
 #=====================================================================
 # Filtering/Excluding Similar POSCAR Files ===========================
+# within the final directory of the POSCAR files =====================
 #=====================================================================
 
-n += 1
-#----------
-print(f' ')
-print(f'========================================')
-print(f'Step {n}: Deleting Similar Lattices ====')
-print(f'========================================')
+#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+# Introduzir o loop de exclusão, até que duas interações retornem exclusão nula
+#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
-files = os.listdir(diret2)
-n_exclusion = 0
+n_exclusion_cont = 0
 
-for i in range(len(files)):
-    #---------------------------------
-    print(f'Progress {i+1}/{len(files)}') 
-    #----------------
-    files_delete = []
-    #------------------------------------
-    if os.path.exists(diret2 + files[i]):
-       #---------------------------------------------------------
-       structure1 = Poscar.from_file(diret2 + files[i]).structure
-       #---------------------------------------------------------
-       poscar = open(diret2 + files[i], "r")
-       #--------------------------------
-       VTemp = poscar.readline().split()
-       if (n_Lattice == 2):
-          Angle_AB1 = float(VTemp[27])
-       if (n_Lattice == 3):
-          temp_angle = VTemp[28].replace('_',' ').split()
-          Angle_AB1 = float(temp_angle[0]) 
-       #---------------------------------------------- 
-       param = poscar.readline(); param = float(param)
-       VTemp = poscar.readline().split();  Ax = float(VTemp[0])*param; Ay = float(VTemp[1])*param;  A = np.array([Ax, Ay]);  mA = np.linalg.norm(A)
-       VTemp = poscar.readline().split();  Bx = float(VTemp[0])*param; By = float(VTemp[1])*param;  B = np.array([Bx, By]);  mB = np.linalg.norm(B)
-       dAB1 = abs(mA - mB)
-       poscar.close()
-       #-------------
-
-       for j in range(len(files)):
-           if (i != j):
-              #------------------------------------
-              if os.path.exists(diret2 + files[j]):
-                 #---------------------------------------------------------
-                 structure2 = Poscar.from_file(diret2 + files[j]).structure
-                 #----------------------------------------------------------------------------------
-                 matcher = StructureMatcher()                   # Creating a StructureMatcher object
-                 if matcher.fit(structure1, structure2):        # Checking if the structures are similar
-                    #-----------------------------------------------------------------------------------
-                    poscar = open(diret2 + files[j], "r")
-                    #--------------------------------
-                    VTemp = poscar.readline().split()
-                    if (n_Lattice == 2):
-                       Angle_AB2 = float(VTemp[27])
-                    if (n_Lattice == 3):
-                       temp_angle = VTemp[28].replace('_',' ').split()
-                       Angle_AB2 = float(temp_angle[0]) 
-                    #----------------------------------------------
-                    param = poscar.readline(); param = float(param)
-                    VTemp = poscar.readline().split();  Ax = float(VTemp[0])*param; Ay = float(VTemp[1])*param;  A = np.array([Ax, Ay]);  mA = np.linalg.norm(A)
-                    VTemp = poscar.readline().split();  Bx = float(VTemp[0])*param; By = float(VTemp[1])*param;  B = np.array([Bx, By]);  mB = np.linalg.norm(B)
-                    dAB2 = abs(mA - mB)
-                    poscar.close()
-                    #----------------
-                    if (dAB1 < dAB2):
-                       if os.path.exists(diret2 + + files[j]):
-                          os.remove(diret2 + files[j])
-                          n_exclusion += 1
-                    #----------------
-                    if (dAB1 > dAB2):
-                       if os.path.exists(diret2 + files[i]):
-                          os.remove(diret2 + files[i])
-                          n_exclusion += 1
-                    #-----------------
-                    if (dAB1 == dAB2):
-                       cont = 1
-                       #--------------------------------------------------------
-                       if ( (Angle_AB1 < 0) and (Angle_AB2 > 0) and cont == 1 ):
-                          if os.path.exists(diret2 + files[i]):
-                             os.remove(diret2 + files[i])
-                             n_exclusion += 1; cont = 0
-                       if ( (Angle_AB1 > 0) and (Angle_AB2 < 0) and cont == 1 ):
-                          if os.path.exists(diret2 + files[j]):
-                             os.remove(diret2 + files[j])
-                             n_exclusion += 1; cont = 0
-                       if (abs(Angle_AB1) < abs(Angle_AB2) and cont == 1):
-                          if os.path.exists(diret2 + files[j]):
+while n_exclusion_cont != 2:
+
+   n += 1
+   #----------
+   print(f' ')
+   print(f'========================================')
+   print(f'Step {n}: Deleting Similar Lattices ====')
+   print(f'========================================')
+
+   files = os.listdir(diret2)
+   n_exclusion = 0
+
+   for i in range(len(files)):
+       #---------------------------------
+       print(f'Progress {i+1}/{len(files)}') 
+       #----------------
+       files_delete = []
+       #------------------------------------
+       if os.path.exists(diret2 + files[i]):
+          #---------------------------------------------------------
+          structure1 = Poscar.from_file(diret2 + files[i]).structure
+          #---------------------------------------------------------
+          poscar = open(diret2 + files[i], "r")
+          #--------------------------------
+          VTemp = poscar.readline().split()
+          if (n_Lattice == 2):
+             Angle_AB1 = float(VTemp[27])
+          if (n_Lattice == 3):
+             temp_angle = VTemp[28].replace('_',' ').split()
+             Angle_AB1 = float(temp_angle[0]) 
+          #---------------------------------------------- 
+          param = poscar.readline(); param = float(param)
+          VTemp = poscar.readline().split();  Ax = float(VTemp[0])*param; Ay = float(VTemp[1])*param;  A = np.array([Ax, Ay]);  mA = np.linalg.norm(A)
+          VTemp = poscar.readline().split();  Bx = float(VTemp[0])*param; By = float(VTemp[1])*param;  B = np.array([Bx, By]);  mB = np.linalg.norm(B)
+          dAB1 = abs(mA - mB)
+          poscar.close()
+          #-------------
+          for j in range(len(files)):
+              if (i != j):
+                 #------------------------------------
+                 if os.path.exists(diret2 + files[j]):
+                    #---------------------------------------------------------
+                    structure2 = Poscar.from_file(diret2 + files[j]).structure
+                    #----------------------------------------------------------------------------------
+                    matcher = StructureMatcher()                   # Creating a StructureMatcher object
+                    if matcher.fit(structure1, structure2):        # Checking if the structures are similar
+                       #-----------------------------------------------------------------------------------
+                       poscar = open(diret2 + files[j], "r")
+                       #--------------------------------
+                       VTemp = poscar.readline().split()
+                       if (n_Lattice == 2):
+                          Angle_AB2 = float(VTemp[27])
+                       if (n_Lattice == 3):
+                          temp_angle = VTemp[28].replace('_',' ').split()
+                          Angle_AB2 = float(temp_angle[0]) 
+                       #----------------------------------------------
+                       param = poscar.readline(); param = float(param)
+                       VTemp = poscar.readline().split();  Ax = float(VTemp[0])*param; Ay = float(VTemp[1])*param;  A = np.array([Ax, Ay]);  mA = np.linalg.norm(A)
+                       VTemp = poscar.readline().split();  Bx = float(VTemp[0])*param; By = float(VTemp[1])*param;  B = np.array([Bx, By]);  mB = np.linalg.norm(B)
+                       dAB2 = abs(mA - mB)
+                       poscar.close()
+                       #----------------
+                       if (dAB1 < dAB2):
+                          if os.path.exists(diret2 + + files[j]):
                              os.remove(diret2 + files[j])
-                             n_exclusion += 1; cont = 0
-                       if (abs(Angle_AB1) > abs(Angle_AB2) and cont == 1):
+                             n_exclusion += 1
+                       #----------------
+                       if (dAB1 > dAB2):
                           if os.path.exists(diret2 + files[i]):
                              os.remove(diret2 + files[i])
-                             n_exclusion += 1; cont = 0
-                       if (Angle_AB1 == Angle_AB2 and cont == 1):
-                          if os.path.exists(diret2 + files[j]):
-                             os.remove(diret2 + files[j])
                              n_exclusion += 1
-
-#-------------------------
-files = os.listdir(diret2)
-n_structures = len(files)
-#------------------------
-
-print(f'============================================')
-print(f'{n_exclusion} cells were filtered/excluded')
-print(f'--------------------------------------------')
-if (n_structures == 1): print(f'Only one cell was found')
-if (n_structures >  1): print(f'{n_structures} cells were found')
-print(f'============================================')
+                       #-----------------
+                       if (dAB1 == dAB2):
+                          cont = 1
+                          #--------------------------------------------------------
+                          if ( (Angle_AB1 < 0) and (Angle_AB2 > 0) and cont == 1 ):
+                             if os.path.exists(diret2 + files[i]):
+                                os.remove(diret2 + files[i])
+                                n_exclusion += 1; cont = 0
+                          if ( (Angle_AB1 > 0) and (Angle_AB2 < 0) and cont == 1 ):
+                             if os.path.exists(diret2 + files[j]):
+                                os.remove(diret2 + files[j])
+                                n_exclusion += 1; cont = 0
+                          if (abs(Angle_AB1) < abs(Angle_AB2) and cont == 1):
+                             if os.path.exists(diret2 + files[j]):
+                                os.remove(diret2 + files[j])
+                                n_exclusion += 1; cont = 0
+                          if (abs(Angle_AB1) > abs(Angle_AB2) and cont == 1):
+                             if os.path.exists(diret2 + files[i]):
+                                os.remove(diret2 + files[i])
+                                n_exclusion += 1; cont = 0
+                          if (Angle_AB1 == Angle_AB2 and cont == 1):
+                             if os.path.exists(diret2 + files[j]):
+                                os.remove(diret2 + files[j])
+                                n_exclusion += 1
+   #-------------------------
+   files = os.listdir(diret2)
+   n_structures = len(files)
+   #------------------------
+
+   print(f'============================================')
+   print(f'{n_exclusion} cells were filtered/excluded')
+   print(f'--------------------------------------------')
+   if (n_structures == 1): print(f'Only one cell was found')
+   if (n_structures >  1): print(f'{n_structures} cells were found')
+   print(f'============================================')
+
+   if (n_exclusion != 0):  n_exclusion_cont  = 0
+   if (n_exclusion == 0):  n_exclusion_cont += 1
+   if (n_structures == 1): n_exclusion_cont  = 2
 
 
 #=====================================================================

{samba_ilum-1.0.0.541 → samba_ilum-1.0.0.543}/samba_ilum/src/cell_check.py

@@ -177,7 +177,7 @@ for name in files0:
         menor_0 = min(angulos_z) if angulos_z else 0.0
         menor_1 = min(angulos_z_na_origem) if angulos_z_na_origem else 0.0
         #--------------------------------------------------------------------
-        # If there is a difference, calculate, correct and save the new file.
+        # If there is a difference, the rotation angle position is fixed ----
         #--------------------------------------------------------------------
         if (menor_0 != menor_1 and menor_0 != 0.0):
             axis_correction += 1

samba_ilum-1.0.0.543/samba_ilum/src/energy_scan.py

@@ -0,0 +1,57 @@
+# SAMBA_ilum Copyright (C) 2025 - Closed source
+
+
+import os
+
+
+try:
+    path_dir, name_dir = os.path.split(os.getcwd())
+except Exception as e:
+    pass
+
+
+#----------------
+converged = False
+#----------------
+try:
+    with open('OUTCAR', 'r') as f:
+        for line in f:
+            if 'reached' in line:
+                converged = True
+                break
+except FileNotFoundError:
+    pass
+except Exception as e:
+    pass
+
+
+if converged:
+    try:
+        with open('OSZICAR', 'r') as f:
+            last_line = f.readlines()[-1].strip()
+            energy_value = last_line.replace('=', ' ').split()[4]
+        #--------------------------------------------------------
+        with open('../energy_scan.txt', "a") as energy_file:
+            energy_file.write(f'{name_dir} {energy_value}\n')
+            #---------------------------------------------
+            temp_name = name_dir.replace('_', ' ').split()
+            #---------------------------------------------
+            if len(temp_name) == 2:
+                t_temp_name_original = temp_name.copy()
+                #---------------------------------------------
+                if temp_name[0] == '0.0': temp_name[0] = '1.0'
+                if temp_name[1] == '0.0': temp_name[1] = '1.0'
+                #---------------------------------------------
+                if temp_name[0] == '1.0' or temp_name[1] == '1.0':
+                    new_name_dir = f'{temp_name[0]}_{temp_name[1]}'
+                    energy_file.write(f'{new_name_dir} {energy_value}\n')
+                #------------------------------------------------------------------------
+                if t_temp_name_original[0] == '0.0' and t_temp_name_original[1] == '0.0':
+                    energy_file.write(f'1.0_0.0 {energy_value}\n')
+                    energy_file.write(f'0.0_1.0 {energy_value}\n')
+    except FileNotFoundError:
+        pass
+    except IndexError:
+        pass
+    except Exception as e:
+        pass

{samba_ilum-1.0.0.541 → samba_ilum-1.0.0.543}/samba_ilum/src/job.py

@@ -130,8 +130,8 @@ job.write(f'                         $vasp_std || echo "Error in VASP in $dir0/$
 job.write(f'                         sleep 10 \n')
 job.write(f'                         python3 contcar_update.py \n')
 job.write(f'                         python3 energy_scan.py \n')
-job.write(f'                         rm -r CHG CHGCAR DOSCAR PROCAR WAVECAR INCAR KPOINTS OSZICAR PCDAT POTCAR REPORT XDATCAR vdw_kernel.bindat EIGENVAL IBZKPT vasprun.xml *.py \n')
-job.write(f'                         bzip2 OUTCAR \n')
+job.write(f'                         rm -r CHG CHGCAR DOSCAR PROCAR WAVECAR INCAR KPOINTS PCDAT POTCAR REPORT XDATCAR vdw_kernel.bindat EIGENVAL IBZKPT vasprun.xml *.py \n')
+job.write(f'                         bzip2 OUTCAR OSZICAR \n')
 job.write(f'                     done \n')
 job.write(f'                     #------------- \n')
 job.write(f'                     cd $dir0/$i/$j \n')
@@ -174,8 +174,8 @@ job.write(f'                                $vasp_std || echo "Error in VASP in
 job.write(f'                                sleep 10 \n')
 job.write(f'                                python3 contcar_update.py \n')
 job.write(f'                                python3 energy_scan.py \n')
-job.write(f'                                rm -r CHG CHGCAR DOSCAR PROCAR WAVECAR INCAR KPOINTS OSZICAR PCDAT POTCAR REPORT XDATCAR vdw_kernel.bindat EIGENVAL IBZKPT vasprun.xml *.py \n')
-job.write(f'                                bzip2 OUTCAR \n')
+job.write(f'                                rm -r CHG CHGCAR DOSCAR PROCAR WAVECAR INCAR KPOINTS PCDAT POTCAR REPORT XDATCAR vdw_kernel.bindat EIGENVAL IBZKPT vasprun.xml *.py \n')
+job.write(f'                                bzip2 OUTCAR OSZICAR \n')
 job.write(f'                                cd $dir0/$i/$j \n')
 job.write(f'                                mv $k OK_$k \n')
 job.write(f'                             fi \n')

samba_ilum-1.0.0.543/samba_ilum/src/poscar_fix_direct_coord.py

@@ -0,0 +1,117 @@
+# SAMBA_ilum Copyright (C) 2025 - Closed source
+
+
+import numpy as np
+import shutil
+import glob
+import sys
+import os
+
+
+#----------------------------------------------------------------------------
+# Function to delete hidden files -------------------------------------------
+#----------------------------------------------------------------------------
+def delete_hidden_files(target_directory):
+    search_files = os.path.join(target_directory, '.*')
+    for item in glob.glob(search_files):
+        try:
+            if os.path.isfile(item): os.remove(item)
+        except OSError: pass
+
+
+#----------------------------------------------------
+# Listing files inside the "dir_poscar" directory ---
+#----------------------------------------------------
+poscar_dir_path = os.path.join(dir_files, dir_poscar)
+#-----------------------------------
+delete_hidden_files(poscar_dir_path)
+files0 = [name for name in os.listdir(poscar_dir_path) if os.path.isfile(os.path.join(poscar_dir_path, name))]
+
+
+#======================================================
+# Checking/Correcting POSCAR files ====================
+#======================================================
+poscar_backup_dir_path = os.path.join(dir_files, 'POSCAR_Backup_temp')
+os.rename(poscar_dir_path, poscar_backup_dir_path)
+os.makedirs(poscar_dir_path, exist_ok=True)
+#------------------------------------------
+files1 = [name for name in os.listdir(poscar_backup_dir_path) if os.path.isfile(os.path.join(poscar_backup_dir_path, name))]
+for file_in in files1:
+    #-------------------
+    type_poscar = 'none'
+    poscar_in  = open(poscar_backup_dir_path + '/' + file_in, "r")
+    for i in range(8): VTemp = poscar_in.readline().split()
+    if (VTemp[0][0] == 'D' or VTemp[0][0] == 'd'): type_poscar = 'direct'
+    if (VTemp[0][0] == 'C' or VTemp[0][0] == 'c'): type_poscar = 'cartesian'
+    if (VTemp[0][0] == 'S' or VTemp[0][0] == 's'): type_poscar = 'selective'
+    poscar_in.close()
+    #----------------
+
+    if (type_poscar != 'cartesian'): shutil.copy2(poscar_backup_dir_path + '/' + file_in, poscar_dir_path + '/' + file_in)
+
+    if (type_poscar == 'cartesian'):
+       #-------------------------------------------------------------
+       poscar_in  = open(poscar_backup_dir_path + '/' + file_in, "r")
+       poscar_out = open(poscar_dir_path + '/' + file_in, "w")
+       #------------------------------------------------------
+       VTemp = poscar_in.readline()
+       poscar_out.write(f'{VTemp}')
+       #------------------------------------------------------------
+       VTemp = poscar_in.readline().split(); param = float(VTemp[0])
+       poscar_out.write(f'1.0 \n')
+       #-----------------------------------
+       VTemp = poscar_in.readline().split()
+       A1x = float(VTemp[0])*param
+       A1y = float(VTemp[1])*param
+       A1z = float(VTemp[2])*param
+       poscar_out.write(f'{A1x:.16} {A1y:.16} {A1z:.16} \n')
+       #----------------------------------------------------
+       VTemp = poscar_in.readline().split()
+       A2x = float(VTemp[0])*param
+       A2y = float(VTemp[1])*param
+       A2z = float(VTemp[2])*param
+       poscar_out.write(f'{A2x:.16} {A2y:.16} {A2z:.16} \n')
+       #----------------------------------------------------
+       VTemp = poscar_in.readline().split()
+       A3x = float(VTemp[0])*param
+       A3y = float(VTemp[1])*param
+       A3z = float(VTemp[2])*param
+       poscar_out.write(f'{A3x:.16} {A3y:.16} {A3z:.16} \n')
+       #----------------------------------------------------
+       VTemp = poscar_in.readline()
+       poscar_out.write(f'{VTemp}')
+       #---------------------------
+       nions = 0
+       VTemp = poscar_in.readline()
+       poscar_out.write(f'{VTemp}')
+       VTemp = VTemp.split()
+       for i in range(len(VTemp)): nions += int(VTemp[i])
+       #-------------------------------------------------
+       VTemp = poscar_in.readline()
+       poscar_out.write(f'Direct \n')
+       #-----------------------------
+       a = np.array([A1x, A1y, A1z])
+       b = np.array([A2x, A2y, A2z])
+       c = np.array([A3x, A3y, A3z])
+       T = np.linalg.inv(np.array([a, b, c]).T)      # Defines the transformation matrix
+       #---------------------------------------
+       for i in range(nions):
+           #-----------------------------------
+           VTemp = poscar_in.readline().split()
+           x = float(VTemp[0])*param
+           y = float(VTemp[1])*param
+           z = float(VTemp[2])*param
+           #------------------------
+           r = np.array([x, y, z])      # Defines the Cartesian position of the atom
+           f = np.dot(T, r)             # Calculates fractional position
+           f = f % 1.0
+           #-------------------------------------------------------
+           poscar_out.write(f'{f[0]:.16} {f[1]:.16} {f[2]:.16} \n')
+       #-----------------------------------------------------------
+       poscar_in.close()   
+       poscar_out.close()
+       #-----------------
+
+#===============================
+os.rmdir(poscar_backup_dir_path)
+#===============================

{samba_ilum-1.0.0.541 → samba_ilum-1.0.0.543}/setup.py

@@ -6,7 +6,7 @@ import json
 
 setup(
     name = "SAMBA_ilum",
-    version = "1.0.0.541",
+    version = "1.0.0.543",
     entry_points={'console_scripts': ['samba_ilum = samba_ilum:main']},
     description = "...",
     author = "Augusto de Lelis Araujo", 

samba_ilum-1.0.0.541/samba_ilum/src/energy_scan.py

@@ -1,35 +0,0 @@
-# SAMBA_ilum Copyright (C) 2025 - Closed source
-
-
-import os
-
-#===================================================
-# Getting the path and name of the current directory:
-path_dir, name_dir = os.path.split(os.getcwd())
-#===============================================
-
-if os.path.isfile('../energy_scan.txt'):  energy = open('../energy_scan.txt', "a")
-else:  energy = open('../energy_scan.txt', "w")
-
-with open('OSZICAR') as file:
-    lines = file.readlines()
-VTemp = lines[-1].strip()
-
-energia = VTemp.replace('=',' ').split()
-
-energy.write(f'{name_dir} {energia[4]} \n')
-
-temp_name = name_dir.replace('_', ' ').split()
-t_temp_name = temp_name
-#------------------------
-if (len(temp_name) == 2):
-   if (temp_name[0] == '0.0'): temp_name[0] = '1.0'
-   if (temp_name[1] == '0.0'): temp_name[1] = '1.0'
-   if (temp_name[0] == '1.0' or temp_name[1] == '1.0'):
-      new_name_dir = str(temp_name[0]) + '_' + str(temp_name[1])
-      energy.write(f'{new_name_dir} {energia[4]} \n')
-   if (t_temp_name[0] == '1.0' and t_temp_name[1] == '1.0'):
-      energy.write(f'1.0_0.0 {energia[4]} \n')
-      energy.write(f'0.0_1.0 {energia[4]} \n')
-
-energy.close()