[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