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

[Bill Paxton]

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
>