Final fixes -- all tests passing again!

After hours of sleuthing and @kecnry's help and keen eye, everything's been fixed up. This makes the PR finally viable.
phoebe-project · Mar 5, 2025 · c3e4f3d · c3e4f3d
1 parent 7ffe0f6
commit c3e4f3d
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 4 deletions.
diff --git a/phoebe/atmospheres/passbands.py b/phoebe/atmospheres/passbands.py
@@ -572,8 +572,9 @@ def load(cls, archive, load_content=True):
                         # store axes:
                         hdul['bb_teffs'] = fits.table_to_hdu(Table({'teffs': ndp.axes[0]}, meta={'extname': 'bb_teffs'}))
                         if 'blackbody:ext' in self.content:
-                            hdul['bb_ebvs'] = fits.table_to_hdu(Table({'ebvs': ndp.axes[1]}, meta={'extname': 'bb_ebvs'}))
-                            hdul['bb_rvs'] = fits.table_to_hdu(Table({'rvs': ndp.axes[2]}, meta={'extname': 'bb_rvs'}))
+                            associated_axes = ndp.table['ext@photon']['associated_axes']
+                            hdul['bb_ebvs'] = fits.table_to_hdu(Table({'ebvs': associated_axes[0]}, meta={'extname': 'bb_ebvs'}))
+                            hdul['bb_rvs'] = fits.table_to_hdu(Table({'rvs': associated_axes[1]}, meta={'extname': 'bb_rvs'}))
 
                         # store tables:
                         hdul.append(fits.ImageHDU(self.ndp['blackbody'].table['inorm@energy']['grid'], name='bbnegrid'))
@@ -625,7 +626,8 @@ def load(cls, archive, load_content=True):
                     if atm.external:
                         continue
 
-                    basic_axes = tuple(np.asarray(hdul[f'{atm.prefix}_{name}'].data[name]) for name in atm.basic_axis_names)
+                    basic_axes = tuple([np.array(list(hdul[f'{atm.prefix}_{name}'].data[name])) for name in atm.basic_axis_names])
+                    # basic_axes = tuple(np.asarray(hdul[f'{atm.prefix}_{name}'].data[name]) for name in atm.basic_axis_names)
                     self.ndp[atm.name] = ndpolator.Ndpolator(basic_axes=basic_axes)
 
                     if f'{atm.name}:Inorm' in self.content:

diff --git a/phoebe/backend/universe.py b/phoebe/backend/universe.py
@@ -338,7 +338,7 @@ def handle_reflection(self,  **kwargs):
         for starref, body in self.items():
             if body.mesh is None:
                 continue
-            abs_normal_intensities = sigma_sb * body.mesh.teffs.for_computations**4 / np.pi  # bolometric intensities
+            abs_normal_intensities = sigma_sb.value * body.mesh.teffs.for_computations**4 / np.pi  # bolometric intensities
             fluxes_intrins_per_body.append(abs_normal_intensities * np.pi)
 
         fluxes_intrins_flat = meshes.pack_column_flat(fluxes_intrins_per_body)