From df1dd30f7ddac253438dfca53f8bc18701b5eee8 Mon Sep 17 00:00:00 2001
From: Ryan Thomas Wollaeger <wollaeger@lanl.gov>
Date: Tue, 7 Jan 2025 18:42:07 -0700
Subject: [PATCH] Roll time unit into mean non-cgs inverse emissivity unit.

---
 singularity-opac/photons/non_cgs_photons.hpp | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/singularity-opac/photons/non_cgs_photons.hpp b/singularity-opac/photons/non_cgs_photons.hpp
index 4d29001..36aff7e 100644
--- a/singularity-opac/photons/non_cgs_photons.hpp
+++ b/singularity-opac/photons/non_cgs_photons.hpp
@@ -248,7 +248,8 @@ class MeanNonCGSUnits {
       : mean_opac_(std::forward<MeanOpac>(mean_opac)), time_unit_(time_unit),
         mass_unit_(mass_unit), length_unit_(length_unit), temp_unit_(temp_unit),
         rho_unit_(mass_unit_ / (length_unit_ * length_unit_ * length_unit_)),
-        inv_emiss_unit_(length_unit_ * time_unit_ * time_unit_ / mass_unit_) {}
+        inv_emiss_unit_(length_unit_ * time_unit_ * time_unit_ * time_unit_ /
+          mass_unit_) {}
 
   auto GetOnDevice() {
     return MeanNonCGSUnits<MeanOpac>(mean_opac_.GetOnDevice(), time_unit_,
@@ -300,7 +301,7 @@ class MeanNonCGSUnits {
                   const int gmode = Rosseland,
                   Real *lambda = nullptr) const {
     const Real J = mean_opac_.Emissivity(rho, temp, gmode);
-    return J * inv_emiss_unit_ * time_unit_;
+    return J * inv_emiss_unit_;
   }
 
   PORTABLE_INLINE_FUNCTION RuntimePhysicalConstants