Show time days

artis-mcrt · Sep 12, 2024 · b33f733 · b33f733
1 parent f474fd8
commit b33f733
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 5 deletions.
diff --git a/gammapkt.cc b/gammapkt.cc
@@ -688,8 +688,7 @@ void update_gamma_dep(const Packet &pkt, const double dist, const int mgi, const
   assert_testmodeonly(std::isfinite(heating_cont));
   atomicadd(globals::dep_estimator_gamma[nonemptymgi], heating_cont);
 
-  const double kappa_E_cont = heating_cont * get_kappa(pkt);
-  atomicadd(globals::estimator_gamma_kappa_integrated, kappa_E_cont);
+  atomicadd(globals::estimator_gamma_kappa_integrated, heating_cont * get_kappa(pkt));
 }
 
 // handle physical pair production event

diff --git a/sn3d.cc b/sn3d.cc
@@ -659,9 +659,10 @@ void normalise_deposition_estimators(int nts) {
   const double mean_en_mev = H * globals::estimator_gamma_nu_cmf_decayspec / MEV;
 
   printout(
-      "timestep %d kappa_eff [cm2/g]: path-integrated %g decay_spec %g deposit_spec %g decay gamma mean E[MeV] %g\n",
-      nts, globals::estimator_gamma_kappa_integrated, globals::estimator_gamma_kappa_decayspec,
-      globals::estimator_gamma_kappa_absorbedspec, mean_en_mev);
+      "timestep %d %g days kappa_eff [cm2/g]: path-integrated %g decay_spec %g deposit_spec %g decay gamma mean E[MeV] "
+      "%g\n",
+      nts, globals::timesteps[nts].mid / DAY, globals::estimator_gamma_kappa_integrated,
+      globals::estimator_gamma_kappa_decayspec, globals::estimator_gamma_kappa_absorbedspec, mean_en_mev);
 
   for (int nonemptymgi = 0; nonemptymgi < grid::get_nonempty_npts_model(); nonemptymgi++) {
     const int mgi = grid::get_mgi_of_nonemptymgi(nonemptymgi);