Closed zingale closed 3 years ago
@dwillcox can we chat about this together?
I've fixed this in pynucastro. See PR #98 for an updated version of subch.
At the moment, this is only on the development branch of pynucastro.
what about the Urca nets?
it looks like this still affects:
ignition_reaclib/C-test ignition_reaclib/URCA-simple ignition_reaclib/C-burn-simple
When do we use ener_gener_rate
outside of the RHS?
I no longer remember. Probably in one of the unit tests
OK. I chose not to implement ener_gener_rate
explicitly in the C++ port of powerlaw, we can explicitly implement that at some point if we need.
I no longer think that this is an issue.
It seems the pynucastro nets put
ener_gener_rate
inactual_network.f90
, which most others have it inactual_rhs.f90
. This makes it hard to use this routine outside of a net. We should always put it in the same place.