Skip to content

Commit

Permalink
Fixed bug to generalize metal used
Browse files Browse the repository at this point in the history
  • Loading branch information
tdprice-858 committed Dec 9, 2024
1 parent 58b33d1 commit 7574d03
Showing 1 changed file with 6 additions and 2 deletions.
8 changes: 6 additions & 2 deletions pynta/postprocessing.py
Original file line number Diff line number Diff line change
Expand Up @@ -569,6 +569,8 @@ def __init__(self,
dT = 10.0 # temperature increment
self.temperature = np.append(temperature, np.arange(T_low, T_high + dT, dT))

self.cat_element = cat_element

def compute_heat_of_formation(self):
'''asdfas'''
# What do we do for other atomic adsorbed species?
Expand Down Expand Up @@ -718,6 +720,7 @@ def _get_vibrational_thermo(self):
return

def fit_NASA(self):
'''I have not yet check the math on all this stuff from Goldsmith group'''
R = self.R
heat_capacity = self.Cp
reference_enthalpy = self.H[0]
Expand Down Expand Up @@ -832,8 +835,7 @@ def fit_NASA(self):
rmg_line += "\t\t Tmin = (%.1F, 'K'),\n" % (298.0)
rmg_line += "\t\t Tmax = (%.1F, 'K'),\n" % (max(temperature))
rmg_line += "\t),\n"
# SHINAE CHANGE TO PT
rmg_line += '\tmetal = "Ag",\n'
rmg_line += f'\tmetal = {self.cat_element},\n'
rmg_line += '\tfacet = "111",\n'

self.thermo_lines = rmg_line
Expand All @@ -843,5 +845,7 @@ def fit_NASA(self):

return





0 comments on commit 7574d03

Please sign in to comment.