Source `apar_m` not weighted by `fapar`

Issue #103 new
David Dickinson created an issue

In get_source_term and related code we weight the potentials with fphi, fapar and fbpar as appropriate. We do not currently weight apar_m in set_source which uses apar directly (this is the dApar/dt term) so this seems inconsistent for fapar/=0 or fapar/=1.

Comments (0)

  1. Log in to comment