[Mesa-users] Log_LHe as a function of shell in star

[Francis Timmes]

yes. thanks bill.

fxt



> On Sep 2, 2020, at 11:15 AM, Bill Paxton <paxton at kitp.ucsb.edu> wrote:
> 
> Hi Frank,
> 
> I think we might want to show have profiles based on the eps_nuc_categories info.
> 
> see report.f90 - where it 1st sets L_by_category for entire model
> 
>         do j=1,num_categories
>            s% L_by_category(j) = &
>               dot_product(s% dm(1:nz), s% eps_nuc_categories(j,1:nz))/Lsun
>         end do
> 
> then it stores some special cases
> 
>            s% power_h_burn = s% L_by_category(ipp) + s% L_by_category(icno)
>            s% power_he_burn = s% L_by_category(i3alf)
> 
> for profiles, we might want something like
> 
> 	eps_h_burn(k) = s% eps_nuc_categories(ipp,k) + s% eps_nuc_categories(icno,k)
> 	eps_he_burn(k) = s% eps_nuc_categories(i3alf,k)
> 	etc
> 
> -b
> 
> 
> 
>> On Sep 2, 2020, at 10:41 AM, Francis Timmes via Mesa-users <mesa-users at lists.mesastar.org> wrote:
>> 
>> oops.
>> 
>> fxt
>> 
>> 
>>> On Sep 2, 2020, at 10:38 AM, Josiah Schwab <jwschwab at ucsc.edu> wrote:
>>> 
>>> from $MESA_DIR/star/defaults/profile_columns.list
>>> 
>>> one can put
>>> 
>>> power_h_burn
>>> power_he_burn
>>> 
>>> in the profiles. externally one can then integrate this
>>> quantities over mass to get the LH and LHe.
>>> 
>>> Pretty sure this is incorrect.  Those values are already the integrals and just live in the profile headers.
>>> 
>>> Josiah 
>> 
>> _______________________________________________
>> mesa-users at lists.mesastar.org
>> https://lists.mesastar.org/mailman/listinfo/mesa-users
>> 
>