On Jul 3, 2010, at 10:01 AM, jingluan wrote:
> Dear Bill,
>
> 1, In the inlist_massive (under work_massive_star), there are statements:
>
> ! at some point, trying to satisfy the residual limits becomes
> hopeless, so you might as
> ! well stop wasting lots of iterations trying. Here's what I did
> after the end of oxygen burning.
> !
> ! ! turn off resid checks early
> ! max_iter_for_resid_tol1 = 3
> ! max_iter_for_resid_tol2 = -1
> ! max_iter_for_resid_tol3 = -1
>
> What is residual limits please?
the "residual" for an equation A = B is the value A - B which would be 0 for a perfect solution.
the newton iteration will not accept a trial solution until the residuals are small enough.
> And where shall one set the value of
> max_iter_for_resid_tol1?
did you try grep to look for it in mesa/star/public?
>
> 2, The nuclear reaction network "basic.net" is used for which range of
> mass? What about "approx21.net"?
It isn't a question of mass range but of stage of burning.
Take a look in mesa/data/net_data/nets to see what isotopes and reactions are in basic.net.
Do the same for approx21.net
Alternatively, in your inlist, change to the net you are interested in and then set
show_net_species_info = .true.
show_net_reactions_info = .true.
>
> 3, Does EZ change nuclear reaction network due to given initial mass please?
EZ only has 1 net that it uses for everything. But I doubt if you are using EZ, right?
mesa/star uses the net you request -- with one exception:
auto_extend_net = .true.
! if true, then automatically extend the basic net as needed
! first adds CO_burn extras and alpha to S32
! later adds rest of alpha chain to Ni56
-B
|