[mesa-users] dxdt_nuc_factor in report.f90

Josiah Schwab jwschwab at berkeley.edu
Mon Sep 5 11:54:59 EDT 2016


Hi Nathan,

> 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? :)

Thanks for digging into the code!  You set a great example for curious
users everywhere.

I agree with you that the if statement is not needed at all.  I will
make that change.

Thanks,
Josiah




More information about the Mesa-users mailing list