Merge pull request #46 from aidancrilly/Hale_TT_reac

Correct MeV to keV conversion
aidancrilly · Sep 18, 2024 · c1a8ad7 · c1a8ad7
2 parents 8bbb153 + 7036239
commit c1a8ad7
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/NeSST/spectral_model.py b/src/NeSST/spectral_model.py
@@ -248,7 +248,7 @@ def matrix_interpolate_gaussian(self,E,vbar,dv):
 # TT_reac_data = np.loadtxt(data_dir + "TT_reac_ENDF.dat")       # sigmav im m^3/s   # From ENDF
 # TT_reac_spline = interpolate_1d(TT_reac_data[:,0],TT_reac_data[:,1],method='linear',bounds_error=False,fill_value=0.0)
 TT_reac_data = np.loadtxt(data_dir + "TT_reac_Hale.dat")       # T in MeV, sigmav im cm^3/s   # From Hale
-TT_reac_spline = interpolate_1d(TT_reac_data[:,0]*1e-3,TT_reac_data[:,1]*1e-6,method='linear',bounds_error=False,fill_value=0.0)
+TT_reac_spline = interpolate_1d(TT_reac_data[:,0]*1e3,TT_reac_data[:,1]*1e-6,method='linear',bounds_error=False,fill_value=0.0)
 
 ########################
 # Primary reactivities #