[mesa-users] dxdt_nuc_factor in report.f90

[Nathan Grin]

Hello,

while playing around i found that when setting
dxdt_nuc_factor = 0 ! shut down composition changes due to nuclear burning
in the controls inlist, consequently s% L_nuc_burn_total, among others, is
set to 0, which i didnt understand because as one shuts off composition
changes, the power due to burning should still be there (and it is, one can
access s% eps_nuc(:) as in the code below)

It traces back to the following piece of code in star/private/report.f90
(mesa-r8845)

         if (s% dxdt_nuc_factor > 0d0) then
            s% power_nuc_burn = dot_product(s% dm(1:nz), s%
eps_nuc(1:nz))/Lsun
            s% power_h_burn = s% L_by_category(ipp) + s% L_by_category(icno)
            s% power_he_burn = s% L_by_category(i3alf)
            s% power_c_burn = s% L_by_category(i_burn_c)
            s% power_nuc_neutrinos = &
               dot_product(s% dm(1:nz),s% eps_nuc_neu_total(1:nz))/Lsun
            s% power_nonnuc_neutrinos = &
               dot_product(s% dm(1:nz),s% non_nuc_neu(1:nz))/Lsun
            s% power_neutrinos = s% power_nuc_neutrinos + s%
power_nonnuc_neutrinos
         else
            s% power_nuc_burn = 0d0
            s% power_h_burn = 0d0
            s% power_he_burn = 0d0
            s% power_c_burn = 0d0
            s% power_nuc_neutrinos = 0d0
            s% power_nonnuc_neutrinos = 0d0
            s% power_neutrinos = 0d0
         end if

         s% total_radiation = s% total_radiation + &
            s% dt*(luminosity + s% power_neutrinos)

         s% L_nuc_burn_total = s% power_nuc_burn

so, i was wondering:
perhaps s% dxdt_nuc_factor > 0d0 should be replaced by s% eps_nuc_factor >
0d0 ?
or maybe the if-statement is not needed, because all  statements in there
should evaluate to 0 as you set eps_nuc_factor=0 anyhow (except for the
neutrino powers maybe?)
or, ofcourse, am i missing something? :)

cheers!
nathan
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20160905/b6f120e9/attachment.html>