@@ -281,10 +281,12 @@ def check_if_melted(self, lmp, filename):
281281 lmp .close ()
282282 # Preserve log file on error
283283 logfile = os .path .join (self .simfolder , "log.lammps" )
284- if os . path . exists ( logfile ) :
284+ try :
285285 os .rename (
286286 logfile , os .path .join (self .simfolder , "melted_error.log.lammps" )
287287 )
288+ except OSError as e :
289+ self .logger .warning (f"Failed to rename log file: { e } " )
288290 raise MeltedError (
289291 "System melted, increase size or reduce temp!\n Solid detection algorithm only works with BCC/FCC/HCP/SC/DIA. Detection algorithm can be turned off by setting:\n tolerance.solid_fraction: 0"
290292 )
@@ -296,10 +298,12 @@ def check_if_solidfied(self, lmp, filename):
296298 lmp .close ()
297299 # Preserve log file on error
298300 logfile = os .path .join (self .simfolder , "log.lammps" )
299- if os . path . exists ( logfile ) :
301+ try :
300302 os .rename (
301303 logfile , os .path .join (self .simfolder , "solidified_error.log.lammps" )
302304 )
305+ except OSError as e :
306+ self .logger .warning (f"Failed to rename log file: { e } " )
303307 raise SolidifiedError ("System solidified, increase temperature" )
304308
305309 def fix_nose_hoover (
@@ -609,13 +613,15 @@ def run_iterative_pressure_convergence(self, lmp):
609613 lmp .close ()
610614 # Preserve log file on error
611615 logfile = os .path .join (self .simfolder , "log.lammps" )
612- if os . path . exists ( logfile ) :
616+ try :
613617 os .rename (
614618 logfile ,
615619 os .path .join (
616620 self .simfolder , "pressure_convergence_error.log.lammps"
617621 ),
618622 )
623+ except OSError as e :
624+ self .logger .warning (f"Failed to rename log file: { e } " )
619625 raise ValueError (
620626 "Pressure did not converge after MD runs, maybe change lattice_constant and try?"
621627 )
@@ -732,13 +738,15 @@ def run_iterative_constrained_pressure_convergence(self, lmp):
732738 lmp .close ()
733739 # Preserve log file on error
734740 logfile = os .path .join (self .simfolder , "log.lammps" )
735- if os . path . exists ( logfile ) :
741+ try :
736742 os .rename (
737743 logfile ,
738744 os .path .join (
739745 self .simfolder , "constrained_pressure_error.log.lammps"
740746 ),
741747 )
748+ except OSError as e :
749+ self .logger .warning (f"Failed to rename log file: { e } " )
742750 raise ValueError ("pressure did not converge" )
743751
744752 def process_pressure (
@@ -1333,7 +1341,7 @@ def temperature_scaling(self, iteration=1):
13331341 lmp .command (f"pair_style { self .calc ._pair_style_with_options [0 ]} " )
13341342
13351343 # read in conf
1336- # conf = os.path.join(self.simfolder, "conf.equilibration. dump")
1344+ # conf = os.path.join(self.simfolder, "conf.dump")
13371345 conf = os .path .join (self .simfolder , "conf.equilibration.data" )
13381346 lmp = ph .read_data (lmp , conf )
13391347
@@ -1421,12 +1429,12 @@ def temperature_scaling(self, iteration=1):
14211429 lmp .command (
14221430 "fix f2 all npt temp %f %f %f %s %f %f %f"
14231431 % (
1424- tf ,
1432+ t0 ,
14251433 t0 ,
14261434 self .calc .md .thermostat_damping [1 ],
14271435 self .iso ,
1428- pf ,
14291436 p0 ,
1437+ pf ,
14301438 self .calc .md .barostat_damping [1 ],
14311439 )
14321440 )
0 commit comments