[Mesa-users] epoch-dependent G
Marco Pereira
ny2292000 at gmail.com
Sun Oct 9 13:29:02 UTC 2022
Dear Francis and Users,
I sought files where the word grav was mentioned. My goal is to just add a
time-dependent G to the evolution of Stars.
I found do_evolve_one_step method in run_star_support8.f90. I suspect that
that is where I should make my modifications.
I also know that grav will affect star winds and star oscillations
frequencies (gyre).
In addition, to star evolution, I recognize the study of binary systems in
the package is also relevant since I did simple modeling of type 1a
Supernovae to get their Absolute Luminosity G-dependence as G^{-3.33}
Attached is the full list of mentions of the word "grav" and the files
where they are located.
My naive expectation was that there would be a nuclear chemistry propagator
(where the nuclear chemistry rate equations would be propagated step by
step taking into consideration density and composition). I thought that
that would be connected to the evolve_one_step method. The other relevant
point would be the calculation of the star winds since they would push away
the feeding hydrogen cloud. That should account for a smaller initial mass
and larger G. That is relevant because in my model the initial mass of the
sun is around 47% of the current mass and the pushed-away gas cloud would
fall back into the star over 4 billion years.
So, my model adds a falling hydrogen cloud to the modeling of the sun.
Any comments would be welcomed. At this exact moment, I am still incapable
of stepping (debugging) through this Fortran program and so I am not quite
sure where it goes. In fact, I have CLion and code::blocks. How do people
debug this program?
Thanks,
Marco
On Tue, Oct 4, 2022 at 7:45 PM Marco Pereira <ny2292000 at gmail.com> wrote:
> The hydrogen density would also be around 10 times higher.
>
> On Tue, Oct 4, 2022 at 7:01 PM Marco Pereira <ny2292000 at gmail.com> wrote:
>
>> Dear Francis and Mesa-Users,
>>
>> In the past, it seems that epoch-dependent G investigated only small
>> variations.
>> Article a) tried to test a variable G but did so while trying to keep G
>> constant. They didn't explore my option.
>> Article b) just proposed that the Sun was more massive but lost mass
>> through solar winds (a diet).
>> So, what I am proposing is distinct and hasn't been done yet, as far as I
>> can tell.
>>
>> a) Can a variable gravitational constant resolve the Faint Young Sun
>> Paradox?
>> Varun Sahni, Yuri Shtanov
>> https://arxiv.org/abs/1405.4369
>>
>>
>> *b) Assessing the massive young Sun hypothesis to solve the warm young
>> Earth puzzleDavid A. Minton, Renu Malhotra*
>> https://arxiv.org/abs/astro-ph/0612321
>>
>> #######################################################
>> I would like to investigate large variations of G.
>>
>> I want to calculate G to be inversely proportional to the 4D radius of
>> the Universe (inverse of the epoch). This means that G would be 1.5 the
>> current G0. That would work if the Sun started as 46% of the current Sun
>> mass and accreted 9% every billion years.
>> Of course, the whole galaxy would also have a radius that was 0.46 the
>> current radius. Angular momentum conservation would have the sun rotating
>> twice as fast around the galaxy. The hydrogen density would also be around
>> 300% higher.
>>
>> So, I would like to increase the rate at which hydrogen is accreted into
>> the sun, have a variable derivative for G, and start with a smaller mass.
>>
>> So, that is what I would like to try. This is the back-of-envelope
>> calculation. We all know that the sun loses mass through solar wind. The
>> sun should gain mass from falling hydrogen.
>> [image: image.png]
>>
>> I will dig deeper into your suggestions. If anyone wants to help me do
>> the calculation, I would be more than grateful.
>>
>> Cheers,
>>
>> Marco Pereira
>>
>>
>> On Tue, Oct 4, 2022 at 5:00 PM Francis Timmes <fxt44 at mac.com> wrote:
>>
>>> hi marco,
>>>
>>> the topic of variable G arises every now and then. one example is
>>>
>>> https://lists.mesastar.org/pipermail/mesa-users/2022-February/013504.html
>>>
>>> and
>>>
>>>
>>> https://docs.mesastar.org/en/release-r22.05.1/using_mesa/extending_mesa.html#using-the-other-hooks
>>>
>>> https://docs.mesastar.org/en/release-r22.05.1/
>>>
>>> seem relevant.
>>>
>>> fxt
>>>
>>>
>>>
>>>
>>>
>>>
>>> > On Oct 4, 2022, at 12:56 PM, Marco Pereira via Mesa-users <
>>> mesa-users at lists.mesastar.org> wrote:
>>> >
>>> > Hi,
>>> >
>>> > I installed MESA-r22.05.1 and its SDK. I finally was able to install
>>> the program and create the Dox folder using doxygen.
>>> > So, I am ready to try to understand the program. That said, what I
>>> want to do might be quite localized.
>>> > I want to add into the integration step the knowledge of the value of
>>> G (Newton's Gravitational Constant) at that time.
>>> >
>>> > Any guidance would be much appreciated. Also, what I have (oxygen
>>> HTML files) doesn't see to be documented enough for me to understand the
>>> program with just that.
>>> >
>>> > Where can I find broader information on what each module and each file
>>> does?
>>> >
>>> > That said, my thirst for knowledge might not be relevant if the
>>> epoch-dependent G can be surgically added.
>>> >
>>> > Thank you in advance for your help,
>>> >
>>> > Marco
>>> > _______________________________________________
>>> > mesa-users at lists.mesastar.org
>>> > https://lists.mesastar.org/mailman/listinfo/mesa-users
>>> >
>>>
>>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20221009/ff4ab30d/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 30991 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20221009/ff4ab30d/attachment.png>
-------------- next part --------------
./atm/private/atm_utils.f90:94: subroutine eval_Teff_g(L, R, M, cgrav, Teff, g)
./atm/private/atm_utils.f90:99: real(dp), intent(in) :: cgrav
./atm/private/atm_utils.f90:103: ! Evaluate the effective temperature and surface gravity
./atm/private/atm_utils.f90:107: g = cgrav * M / (R*R)
./atm/private/atm_utils.f90:114: T, P, rho, tau, kap, L, M, R, cgrav) result (gradr)
./atm/private/atm_utils.f90:126: real(dp), intent(in) :: cgrav
./atm/private/atm_utils.f90:139: gradr = P*kap*L / (16._dp*pi*clight*M*cgrav*Prad)
./atm/private/atm_utils.f90:145: s = (2._dp*crad*T*T*T*SQRT(R))/(3._dp*cgrav*M*rho)*pow(L/(8._dp*pi*boltz_sigma), 0.25_dp)
./atm/private/atm_utils.f90:151: dilution_factor = (1._dp + f*s*(4._dp*pi*cgrav*clight*M)/(kap*L))/(1._dp + f*s)
./atm/private/atm_irradiated.f90:51: L, R, M, cgrav, T_eq, P_surf, kap_guess, kap_v, gamma, &
./atm/private/atm_irradiated.f90:64: real(dp), intent(in) :: cgrav
./atm/private/atm_irradiated.f90:112: ! Evaluate the 'interior' temperature & gravity
./atm/private/atm_irradiated.f90:114: call eval_Teff_g(L, R, M, cgrav, T_int, g)
./atm/private/atm_t_tau_varying.f90:53: tau_surf, L, R, M, cgrav, &
./atm/private/atm_t_tau_varying.f90:67: real(dp), intent(in) :: cgrav
./atm/private/atm_t_tau_varying.f90:113: ! Evaluate the gravity
./atm/private/atm_t_tau_varying.f90:115: g = cgrav*M/(R*R)
./atm/private/atm_t_tau_varying.f90:518: tau_surf, L, R, Teff, M, cgrav, lnP_surf, tau_outer, &
./atm/private/atm_t_tau_varying.f90:532: real(dp), intent(in) :: cgrav
./atm/private/atm_t_tau_varying.f90:586: ! Evaluate the gravity
./atm/private/atm_t_tau_varying.f90:588: g = cgrav*M/(R*R)
./atm/private/atm_t_tau_varying.f90:828: gradr = eval_Paczynski_gradr(exp(lnT), exp(lnP), exp(lnRho), tau, kap, L, M, R, cgrav)
./atm/private/atm_t_tau_uniform.f90:54: tau_surf, L, R, M, cgrav, kap_guess, Pextra_factor, &
./atm/private/atm_t_tau_uniform.f90:68: real(dp), intent(in) :: cgrav
./atm/private/atm_t_tau_uniform.f90:124: ! Evaluate the gravity
./atm/private/atm_t_tau_uniform.f90:126: g = cgrav*M/(R*R)
./atm/private/atm_t_tau_uniform.f90:134: tau_surf, Teff, g, L, M, cgrav, &
./atm/private/atm_t_tau_uniform.f90:194: tau_surf, Teff, g, L, M, cgrav, &
./atm/private/atm_t_tau_uniform.f90:264: tau_surf, L, R, Teff, M, cgrav, kap, Pextra_factor, tau_outer, &
./atm/private/atm_t_tau_uniform.f90:278: real(dp), intent(in) :: cgrav
./atm/private/atm_t_tau_uniform.f90:328: ! Evaluate the gravity
./atm/private/atm_t_tau_uniform.f90:330: g = cgrav*M/(R*R)
./atm/private/atm_t_tau_uniform.f90:538: tau, Teff, g, L, M, cgrav, &
./atm/private/atm_t_tau_uniform.f90:561: gradr = eval_Paczynski_gradr(exp(lnT), exp(lnP), exp(lnRho), tau, kap, L, M, R, cgrav)
./atm/private/atm_t_tau_uniform.f90:594: tau, Teff, g, L, M, cgrav, &
./atm/private/atm_t_tau_uniform.f90:607: real(dp), intent(in) :: cgrav
./atm/private/atm_t_tau_uniform.f90:653: Pextra = Pextra_factor*(kap/tau)*(L/M)/(6._dp*pi*clight*cgrav)
./atm/private/atm_table.f90:52: L, R, M, cgrav, id, Z, skip_partials, &
./atm/private/atm_table.f90:66: real(dp), intent(in) :: cgrav
./atm/private/atm_table.f90:122: ! Evaluate the gravity
./atm/private/atm_table.f90:124: g = cgrav*M/(R*R)
./atm/private/atm_table.f90:216: L, Teff, R, M, cgrav, id, alfa, beta, ierr)
./atm/private/atm_table.f90:226: real(dp), intent(in) :: cgrav
./atm/private/atm_table.f90:280: ! Evaluate the gravity
./atm/private/atm_table.f90:282: g = cgrav*M/(R*R)
./stella/src/stl/phys_constants.f90:32: ! *** Newtonian constant of gravitation G [CGS] ***
./stella/res/stella_extras.f90:354: log_g = log10(standard_cgrav*mphot*Msun/(rphot*rphot))
./star/public/star_lib.f90:689: subroutine star_build_atm(s, L, R, Teff, M, cgrav, ierr)
./star/public/star_lib.f90:693: real(dp), intent(in) :: L, R, Teff, M, cgrav
./star/public/star_lib.f90:695: call build_atm(s, L, R, Teff, M, cgrav, ierr)
./star/public/star_lib.f90:2220: id, tau_surf, L, R, M, cgrav, skip_partials, Teff, &
./star/public/star_lib.f90:2226: real(dp), intent(in) :: tau_surf, L, R, M, cgrav
./star/public/star_lib.f90:2238: s, tau_surf, L, R, M, cgrav, skip_partials, &
./star/public/star_lib.f90:2314: total_internal_energy, total_gravitational_energy, &
./star/public/star_lib.f90:2320: total_internal_energy, total_gravitational_energy, &
./star/public/star_lib.f90:2328: total_internal_energy, total_gravitational_energy, &
./star/public/star_lib.f90:3063: iso, XH1, cgrav, m, gradL_composition_term, mixing_length_alpha, &
./star/public/star_lib.f90:3071: XH1, cgrav, m, gradL_composition_term, mixing_length_alpha
./star/public/star_lib.f90:3080: iso, XH1, cgrav, m, gradL_composition_term, mixing_length_alpha, &
./star/public/star_lib.f90:3087: iso, XH1, cgrav, m, gradL_composition_term, mixing_length_alpha, &
./star/public/star_lib.f90:3100: XH1, cgrav, m, gradL_composition_term, &
./star/public/star_lib.f90:3113: iso, XH1, cgrav, m, gradL_composition_term, mixing_length_alpha, &
./star/job/run_star_support.f90:703: tmp = abs(1d0 + s% total_gravitational_energy_end/s% virial_thm_P_avg)
./star/job/run_star_support.f90:707: abs(s% total_gravitational_energy_end), s% virial_thm_P_avg
./star/job/run_star_support.f90:1363: log_g = safe_log10(s% grav(1))
./star/job/run_star_support.f90:1519: write(*,1) 'g', s% cgrav(1)*s% mstar/(s% r(1)*s% r(1))
./star/job/run_star_support.f90:2780: if (s% eps_grav_factor /= 1) &
./star/job/run_star_support.f90:2781: write(*,1) 'eps_grav_factor', s% eps_grav_factor
./star/private/brunt.f90:146: ! correct for difference between gravitational mass density and baryonic mass density (rho)
./star/private/brunt.f90:168: f = pow2(s% grav(k))*rho_P_chiT_chiRho_face(k)
./star/private/brunt.f90:311: dlnP_dm = -s% cgrav(k)*s% m(k)/(pi4*pow4(s% r(k))*Ppoint)
./star/private/brunt.f90:317: ! add term accounting for the composition-related gradient in gravitational mass
./star/private/eps_grav.f90:26: module eps_grav
./star/private/eps_grav.f90:37: public :: eval_eps_grav_and_partials, zero_eps_grav_and_partials
./star/private/eps_grav.f90:42: subroutine eval_eps_grav_and_partials(s, k, ierr)
./star/private/eps_grav.f90:46: type(auto_diff_real_star_order1) :: eps_grav
./star/private/eps_grav.f90:55: call zero_eps_grav_and_partials(s, k)
./star/private/eps_grav.f90:58: ! if include_composition_in_eps_grav is true
./star/private/eps_grav.f90:59: ! then these are set in the call to eval_eps_grav_composition
./star/private/eps_grav.f90:60: s% d_eps_grav_dx(:,k) = 0
./star/private/eps_grav.f90:62: call eval1_eps_grav_and_partials(s, k, eps_grav, ierr)
./star/private/eps_grav.f90:65: s% eps_grav_ad(k) = eps_grav
./star/private/eps_grav.f90:67: if (s% use_other_eps_grav) then
./star/private/eps_grav.f90:69: call s% other_eps_grav(s% id, k, s% dt, ierr)
./star/private/eps_grav.f90:74: s% eps_grav_ad(k) = s% eps_grav_factor * s% eps_grav_ad(k)
./star/private/eps_grav.f90:76: end subroutine eval_eps_grav_and_partials
./star/private/eps_grav.f90:79: subroutine eval1_eps_grav_and_partials(s, k, eps_grav, ierr)
./star/private/eps_grav.f90:82: type(auto_diff_real_star_order1), intent(out) :: eps_grav
./star/private/eps_grav.f90:85: type(auto_diff_real_star_order1) :: eps_grav_lnS, eps_grav_std
./star/private/eps_grav.f90:94: if (using_PC .and. s% gam_start(k) >= s% Gamma_lnS_eps_grav_full_on) then
./star/private/eps_grav.f90:95: call do_lnS_eps_grav(s, k, eps_grav, ierr)
./star/private/eps_grav.f90:96: else if (using_PC .and. s% gam_start(k) > s% Gamma_lnS_eps_grav_full_off) then
./star/private/eps_grav.f90:98: alfa = (Gamma - s% Gamma_lnS_eps_grav_full_off) / &
./star/private/eps_grav.f90:99: (s% Gamma_lnS_eps_grav_full_on - s% Gamma_lnS_eps_grav_full_off)
./star/private/eps_grav.f90:100: call do_lnS_eps_grav(s, k, eps_grav_lnS, ierr)
./star/private/eps_grav.f90:102: call do_std_eps_grav(s, k, eps_grav_std, ierr)
./star/private/eps_grav.f90:107: eps_grav = alfa * eps_grav_lnS + (1d0 - alfa) * eps_grav_std
./star/private/eps_grav.f90:109: call do_std_eps_grav(s, k, eps_grav, ierr)
./star/private/eps_grav.f90:113: if (ierr /= 0 .or. is_bad(eps_grav% val)) then
./star/private/eps_grav.f90:115: s% retry_message = 'failed in eval_eps_grav_and_partials'
./star/private/eps_grav.f90:118: 'failed in eval_eps_grav_and_partials', k, eps_grav% val
./star/private/eps_grav.f90:121: call mesa_error(__FILE__,__LINE__,'eval1_eps_grav_and_partials')
./star/private/eps_grav.f90:126: end subroutine eval1_eps_grav_and_partials
./star/private/eps_grav.f90:130: subroutine do_std_eps_grav(s, k, eps_grav, ierr)
./star/private/eps_grav.f90:135: type(auto_diff_real_star_order1), intent(out) :: eps_grav
./star/private/eps_grav.f90:139: latent_ddlnT, latent_ddlnRho, eps_grav_start, eps_grav_composition_term
./star/private/eps_grav.f90:152: if (s% use_time_centered_eps_grav) then
./star/private/eps_grav.f90:170: eps_grav = -T*Cp * ((1d0 - grada*chiT)*dlnT_dt - grada*chiRho*dlnd_dt)
./star/private/eps_grav.f90:176: eps_grav = eps_grav - (dlnd_dt * latent_ddlnRho + dlnT_dt * latent_ddlnT)
./star/private/eps_grav.f90:180: if (s% use_time_centered_eps_grav) then
./star/private/eps_grav.f90:183: eps_grav_start = -s% T_start(k)*s% cp_start(k) * ((1d0 - s% grada_start(k)*s% chiT_start(k))*dlnT_dt - s% grada_start(k)*s% chiRho_start(k)*dlnd_dt)
./star/private/eps_grav.f90:186: eps_grav_start = eps_grav_start - (dlnd_dt * s% latent_ddlnRho_start(k) + dlnT_dt * s% latent_ddlnT_start(k))
./star/private/eps_grav.f90:188: eps_grav = theta * eps_grav + (1d0-theta) * eps_grav_start
./star/private/eps_grav.f90:193: if (s% include_composition_in_eps_grav) then
./star/private/eps_grav.f90:194: call eval_eps_grav_composition(s, k, eps_grav_composition_term, ierr)
./star/private/eps_grav.f90:196: eps_grav = eps_grav + eps_grav_composition_term
./star/private/eps_grav.f90:199: if (is_bad(eps_grav% val)) then
./star/private/eps_grav.f90:201: s% retry_message = 'do_lnd_eps_grav -- bad value for eps_grav'
./star/private/eps_grav.f90:203: write(*,2) 'do_lnd_eps_grav -- bad value for eps_grav', k, eps_grav% val
./star/private/eps_grav.f90:204: if (s% stop_for_bad_nums) call mesa_error(__FILE__,__LINE__,'do_lnd_eps_grav')
./star/private/eps_grav.f90:212: write(*,*) 'do_std_eps_grav chiT', s% solver_test_partials_var
./star/private/eps_grav.f90:215: end subroutine do_std_eps_grav
./star/private/eps_grav.f90:218: subroutine do_lnS_eps_grav(s, k, eps_grav, ierr)
./star/private/eps_grav.f90:222: type(auto_diff_real_star_order1), intent(out) :: eps_grav
./star/private/eps_grav.f90:226: type(auto_diff_real_star_order1) :: entropy, T, eps_grav_composition_term
./star/private/eps_grav.f90:235: eps_grav = -T*(entropy - entropy_start)/s% dt
./star/private/eps_grav.f90:250: ! if (s% include_composition_in_eps_grav) then
./star/private/eps_grav.f90:251: ! call eval_eps_grav_composition(s, k, eps_grav_composition_term, ierr)
./star/private/eps_grav.f90:253: ! eps_grav = eps_grav + eps_grav_composition_term
./star/private/eps_grav.f90:256: if (is_bad(eps_grav% val)) then
./star/private/eps_grav.f90:258: s% retry_message = 'do_lnS_eps_grav -- bad value for eps_grav'
./star/private/eps_grav.f90:260: write(*,2) 'do_lnS_eps_grav -- bad value for eps_grav', k, eps_grav% val
./star/private/eps_grav.f90:261: if (s% stop_for_bad_nums) call mesa_error(__FILE__,__LINE__,'do_lnS_eps_grav')
./star/private/eps_grav.f90:265: end subroutine do_lnS_eps_grav
./star/private/eps_grav.f90:268: subroutine eval_eps_grav_composition(s, k, eps_grav_composition_term, ierr)
./star/private/eps_grav.f90:275: type(auto_diff_real_star_order1), intent(out) :: eps_grav_composition_term
./star/private/eps_grav.f90:292: eps_grav_composition_term = 0
./star/private/eps_grav.f90:294: if (s% use_time_centered_eps_grav) then
./star/private/eps_grav.f90:307: ! eps_grav = -de/dt + (P/rho) * dlnd/dt
./star/private/eps_grav.f90:309: s% d_eps_grav_dx(j,k) = -s% energy(k) * s% dlnE_dxa_for_partials(j,k)/s% dt + &
./star/private/eps_grav.f90:332: if (s% use_time_centered_eps_grav) then
./star/private/eps_grav.f90:359: ! s% d_eps_grav_dx(j,k) = theta * s% d_eps_grav_dx(j,k) + (1d0 - theta) * &
./star/private/eps_grav.f90:366: call wrap(eps_grav_composition_term, -de/s% dt, &
./star/private/eps_grav.f90:380: s% eps_grav_composition_term(k) = eps_grav_composition_term% val
./star/private/eps_grav.f90:392: if (is_bad(eps_grav_composition_term% val)) then
./star/private/eps_grav.f90:395: write(*,2) 'eps_grav_composition_term', k, eps_grav_composition_term% val
./star/private/eps_grav.f90:396: !call mesa_error(__FILE__,__LINE__,'eval_eps_grav_composition')
./star/private/eps_grav.f90:399: write(*,2) 'include_composition_in_eps_grav -- bad value for eps_grav_composition_term', k, eps_grav_composition_term% val
./star/private/eps_grav.f90:400: call mesa_error(__FILE__,__LINE__,'eval_eps_grav_composition')
./star/private/eps_grav.f90:405: end subroutine eval_eps_grav_composition
./star/private/eps_grav.f90:408: subroutine zero_eps_grav_and_partials(s, k)
./star/private/eps_grav.f90:411: s% eps_grav_ad(k) = 0
./star/private/eps_grav.f90:412: end subroutine zero_eps_grav_and_partials
./star/private/eps_grav.f90:414: end module eps_grav
./star/private/conv_premix.f90:1011: z_s = MIN(zi%vc_b*zi%vc_b/(2._dp*s%grav(kf)*(1._dp - mu_i/mu_e)), s%mlt_mixing_length(kf))
./star/private/conv_premix.f90:1036: z_s = MIN(zi%vc_t*zi%vc_t/(2._dp*s%grav(kf)*(1._dp - mu_e/mu_i)), s%mlt_mixing_length(kf))
./star/private/star_profile_def.f90:64: integer, parameter :: p_grav = p_log_g + 1
./star/private/star_profile_def.f90:65: integer, parameter :: p_r_div_g = p_grav + 1
./star/private/star_profile_def.f90:115: integer, parameter :: p_m_grav = p_dm_bar + 1
./star/private/star_profile_def.f90:116: integer, parameter :: p_mass_correction_factor = p_m_grav + 1
./star/private/star_profile_def.f90:117: integer, parameter :: p_m_grav_div_m_baryonic = p_mass_correction_factor + 1
./star/private/star_profile_def.f90:119: integer, parameter :: p_xr = p_m_grav_div_m_baryonic + 1
./star/private/star_profile_def.f90:375: integer, parameter :: p_cgrav_factor = p_alpha_mlt + 1
./star/private/star_profile_def.f90:376: integer, parameter :: p_log_extra_L = p_cgrav_factor + 1
./star/private/star_profile_def.f90:380: integer, parameter :: p_extra_grav = p_extra_omegadot + 1
./star/private/star_profile_def.f90:381: integer, parameter :: p_extra_heat = p_extra_grav + 1
./star/private/star_profile_def.f90:386: integer, parameter :: p_dvdt_grav = p_d_v_div_r_dr + 1
./star/private/star_profile_def.f90:387: integer, parameter :: p_dvdt_dPdm = p_dvdt_grav + 1
./star/private/star_profile_def.f90:397: integer, parameter :: p_dPdr_div_grav = p_gradP_div_rho + 1
./star/private/star_profile_def.f90:398: integer, parameter :: p_env_eps_grav = p_dPdr_div_grav + 1
./star/private/star_profile_def.f90:400: integer, parameter :: p_log_abs_eps_grav_dm_div_L = p_env_eps_grav + 1
./star/private/star_profile_def.f90:401: integer, parameter :: p_eps_grav_composition_term = p_log_abs_eps_grav_dm_div_L + 1
./star/private/star_profile_def.f90:403: integer, parameter :: p_eps_grav_plus_eps_mdot = p_eps_grav_composition_term + 1
./star/private/star_profile_def.f90:404: integer, parameter :: p_ergs_eps_grav_plus_eps_mdot = p_eps_grav_plus_eps_mdot + 1
./star/private/star_profile_def.f90:406: integer, parameter :: p_ergs_mdot = p_ergs_eps_grav_plus_eps_mdot + 1
./star/private/star_profile_def.f90:409: integer, parameter :: p_dm_eps_grav = p_eps_mdot + 1
./star/private/star_profile_def.f90:411: integer, parameter :: p_log_xm_div_delta_m = p_dm_eps_grav + 1
./star/private/star_profile_def.f90:413: integer, parameter :: p_eps_grav = p_xm_div_delta_m + 1
./star/private/star_profile_def.f90:414: integer, parameter :: p_signed_log_eps_grav = p_eps_grav + 1
./star/private/star_profile_def.f90:415: integer, parameter :: p_log_conv_L_div_L = p_signed_log_eps_grav + 1
./star/private/star_profile_def.f90:755: profile_column_name(p_grav) = 'grav'
./star/private/star_profile_def.f90:817: profile_column_name(p_m_grav) = 'm_grav'
./star/private/star_profile_def.f90:819: profile_column_name(p_m_grav_div_m_baryonic) = 'm_grav_div_m_baryonic'
./star/private/star_profile_def.f90:1066: profile_column_name(p_extra_grav) = 'extra_grav'
./star/private/star_profile_def.f90:1069: profile_column_name(p_cgrav_factor) = 'cgrav_factor'
./star/private/star_profile_def.f90:1074: profile_column_name(p_dvdt_grav) = 'dvdt_grav'
./star/private/star_profile_def.f90:1083: profile_column_name(p_dPdr_div_grav) = 'dPdr_div_grav'
./star/private/star_profile_def.f90:1085: profile_column_name(p_log_abs_eps_grav_dm_div_L) = 'log_abs_eps_grav_dm_div_L'
./star/private/star_profile_def.f90:1086: profile_column_name(p_eps_grav_composition_term) = 'eps_grav_composition_term'
./star/private/star_profile_def.f90:1087: profile_column_name(p_dm_eps_grav) = 'dm_eps_grav'
./star/private/star_profile_def.f90:1089: profile_column_name(p_eps_grav_plus_eps_mdot) = 'eps_grav_plus_eps_mdot'
./star/private/star_profile_def.f90:1090: profile_column_name(p_ergs_eps_grav_plus_eps_mdot) = 'ergs_eps_grav_plus_eps_mdot'
./star/private/star_profile_def.f90:1096: profile_column_name(p_eps_grav) = 'eps_grav'
./star/private/star_profile_def.f90:1097: profile_column_name(p_env_eps_grav) = 'env_eps_grav'
./star/private/star_profile_def.f90:1098: profile_column_name(p_signed_log_eps_grav) = 'signed_log_eps_grav'
./star/private/mix_info.f90:1132: Hp = s% Peos(ktop)/(s% rho(ktop)*s% grav(ktop))
./star/private/mix_info.f90:1194: Hp = s% Peos(ktop)/(s% rho(ktop)*s% grav(ktop))
./star/private/mix_info.f90:2399: if (s% cgrav(k) <= 0) then
./star/private/mix_info.f90:2402: Hp = P_face/(rho_face*s% cgrav(k)* &
./star/private/ctrls_io.f90:148: fix_eps_grav_transition_to_grid, make_gradr_sticky_in_solver_iters, min_logT_for_make_gradr_sticky_in_solver_iters, &
./star/private/ctrls_io.f90:350: use_time_centered_eps_grav, &
./star/private/ctrls_io.f90:351: use_mass_corrections, use_gravity_rotation_correction, eps_grav_factor, eps_mdot_factor, &
./star/private/ctrls_io.f90:352: include_composition_in_eps_grav, no_dedt_form_during_relax, &
./star/private/ctrls_io.f90:354: max_num_surf_revisions, Gamma_lnS_eps_grav_full_off, Gamma_lnS_eps_grav_full_on, &
./star/private/ctrls_io.f90:532: use_other_build_initial_model, use_other_cgrav, use_other_energy_implicit, use_other_momentum, &
./star/private/ctrls_io.f90:533: use_other_energy, use_other_mesh_functions, use_other_eps_grav, use_other_gradr_factor, &
./star/private/ctrls_io.f90:634: if (.not. (trim(s% energy_eqn_option) == 'dedt' .or. trim(s% energy_eqn_option) == 'eps_grav')) then
./star/private/ctrls_io.f90:637: write(*,*) "Available options are 'dedt' or 'eps_grav'"
./star/private/ctrls_io.f90:1083: s% fix_eps_grav_transition_to_grid = fix_eps_grav_transition_to_grid
./star/private/ctrls_io.f90:1871: s% use_time_centered_eps_grav = use_time_centered_eps_grav
./star/private/ctrls_io.f90:1875: s% use_gravity_rotation_correction = use_gravity_rotation_correction
./star/private/ctrls_io.f90:1876: s% eps_grav_factor = eps_grav_factor
./star/private/ctrls_io.f90:1878: s% include_composition_in_eps_grav = include_composition_in_eps_grav
./star/private/ctrls_io.f90:1881: s% Gamma_lnS_eps_grav_full_off = Gamma_lnS_eps_grav_full_off
./star/private/ctrls_io.f90:1882: s% Gamma_lnS_eps_grav_full_on = Gamma_lnS_eps_grav_full_on
./star/private/ctrls_io.f90:2440: s% use_other_cgrav = use_other_cgrav
./star/private/ctrls_io.f90:2449: s% use_other_eps_grav = use_other_eps_grav
./star/private/ctrls_io.f90:2485: s% surface_accel_div_grav_limit = surface_accel_div_grav_limit
./star/private/ctrls_io.f90:2773: fix_eps_grav_transition_to_grid = s% fix_eps_grav_transition_to_grid
./star/private/ctrls_io.f90:3548: use_time_centered_eps_grav = s% use_time_centered_eps_grav
./star/private/ctrls_io.f90:3552: use_gravity_rotation_correction = s% use_gravity_rotation_correction
./star/private/ctrls_io.f90:3553: eps_grav_factor = s% eps_grav_factor
./star/private/ctrls_io.f90:3555: include_composition_in_eps_grav = s% include_composition_in_eps_grav
./star/private/ctrls_io.f90:3558: Gamma_lnS_eps_grav_full_off = s% Gamma_lnS_eps_grav_full_off
./star/private/ctrls_io.f90:3559: Gamma_lnS_eps_grav_full_on = s% Gamma_lnS_eps_grav_full_on
./star/private/ctrls_io.f90:4117: use_other_cgrav = s% use_other_cgrav
./star/private/ctrls_io.f90:4126: use_other_eps_grav = s% use_other_eps_grav
./star/private/ctrls_io.f90:4162: surface_accel_div_grav_limit = s% surface_accel_div_grav_limit
./star/private/gravity_darkening.f90:2:!!% This module contains code to interpolate gravity darkening coefficients.
./star/private/gravity_darkening.f90:25:module gravity_darkening
./star/private/gravity_darkening.f90:47: public :: gravity_darkening_Teff_coeff, gravity_darkening_L_coeff
./star/private/gravity_darkening.f90:68: coefficient_filename = trim(mesa_data_dir) // '/star_data/gravity_darkening_coefficients.data'
./star/private/gravity_darkening.f90:117: function gravity_darkening_Teff_coeff(omega,inclination) result(Teff_coeff)
./star/private/gravity_darkening.f90:121: end function gravity_darkening_Teff_coeff
./star/private/gravity_darkening.f90:124: function gravity_darkening_L_coeff(omega,inclination) result(L_coeff)
./star/private/gravity_darkening.f90:128: end function gravity_darkening_L_coeff
./star/private/gravity_darkening.f90:131:end module gravity_darkening
./star/private/hydro_vars.f90:36: public :: set_vars_if_needed, set_vars, set_final_vars, update_vars, set_cgrav, &
./star/private/hydro_vars.f90:67: skip_m_grav_and_grav = .false., &
./star/private/hydro_vars.f90:83: s, skip_m_grav_and_grav, &
./star/private/hydro_vars.f90:99: skip_m_grav_and_grav, &
./star/private/hydro_vars.f90:107: skip_other_cgrav, &
./star/private/hydro_vars.f90:121: skip_other_cgrav = .false.
./star/private/hydro_vars.f90:123: skip_m_grav_and_grav = .false.
./star/private/hydro_vars.f90:147: skip_micro_vars, skip_m_grav_and_grav, skip_eos, skip_net, skip_neu, &
./star/private/hydro_vars.f90:148: skip_kap, skip_grads, skip_rotation, skip_brunt, skip_other_cgrav, &
./star/private/hydro_vars.f90:167: s, skip_m_grav_and_grav, &
./star/private/hydro_vars.f90:174: skip_m_grav_and_grav, &
./star/private/hydro_vars.f90:183: logical, parameter :: skip_other_cgrav = .false.
./star/private/hydro_vars.f90:193: skip_m_grav_and_grav, skip_net, skip_neu, skip_kap, &
./star/private/hydro_vars.f90:194: skip_grads, skip_rotation, skip_brunt, skip_other_cgrav, &
./star/private/hydro_vars.f90:208: skip_m_grav_and_grav, skip_net, skip_neu, skip_kap, &
./star/private/hydro_vars.f90:209: skip_grads, skip_rotation, skip_brunt, skip_other_cgrav, &
./star/private/hydro_vars.f90:216: skip_m_grav_and_grav, skip_net, skip_neu, skip_kap, skip_grads, &
./star/private/hydro_vars.f90:217: skip_rotation, skip_brunt, skip_other_cgrav, &
./star/private/hydro_vars.f90:241: skip_micro_vars, skip_m_grav_and_grav, skip_eos, skip_net, skip_neu, &
./star/private/hydro_vars.f90:242: skip_kap, skip_grads, skip_rotation, skip_brunt, skip_other_cgrav, &
./star/private/hydro_vars.f90:459: s, nzlo, nzhi, skip_basic_vars, skip_micro_vars, skip_m_grav_and_grav, &
./star/private/hydro_vars.f90:461: skip_brunt, skip_other_cgrav, &
./star/private/hydro_vars.f90:468: set_m_grav_and_grav, set_scale_height, get_tau, &
./star/private/hydro_vars.f90:478: skip_basic_vars, skip_micro_vars, skip_m_grav_and_grav, &
./star/private/hydro_vars.f90:480: skip_grads, skip_rotation, skip_other_cgrav, &
./star/private/hydro_vars.f90:510: if (.not. skip_other_cgrav) call set_cgrav(s, ierr)
./star/private/hydro_vars.f90:515: if (.not. skip_m_grav_and_grav) then
./star/private/hydro_vars.f90:516: ! don't change m_grav or grav during solver iteratons
./star/private/hydro_vars.f90:517: if (dbg) write(*,*) 'call set_m_grav_and_grav'
./star/private/hydro_vars.f90:518: call set_m_grav_and_grav(s)
./star/private/hydro_vars.f90:738: subroutine set_cgrav(s, ierr)
./star/private/hydro_vars.f90:741: if (s% use_other_cgrav) then
./star/private/hydro_vars.f90:742: call s% other_cgrav(s% id, ierr)
./star/private/hydro_vars.f90:745: write(*,*) 'other_cgrav returned ierr', ierr
./star/private/hydro_vars.f90:749: s% cgrav(1:s% nz) = standard_cgrav
./star/private/hydro_vars.f90:751: end subroutine set_cgrav
./star/private/hydro_vars.f90:835: s, tau_surf, L_surf, R_surf, s% m(1), s% cgrav(1), skip_partials, &
./star/private/hydro_vars.f90:846: write(*,1) 's% cgrav(1)', s% cgrav(1)
./star/private/star_private_def.f90:201: termination_code_str(t_surface_accel_div_grav_limit) = 'surface_accel_div_grav_limit'
./star/private/hydro_riemann.f90:47: ! "Well-balanced schemes for the Euler equations with gravitation",
./star/private/hydro_riemann.f90:93: geometry_source_ad, gravity_source_ad, &
./star/private/hydro_riemann.f90:130: call setup_gravity_source
./star/private/hydro_riemann.f90:134: geometry_source_ad + gravity_source_ad + diffusion_source_ad
./star/private/hydro_riemann.f90:206: subroutine setup_gravity_source
./star/private/hydro_riemann.f90:209: ! left 1/2 of dm gets gravity force at left face
./star/private/hydro_riemann.f90:210: ! right 1/2 of dm gets gravity force at right face.
./star/private/hydro_riemann.f90:211: ! this form is to match the gravity force equilibrium reconstruction.
./star/private/hydro_riemann.f90:227: gravity_source_ad = gsL + gsR ! total gravitational force on cell
./star/private/hydro_riemann.f90:230: end subroutine setup_gravity_source
./star/private/hydro_riemann.f90:286: real(dp) :: cgrav
./star/private/hydro_riemann.f90:287: cgrav = s% cgrav(k)
./star/private/hydro_riemann.f90:288: G = cgrav
./star/private/hydro_riemann.f90:289: if (s% rotation_flag .and. s% use_gravity_rotation_correction) &
./star/private/hydro_riemann.f90:305: gamma1L_ad, gamma1R_ad, csL_ad, csR_ad, G_ad, dPdm_grav_ad, &
./star/private/hydro_riemann.f90:348: ! change PR and PL for gravity
./star/private/hydro_riemann.f90:351: dPdm_grav_ad = -G_ad*s% m_grav(k)/(pow2(r_ad)*A_ad) ! cm^-1 s^-2
./star/private/hydro_riemann.f90:354: PL_ad = PL_ad + delta_m*dPdm_grav_ad
./star/private/hydro_riemann.f90:357: PR_ad = PR_ad + delta_m*dPdm_grav_ad
./star/private/adjust_mass.f90:93: s% total_gravitational_energy_before_adjust_mass, &
./star/private/adjust_mass.f90:1261: w_div_wcrit_roche = w_div_w_roche_jrot(r00,s% m(k),s% j_rot(k),s% cgrav(k), &
./star/private/adjust_mass.f90:1264: w_div_wcrit_roche = w_div_w_roche_omega(r00,s% m(k),s% omega(k),s% cgrav(k), &
./star/private/pulse_gyre.f90:98: call build_atm(s, s%L(1), s%r(1), s%Teff, s%m_grav(1), s%cgrav(1), ierr)
./star/private/pulse_gyre.f90:129: m_outer = s%m_grav(1)
./star/private/pulse_gyre.f90:203: real(dp) :: grav
./star/private/pulse_gyre.f90:229: m = s%m_grav(1) !+ s%atm_structure(atm_delta_m,k)
./star/private/pulse_gyre.f90:239: grav = s%cgrav(1)*m/(r*r)
./star/private/pulse_gyre.f90:240: N2 = grav*grav*(rho/P)*delta*(nabla_ad - nabla)
./star/private/pulse_gyre.f90:295: m = s%m_grav(k)
./star/private/pulse_gyre.f90:304: N2 = eval_face_A_ast(s, k, k_a, k_b)*s%grav(k)/s%r(k)
./star/private/pulse_gyre.f90:372: d2P_dr2_c = -four_thirds*pi*s% cgrav(s% nz)*rho**2
./star/private/turb_support.f90:65: iso, XH1, cgrav, m, gradL_composition_term, mixing_length_alpha, &
./star/private/turb_support.f90:71: XH1, cgrav, m, gradL_composition_term, mixing_length_alpha
./star/private/turb_support.f90:95: iso, XH1, cgrav, m, gradL_composition_term, mixing_length_alpha, &
./star/private/turb_support.f90:121: real(dp) :: cgrav, m, XH1, gradL_old, grada_face_old
./star/private/turb_support.f90:127: cgrav = s% cgrav(k)
./star/private/turb_support.f90:128: m = s% m_grav(k)
./star/private/turb_support.f90:145: iso, XH1, cgrav, m, gradL_composition_term, mixing_length_alpha, &
./star/private/turb_support.f90:151: iso, XH1, cgrav, m, gradL_composition_term, mixing_length_alpha, &
./star/private/turb_support.f90:161: iso, XH1, cgrav, m, gradL_composition_term, mixing_length_alpha, &
./star/private/turb_support.f90:172: XH1, cgrav, m, gradL_composition_term, &
./star/private/turb_support.f90:178: type(auto_diff_real_star_order1) :: Pr, Pg, grav, Lambda, gradL, beta
./star/private/turb_support.f90:196: grav = cgrav*m/pow2(r)
./star/private/turb_support.f90:214: .or. m < 1d-10 .or. r%val < 1d-10 .or. cgrav < 1d-10) return
./star/private/turb_support.f90:255: conv_vel_start, mixing_length_alpha, s% alpha_TDC_DAMP, s%alpha_TDC_DAMPR, s%alpha_TDC_PtdVdt, s%dt, cgrav, m, report, &
./star/private/turb_support.f90:267: chiT, chiRho, Cp, grav, Lambda, rho, P, T, opacity, &
./star/private/turb_support.f90:323: chiT, chiRho, Cp, grav, Lambda, rho, P, T, opacity, &
./star/private/turb_support.f90:348: s% semiconvection_option, cgrav, Cp, gradr, grada, gradL, &
./star/private/pulse_fgong.f90:135: call build_atm(s, s%L(1), s%r(1), s%Teff, s%m_grav(1), s%cgrav(1), ierr)
./star/private/pulse_fgong.f90:167: m_outer = s%m_grav(1)
./star/private/pulse_fgong.f90:185: d2P_dr2_c = -four_thirds*pi*s% cgrav(s% nz)*rho_c**2
./star/private/pulse_fgong.f90:191: global_data(15) = standard_cgrav
./star/private/pulse_fgong.f90:256: real(dp) :: grav
./star/private/pulse_fgong.f90:281: eps_grav => point_data(19,j), &
./star/private/pulse_fgong.f90:297: lnq = log(s%m_grav(1)/m_outer)
./star/private/pulse_fgong.f90:311: grav = s%cgrav(1)*s%m_grav(1)/(r*r)
./star/private/pulse_fgong.f90:312: N2 = grav*grav*(rho/P)*delta*(nabla_ad - s%atm_structure(atm_gradT,k))
./star/private/pulse_fgong.f90:313: A_ast = N2*r/grav
./star/private/pulse_fgong.f90:320: eps_grav = 0d0
./star/private/pulse_fgong.f90:375: eps_grav => point_data(19,j), &
./star/private/pulse_fgong.f90:391: lnq = log(s%m_grav(k)/m_outer)
./star/private/pulse_fgong.f90:402: eps = eval_face(s%dq, s%eps_nuc, k, k_a, k_b) + eval_face(s%dq, s%eps_grav_ad%val, k, k_a, k_b)
./star/private/pulse_fgong.f90:418: eps_grav = eval_face(s%dq, s%eps_grav_ad%val, k, k_a, k_b)
./star/private/pulse_fgong.f90:472: eps_grav => point_data(19,j), &
./star/private/pulse_fgong.f90:495: eps = eval_center(s%rmid, s%eps_nuc, k_a, k_b) + eval_center(s%rmid, s%eps_grav_ad%val, k_a, k_b)
./star/private/pulse_fgong.f90:507: eps_grav = eval_center(s%rmid, s%eps_grav_ad%val, k_a, k_b)
./star/private/struct_burn_mix.f90:339: s% m_grav_start(k) = s% m_grav(k)
./star/private/struct_burn_mix.f90:360: s% total_gravitational_energy_start, &
./star/private/adjust_mesh_split_merge.f90:176: w_div_w_roche_jrot(r00,s% m(k),s% j_rot(k),s% cgrav(k), &
./star/private/adjust_mesh_split_merge.f90:654: s% cgrav(im) = s% cgrav(i0)
./star/private/adjust_mesh_split_merge.f90:734: PE = -s% cgrav(k)*mC*dm/rC
./star/private/adjust_mesh_split_merge.f90:766: if (s% cgrav(k) == 0) then
./star/private/adjust_mesh_split_merge.f90:776: PE = -s% cgrav(k)*mC*dm/rC
./star/private/adjust_mesh_split_merge.f90:1037: s% cgrav(jp) = s% cgrav(j)
./star/private/adjust_mesh_split_merge.f90:1241: s% cgrav(ip) = s% cgrav(i)
./star/private/adjust_mesh_split_merge.f90:1432: total_PE = total_PE - s% cgrav(k)*mC*s% dm(k)/rC
./star/private/star_history_def.f90:80: integer, parameter :: h_star_gravitational_mass = h_kh_div_mdot_timescales + 1
./star/private/star_history_def.f90:81: integer, parameter :: h_star_mass_grav_div_mass = h_star_gravitational_mass + 1
./star/private/star_history_def.f90:83: integer, parameter :: h_conv_mx1_top = h_star_mass_grav_div_mass + 1
./star/private/star_history_def.f90:200: integer, parameter :: h_total_gravitational_energy_after_adjust_mass = &
./star/private/star_history_def.f90:203: h_total_gravitational_energy_after_adjust_mass + 1
./star/private/star_history_def.f90:212: integer, parameter :: h_total_gravitational_energy = h_total_internal_energy + 1
./star/private/star_history_def.f90:213: integer, parameter :: h_total_turbulent_energy = h_total_gravitational_energy + 1
./star/private/star_history_def.f90:220: integer, parameter :: h_log_total_gravitational_energy = h_log_total_internal_energy + 1
./star/private/star_history_def.f90:221: integer, parameter :: h_log_total_turbulent_energy = h_log_total_gravitational_energy + 1
./star/private/star_history_def.f90:232: integer, parameter :: h_total_eps_grav = h_virial_thm_rel_err + 1
./star/private/star_history_def.f90:233: integer, parameter :: h_work_outward_at_surface = h_total_eps_grav + 1
./star/private/star_history_def.f90:300: integer, parameter :: h_gravity = h_effective_T + 1
./star/private/star_history_def.f90:301: integer, parameter :: h_log_g = h_gravity + 1
./star/private/star_history_def.f90:329: integer, parameter :: h_center_eps_grav = h_center_non_nuc_neu + 1
./star/private/star_history_def.f90:330: integer, parameter :: h_center_degeneracy = h_center_eps_grav + 1
./star/private/star_history_def.f90:346: integer, parameter :: h_log_abs_Lgrav = h_dt_div_min_tau_conv + 1
./star/private/star_history_def.f90:347: integer, parameter :: h_eps_grav_integral = h_log_abs_Lgrav + 1
./star/private/star_history_def.f90:348: integer, parameter :: h_log_extra_L = h_eps_grav_integral + 1
./star/private/star_history_def.f90:433: integer, parameter :: h_cz_grav = h_cz_scale_height + 1
./star/private/star_history_def.f90:434: integer, parameter :: h_cz_log_eps_nuc = h_cz_grav + 1
./star/private/star_history_def.f90:457: integer, parameter :: h_cz_top_grav = h_cz_top_scale_height + 1
./star/private/star_history_def.f90:458: integer, parameter :: h_cz_top_log_eps_nuc = h_cz_top_grav + 1
./star/private/star_history_def.f90:636: integer, parameter :: h_grav_dark_L_polar = h_diffusion_solver_iters + 1
./star/private/star_history_def.f90:637: integer, parameter :: h_grav_dark_Teff_polar = h_grav_dark_L_polar + 1
./star/private/star_history_def.f90:638: integer, parameter :: h_grav_dark_L_equatorial = h_grav_dark_Teff_polar + 1
./star/private/star_history_def.f90:639: integer, parameter :: h_grav_dark_Teff_equatorial = h_grav_dark_L_equatorial + 1
./star/private/star_history_def.f90:641: integer, parameter :: h_apsidal_constant_k2 = h_grav_dark_Teff_equatorial + 1
./star/private/star_history_def.f90:731: history_column_name(h_star_gravitational_mass) = 'star_gravitational_mass'
./star/private/star_history_def.f90:732: history_column_name(h_star_mass_grav_div_mass) = 'star_mass_grav_div_mass'
./star/private/star_history_def.f90:840: history_column_name(h_total_gravitational_energy_after_adjust_mass) = &
./star/private/star_history_def.f90:841: 'total_gravitational_energy_after_adjust_mass'
./star/private/star_history_def.f90:850: history_column_name(h_total_gravitational_energy) = 'total_gravitational_energy'
./star/private/star_history_def.f90:858: history_column_name(h_log_total_gravitational_energy) = 'log_total_gravitational_energy'
./star/private/star_history_def.f90:868: history_column_name(h_total_eps_grav) = 'total_eps_grav'
./star/private/star_history_def.f90:912: history_column_name(h_gravity) = 'gravity'
./star/private/star_history_def.f90:983: history_column_name(h_center_eps_grav) = 'center_eps_grav'
./star/private/star_history_def.f90:1059: history_column_name(h_log_abs_Lgrav) = 'log_abs_Lgrav'
./star/private/star_history_def.f90:1060: history_column_name(h_eps_grav_integral) = 'eps_grav_integral'
./star/private/star_history_def.f90:1138: history_column_name(h_cz_grav) = 'cz_grav'
./star/private/star_history_def.f90:1162: history_column_name(h_cz_top_grav) = 'cz_top_grav'
./star/private/star_history_def.f90:1305: history_column_name(h_grav_dark_L_polar) = 'grav_dark_L_polar'
./star/private/star_history_def.f90:1306: history_column_name(h_grav_dark_Teff_polar) = 'grav_dark_Teff_polar'
./star/private/star_history_def.f90:1308: history_column_name(h_grav_dark_L_equatorial) = 'grav_dark_L_equatorial'
./star/private/star_history_def.f90:1309: history_column_name(h_grav_dark_Teff_equatorial) = 'grav_dark_Teff_equatorial'
./star/private/hydro_eqns.f90:64: use eps_grav, only: zero_eps_grav_and_partials
./star/private/hydro_eqns.f90:126: if (trim(s% energy_eqn_option) == 'eps_grav') then
./star/private/hydro_eqns.f90:127: s% eps_grav_form_for_energy_eqn = .true.
./star/private/hydro_eqns.f90:129: s% eps_grav_form_for_energy_eqn = .false.
./star/private/hydro_eqns.f90:237: call zero_eps_grav_and_partials(s, k)
./star/private/hydro_eqns.f90:638: jrot_ratio = jrot00/sqrt(s% cgrav(k)*s% m(k)*r00)
./star/private/hydro_eqns.f90:640: jrot_ratio = s% j_rot(k)/sqrt(s% cgrav(k)*s% m(k)*r00)
./star/private/hydro_eqns.f90:684: scale = 1d6*max(1d2*sqrt(s% cgrav(k)*s% m(k)*s%r_start(k))/s%dt,&
./star/private/hydro_eqns.f90:853: dP0 = s% cgrav(1)*s% m_grav(1)*s% dm(1)/(8*pi*pow4(r))
./star/private/diffusion_support.f90:141: s% grav(k), dlnne_dr_face, &
./star/private/diffusion_support.f90:158: amu*(s% grav(k))*e_ap(k) ! + e_ar(k) should be the radiation term, skipping for now
./star/private/diffusion_support.f90:166: g_field_face(k) = s% grav(k)
./star/private/diffusion_support.f90:168: ! Electric and gravitational field from Thoul solve.
./star/private/diffusion_support.f90:255: dlnP_dr_face, dlnT_dr_face, dlnRho_dr_face, grav, dlnne_dr_face, &
./star/private/diffusion_support.f90:273: dlnRho_dr_face, grav, dlnne_dr_face
./star/private/diffusion_support.f90:312: grav, dlnne_dr_face, X_face, &
./star/private/diffusion_support.f90:630: dlnP_dr, dlnT_dr, dlnRho_dr, grav, dlnne_dr, X_face, &
./star/private/diffusion_support.f90:639: real(dp), intent(in) :: dlnP_dr, dlnT_dr, dlnRho_dr, grav, dlnne_dr
./star/private/diffusion_support.f90:684: vlnP(i) = AP(i)*amu*grav*diffusion_factor(i)*limit_coeff
./star/private/diffusion_support.f90:806:! the gravitational force g.
./star/private/diffusion_support.f90:823: ! For I=N, we get the gravitational force g
./star/private/diffusion_support.f90:1110: ! Right now this is for thermal diffusion off, assuming gravity
./star/private/diffusion_support.f90:1156: ! gravity, k_b*T, and number density gradients to provide the
./star/private/diffusion_support.f90:1238: ! and making up for losing that equation by treating grav as known.
./star/private/diffusion_support.f90:1283: ! Assuming gravity is a known, so there are m uknown
./star/private/diffusion_support.f90:1395: ! gravity, k_b*T, and number density gradients to provide the
./star/private/hydro_rotation.f90:45: real(dp) function w_div_w_roche_omega(rphi,Mphi,omega,cgrav, max_w, max_w2, w_div_wc_flag) result(w_roche)
./star/private/hydro_rotation.f90:46: real(dp), intent(in) :: rphi,Mphi,omega,cgrav, max_w, max_w2
./star/private/hydro_rotation.f90:55: dimless_rphi = rphi*pow(abs(omega), two_thirds)/pow(cgrav*Mphi,one_third)
./star/private/hydro_rotation.f90:101: real(dp) function w_div_w_roche_jrot(rphi,Mphi,jrot,cgrav, max_w, max_w2, w_div_wc_flag) result(w_roche)
./star/private/hydro_rotation.f90:102: real(dp), intent(in) :: rphi,Mphi,jrot,cgrav, max_w, max_w2
./star/private/hydro_rotation.f90:112: dimless_factor = abs(jrot)/sqrt(cgrav*Mphi*rphi)
./star/private/hydro_rotation.f90:226: w_div_w_roche_jrot(s% r(k),s% m(k),s% j_rot(k),s% cgrav(k), &
./star/private/hydro_rotation.f90:309: w_div_w_roche_jrot(get_r_from_xh(s,k),s% m(k),s% j_rot(k),s% cgrav(k), &
./star/private/hydro_rotation.f90:325: w_div_w_roche_omega(get_r_from_xh(s,k),s% m(k),s% omega(k),s% cgrav(k), &
./star/private/hydro_rotation.f90:512: cgrav, kap, mmid, Lmid, rmid, logT_sum, logRho_sum
./star/private/hydro_rotation.f90:547: mmid = 0.5d0*(s% m_grav(k) + s% m_grav(k+1))
./star/private/hydro_rotation.f90:549: cgrav = 0.5d0*(s% cgrav(k) + s% cgrav(k+1))
./star/private/hydro_rotation.f90:618: cgrav = 0.5d0*(s% cgrav(k) + s% cgrav(k+1))
./star/private/hydro_rotation.f90:619: mmid = 0.5d0*(s% m_grav(k) + s% m_grav(k+1))
./star/private/hydro_rotation.f90:626: omega_crit = sqrt(gamma_factor*cgrav*mmid/pow3(rmid))
./star/private/rsp_def.f90:150: use const_def, only: standard_cgrav, boltz_sigma, &
./star/private/rsp_def.f90:158: G=standard_cgrav
./star/private/pulse_saio.f90:93: call build_atm(s, s%L(1), s%r(1), s%Teff, s%m_grav(1), s%cgrav(1), ierr)
./star/private/pulse_saio.f90:120: m_outer = s%m_grav(1)
./star/private/pulse_saio.f90:212: m = s% m_grav(1)
./star/private/pulse_saio.f90:232: Lrad = (16*pi*clight*crad*s%cgrav(k)*m*T*T*T*T*nabla)/(3*P*kap)
./star/private/pulse_saio.f90:280: m = s%m_grav(k)
./star/private/pulse_saio.f90:288: eps = eval_face(s%dq, s%eps_nuc, k, k_a, k_b) + eval_face(s%dq, s%eps_grav_ad%val, k, k_a, k_b)
./star/private/pulse_saio.f90:303: Lrad = (16*pi*clight*crad*s%cgrav(k)*m*T*T*T*T*nabla)/(3*P*kap)
./star/private/profile_getval.f90:523: val = safe_log10(s% grav(k))
./star/private/profile_getval.f90:524: case (p_grav)
./star/private/profile_getval.f90:525: val = s% grav(k)
./star/private/profile_getval.f90:527: val = s% grav(k)/s% r(k)
./star/private/profile_getval.f90:529: val = s% r(k)/s% grav(k)
./star/private/profile_getval.f90:530: case (p_signed_log_eps_grav)
./star/private/profile_getval.f90:531: val = s% eps_grav_ad(k)% val
./star/private/profile_getval.f90:541: val = s% eps_nuc(k) - s% non_nuc_neu(k) + s% eps_grav_ad(k)% val
./star/private/profile_getval.f90:648: case (p_m_grav)
./star/private/profile_getval.f90:649: val = s% m_grav(k)/Msun
./star/private/profile_getval.f90:650: case (p_m_grav_div_m_baryonic)
./star/private/profile_getval.f90:651: val = s% m_grav(k)/s% m(k)
./star/private/profile_getval.f90:1013: val = safe_log10(sqrt(2*s% cgrav(k)*s% m(k)/(s% r(k))))
./star/private/profile_getval.f90:1020: val = val/sqrt(2*s% cgrav(k)*s% m(k)/(s% r(k)))
./star/private/profile_getval.f90:1027: val = val/sqrt(2d0*s% cgrav(k)*s% m(k)/(s% r(k)))
./star/private/profile_getval.f90:1128: case (p_cgrav_factor)
./star/private/profile_getval.f90:1129: val = s% cgrav(k)/standard_cgrav
./star/private/profile_getval.f90:1139: case (p_extra_grav)
./star/private/profile_getval.f90:1140: val = s% extra_grav(k)%val
./star/private/profile_getval.f90:1147: case (p_log_abs_eps_grav_dm_div_L)
./star/private/profile_getval.f90:1149: abs(s% eps_grav_ad(k)% val)*s% dm(k)/max(1d0,abs(s% L(k))))
./star/private/profile_getval.f90:1151: case (p_eps_grav_composition_term)
./star/private/profile_getval.f90:1152: if (s% include_composition_in_eps_grav) &
./star/private/profile_getval.f90:1153: val = s% eps_grav_composition_term(k)
./star/private/profile_getval.f90:1155: case (p_eps_grav_plus_eps_mdot)
./star/private/profile_getval.f90:1156: val = s% eps_grav_ad(k)% val + s% eps_mdot(k)
./star/private/profile_getval.f90:1157: case (p_ergs_eps_grav_plus_eps_mdot)
./star/private/profile_getval.f90:1158: val = (s% eps_grav_ad(k)% val + s% eps_mdot(k))*s% dm(k)*s% dt
./star/private/profile_getval.f90:1229: case (p_dPdr_div_grav)
./star/private/profile_getval.f90:1230: if (k > 1 .and. k < nz .and. s% cgrav(k) > 0d0 .and. s% RTI_flag) then
./star/private/profile_getval.f90:1241: case (p_dvdt_grav)
./star/private/profile_getval.f90:1242: val = -s% cgrav(k)*s% m(k)/(s% r(k)*s% r(k))
./star/private/profile_getval.f90:1246: case (p_dm_eps_grav)
./star/private/profile_getval.f90:1247: val = s% eps_grav_ad(k)% val*s% dm(k)
./star/private/profile_getval.f90:1248: case (p_eps_grav)
./star/private/profile_getval.f90:1249: val = s% eps_grav_ad(k)% val
./star/private/profile_getval.f90:1256: case (p_env_eps_grav)
./star/private/profile_getval.f90:1257: val = -s% gradT_sub_grada(k)*s% grav(k)*s% mstar_dot*s% Cp(k)*s% T(k) / &
./star/private/profile_getval.f90:1458: s% Cp(k)*s% T(k)/(s% Peos(k)/(s% rho(k)*s% grav(k))/s% csound(k)))
./star/private/profile_getval.f90:1460: val = safe_log10(s% Peos(k)/(s% rho(k)*s% grav(k))/s% csound(k))
./star/private/profile_getval.f90:1462: val = s% Peos(k)/(s% rho(k)*s% grav(k))/Rsun
./star/private/profile_getval.f90:1464: val = s% Peos(k)/(s% rho(k)*s% grav(k))
./star/private/profile_getval.f90:2138: if (s% calculate_Brunt_N2) val = s% brunt_N2(k)*s% r(k)/s% grav(k)
./star/private/profile_getval.f90:2141: if (s% calculate_Brunt_N2) val = s% brunt_N2(k)*s% r(k)/s% grav(k)/x/x
./star/private/profile_getval.f90:2144: safe_log10(s% brunt_N2(k)/(3*s% cgrav(1)*s% m_grav(1)/pow3(s% r(1))))
./star/private/profile_getval.f90:2147: s% brunt_N2(k)/(3*s% cgrav(1)*s% m_grav(1)/pow3(s% r(1)))
./star/private/profile_getval.f90:2150: sqrt(max(0d0,s% brunt_N2(k))/(3*s% cgrav(1)*s% m_grav(1)/pow3(s% r(1))))
./star/private/create_initial_model.f90:160: G = standard_cgrav
./star/private/create_initial_model.f90:270: real(dp) :: kap,tau,grav
./star/private/create_initial_model.f90:350: grav = G*y(2)/(y(1)*y(1))
./star/private/create_initial_model.f90:351: tau = kap * P / grav
./star/private/hydro_temperature.f90:303: use hydro_momentum, only: expected_HSE_grav_term
./star/private/hydro_temperature.f90:310: type(auto_diff_real_star_order1) :: grav, area, P00, Pm1
./star/private/hydro_temperature.f90:318: ! for rotation, multiply gravity by factor fp. MESA 2, eqn 22.
./star/private/hydro_temperature.f90:320: call expected_HSE_grav_term(s, k, grav, area, ierr)
./star/private/hydro_temperature.f90:335: dlnPdm_qhse = grav/(area*Ppoint) ! note that expected_HSE_grav_term is negative
./star/private/hydro_rsp2_support.f90:110: skip_m_grav_and_grav = .false., &
./star/private/hydro_rsp2_support.f90:117: skip_other_cgrav = .true., &
./star/private/hydro_rsp2_support.f90:127: skip_micro_vars, skip_m_grav_and_grav, skip_eos, skip_net, skip_neu, &
./star/private/hydro_rsp2_support.f90:128: skip_kap, skip_grads, skip_rotation, skip_brunt, skip_other_cgrav, &
./star/private/hydro_rsp2_support.f90:146: s, s% tau_factor*s% tau_base, s% L(1), s% r(1), s% m(1), s% cgrav(1), skip_partials, &
./star/private/hydro_rsp2_support.f90:368: P_00 = P_m1 + s% cgrav(k)*s% m(k)*dm_face/(4d0*pi*pow4(s% r(k)))
./star/private/hydro_rsp2_support.f90:440: Hp_face = pow2(r_00)*Peos_div_rho/(s% cgrav(k)*s% m(k))
./star/private/hydro_rsp2_support.f90:446: Hp_face = pow2(r_00)*Peos_div_rho/(s% cgrav(k)*s% m(k))
./star/private/remove_shells.f90:288: pe = -s% cgrav(k)*m_cntr*s% dm(k)/rC
./star/private/remove_shells.f90:318: pe = -s% cgrav(k)*m_cntr*s% dm(k)/rC
./star/private/remove_shells.f90:855: vesc = sqrt(2*s% cgrav(k)*s% m(k)/(s% r(k)))
./star/private/overshoot_utils.f90:188: Hp = P/(rho*s%cgrav(k_m)* &
./star/private/eps_mdot.f90:491: s% total_gravitational_energy_after_adjust_mass = 0d0
./star/private/eps_mdot.f90:623: ! In the absence of composition changes this is just dm_new * (change in gravitational potential),
./star/private/eps_mdot.f90:665: s% total_gravitational_energy_after_adjust_mass, &
./star/private/pulse_utils.f90:290: A_ast_face = s%brunt_N2(k)*s%r(k)/s%grav(k)
./star/private/pulse_utils.f90:296: A_ast_1 = s%brunt_N2(k_a+1)*s%r(k_a+1)/s%grav(k_a+1)
./star/private/pulse_utils.f90:297: A_ast_2 = s%brunt_N2(k_a+2)*s%r(k_a+2)/s%grav(k_a+2)
./star/private/pulse_utils.f90:303: A_ast_1 = s%brunt_N2(k_b-1)*s%r(k_b-1)/s%grav(k_b)
./star/private/pulse_utils.f90:304: A_ast_2 = s%brunt_N2(k_b)*s%r(k_b)/s%grav(k_b)
./star/private/pulse_utils.f90:310: A_ast_face = s%brunt_N2(k)*s%r(k)/s%grav(k)
./star/private/hydro_momentum.f90:39: expected_HSE_grav_term
./star/private/hydro_momentum.f90:108: other_ad, dm_div_A_ad, grav_ad, area_ad, dPtot_ad, d_mlt_Pturb_ad, &
./star/private/hydro_momentum.f90:109: iPtotavg_ad, other_dm_div_A_ad, grav_dm_div_A_ad, &
./star/private/hydro_momentum.f90:121:! dv/dt = - G*m/r^2 - (dPtot_ad + d_mlt_Pturb_ad)*area/dm + extra_grav + Uq + RTI_diffusion + RTI_kick
./star/private/hydro_momentum.f90:123:! grav_ad = expected_HSE_grav_term = -G*m/r^2 with possible modifications for rotation
./star/private/hydro_momentum.f90:124:! other_ad = expected_non_HSE_term = extra_grav - dv/dt + Uq
./star/private/hydro_momentum.f90:125:! extra_grav is from the other_momentum hook
./star/private/hydro_momentum.f90:133:! 0 = extra_grav - dv/dt + Uq - G*m/r^2 - RTI_diffusion - RTI_kick - (dPtot_ad + d_mlt_Pturb_ad)*area/dm
./star/private/hydro_momentum.f90:134:! 0 = other + grav - RTI_terms - (dPtot_ad + d_mlt_Pturb_ad)*area/dm
./star/private/hydro_momentum.f90:135:! 0 = (other + grav - RTI_terms)*dm/area - dPtot_ad - d_mlt_Pturb_ad
./star/private/hydro_momentum.f90:136:! 0 = other_dm_div_A_ad + grav_dm_div_A_ad - dPtot_ad - d_mlt_Pturb_ad + RTI_terms_dm_div_A_ad
./star/private/hydro_momentum.f90:138: call setup_HSE(dm_div_A, ierr); if (ierr /= 0) return ! grav_ad and dm_div_A_ad
./star/private/hydro_momentum.f90:139: call setup_non_HSE(ierr); if (ierr /= 0) return ! other = s% extra_grav(k) - dv_dt
./star/private/hydro_momentum.f90:145: grav_dm_div_A_ad = grav_ad*dm_div_A_ad
./star/private/hydro_momentum.f90:150: other_dm_div_A_ad + grav_dm_div_A_ad - dPtot_ad - d_mlt_Pturb_ad + RTI_terms_dm_div_A_ad
./star/private/hydro_momentum.f90:189: ! rho represents the inertial (gravitational) mass density.
./star/private/hydro_momentum.f90:213: call expected_HSE_grav_term(s, k, grav_ad, area_ad, ierr)
./star/private/hydro_momentum.f90:224: ! other = extra_grav - dv/dt
./star/private/hydro_momentum.f90:343: subroutine expected_HSE_grav_term(s, k, grav, area, ierr)
./star/private/hydro_momentum.f90:347: type(auto_diff_real_star_order1), intent(out) :: area, grav
./star/private/hydro_momentum.f90:359: if (s% rotation_flag .and. s% use_gravity_rotation_correction) then
./star/private/hydro_momentum.f90:360: grav = -s% cgrav(k)*s% m_grav(k)*inv_R2*s% fp_rot(k)
./star/private/hydro_momentum.f90:362: grav = -s% cgrav(k)*s% m_grav(k)*inv_R2
./star/private/hydro_momentum.f90:372: write(*,*) 'expected_HSE_grav_term', s% solver_test_partials_var
./star/private/hydro_momentum.f90:375: end subroutine expected_HSE_grav_term
./star/private/hydro_momentum.f90:378: ! other = s% extra_grav(k) - s% dv_dt(k)
./star/private/hydro_momentum.f90:400: extra_ad = s% extra_grav(k)
./star/private/do_one_utils.f90:94: write(*,'(a)') " 'lg_gsurf' is log10 surface gravity,"
./star/private/do_one_utils.f90:579: real(dp) :: log_surface_gravity, log_surface_temperature, log_surface_density, &
./star/private/do_one_utils.f90:641: if (s% cgrav(1) <= 0d0) then
./star/private/do_one_utils.f90:644: v_surf_div_v_esc = v_surf/sqrt(2*s% cgrav(1)*s% m(1)/(s% photosphere_r*Rsun))
./star/private/do_one_utils.f90:651: log_surface_gravity = safe_log10(s%grav(1))
./star/private/do_one_utils.f90:772: sqrt(s% cgrav(k_omega)*s% m(k_omega)/ pow3(s% r_equatorial(k_omega)))
./star/private/do_one_utils.f90:1058: else if (log_surface_gravity <= s% log_g_lower_limit) then
./star/private/do_one_utils.f90:1059: call compare_to_target('log_surface_gravity <= log_g_lower_limit', &
./star/private/do_one_utils.f90:1060: log_surface_gravity, s% log_g_lower_limit, t_log_g_lower_limit)
./star/private/do_one_utils.f90:1062: else if (log_surface_gravity >= s% log_g_upper_limit) then
./star/private/do_one_utils.f90:1063: call compare_to_target('log_surface_gravity >= log_g_upper_limit', &
./star/private/do_one_utils.f90:1064: log_surface_gravity, s% log_g_upper_limit, t_log_g_upper_limit)
./star/private/do_one_utils.f90:1483: real(dp) :: v, surf_grav, power_he_burn, power_z_burn, &
./star/private/evolve.f90:232: call set_to_NaN(s% total_eps_grav)
./star/private/evolve.f90:235: call set_to_NaN(s% total_gravitational_energy)
./star/private/evolve.f90:236: call set_to_NaN(s% total_gravitational_energy_after_adjust_mass)
./star/private/evolve.f90:237: call set_to_NaN(s% total_gravitational_energy_before_adjust_mass)
./star/private/evolve.f90:238: call set_to_NaN(s% total_gravitational_energy_end)
./star/private/evolve.f90:239: call set_to_NaN(s% total_gravitational_energy_initial)
./star/private/evolve.f90:240: call set_to_NaN(s% total_gravitational_energy_old)
./star/private/evolve.f90:241: call set_to_NaN(s% total_gravitational_energy_start)
./star/private/evolve.f90:290: use hydro_vars, only: set_vars_if_needed, set_vars, set_cgrav
./star/private/evolve.f90:344: call set_cgrav(s, ierr)
./star/private/evolve.f90:345: if (failed('set_cgrav')) return
./star/private/evolve.f90:424: s% total_gravitational_energy_old, &
./star/private/evolve.f90:435: s% total_gravitational_energy_old, &
./star/private/evolve.f90:447: s% total_gravitational_energy_initial = &
./star/private/evolve.f90:448: s% total_gravitational_energy_old
./star/private/evolve.f90:513: use hydro_vars, only: set_vars_if_needed, set_vars, set_final_vars, set_cgrav
./star/private/evolve.f90:590: call set_cgrav(s, ierr)
./star/private/evolve.f90:591: if (failed('set_cgrav')) return
./star/private/evolve.f90:848: kh_timescale = eval_kh_timescale(s% cgrav(1), s% mstar, r_phot, s% L(1))
./star/private/evolve.f90:1117: diff_total_gravitational_energy, diff_total_internal_energy, diff_total_kinetic_energy, &
./star/private/evolve.f90:1122: total_energy_from_fixed_m_grav
./star/private/evolve.f90:1153: s% total_gravitational_energy_end, &
./star/private/evolve.f90:1160: s% total_gravitational_energy_old = s% total_gravitational_energy_end
./star/private/evolve.f90:1170: s% total_gravitational_energy_end, &
./star/private/evolve.f90:1189: s% total_eps_grav = dt*dot_product(s% dm(1:nz), s% eps_grav_ad(1:nz)% val)
./star/private/evolve.f90:1191: if (s% u_flag .and. s% total_eps_grav /= 0d0) then
./star/private/evolve.f90:1192: write(*,2) 'u_flag energy accounting ignores total_eps_grav', s% model_number, s% total_eps_grav
./star/private/evolve.f90:1193: s% total_eps_grav = 0
./star/private/evolve.f90:1310: ! during the newton iterations, m_grav remains fixed,
./star/private/evolve.f90:1313: ! -sum dm G Delta(m_grav) / r
./star/private/evolve.f90:1315: ! after the iterations, we update m_grav before doing the energy accounting.
./star/private/evolve.f90:1318: total_energy_from_fixed_m_grav = &
./star/private/evolve.f90:1319: -dot_product(s% cgrav(1:nz)/s%r(1:nz)*(s%m_grav(1:nz)-s%m_grav_start(1:nz)), s% dm(1:nz))
./star/private/evolve.f90:1321: total_energy_from_fixed_m_grav = 0
./star/private/evolve.f90:1377: s% total_gravitational_energy = s% total_gravitational_energy_end
./star/private/evolve.f90:1385: write(*,2) 'INFO: use_mass_corrections incurred rel_E_err', s% model_number, -total_energy_from_fixed_m_grav/s% total_energy
./star/private/evolve.f90:1511: diff_total_gravitational_energy = &
./star/private/evolve.f90:1512: s% total_gravitational_energy_end - s% total_gravitational_energy_start
./star/private/evolve.f90:1524: (sum_cell_dpe - diff_total_gravitational_energy)/s% total_energy, &
./star/private/evolve.f90:1525: sum_cell_dpe, diff_total_gravitational_energy
./star/private/pulse_osc.f90:138: call build_atm(s, s%L(1), s%r(1), s%Teff, s%m_grav(1), s%cgrav(1), ierr)
./star/private/pulse_osc.f90:170: m_outer = s%m_grav(1)
./star/private/pulse_osc.f90:191: d2P_dr2_c = -four_thirds*pi*s% cgrav(s% nz)*rho_c**2
./star/private/pulse_osc.f90:269: real(dp) :: grav
./star/private/pulse_osc.f90:314: lnq = log(s%m_grav(1)/m_outer)
./star/private/pulse_osc.f90:328: grav = s%cgrav(1)*s%m_grav(1)/(r*r)
./star/private/pulse_osc.f90:329: N2 = grav*grav*(rho/P)*delta*(nabla_ad - nabla)
./star/private/pulse_osc.f90:330: A_ast = N2*r/grav
./star/private/pulse_osc.f90:417: lnq = log(s%m_grav(k)/m_outer)
./star/private/pulse_osc.f90:428: eps = eval_face(s%dq, s%eps_nuc, k, k_a, k_b) + eval_face(s%dq, s%eps_grav_ad%val, k, k_a, k_b)
./star/private/pulse_osc.f90:525: eps = eval_center(s%rmid, s%eps_nuc, k_a, k_b) + eval_center(s%rmid, s%eps_grav_ad%val, k_a, k_b)
./star/private/rotation_mix_info.f90:64: r, m, L, j_rot, gradT, grada, grav, visc, Ri
./star/private/rotation_mix_info.f90:191: rho, T, r, L, omega, Cp, abar, zbar, delta, grav, &
./star/private/rotation_mix_info.f90:466: grav => s% grav
./star/private/rotation_mix_info.f90:533: scale_height(i) = P(i)*r(i)*r(i)/(s% cgrav(i)*m(i)*rho(i))
./star/private/rotation_mix_info.f90:534: scale_height2 = sqrt(P(i)/s% cgrav(i))/rho(i)
./star/private/rotation_mix_info.f90:595: N2(i) = -grav(i)*(1/gamma1(i) - dlnRho_dlnP)/scale_height(i)
./star/private/rotation_mix_info.f90:596: N2_mu(i) = -grav(i)/scale_height(i)*(1/chiRho(i) - delta(i)*dlnT_dlnP - dlnRho_dlnP)
./star/private/rotation_mix_info.f90:619: t_dyn(i) = sqrt(r(i)*r(i)*r(i)/(s% cgrav(i)*m(i)))
./star/private/rotation_mix_info.f90:620: t_kh(i) = s% cgrav(i)*m(i)*m(i)/(r(i)*max(1d0,L(i)+L_neu(i)))
./star/private/rotation_mix_info.f90:625: ri0 = (rho(i)*delta(i)/P(i))*pow2(dlnR_domega(i)*grav(i))
./star/private/rotation_mix_info.f90:677: denom = (-gradT_sub_grada(i)*delta(i)*pow2(s% cgrav(i)*m(i)))
./star/private/rotation_mix_info.f90:800: D = grav(i)/rho(i)*(dRho_dr_ad(i)-dRho_dr(i))+dr2omega(i)/(r(i)*r(i)*r(i))
./star/private/rotation_mix_info.f90:826: D = pow2(height*s% D_SH(k)*r(k)/grav(k))/t_dyn(k)
./star/private/rotation_mix_info.f90:1147: rho, T, r, L, omega, Cp, abar, zbar, delta, grav, &
./star/private/rotation_mix_info.f90:1153: rho, T, r, L, omega, Cp, abar, zbar, delta, grav, &
./star/private/rotation_mix_info.f90:1322: pow(grav(k)*delta(k)*scale_height(k)*MAX(0.D0,L(k))/ &
./star/private/relax.f90:311: logical :: do_element_diffusion, include_composition_in_eps_grav
./star/private/relax.f90:354: include_composition_in_eps_grav = s% include_composition_in_eps_grav
./star/private/relax.f90:355: s% include_composition_in_eps_grav = .false. ! don't need energetic effects of artificial changes
./star/private/relax.f90:367: s% include_composition_in_eps_grav = include_composition_in_eps_grav
./star/private/init.f90:434: use other_cgrav, only: default_other_cgrav
./star/private/init.f90:441: use other_eps_grav, only: null_other_eps_grav
./star/private/init.f90:591: s% other_cgrav => default_other_cgrav
./star/private/init.f90:605: s% other_eps_grav => null_other_eps_grav
./star/private/init.f90:689: s% total_gravitational_energy = 0d0
./star/private/init.f90:1032: s% total_gravitational_energy, &
./star/private/init.f90:1043: s% total_gravitational_energy, &
./star/private/pre_ms_model.f90:55: lgL, eps_grav, lnd, dlnd, lnd1, lnd3, y1, y3, epsx, epsy, &
./star/private/pre_ms_model.f90:180: ! use uniform eps_grav to give that luminosity
./star/private/pre_ms_model.f90:181: eps_grav = exp10(lgL)*Lsun/mstar
./star/private/pre_ms_model.f90:188: write(*,1) 'eps_grav', eps_grav
./star/private/pre_ms_model.f90:200: rpar(i+1) = eps_grav; i = i+1
./star/private/pre_ms_model.f90:306: real(dp) :: rho_c, T_c, eps_grav, x, z, abar, zbar, d_log10_P
./star/private/pre_ms_model.f90:337: eps_grav = rpar(i+1); i = i+1
./star/private/pre_ms_model.f90:355: s, T_c, rho_c, d_log10_P, eps_grav, &
./star/private/pre_ms_model.f90:380: s, T_c, rho_c, d_log10_P_in, eps_grav_in, &
./star/private/pre_ms_model.f90:392: T_c, rho_c, d_log10_P_in, eps_grav_in, &
./star/private/pre_ms_model.f90:406: res(num_eos_basic_results), eps_grav, P_c, logP, m, &
./star/private/pre_ms_model.f90:411: cgrav, r, rmid, rho, logRho, T, lnT, L, P, P0, dm, m0, L0, r0, lnT0, T0, &
./star/private/pre_ms_model.f90:432: cgrav = standard_cgrav
./star/private/pre_ms_model.f90:434: eps_grav = eps_grav_in
./star/private/pre_ms_model.f90:435: if (dbg) write(*,1) 'eps_grav', eps_grav
./star/private/pre_ms_model.f90:459: P = P_c - 3*cgrav/(8*pi)*pow(pi4*rho_c/3,4d0/3d0)*pow(m,two_thirds)
./star/private/pre_ms_model.f90:467: pow(pi/6,one_third)*cgrav*grada*pow(rho_c*rho_c*m,two_thirds)/P_c
./star/private/pre_ms_model.f90:486: L = eps_grav*m ! L at nz
./star/private/pre_ms_model.f90:499: eps_grav = 1.1d0*eps_grav
./star/private/pre_ms_model.f90:558: dm = -pi4*pow4(rmid)*(P-P0)/(cgrav*mmid)
./star/private/pre_ms_model.f90:564: L = L0 + dm*eps_grav ! luminosity at point k
./star/private/pre_ms_model.f90:701: opacity, grav, scale_height, scale_height2, gradr, cgrav
./star/private/pre_ms_model.f90:726: cgrav = standard_cgrav
./star/private/pre_ms_model.f90:727: grav = cgrav*m/pow2(r)
./star/private/pre_ms_model.f90:728: scale_height = P/(grav*rho) ! this assumes HSE
./star/private/pre_ms_model.f90:730: scale_height2 = sqrt(P/cgrav)/rho
./star/private/pre_ms_model.f90:735: gradr = P*opacity*L/(16d0*pi*clight*m*cgrav*crad*pow4(T)/3d0)
./star/private/pre_ms_model.f90:739: s% net_iso(ih1), x, standard_cgrav, m, gradL_composition_term, s% mixing_length_alpha, &
./star/private/adjust_xyz.f90:115: call realloc(s% d_eps_grav_dX); if (ierr /= 0) return
./star/private/atm_support.f90:53: s, tau_surf, L, R, M, cgrav, skip_partials, &
./star/private/atm_support.f90:60: real(dp), intent(in) :: tau_surf, L, R, M, cgrav
./star/private/atm_support.f90:95: s, tau_surf, L, R, M, cgrav, &
./star/private/atm_support.f90:105: s, skip_partials, L, R, M, cgrav, &
./star/private/atm_support.f90:114: s, s% atm_irradiated_opacity, skip_partials, L, R, M, cgrav, &
./star/private/atm_support.f90:133: s, tau_surf, L, R, M, cgrav, skip_partials, &
./star/private/atm_support.f90:140: real(dp), intent(in) :: tau_surf, L, R, M, cgrav
./star/private/atm_support.f90:160: s, tau_surf, L, R, M, cgrav, 'Eddington', 'iterated', skip_partials, &
./star/private/atm_support.f90:246: s, tau_surf, L, R, M, cgrav, T_tau_relation, T_tau_opacity, skip_partials, &
./star/private/atm_support.f90:258: real(dp), intent(in) :: tau_surf, L, R, M, cgrav
./star/private/atm_support.f90:311: tau_surf, L, R, M, cgrav, s% opacity(1), s% Pextra_factor, &
./star/private/atm_support.f90:335: tau_surf, L, R, M, cgrav, kap_guess, s% Pextra_factor, &
./star/private/atm_support.f90:353: tau_surf, L, R, M, cgrav, &
./star/private/atm_support.f90:465: s, skip_partials, L, R, M, cgrav, &
./star/private/atm_support.f90:479: real(dp), intent(in) :: L, R, M, cgrav
./star/private/atm_support.f90:558: L, Teff, R, M, cgrav, table_id, alfa, beta, ierr)
./star/private/atm_support.f90:610: L, R, M, cgrav, table_id, Z, skip_partials, &
./star/private/atm_support.f90:647: s, s% tau_base, L, R, M, cgrav, &
./star/private/atm_support.f90:779: s, irradiated_opacity, skip_partials, L, R, M, cgrav, &
./star/private/atm_support.f90:791: real(dp), intent(in) :: L, R, M, cgrav
./star/private/atm_support.f90:837: L, R, M, cgrav, s% atm_irradiated_T_eq, s% atm_irradiated_P_surf, &
./star/private/atm_support.f90:855: L, R, M, cgrav, s% atm_irradiated_T_eq, s% atm_irradiated_P_surf, &
./star/private/atm_support.f90:1051: s, L, R, Teff, M, cgrav, ierr)
./star/private/atm_support.f90:1054: real(dp), intent(in) :: L, R, Teff, M, cgrav
./star/private/atm_support.f90:1065: s, s% tau_factor*s% tau_base, L, R, Teff, M, cgrav, &
./star/private/atm_support.f90:1084: s, tau_surf, L, R, Teff, M, cgrav, T_tau_relation, T_tau_opacity, &
./star/private/atm_support.f90:1093: real(dp), intent(in) :: tau_surf, L, R, Teff, M, cgrav
./star/private/atm_support.f90:1118: s, tau_surf, L, R, M, cgrav, T_tau_relation, T_tau_opacity, .TRUE., &
./star/private/atm_support.f90:1146: tau_surf, L, R, Teff, M, cgrav, kap, s% Pextra_factor, s% atm_build_tau_outer, &
./star/private/atm_support.f90:1160: tau_surf, L, R, Teff, M, cgrav, lnP_surf, s% atm_build_tau_outer, &
./star/private/solver_support.f90:61: s% x_scale(i,k) = 10d0*sqrt(s% cgrav(k)*s% m(k)*s% r_start(k))
./star/private/solver_support.f90:843: skip_m_grav_and_grav = .true., &
./star/private/solver_support.f90:847: skip_other_cgrav = .true.
./star/private/solver_support.f90:1090: ! note: m_grav is held constant during solver iterations
./star/private/solver_support.f90:1091: s% grav(k) = s% cgrav(k)*s% m_grav(k)/(s% r(k)*s% r(k))
./star/private/solver_support.f90:1105: skip_m_grav_and_grav, skip_eos, skip_net, skip_neu, skip_kap, &
./star/private/solver_support.f90:1106: skip_grads, skip_rotation, skip_brunt, skip_other_cgrav, &
./star/private/solver_support.f90:1515: write(*,2) 'eps_grav', k, s% eps_grav_ad(k)% val
./star/private/hydro_rsp2.f90:257: Hp_face = pow2(r_00)*Peos_div_rho/(s% cgrav(k)*s% m(k))
./star/private/hydro_rsp2.f90:263: Hp_face = pow2(r_00)*Peos_div_rho/(s% cgrav(k)*s% m(k))
./star/private/hydro_rsp2.f90:268: ! Hp_face%val, pow2(r_00%val)/(s% cgrav(k)*s% m(k)), &
./star/private/hydro_rsp2.f90:276: alt_Hp_face = sqrt(Peos_face/s% cgrav(k))/d_face
./star/private/hydro_rsp2.f90:634: real(dp) :: mmid, cgrav_mid
./star/private/hydro_rsp2.f90:647: cgrav_mid = 0.5d0*(s% cgrav(k) + s% cgrav(k+1))
./star/private/hydro_rsp2.f90:651: cgrav_mid = s% cgrav(k)
./star/private/hydro_rsp2.f90:653: Hp_cell = pow2(rmid)*Peos_00/(d_00*cgrav_mid*mmid)
./star/private/report.f90:50: vesc = sqrt(2*s% cgrav(k)*s% m(k)/(s% r(k)))
./star/private/report.f90:244: s% kh_timescale = eval_kh_timescale(s% cgrav(1), mstar, radius, luminosity)/secyer
./star/private/report.f90:249: write(*,1) 's% cgrav(1)', s% cgrav(1)
./star/private/report.f90:262: if (s% cgrav(1) <= 0d0) then
./star/private/report.f90:265: s% dynamic_timescale = 2*pi*sqrt(radius*radius*radius/(s% cgrav(1)*mstar))
./star/private/report.f90:381: 0.25d6/pi*s% grav(1)*sqrt(s% gamma1(1)*s% rho(1)/s% Peos(1))
./star/private/report.f90:1319: bdy_g = s% grav(1)
./star/private/report.f90:1355: bdy_g = interp2(s% grav(k), s% grav(k+1))
./star/private/winds.f90:429: kh_timescale = eval_kh_timescale(s% cgrav(1), M1, R1, L1)
./star/private/winds.f90:630: vesc2 = 2d0 * s% cgrav(1)*M/R
./star/private/winds.f90:654: scale_height = s% Peos(1) / (s% cgrav(1)*s% m(1)*s% rho(1) / (s% r(1)**2)) / Rsun
./star/private/read_model.f90:79: use star_utils, only: set_m_and_dm, set_m_grav_and_grav, set_dm_bar, &
./star/private/read_model.f90:103: call set_m_grav_and_grav(s)
./star/private/read_model.f90:144: call fill_ad_with_zeros(s% eps_grav_ad,1,-1)
./star/private/hydro_energy.f90:68: use eps_grav, only: eval_eps_grav_and_partials
./star/private/hydro_energy.f90:79: dwork_dm_ad, eps_grav_ad, dke_dt_ad, dpe_dt_ad, de_dt_ad
./star/private/hydro_energy.f90:85: logical :: test_partials, doing_op_split_burn, eps_grav_form
./star/private/hydro_energy.f90:95: call setup_eps_grav(ierr); if (ierr /= 0) return ! do this first - it sets eps_grav_form
./star/private/hydro_energy.f90:110: if (eps_grav_form) then ! for this case, dwork_dm doesn't include work by P since that is in eps_grav
./star/private/hydro_energy.f90:112: others_ad - d_turbulent_energy_dt_ad - dwork_dm_ad + eps_grav_ad
./star/private/hydro_energy.f90:136: if (eps_grav_form) write(*,*) 'eps_grav_form', eps_grav_form
./star/private/hydro_energy.f90:177: skip_P = eps_grav_form
./star/private/hydro_energy.f90:337: subroutine setup_eps_grav(ierr)
./star/private/hydro_energy.f90:342: if (s% u_flag) then ! for now, assume u_flag means no eps_grav
./star/private/hydro_energy.f90:343: eps_grav_form = .false.
./star/private/hydro_energy.f90:348: eps_grav_form = s% eps_grav_form_for_energy_eqn
./star/private/hydro_energy.f90:350: if (.not. eps_grav_form) then ! check if want it true
./star/private/hydro_energy.f90:351: if (s% doing_relax .and. s% no_dedt_form_during_relax) eps_grav_form = .true.
./star/private/hydro_energy.f90:354: if (eps_grav_form) then
./star/private/hydro_energy.f90:356: call mesa_error(__FILE__,__LINE__,'cannot use eps_grav with et yet. fix energy eqn.')
./star/private/hydro_energy.f90:358: call eval_eps_grav_and_partials(s, k, ierr) ! get eps_grav info
./star/private/hydro_energy.f90:360: if (s% report_ierr) write(*,2) 'failed in eval_eps_grav_and_partials', k
./star/private/hydro_energy.f90:363: eps_grav_ad = s% eps_grav_ad(k)
./star/private/hydro_energy.f90:366: end subroutine setup_eps_grav
./star/private/hydro_energy.f90:381: if (.not. eps_grav_form) then
./star/private/hydro_energy.f90:439: if (.not. eps_grav_form) then
./star/private/hydro_energy.f90:448: dequ = scal*s% d_eps_grav_dx(j,k)
./star/private/hydro_energy.f90:449: if (checking) call check_dequ(dequ,'d_eps_grav_dx')
./star/private/mesh_adjust.f90:64: use star_utils, only: set_m_grav_and_grav
./star/private/mesh_adjust.f90:325: ! ensure the things we need to calculate m_grav and grav are up to date
./star/private/mesh_adjust.f90:342: call set_m_grav_and_grav(s)
./star/private/mesh_adjust.f90:2112: w_div_w_roche_jrot(r00,s% m(k),s% j_rot(k),s% cgrav(k), &
./star/private/element_diffusion.f90:415: write(*,*) "g_field_element_diffusion/grav = ", s% g_field_element_diffusion(k) / s% grav(k)
./star/private/element_diffusion.f90:567: real(dp) :: grav, area, P_face
./star/private/element_diffusion.f90:574: grav = -s% cgrav(k)*s% m(k)/s% r(k)**2
./star/private/element_diffusion.f90:577: dlnPdm(k) = grav/(area*P_face) ! estimate based on QHSE
./star/private/rsp_step.f90:552: ! EGRV = EGRV - s% cgrav(k) * (s%m(k)-0.5d0*s%dm(k))*s%dm(k)/(0.5d0*(s%r(k)+s%r(k+1)))
./star/private/rsp_step.f90:554: ! EGRV = EGRV - s% cgrav(k) * (s%m(k)-0.5d0*s%dm(k))*s%dm(k)/(0.5d0*(s%r(k)+s%r_center))
./star/private/rsp_step.f90:1162: POM = (s% r(k)**2)/(2.d0*s% cgrav(k)*s% m(k))
./star/private/rsp_step.f90:1194: POM = (s% r(k)**2)/(s% cgrav(k)*s% M(k))
./star/private/rsp_step.f90:2704: grav, dXP_dm, Uq1, A_dm, R_00, dv_dr, dvdt_factor, &
./star/private/rsp_step.f90:2760: grav = - s% cgrav(k)*s% m(k)/(s% r(k)*s% r_start(k))
./star/private/rsp_step.f90:2869: + area*dXP_dm - grav - s% Uq(k) - Fr_term
./star/private/rsp_step.f90:2903: + grav/s% r(k) &
./star/private/rsp_step.f90:2990: write(*,2) 'grav', k, grav
./star/private/star_utils.f90:259: subroutine set_m_grav_and_grav(s) ! using mass_corrections
./star/private/star_utils.f90:267: s% m_grav(k) = s% m(k)
./star/private/star_utils.f90:270: s% m_grav(nz) = &
./star/private/star_utils.f90:273: s% m_grav(k) = &
./star/private/star_utils.f90:274: s% m_grav(k+1) + s% dm(k)*s% mass_correction(k)
./star/private/star_utils.f90:278: ! We need to call set_m_grav_and_grav during model loading before we have set all vars
./star/private/star_utils.f90:280: s% grav(k) = s% cgrav(k)*s% m_grav(k)/(r*r)
./star/private/star_utils.f90:282: end subroutine set_m_grav_and_grav
./star/private/star_utils.f90:985: logg = safe_log10(s% cgrav(1)*m/(r*r))
./star/private/star_utils.f90:1020: logg = safe_log10(s% cgrav(k_phot)*m/(r*r))
./star/private/star_utils.f90:1038: logg = safe_log10(s% cgrav(k_phot)*m/(r*r))
./star/private/star_utils.f90:1558: if (s% cgrav(k) == 0) then
./star/private/star_utils.f90:1575: Hp = Peos_face/(rho_face*s% grav(k))
./star/private/star_utils.f90:1576: if (s% cgrav(k) <= 0) then
./star/private/star_utils.f90:1579: alt_Hp = sqrt(Peos_face / s% cgrav(k)) / rho_face
./star/private/star_utils.f90:1595: if (s% cgrav(k) <= 0d0) then
./star/private/star_utils.f90:1600: g = s% cgrav(k)*s% m_grav(k)/(s% r(k)*s% r(k))
./star/private/star_utils.f90:1603: Pextra_factor*(s% L(k)/s% m_grav(k))/(6d0*pi*clight*s% cgrav(k)))
./star/private/star_utils.f90:1616: if (s% cgrav(1) <= 0d0) return
./star/private/star_utils.f90:1622: Ledd_sum = s% dq(1)*pi4*clight*s% cgrav(1)*s% m_grav(1)/s% opacity(1)
./star/private/star_utils.f90:1629: s% dq(k)*pi4*clight*s% cgrav(1)*s% m_grav(1)/s% opacity(k)
./star/private/star_utils.f90:1781: vesc = sqrt(2d0*s% cgrav(k)*s% m(k)/s% r(k))
./star/private/star_utils.f90:1910: Ledd = pi4*clight*s% cgrav(j)*s% m_grav(j)/kap_face
./star/private/star_utils.f90:1959: if (s% cgrav(k) <= 0) then
./star/private/star_utils.f90:1972: -(area*area*crad*(del_T4/del_m)/3)/(pi4*s% cgrav(j)*s% m_grav(j))
./star/private/star_utils.f90:2076: real(qp) :: qhalf, rp1, r00, mp1, m00, Gp1, G00, gravp1, grav00, Mbar
./star/private/star_utils.f90:2077: real(dp) :: d_grav00_dlnR00, d_gravp1_dlnRp1
./star/private/star_utils.f90:2088: Gp1 = s% cgrav(s% nz)
./star/private/star_utils.f90:2091: mp1 = s% m_grav(k+1)
./star/private/star_utils.f90:2092: Gp1 = s% cgrav(k+1)
./star/private/star_utils.f90:2095: gravp1 = 0d0
./star/private/star_utils.f90:2096: d_gravp1_dlnRp1 = 0d0
./star/private/star_utils.f90:2098: gravp1 = -Gp1*mp1/rp1
./star/private/star_utils.f90:2099: d_gravp1_dlnRp1 = -gravp1
./star/private/star_utils.f90:2102: m00 = s% m_grav(k)
./star/private/star_utils.f90:2103: G00 = s% cgrav(k)
./star/private/star_utils.f90:2104: grav00 = -G00*m00/r00
./star/private/star_utils.f90:2105: d_grav00_dlnR00 = -grav00
./star/private/star_utils.f90:2106: cell_specific_PE_qp = qhalf*Mbar*(gravp1 + grav00)
./star/private/star_utils.f90:2107: d_dlnR00 = qhalf*Mbar*d_grav00_dlnR00
./star/private/star_utils.f90:2108: d_dlnRp1 = qhalf*Mbar*d_gravp1_dlnRp1
./star/private/star_utils.f90:2111: write(*,2) 'gravp1', k, gravp1
./star/private/star_utils.f90:2112: write(*,2) 'grav00', k, grav00
./star/private/star_utils.f90:2130: real(qp) :: qhalf, rp1, r00, mp1, m00, Gp1, G00, gravp1, grav00, Mbar
./star/private/star_utils.f90:2141: Gp1 = s% cgrav(s% nz)
./star/private/star_utils.f90:2144: mp1 = s% m_grav_start(k+1)
./star/private/star_utils.f90:2145: Gp1 = s% cgrav(k+1)
./star/private/star_utils.f90:2148: gravp1 = 0d0
./star/private/star_utils.f90:2150: gravp1 = -Gp1*mp1/rp1
./star/private/star_utils.f90:2153: m00 = s% m_grav_start(k)
./star/private/star_utils.f90:2154: G00 = s% cgrav(k)
./star/private/star_utils.f90:2155: grav00 = -G00*m00/r00
./star/private/star_utils.f90:2156: cell_start_specific_PE_qp = qhalf*Mbar*(gravp1 + grav00)
./star/private/star_utils.f90:2159: write(*,2) 'gravp1', k, gravp1
./star/private/star_utils.f90:2160: write(*,2) 'grav00', k, grav00
./star/private/star_utils.f90:2281: total_internal_energy, total_gravitational_energy, &
./star/private/star_utils.f90:2290: total_internal_energy, total_gravitational_energy, &
./star/private/star_utils.f90:2299: total_gravitational_energy = 0d0
./star/private/star_utils.f90:2323: total_gravitational_energy = total_gravitational_energy + cell1
./star/private/star_utils.f90:2345: sum_total = total_internal_energy + total_gravitational_energy + &
./star/private/star_utils.f90:2398: total_internal_energy, total_gravitational_energy, &
./star/private/star_utils.f90:2406: total_internal_energy, total_gravitational_energy, &
./star/private/star_utils.f90:2413: total_internal_energy, total_gravitational_energy, &
./star/private/star_utils.f90:2418: total_internal_energy, total_gravitational_energy, &
./star/private/star_utils.f90:2426: total_internal_energy, total_gravitational_energy, &
./star/private/star_utils.f90:2434: total_internal_energy, total_gravitational_energy, &
./star/private/star_utils.f90:2439: total_internal_energy, total_gravitational_energy, &
./star/private/star_utils.f90:2445: total_internal_energy, total_gravitational_energy, &
./star/private/star_utils.f90:2456: total_internal_energy, total_gravitational_energy, &
./star/private/star_utils.f90:2464: total_internal_energy, total_gravitational_energy, &
./star/private/star_utils.f90:2475: total_internal_energy, total_gravitational_energy, &
./star/private/star_utils.f90:2483: total_internal_energy, total_gravitational_energy, &
./star/private/star_utils.f90:2600: Ledd = pi4*clight*s% cgrav(k)*s% m_grav(k)/s% opacity(k)
./star/private/star_utils.f90:2603: omega_crit = sqrt(gamma_factor*s% cgrav(k)*s% m_grav(k)/pow3(rmid))
./star/private/star_utils.f90:3569: f = pow2(s% grav(k))*rho_face/Peos_face*chiT_face/chiRho_face
./star/private/star_utils.f90:3617: tau_qhse = abs_dv/(s% cgrav(k)*s% m_grav(k)/pow2(s% r(k)))
./star/private/star_utils.f90:3782: gradr = P*opacity*L/(16d0*pi*clight*s% m_grav(k)*s% cgrav(k)*Pr)
./star/private/star_utils.f90:3790: type(auto_diff_real_star_order1) :: grav, scale_height2, P, rho
./star/private/star_utils.f90:3793: G = s% cgrav(k)
./star/private/star_utils.f90:3794: grav = G*s% m_grav(k)/pow2(wrap_r_00(s,k))
./star/private/star_utils.f90:3797: scale_height = P/(grav*rho) ! this assumes HSE
./star/private/star_utils.f90:3812: real(dp) :: G, grav, scale_height2, P, rho
./star/private/star_utils.f90:3814: G = s% cgrav(k)
./star/private/star_utils.f90:3815: grav = G*s% m_grav(k)/pow2(s% r(k))
./star/private/star_utils.f90:3820: scale_height = P/(grav*rho) ! this assumes HSE
./star/private/pulse_cafein.f90:97: call build_atm(s, s%L(1), s%r(1), s%Teff, s%m_grav(1), s%cgrav(1), ierr)
./star/private/pulse_cafein.f90:127: M_star = s%m_grav(1)
./star/private/pulse_cafein.f90:134: s_units = SQRT(pow3(R_star)/(s%cgrav(1)*M_star))
./star/private/pulse_cafein.f90:137: erg_units = s%cgrav(1)*M_star**2/R_star
./star/private/pulse_cafein.f90:138: K_units = s%cgrav(1)*M_star**2/(R_star*cgas)
./star/private/pulse_cafein.f90:343: m = s%m_grav(1) !+ s%atm_structure(atm_delta_m,k)
./star/private/pulse_cafein.f90:374: g = s%cgrav(1)*m/(r*r)
./star/private/pulse_cafein.f90:389: c_4 = pi4*r**3*rho*T*c_P/l_rad(j)*SQRT(s%cgrav(1)*M_star/R_star**3)
./star/private/pulse_cafein.f90:454: m = s%m_grav(k)
./star/private/pulse_cafein.f90:496: g = s%grav(k)
./star/private/pulse_cafein.f90:511: c_4 = pi4*r**3*rho*T*c_P/l_rad(j)*SQRT(s%cgrav(1)*M_star/R_star**3)
./star/private/pulse_cafein.f90:536: real(dp) :: cgrav
./star/private/pulse_cafein.f90:626: cgrav = eval_center(s%r, s%cgrav, 1, s%nz)
./star/private/pulse_cafein.f90:630: c_3 = 9d0*eps*kap*P/(16*pi*crad*clight*T**4*nabla*cgrav)
./star/private/pulse_cafein.f90:631: c_4 = 9d0*T*c_P*kap*P/(16*pi*crad*clight*T**4*nabla*cgrav)*SQRT(s%cgrav(s%nz)*M_star/R_star**3)
./star/private/rsp.f90:160: normalize_dqs, set_qs, set_m_and_dm, set_dm_bar, set_m_grav_and_grav, &
./star/private/rsp.f90:203: call set_m_grav_and_grav(s)
./star/private/rsp.f90:438: (s% Peos(k-1) - s% Peos(k))/(-s% cgrav(k)*s% m(k)*s% dm_bar(k)/(4d0*pi*s% r(k)**4)) - 1d0
./star/private/rsp.f90:958: total_internal_energy, total_gravitational_energy, &
./star/private/rsp.f90:963: total_internal_energy, total_gravitational_energy, &
./star/private/rsp.f90:969: total_gravitational_energy = EGRV
./star/private/history.f90:1181: use gravity_darkening
./star/private/history.f90:1326: !val = get_lum_band_by_id(i,safe_log10(s% T(k)),safe_log10(s% grav(k)),m_div_h,s% L(k)/lsun, ierr)
./star/private/history.f90:1338: !val = get_abs_mag_by_id(i,safe_log10(s% T(k)),safe_log10(s% grav(k)),m_div_h,s% L(k)/lsun, ierr)
./star/private/history.f90:1349: !val = get_bc_by_id(i,safe_log10(s% T(k)),safe_log10(s% grav(k)),m_div_h, ierr)
./star/private/history.f90:1479: case(h_star_gravitational_mass)
./star/private/history.f90:1480: val = s% m_grav(1)/Msun
./star/private/history.f90:1481: case(h_star_mass_grav_div_mass)
./star/private/history.f90:1482: val = s% m_grav(1)/s% m(1)
./star/private/history.f90:1770: case(h_gravity)
./star/private/history.f90:1771: val = s% grav(1)
./star/private/history.f90:1773: val = safe_log10(s% grav(1))
./star/private/history.f90:1816: case(h_total_gravitational_energy_after_adjust_mass)
./star/private/history.f90:1817: val = s% total_gravitational_energy_after_adjust_mass
./star/private/history.f90:1829: case(h_total_gravitational_energy)
./star/private/history.f90:1830: val = s% total_gravitational_energy_end
./star/private/history.f90:1844: case(h_log_total_gravitational_energy)
./star/private/history.f90:1845: val = safe_log10(abs(s% total_gravitational_energy_end))
./star/private/history.f90:1885: val = (val + s% total_gravitational_energy_end)/val
./star/private/history.f90:1887: case(h_total_eps_grav)
./star/private/history.f90:1888: val = s% total_eps_grav
./star/private/history.f90:1953: val = sqrt(2*s% cgrav(1)*s% m(1)/(s% r(1)))
./star/private/history.f90:1955: val = v_surf/sqrt(2*s% cgrav(1)*s% m(1)/(s% r(1)))
./star/private/history.f90:1957: val = v_surf/sqrt(2*s% cgrav(1)*s% m(1)/(s% r(1)))
./star/private/history.f90:2075: case(h_center_eps_grav)
./star/private/history.f90:2076: val = center_value(s, s% eps_grav_ad(1:nz)% val)
./star/private/history.f90:2294: case(h_eps_grav_integral)
./star/private/history.f90:2295: val = dot_product(s% dm(1:nz), s% eps_grav_ad(1:nz)% val)/Lsun
./star/private/history.f90:2300: case(h_log_abs_Lgrav)
./star/private/history.f90:2301: val = safe_log10(abs(dot_product(s% dm(1:nz), s% eps_grav_ad(1:nz)%val)/Lsun))
./star/private/history.f90:2539: case(h_cz_grav)
./star/private/history.f90:2542: val = s% grav(k)
./star/private/history.f90:2665: case(h_cz_top_grav)
./star/private/history.f90:2668: val = s% grav(k)
./star/private/history.f90:2864: val = s% total_gravitational_energy_end
./star/private/history.f90:2866: val = safe_log10(abs(s% total_gravitational_energy_end))
./star/private/history.f90:2900: case(h_grav_dark_L_polar) ! pole is at inclination = 0
./star/private/history.f90:2902: w_div_w_Kep = if_rot(s% omega(1)*sqrt(pow3(s% r_equatorial(1))/(s% cgrav(1)*s% m(1))))
./star/private/history.f90:2903: val = gravity_darkening_L_coeff(w_div_w_Kep, 0.0d0)*s% L_surf
./star/private/history.f90:2907: case(h_grav_dark_Teff_polar)
./star/private/history.f90:2909: w_div_w_Kep = if_rot(s% omega(1)*sqrt(pow3(s% r_equatorial(1))/(s% cgrav(1)*s% m(1))))
./star/private/history.f90:2910: val = gravity_darkening_Teff_coeff(w_div_w_Kep, 0.0d0)*s% Teff
./star/private/history.f90:2914: case(h_grav_dark_L_equatorial) ! equator is at inclination = pi/2
./star/private/history.f90:2916: w_div_w_Kep = if_rot(s% omega(1)*sqrt(pow3(s% r_equatorial(1))/(s% cgrav(1)*s% m(1))))
./star/private/history.f90:2917: val = gravity_darkening_L_coeff(w_div_w_Kep, 0.5d0*pi) * s% L_surf
./star/private/history.f90:2921: case(h_grav_dark_Teff_equatorial)
./star/private/history.f90:2923: w_div_w_Kep = if_rot(s% omega(1)*sqrt(pow3(s% r_equatorial(1))/(s% cgrav(1)*s% m(1))))
./star/private/history.f90:2924: val = gravity_darkening_Teff_coeff(w_div_w_Kep, 0.5d0*pi) * s% Teff
./star/private/alloc.f90:467: use const_def, only: standard_cgrav
./star/private/alloc.f90:650: call do1(s% m_grav, c% m_grav)
./star/private/alloc.f90:651: if (failed('m_grav')) exit
./star/private/alloc.f90:652: call do1(s% m_grav_start, c% m_grav_start)
./star/private/alloc.f90:653: if (failed('m_grav_start')) exit
./star/private/alloc.f90:764: call do1(s% grav, c% grav)
./star/private/alloc.f90:765: if (failed('grav')) exit
./star/private/alloc.f90:937: call do1_ad(s% extra_grav, c% extra_grav)
./star/private/alloc.f90:938: if (failed('extra_grav')) exit
./star/private/alloc.f90:958: call do1(s% cgrav, c% cgrav)
./star/private/alloc.f90:959: if (failed('cgrav')) exit
./star/private/alloc.f90:962: s% cgrav(1:nz) = standard_cgrav
./star/private/alloc.f90:982: call do1_ad(s% eps_grav_ad, c% eps_grav_ad)
./star/private/alloc.f90:983: if (failed('eps_grav_ad')) exit
./star/private/alloc.f90:984: call do2(s% d_eps_grav_dx, c% d_eps_grav_dx, species, null_str)
./star/private/alloc.f90:985: if (failed('d_eps_grav_dx')) exit
./star/private/alloc.f90:987: call do1(s% eps_grav_composition_term, c% eps_grav_composition_term)
./star/private/alloc.f90:988: if (failed('eps_grav_composition_term')) exit
./star/private/rsp_eval_eos_and_kap.f90:834: s, tau, L, R, M, s% cgrav(1), .TRUE., &
./star/other/other_mlt_results.f90:43: iso, XH1, cgrav, m, gradL_composition_term, mixing_length_alpha, &
./star/other/other_mlt_results.f90:55: XH1, cgrav, m, gradL_composition_term, &
./star/other/other_momentum_implicit.f90:49: ! s% extra_grav(1:s%nz) = 0 ! this is added to -G*m/r^2
./star/other/other_momentum_implicit.f90:50: ! note that extra_grav is type(auto_diff_real_star_order1) so includes partials.
./star/other/other_cgrav.f90:26: module other_cgrav
./star/other/other_cgrav.f90:29: ! control name: use_other_cgrav = .true.
./star/other/other_cgrav.f90:30: ! procedure pointer: s% other_cgrav => my_routine
./star/other/other_cgrav.f90:32: ! note that other_cgrav only changes G in the stellar structure
./star/other/other_cgrav.f90:33: ! the binary module is unaffected by changes in cgrav
./star/other/other_cgrav.f90:43: subroutine default_other_cgrav(id, ierr)
./star/other/other_cgrav.f90:44: use const_def, only: standard_cgrav
./star/other/other_cgrav.f90:52: s% cgrav(:) = standard_cgrav
./star/other/other_cgrav.f90:53: end subroutine default_other_cgrav
./star/other/other_cgrav.f90:55: end module other_cgrav
./star/other/other_eps_grav.f90:26: module other_eps_grav
./star/other/other_eps_grav.f90:29: ! control name: use_other_eps_grav = .true.
./star/other/other_eps_grav.f90:30: ! procedure pointer: s% other_eps_grav => my_routine
./star/other/other_eps_grav.f90:41: subroutine null_other_eps_grav(id, k, dt, ierr)
./star/other/other_eps_grav.f90:47: ! NOTE: this is called after 1st doing the standard eps_grav calculation.
./star/other/other_eps_grav.f90:50: ! NOTE: need to set the auto_diff variable s% eps_grav_ad(k)
./star/other/other_eps_grav.f90:55: end subroutine null_other_eps_grav
./star/other/other_eps_grav.f90:58: end module other_eps_grav
./star/other/other_momentum.f90:50: !s% extra_grav(k) = 0d0 ! this is added to -G*m/r^2
./star/test_suite/accreted_material_j/src/run_star_extras.f90:72: s% x_ctrl(1)*sqrt(s% cgrav(1) * s% mstar * s% photosphere_r*Rsun)
./star/test_suite/magnetic_braking/src/run_star_extras.f90:108: vinf = 1.92d0 * sqrt(2.0d0 * standard_cgrav * s% mstar / (s% photosphere_r * Rsun))
./star/test_suite/carbon_kh/src/run_star_extras.f90:214: names(4) = 'rse_eps_grav'
./star/test_suite/carbon_kh/src/run_star_extras.f90:245: ! difference between dedt and eps_grav - style estimate
./star/test_suite/gyre_in_mesa_ms/src/run_star_extras.f90:88: call gyre_set_constant('G_GRAVITY', standard_cgrav)
./star/test_suite/12M_pre_ms_to_core_collapse/src/run_star_extras.f90:130: call gyre_set_constant('G_GRAVITY', standard_cgrav)
./star/test_suite/ns_c/src/run_star_extras.f90:39: include "other_cgrav.inc"
./star/test_suite/ns_c/src/run_star_extras.f90:56: s% other_cgrav => my_other_cgrav
./star/test_suite/ppisn/src/run_star_extras.f90:461: ! if (s% xh_start(s% i_u,1)>sqrt(2*s% cgrav(1)*s% m(1)/exp(s% lnR_start(1)))) then
./star/test_suite/ppisn/src/run_star_extras.f90:472: if (velocity>sqrt(2*s% cgrav(1)*s% m(1)/radius)) then
./star/test_suite/ppisn/src/run_star_extras.f90:635: names(h_Omega_ejecta) = "Omega_ejecta" ! gravitational energy of layers above vesc
./star/test_suite/ppisn/src/run_star_extras.f90:689: v_esc = sqrt(2*s% cgrav(k)*s% m(k)/(s% r(k)))
./star/test_suite/ppisn/src/run_star_extras.f90:728: vals(h_Omega_ejecta) = vals(h_Omega_ejecta) - s% dm_bar(k0)*s% cgrav(k0)*s% m(k0)/s% r(k0)
./star/test_suite/ppisn/src/run_star_extras.f90:774: names(3) = "specific_grav_e"
./star/test_suite/ppisn/src/run_star_extras.f90:782: vals(k,1) = sqrt(2*s% cgrav(k)*s% m(k)/(s% r(k)))
./star/test_suite/ppisn/src/run_star_extras.f90:784: vals(k,3) = -s% cgrav(k)*s% m(k)/s% r(k)
./star/test_suite/ppisn/src/run_star_extras.f90:792: vals(k,8) = J_inside*clight/(pow2(s% m(k))*standard_cgrav)
./star/test_suite/ppisn/src/run_star_extras.f90:803: vals(k,1) = J_inside*clight/(pow2(s% m(k))*standard_cgrav)
./star/test_suite/ppisn/src/run_star_extras.f90:875: v_esc = sqrt(2*s% cgrav(k)*s% m(k)/(s% r(k)))
./star/test_suite/ppisn/src/run_star_extras.f90:899: avg_v_div_vesc = avg_v_div_vesc + s% dm(k0)*s% u(k0)/sqrt(2*s% cgrav(k0)*s% m(k0)/(s% r(k0)))
./star/test_suite/ppisn/src/run_star_extras.f90:901: 0.5d0*s% dm(k0)*s% u(k0)*s% u(k0) - s% dm(k0)*s% cgrav(k0)*s% m(k0)/s% r(k0) &
./star/test_suite/ppisn/src/run_star_extras.f90:913: 0.5d0*s% u(k0)*s% u(k0) - s% cgrav(k0)*s% m(k0)/s% r(k0) &
./star/test_suite/ppisn/src/run_star_extras.f90:921: s% dm(k0)*s% u(k0)/sqrt(2*s% cgrav(k0)*s% m(k0)/(s% r(k0)))
./star/test_suite/ppisn/src/run_star_extras.f90:1004: v_esc = sqrt(2*s% cgrav(k0)*s% m(k0)/(s% r(k0)))
./star/test_suite/ppisn/src/run_star_extras.f90:1201: tdyn = 1d0/sqrt(s% m(k)/(4d0/3d0*pi*pow3(s% r(k)))*standard_cgrav)
./star/test_suite/ppisn/src/run_star_extras.f90:1249: v_esc = sqrt(2*s% cgrav(k)*s% m(k)/(s% r(k)))
./star/test_suite/gyre_in_mesa_spb/src/run_star_extras.f90:88: call gyre_set_constant('G_GRAVITY', standard_cgrav)
./star/test_suite/ns_he/src/run_star_extras.f90:39: include "other_cgrav.inc"
./star/test_suite/ns_he/src/run_star_extras.f90:56: s% other_cgrav => my_other_cgrav
./star/test_suite/ns_h/src/run_star_extras.f90:39: include "other_cgrav.inc"
./star/test_suite/ns_h/src/run_star_extras.f90:56: s% other_cgrav => my_other_cgrav
./star/test_suite/extended_convective_penetration/src/run_star_extras.f90:451: Hp = P/(rho*s%cgrav(k)* (s%M_center + s%xmstar*s%q(k))/(r*r))
./star/test_suite/zams_to_cc_80/src/run_star_extras.f90:132: call gyre_set_constant('G_GRAVITY', standard_cgrav)
./star/test_suite/gyre_in_mesa_wd/src/run_star_extras.f90:88: call gyre_set_constant('G_GRAVITY', standard_cgrav)
./star/test_suite/ccsn_IIp/src/run_star_extras.f90:713: !windv = sqrt(2*s% cgrav(kk)*s% m(kk)/r0) ! escape velocity
./star/test_suite/ccsn_IIp/src/run_star_extras.f90:760: s% cgrav(k+1)*s% m(k+1)*(s% dm(k+1)+s% dm(k))/(8*pi*r*r*r*r)
./star/test_suite/ccsn_IIp/src/run_star_extras.f90:1167: s% v_center - s% cgrav(s% nz)*s% M_center/(s% R_center*s% R_center)
./star/test_suite/make_env/src/run_star_extras.f90:183: kap, dlnkap_dlnd, dlnkap_dlnT, grav, rho_00, vol, &
./star/test_suite/make_env/src/run_star_extras.f90:222: s% m_grav(1) = s% mstar
./star/test_suite/make_env/src/run_star_extras.f90:223: s% cgrav(1) = standard_cgrav
./star/test_suite/make_env/src/run_star_extras.f90:285: s% m_grav(k) = s% m(k)
./star/test_suite/make_env/src/run_star_extras.f90:381: s% id, tau_surf, s% L(1), s% r(1), s% m(1), s% cgrav(1), skip_partials, &
./star/test_suite/make_env/src/run_star_extras.f90:450: s% cgrav(k) = standard_cgrav
./star/test_suite/make_env/src/run_star_extras.f90:451: grav = -s% cgrav(k)*s% m(k)/r_00**2
./star/test_suite/make_env/src/run_star_extras.f90:460: P_00 = P_m1 - grav*dm_face/area
./star/test_suite/custom_colors/src/run_star_extras.f90:178: safe_log10(s% grav(1)),&
./star/test_suite/gyre_in_mesa_rsg/src/run_star_extras.f90:88: call gyre_set_constant('G_GRAVITY', standard_cgrav)
./star/test_suite/rsp_gyre/src/run_star_extras.f90:102: call gyre_set_constant('G_GRAVITY', standard_cgrav)
./star/test_suite/gyre_in_mesa_envelope/src/run_star_extras.f90:88: call gyre_set_constant('G_GRAVITY', standard_cgrav)
./star/test_suite/simplex_solar_calibration/src/simplex_search_run_support.f90:1184: logg = log10(s% grav(1))
./star/test_suite/gyre_in_mesa_bcep/src/run_star_extras.f90:88: call gyre_set_constant('G_GRAVITY', standard_cgrav)
./turb/test/src/test_turb.f90:26: type(auto_diff_real_star_order1) :: chiT, chiRho, Cp, grav, Lambda, rho, P, T, opacity, gradr, grada, gradL
./turb/test/src/test_turb.f90:45: grav = 1d4
./turb/test/src/test_turb.f90:49: Lambda = P / (rho * grav)
./turb/test/src/test_turb.f90:55: gradr = 3d0 * P * opacity * L / (64 * pi * boltz_sigma * pow4(T) * grav * pow2(r))
./turb/test/src/test_turb.f90:58: chiT, chiRho, Cp, grav, Lambda, rho, P, T, opacity, &
./turb/test/src/test_turb.f90:65: gradr = 3d0 * P * opacity * L / (64 * pi * boltz_sigma * pow4(T) * grav * pow2(r))
./turb/test/src/test_turb.f90:68: chiT, chiRho, Cp, grav, Lambda, rho, P, T, opacity, &
./turb/test/src/test_turb.f90:79: real(dp) :: mixing_length_alpha, conv_vel_start, alpha_TDC_DAMP, alpha_TDC_DAMPR, alpha_TDC_PtdVdt, dt, cgrav, m, scale
./turb/test/src/test_turb.f90:81: r, L, T, P, opacity, rho, dV, chiRho, chiT, Cp, gradr, grada, scale_height, gradL, grav, Lambda
./turb/test/src/test_turb.f90:102: cgrav = standard_cgrav
./turb/test/src/test_turb.f90:103: grav = m * cgrav / pow2(r)
./turb/test/src/test_turb.f90:106: scale_height = P / (rho * grav)
./turb/test/src/test_turb.f90:113: gradr = 3d0 * P * opacity * L / (64 * pi * boltz_sigma * pow4(T) * grav * pow2(r))
./turb/test/src/test_turb.f90:117: gradr = 3d0 * P * opacity * L / (64 * pi * boltz_sigma * pow4(T) * grav * pow2(r))
./turb/test/src/test_turb.f90:138: conv_vel_start, mixing_length_alpha, alpha_TDC_DAMP, alpha_TDC_DAMPR, alpha_TDC_PtdVdt, dt, cgrav, m, report, &
./turb/test/src/test_turb.f90:145: chiT, chiRho, Cp, grav, Lambda, rho, P, T, opacity, &
./turb/test/src/test_turb.f90:155: real(dp) :: mixing_length_alpha, conv_vel_start, alpha_TDC_DAMP, alpha_TDC_DAMPR, alpha_TDC_PtdVdt, dt, cgrav, m, scale
./turb/test/src/test_turb.f90:173: cgrav=6.6743000000000004d-8
./turb/test/src/test_turb.f90:190: gradr = 3d0 * P * opacity * L / (64 * pi * boltz_sigma * pow4(T) * cgrav * m)
./turb/test/src/test_turb.f90:198: conv_vel_start, mixing_length_alpha, alpha_TDC_DAMP, alpha_TDC_DAMPR, alpha_TDC_PtdVdt, dt, cgrav, m, report, &
./turb/public/turb.f90:68: !! @param cgrav gravitational constant (erg*cm/g^2).
./turb/public/turb.f90:93: conv_vel_start, mixing_length_alpha, alpha_TDC_DAMP, alpha_TDC_DAMPR, alpha_TDC_PtdVdt, dt, cgrav, m, report, &
./turb/public/turb.f90:98: real(dp), intent(in) :: conv_vel_start, mixing_length_alpha, alpha_TDC_DAMP, alpha_TDC_DAMPR, alpha_TDC_PtdVdt, dt, cgrav, m, scale
./turb/public/turb.f90:105: type(auto_diff_real_star_order1) :: L, grav, Lambda, Gamma
./turb/public/turb.f90:114: grav = cgrav * m / pow2(r)
./turb/public/turb.f90:115: L = 64 * pi * boltz_sigma * pow4(T) * grav * pow2(r) * gradr / (3d0 * P * opacity)
./turb/public/turb.f90:118: chiT, chiRho, Cp, grav, Lambda, rho, P, T, opacity, &
./turb/public/turb.f90:133: info%L0 = convert((16d0*pi*crad*clight/3d0)*cgrav*m*pow4(T)/(P*opacity)) ! assumes QHSE for dP/dm
./turb/public/turb.f90:171: !! @param cgrav gravitational constant (erg*cm/g^2).
./turb/public/turb.f90:184: semiconvection_option, cgrav, Cp, gradr, grada, gradL, &
./turb/public/turb.f90:191: real(dp), intent(in) :: alpha_semiconvection, cgrav, gradL_composition_term, m
./turb/public/turb.f90:196: semiconvection_option, cgrav, Cp, gradr, grada, gradL, &
./turb/public/turb.f90:210: !! @param grav The acceleration due to gravity (cm/s^2).
./turb/public/turb.f90:226: chiT, chiRho, Cp, grav, Lambda, rho, P, T, opacity, &
./turb/public/turb.f90:230: type(auto_diff_real_star_order1), intent(in) :: chiT, chiRho, Cp, grav, Lambda, rho, P, T, opacity, gradr, grada, gradL
./turb/public/turb.f90:237: chiT, chiRho, Cp, grav, Lambda, rho, P, T, opacity, &
./turb/private/semiconvection.f90:49: !! @param cgrav gravitational constant (erg*cm/g^2).
./turb/private/semiconvection.f90:62: semiconvection_option, cgrav, Cp, gradr, grada, gradL, &
./turb/private/semiconvection.f90:72: real(dp), intent(in) :: alpha_semiconvection, cgrav, gradL_composition_term, m
./turb/private/semiconvection.f90:105: LG = (16d0*pi*clight*m*cgrav*Pr)/(P*opacity)
./turb/private/mlt.f90:45: !! @param grav The acceleration due to gravity (cm/s^2).
./turb/private/mlt.f90:61: chiT, chiRho, Cp, grav, Lambda, rho, P, T, opacity, &
./turb/private/mlt.f90:68: type(auto_diff_real_star_order1), intent(in) :: chiT, chiRho, Cp, grav, Lambda, rho, P, T, opacity, gradr, grada, gradL
./turb/private/mlt.f90:87: Cp*grav*pow2(Lambda)*rho*(sqrt(Q*rho/(2d0*P)))/9d0 ! erg / (K cm sec)
./colors/public/colors_lib.f90:128: real(dp), intent(in) :: log_g ! log_10 of surface gravity
./colors/public/colors_lib.f90:172: real(dp), intent(in) :: log_g ! log_10 of surface gravity
./colors/public/colors_lib.f90:261: real(dp), intent(in) :: log_g ! log_10 of surface gravity
./colors/public/colors_lib.f90:282: real(dp), intent(in) :: log_g ! log_10 of surface gravity
./colors/public/colors_lib.f90:360: real(dp), intent(in) :: log_g ! log_10 of surface gravity
./colors/public/colors_lib.f90:389: real(dp), intent(in) :: log_g ! log_10 of surface gravity
./const/public/const_def.f90:117: real(dp), parameter :: standard_cgrav = 6.67430d-8 ! gravitational constant (g^-1 cm^3 s^-2)
./const/public/const_def.f90:120: ! standard gravitational parameters = G*M, units cm^3 s^-2
./const/public/const_def.f90:127: real(dp), parameter :: Msun = mu_sun / standard_cgrav ! solar mass (g); gravitational mass, not baryonic
./const/public/const_def.f90:134: real(dp), parameter :: m_earth = mu_earth/standard_cgrav! earth mass (g)
./const/public/const_def.f90:138: real(dp), parameter :: m_jupiter = mu_jupiter/standard_cgrav ! jupiter mass (g)
./astero/private/gyre_support.f90:73: call gyre_set_constant('G_GRAVITY', standard_cgrav)
./astero/private/adipls_support.f90:234: common/ccgrav/ cgrav
./astero/private/adipls_support.f90:235: real(dp) :: cgrav
./astero/private/adipls_support.f90:290: cgrav, nn, iconst, ivar, ivers, global_data, point_data, data, aa, nn, ierr)
./astero/private/adipls_support.f90:666: subroutine read_and_store(iriche, iturpr, cgrav)
./astero/private/adipls_support.f90:668: real(dp), intent(in) :: cgrav
./astero/private/adipls_support.f90:685: cgrav, nn, iconst, ivar, ivers, glob, var, data, aa, nn, ierr)
./astero/private/adipls_support.f90:818: cgrav, nn_in, iconst, ivar, ivers, glob, var, data, aa, nn, ierr)
./astero/private/adipls_support.f90:820: real(dp), intent(in) :: cgrav
./astero/private/adipls_support.f90:885: aa(3,n)=cgrav*glob(1)*q(n)*var(5,n)/(var(10,n)*var(4,n)*var(1,n))
./astero/private/adipls_support.f90:909: data(5)=four_thirds_pi*cgrav*(var(5,1)*glob(2))**2/(var(4,1)*var(10,1))
./astero/private/adipls_support.f90:1033: 'dxdt_nuc_h1', 'z', 'dr_to_surf', 'eps_grav', 'xhe3', &
./astero/private/extras_support.f90:214: logg = log10(s% grav(1))
./astero/private/astero_support.f90:107: G = standard_cgrav
./astero/private/astero_support.f90:108: M = s% m_grav(1)
./astero/test_suite/astero_adipls/src/run_star_extras.f90:246: G = s% cgrav(1)
./binary/private/binary_mdot.f90:465: mdot_edd = pi4*standard_cgrav*b% m(b% a_i) &
./binary/private/binary_mdot.f90:468: mdot_edd = pi4*standard_cgrav*b% m(b% a_i)&
./binary/private/binary_mdot.f90:683: real(dp) :: rho_exponent, F1, q, rho, p, grav, hp, v_th, rl3, q_temp
./binary/private/binary_mdot.f90:691: grav = standard_cgrav*b% m(b% d_i)/pow2(b% r(b% d_i)) ! local gravitational acceleration
./binary/private/binary_mdot.f90:692: hp = p/(grav*rho) ! pressure scale height
./binary/private/binary_mdot.f90:702: rl3/(standard_cgrav*b% m(b% d_i)) * rho * F1
./binary/private/binary_mdot.f90:728: real(dp) :: F1, F3, G1, dP, q, rho, p, grav, hp, v_th, rl3, q_temp
./binary/private/binary_mdot.f90:756: mdot_thick = -2.0D0*pi*F1*rl_d*rl_d*rl_d/(standard_cgrav*b% m(b% d_i))*mdot_thick
./binary/private/binary_mdot.f90:790: real(dp) :: q, rho, p, grav, hp, v_th
./binary/private/binary_mdot.f90:800: grav = standard_cgrav*b% m(b% d_i)/pow2(b% r(b% d_i)) ! local gravitational acceleration
./binary/private/binary_mdot.f90:801: hp = p/(grav*rho) ! pressure scale height
./binary/private/binary_mdot.f90:805: G = standard_cgrav
./binary/private/binary_mdot.f90:834: real(dp) :: rho_exponent, F1, q, q_temp, rho, p, grav, hp, v_th, dm
./binary/private/binary_mdot.f90:843: grav = standard_cgrav*b% m(b% d_i)/pow2(b% r(b% d_i)) ! local gravitational acceleration
./binary/private/binary_mdot.f90:844: hp = p/(grav*rho) ! pressure scale height
./binary/private/binary_mdot.f90:859: (standard_cgrav*b% m(b% d_i)) * rho * F1
./binary/private/binary_mdot.f90:972: sqrt(standard_cgrav * b% m(b% a_i) * b% r(b% a_i))
./binary/private/binary_mdot.f90:976: sqrt(standard_cgrav * b% m(b% a_i) * 1.7d0*min_r)
./binary/private/binary_mdot.f90:979: sqrt(standard_cgrav * b% m(b% a_i) * b% r(b% a_i))
./binary/private/binary_ce.f90:169: b% CE_Omega_out(1) = - standard_cgrav*s% m(1)*s% dm_bar(1)/s% r(1)
./binary/private/binary_ce.f90:172: b% CE_Omega_out(4*k-3) = b% CE_Omega_out(4*(k-1)-3) - standard_cgrav*s% m(k)*s% dm_bar(k)/s% r(k)
./binary/private/binary_ce.f90:175: b% CE_Omega_in(4*s% nz-3) = - standard_cgrav*s% m(s% nz)*s% dm_bar(s% nz)/s% r(s% nz)
./binary/private/binary_ce.f90:178: b% CE_Omega_in(4*k-3) = b% CE_Omega_in(4*(k+1)-3) - standard_cgrav*s% m(k)*s% dm_bar(k)/s% r(k)
./binary/private/binary_ce.f90:272: Ecore = Ecore + s% energy(k)*s% dm(k) - standard_cgrav*s% m(k)*s% dm_bar(k)/s% r(k)
./binary/private/binary_ce.f90:282: lambda = -(standard_cgrav*b% CE_initial_Mdonor*(b% CE_initial_Mdonor - s% m(1))) &
./binary/private/binary_ce.f90:286: Ebind = -(standard_cgrav*b% CE_initial_Mdonor*(b% CE_initial_Mdonor - s% m(1))) &
./binary/private/binary_ce.f90:298: initial_Eorb = -standard_cgrav*b% CE_initial_Mdonor*b% CE_initial_Maccretor/(2*b% CE_initial_separation)
./binary/private/binary_ce.f90:300: separation = -b% CE_alpha*standard_cgrav*s% m(1)*b% CE_initial_Maccretor &
./binary/private/binary_jdot.f90:44: ! calculate jdot from gravitational wave radiation
./binary/private/binary_jdot.f90:105: real(dp) :: bs4, clight5, cgrav3
./binary/private/binary_jdot.f90:114: cgrav3 = standard_cgrav*standard_cgrav*standard_cgrav
./binary/private/binary_jdot.f90:115: b% jdot_gr = -32d0 * cgrav3 * b% m(b% a_i) * b% m(b% d_i) * (b% m(b% a_i) + b% m(b% d_i)) / &
./binary/private/binary_jdot.f90:140: sqrt(standard_cgrav * (b% m(1) + b% m(2)) * b% separation)
./binary/private/binary_evolve.f90:355: (b% m(1)+b% m(2)) / standard_cgrav * 1 / (1 - pow2(b% eccentricity))
./binary/private/binary_evolve.f90:360: (standard_cgrav*(b% m(1)+b% m(2))))
./binary/private/binary_tides.f90:147: G = standard_cgrav
./binary/private/binary_utils.f90:164: pow((standard_cgrav*(b% m(1)+b% m(2)))*pow2(b% period/(2*pi)),one_third)
./binary/private/binary_utils.f90:182: (standard_cgrav*(b% m(1)+b% m(2))))
./binary/private/binary_utils.f90:196: b% angular_momentum_j = b% m(1) * b% m(2) * sqrt(standard_cgrav *&
./binary/private/binary_do_one_utils.f90:90: write(*,'(a)') " 'dot_Jgr', time derivative of Jorb due to gravitational waves (g*cm^2/s^2),"
./binary/private/binary_do_one_utils.f90:162: Eorb = -standard_cgrav * b% m(1) * b% m(2) / (2*b% separation)
./binary/private/binary_wind.f90:118: v_orb = sqrt(standard_cgrav * b% m(s_i) / b% separation) !cm/s
./binary/private/binary_wind.f90:121: v_wind = sqrt( 2d0 / 8d0 * standard_cgrav * b% m(s_i) / b% r(s_i) )
./binary/private/binary_wind.f90:126: pow2(standard_cgrav * b% m(3-s_i) / pow2(v_wind)) *&
./binary/test_suite/jdot_gr_check/src/run_binary_extras.f90:117: beta = 64d0/5d0 * pow3(standard_cgrav) * b% m1 * b% m2 * (b% m1 + b% m2) * pow3(Msun) / pow5(clight)
-------------- next part --------------
./atm/private/atm_utils.f90
./atm/private/atm_irradiated.f90
./atm/private/atm_t_tau_varying.f90
./atm/private/atm_t_tau_uniform.f90
./atm/private/atm_table.f90
./stella/src/stl/phys_constants.f90
./stella/res/stella_extras.f90
./star/public/star_lib.f90
./star/job/run_star_support.f90
./star/private/brunt.f90
./star/private/eps_grav.f90
./star/private/conv_premix.f90
./star/private/star_profile_def.f90
./star/private/mix_info.f90
./star/private/ctrls_io.f90
./star/private/gravity_darkening.f90
./star/private/hydro_vars.f90
./star/private/star_private_def.f90
./star/private/hydro_riemann.f90
./star/private/adjust_mass.f90
./star/private/pulse_gyre.f90
./star/private/turb_support.f90
./star/private/pulse_fgong.f90
./star/private/struct_burn_mix.f90
./star/private/adjust_mesh_split_merge.f90
./star/private/star_history_def.f90
./star/private/hydro_eqns.f90
./star/private/diffusion_support.f90
./star/private/hydro_rotation.f90
./star/private/rsp_def.f90
./star/private/pulse_saio.f90
./star/private/profile_getval.f90
./star/private/create_initial_model.f90
./star/private/hydro_temperature.f90
./star/private/hydro_rsp2_support.f90
./star/private/remove_shells.f90
./star/private/overshoot_utils.f90
./star/private/eps_mdot.f90
./star/private/pulse_utils.f90
./star/private/hydro_momentum.f90
./star/private/do_one_utils.f90
./star/private/evolve.f90
./star/private/pulse_osc.f90
./star/private/rotation_mix_info.f90
./star/private/relax.f90
./star/private/init.f90
./star/private/pre_ms_model.f90
./star/private/adjust_xyz.f90
./star/private/atm_support.f90
./star/private/solver_support.f90
./star/private/hydro_rsp2.f90
./star/private/report.f90
./star/private/winds.f90
./star/private/read_model.f90
./star/private/hydro_energy.f90
./star/private/mesh_adjust.f90
./star/private/element_diffusion.f90
./star/private/rsp_step.f90
./star/private/star_utils.f90
./star/private/pulse_cafein.f90
./star/private/rsp.f90
./star/private/history.f90
./star/private/alloc.f90
./star/private/rsp_eval_eos_and_kap.f90
./star/other/other_mlt_results.f90
./star/other/other_momentum_implicit.f90
./star/other/other_cgrav.f90
./star/other/other_eps_grav.f90
./star/other/other_momentum.f90
./star/test_suite/accreted_material_j/src/run_star_extras.f90
./star/test_suite/magnetic_braking/src/run_star_extras.f90
./star/test_suite/carbon_kh/src/run_star_extras.f90
./star/test_suite/gyre_in_mesa_ms/src/run_star_extras.f90
./star/test_suite/12M_pre_ms_to_core_collapse/src/run_star_extras.f90
./star/test_suite/ns_c/src/run_star_extras.f90
./star/test_suite/ppisn/src/run_star_extras.f90
./star/test_suite/gyre_in_mesa_spb/src/run_star_extras.f90
./star/test_suite/ns_he/src/run_star_extras.f90
./star/test_suite/ns_h/src/run_star_extras.f90
./star/test_suite/extended_convective_penetration/src/run_star_extras.f90
./star/test_suite/zams_to_cc_80/src/run_star_extras.f90
./star/test_suite/gyre_in_mesa_wd/src/run_star_extras.f90
./star/test_suite/ccsn_IIp/src/run_star_extras.f90
./star/test_suite/make_env/src/run_star_extras.f90
./star/test_suite/custom_colors/src/run_star_extras.f90
./star/test_suite/gyre_in_mesa_rsg/src/run_star_extras.f90
./star/test_suite/rsp_gyre/src/run_star_extras.f90
./star/test_suite/gyre_in_mesa_envelope/src/run_star_extras.f90
./star/test_suite/simplex_solar_calibration/src/simplex_search_run_support.f90
./star/test_suite/gyre_in_mesa_bcep/src/run_star_extras.f90
./turb/test/src/test_turb.f90
./turb/public/turb.f90
./turb/private/semiconvection.f90
./turb/private/mlt.f90
./colors/public/colors_lib.f90
./const/public/const_def.f90
./astero/private/gyre_support.f90
./astero/private/adipls_support.f90
./astero/private/extras_support.f90
./astero/private/astero_support.f90
./astero/test_suite/astero_adipls/src/run_star_extras.f90
./binary/private/binary_mdot.f90
./binary/private/binary_ce.f90
./binary/private/binary_jdot.f90
./binary/private/binary_evolve.f90
./binary/private/binary_tides.f90
./binary/private/binary_utils.f90
./binary/private/binary_do_one_utils.f90
./binary/private/binary_wind.f90
./binary/test_suite/jdot_gr_check/src/run_binary_extras.f90
More information about the Mesa-users
mailing list