diff --git a/docs/source/changelog.rst b/docs/source/changelog.rst index 34fc942a7..f3cf8c3cf 100644 --- a/docs/source/changelog.rst +++ b/docs/source/changelog.rst @@ -28,6 +28,9 @@ For a more in-depth look at the new build system, see :doc:`developing/build-sys New Features ------------ +Diffusive overshooting (overmixing) prescriptions now support a +``step+exponential`` option, see :ref:`reference/controls:overshoot_scheme`. + MESA no longer stops when reactions for which special rates are set are not in the nuclear network, only a warning is printed. This is intended to make it easier to test various network sizes without having to also change the list of special reactions. .. _Bug Fixes main: @@ -35,6 +38,10 @@ MESA no longer stops when reactions for which special rates are set are not in t Bug Fixes --------- +Fixed a small bug where diffusive overshooting (overmixing) routines did not +respect changes to the mixing length set by ``other_alpha_mlt``, and used the +``mixing_length_alpha`` instead. See `gh-1003 `_. + The parameter ``report_max_infall_inside_fe_core`` was ignored in versions r25.12.1 and r26.4.1 and always had it's default value. See `gh-981 https://github.com/MESAHub/mesa/pull/981`_. ``fe_core_infall_limit`` now obeys ``when_to_stop_rtol`` and ``when_to_stop_atol`` again (broken since r11532). diff --git a/docs/source/known_bugs.rst b/docs/source/known_bugs.rst index 64abefa1b..a285bb1e3 100644 --- a/docs/source/known_bugs.rst +++ b/docs/source/known_bugs.rst @@ -13,6 +13,17 @@ issue, but it may not be complete. r26.4.1 ======= +.. _overshoot_other_alpha_mlt_bug: + +Diffusive overshooting: ``other_alpha_mlt`` ignored +--------------------------------------------------- + +Diffusive overshooting (overmixing) routines did not respect changes to the +mixing length set by ``other_alpha_mlt``, and used the ``mixing_length_alpha`` +instead. + +This is fixed in `gh-1003 `_. + .. _report_max_infall_inside_fe_core_bug: Controls: ``report_max_infall_inside_fe_core`` is ignored diff --git a/star/defaults/controls.defaults b/star/defaults/controls.defaults index 357fe8631..3f15424e3 100644 --- a/star/defaults/controls.defaults +++ b/star/defaults/controls.defaults @@ -2681,8 +2681,8 @@ ! overshooting ! ____________ - ! There are two schemes implemented in MESA to treat overshooting: - ! a step overshoot scheme and an exponential scheme. + ! There are three schemes implemented in MESA to treat overshooting: + ! a step overshoot scheme, an exponential scheme, and a step+exponential scheme. ! Parameters for exponential diffusive overshoot are described in the paper by Falk Herwig, ! "The evolution of AGB stars with convective overshoot", A&A, 360, 952-968 (2000). @@ -2699,7 +2699,7 @@ ! :: - overshoot_scheme(:) = '' ! ``exponential``, ``step``, ``other`` + overshoot_scheme(:) = '' ! ``exponential``, ``step``, ``step+exponential``, ``other`` overshoot_zone_type(:) = '' ! ``burn_H``, ``burn_He``, ``burn_Z``, ``nonburn``, ``any`` overshoot_zone_loc(:) = '' ! ``core``, ``shell``, ``any`` overshoot_bdy_loc(:) = '' ! ``bottom``, ``top``, ``any`` @@ -2711,11 +2711,12 @@ ! :: overshoot_f(:) = 0d0 + overshoot_f2(:) = 0d0 overshoot_f0(:) = 0d0 ! The switch from convective mixing to overshooting happens at a distance f0*Hp into the convection zone ! from the estimated location where ``grad_ad = grad_rad``, where Hp is the pressure scale height at that location. - ! A value <= 0 for f0 is a mistake -- you are required to set f0 as well as f. + ! A value <= 0 for f0 is a mistake -- you are required to set f0 as well as f (and f2 if using step+exponential). ! take a look at the following from an email concerning this: ! Overshooting works by taking the diffusion mixing coefficient at the edge ! of the convection zone and extending it beyond the zone. But -- and here's the issue -- @@ -2736,6 +2737,25 @@ ! These are arrays of size ``NUM_OVERSHOOT_PARAM_SETS`` which is defined in ! ``star_data/public/star_data_def.inc`` (currently 16) + ! If you use the 'step+exponential' scheme, the step is placed first and the + ! exponential component follows it. In this scheme, ``overshoot_f`` sets + ! the width of the step part and ``overshoot_f2`` sets the length scale of + ! the additional exponential tail. The parameterizations for both parts are + ! still evaluated at the convective boundary (as described above); the + ! exponential part is displaced by the extent of the step part. + ! This approach is motivated by: + ! Michielsen et al. 2021 (`A&A 650, A175 `_) + ! Anders et al. 2022a (`ApJ 926, 169 `_) + ! Anders et al. 2022b (`ApJL 928, L10 `_) + ! Johnston et al. 2024 (`ApJ 964, 170 `_) + + ! The step overshoot part models convective entrainment, where a convective + ! zone "eats up" a Ledoux-stable region (typically above a convective core). + ! The exponential overshoot part then models actual overshooting where inertia + ! carries material beyond the Schwarzschild boundary. + ! Note however that no temperature gradients are modified, MESA's implementation of overshoot + ! only extends mixing regions (through the diffusion coefficients, ``D_mix``) + ! :: overshoot_D0(:) = 0d0 diff --git a/star/private/ctrls_io.f90 b/star/private/ctrls_io.f90 index 3973e3970..7b6cbff92 100644 --- a/star/private/ctrls_io.f90 +++ b/star/private/ctrls_io.f90 @@ -147,7 +147,7 @@ module ctrls_io predictive_bdy_q_min, predictive_bdy_q_max, T_mix_limit, RSP_report_undercorrections, & do_conv_premix, conv_premix_avoid_increase, conv_premix_time_factor, & conv_premix_fix_pgas, conv_premix_dump_snapshots, do_premix_heating, & - overshoot_f, overshoot_f0, overshoot_D0, RSP_Qvisc_linear, dq_D_mix_zero_at_H_He_crossover, & + overshoot_f, overshoot_f2, overshoot_f0, overshoot_D0, RSP_Qvisc_linear, dq_D_mix_zero_at_H_He_crossover, & overshoot_Delta0, overshoot_mass_full_on, overshoot_mass_full_off, dq_D_mix_zero_at_H_C_crossover, & overshoot_scheme, overshoot_zone_type, overshoot_zone_loc, RSP_Qvisc_quadratic, & overshoot_bdy_loc, overshoot_D_min, overshoot_brunt_B_max, mlt_gradT_fraction, max_conv_vel_div_csound, & @@ -1125,6 +1125,7 @@ subroutine store_controls(s, ierr) s% do_premix_heating = do_premix_heating s% overshoot_f = overshoot_f + s% overshoot_f2 = overshoot_f2 s% overshoot_f0 = overshoot_f0 s% overshoot_D0 = overshoot_D0 s% overshoot_Delta0 = overshoot_Delta0 @@ -2848,6 +2849,7 @@ subroutine set_controls_for_writing(s, ierr) do_premix_heating = s% do_premix_heating overshoot_f = s% overshoot_f + overshoot_f2 = s% overshoot_f2 overshoot_f0 = s% overshoot_f0 overshoot_D0 = s% overshoot_D0 overshoot_Delta0 = s% overshoot_Delta0 diff --git a/star/private/overshoot.f90 b/star/private/overshoot.f90 index 1643ca4b6..6b9e01118 100644 --- a/star/private/overshoot.f90 +++ b/star/private/overshoot.f90 @@ -172,7 +172,7 @@ subroutine add_overshooting (s, ierr) select case (s%overshoot_scheme(j)) case ('exponential') call eval_overshoot_exp(s, i, j, k_a, k_b, D, vc, ierr) - case ('step') + case ('step', 'step+exponential') call eval_overshoot_step(s, i, j, k_a, k_b, D, vc, ierr) case ('other') call s% other_overshooting_scheme(s% id, i, j, k_a, k_b, D, vc, ierr) diff --git a/star/private/overshoot_step.f90 b/star/private/overshoot_step.f90 index d003b36ef..563b1bd9a 100644 --- a/star/private/overshoot_step.f90 +++ b/star/private/overshoot_step.f90 @@ -40,7 +40,7 @@ subroutine eval_overshoot_step (s, i, j, k_a, k_b, D, vc, ierr) real(dp), intent(out) :: vc(:) integer, intent(out) :: ierr - real(dp) :: f + real(dp) :: f, f2 real(dp) :: f0 real(dp) :: D0 real(dp) :: Delta0 @@ -60,19 +60,28 @@ subroutine eval_overshoot_step (s, i, j, k_a, k_b, D, vc, ierr) ! Evaluate the overshoot diffusion coefficient D(k_a:k_b) and ! mixing velocity vc(k_a:k_b) at the i'th convective boundary, ! using the j'th set of overshoot parameters. The overshoot - ! follows a simple step scheme + ! follows a simple step scheme, with an optional exponential tail. ierr = 0 ! Extract parameters f = s%overshoot_f(j) + f2 = 0._dp f0 = s%overshoot_f0(j) D0 = s%overshoot_D0(j) Delta0 = s%overshoot_Delta0(j) - if (f <= 0._dp .OR. f0 <= 0._dp) then + if (s%overshoot_scheme(j) == 'step+exponential') then + f2 = s%overshoot_f2(j) + if (f <= 0._dp .OR. f0 <= 0._dp .OR. f2 <= 0._dp) then + write(*,*) 'ERROR: for step+exp overshooting, must set f, f2 and f0 > 0' + write(*,*) 'see description of overshooting in star/defaults/control.defaults' + ierr = -1 + return + end if + else if (f <= 0._dp .OR. f0 <= 0._dp) then write(*,*) 'ERROR: for step overshooting, must set f and f0 > 0' write(*,*) 'see description of overshooting in star/defaults/control.defaults' ierr = -1 @@ -120,8 +129,6 @@ subroutine eval_overshoot_step (s, i, j, k_a, k_b, D, vc, ierr) face_loop : do k = k_a, k_b, dk - ! Evaluate the step factor - r = s%r(k) if (outward) then @@ -132,6 +139,8 @@ subroutine eval_overshoot_step (s, i, j, k_a, k_b, D, vc, ierr) if (dr < f*Hp_cb) then factor = 1._dp + else if (f2 > 0._dp) then + factor = exp(-2._dp*(dr - f*Hp_cb)/(f2*Hp_cb)) else factor = 0._dp end if @@ -144,6 +153,7 @@ subroutine eval_overshoot_step (s, i, j, k_a, k_b, D, vc, ierr) else vc(k) = 0d0 end if + ! Check for early overshoot completion if (D(k) < s%overshoot_D_min) then diff --git a/star/private/overshoot_utils.f90 b/star/private/overshoot_utils.f90 index a431a11f1..663722b26 100644 --- a/star/private/overshoot_utils.f90 +++ b/star/private/overshoot_utils.f90 @@ -216,7 +216,7 @@ subroutine eval_conv_bdy_Hp (s, i, Hp, ierr) if (s%overshoot_alpha > 0d0) then if (s%overshoot_alpha*Hp > dr) Hp = dr/s%overshoot_alpha else - if (s%alpha_mlt(k)*Hp > dr) Hp = dr/s%mixing_length_alpha + if (s%alpha_mlt(k)*Hp > dr) Hp = dr/s%alpha_mlt(k) end if end if diff --git a/star/test_suite/20M_pre_ms_to_core_collapse/inlist_common b/star/test_suite/20M_pre_ms_to_core_collapse/inlist_common index 7bc607728..86bc7b973 100644 --- a/star/test_suite/20M_pre_ms_to_core_collapse/inlist_common +++ b/star/test_suite/20M_pre_ms_to_core_collapse/inlist_common @@ -12,8 +12,7 @@ !change_net_net = .true. !new_net_name = 'mesa_206.net' !dump_missing_metals_into_heaviest = .false. - - + ! if using rotation change_D_omega_flag = .true. new_D_omega_flag = .true. @@ -41,7 +40,8 @@ show_timestep_limit_counts_when_terminate = .true. ! pgstar is recommended for diagnosing issues and understanding the evolution - pgstar_flag = .false. + pgstar_flag = .true. + disable_pgstar_during_relax_flag = .false. save_pgstar_files_when_terminate = .true. / ! end of star_job namelist @@ -144,15 +144,16 @@ ! we use step overshooting in H core - overshoot_scheme(1) = 'step' + overshoot_scheme(1) = 'step+exponential' overshoot_zone_type(1) = 'burn_H' overshoot_zone_loc(1) = 'core' overshoot_bdy_loc(1) = 'top' - overshoot_f(1) = 0.345 ! for M>10 - !overshoot_f(1) = 0.21 ! For M<10 + overshoot_f(1) = 0.335 ! for M>10 ! step part, slightly less than brott, + overshoot_f2(1) = 0.01 ! for M>10 ! exponential part, replace outer part of brott with this + !overshoot_f(1) = 0.21 ! For M<10 overshoot_f0(1) = 0.01 - ! exponential in the H core + ! exponential in the He core overshoot_scheme(2) = 'exponential' overshoot_zone_type(2) = 'burn_He' overshoot_zone_loc(2) = 'core' diff --git a/star/test_suite/20M_pre_ms_to_core_collapse/inlist_make_late_pre_zams b/star/test_suite/20M_pre_ms_to_core_collapse/inlist_make_late_pre_zams index 89ed7ec43..f8ed6c7c9 100644 --- a/star/test_suite/20M_pre_ms_to_core_collapse/inlist_make_late_pre_zams +++ b/star/test_suite/20M_pre_ms_to_core_collapse/inlist_make_late_pre_zams @@ -11,8 +11,6 @@ save_model_filename = 'late_pre_zams.mod' required_termination_code_string = 'max_age' - !pgstar_flag = .true. - / ! end of star_job namelist &eos diff --git a/star/test_suite/20M_pre_ms_to_core_collapse/inlist_pgstar b/star/test_suite/20M_pre_ms_to_core_collapse/inlist_pgstar index fc95803e7..afb609442 100644 --- a/star/test_suite/20M_pre_ms_to_core_collapse/inlist_pgstar +++ b/star/test_suite/20M_pre_ms_to_core_collapse/inlist_pgstar @@ -17,7 +17,9 @@ pgstar_show_age = .false. Grid1_win_flag = .true. Grid1_win_width = 12 +Grid1_win_width = 17 Grid1_win_aspect_ratio = 0.666 +Grid1_win_aspect_ratio = 0.53 ! file output Grid1_file_flag = .true. diff --git a/star/test_suite/20M_pre_ms_to_core_collapse/inlist_remove_envelope b/star/test_suite/20M_pre_ms_to_core_collapse/inlist_remove_envelope index 6d787b9fe..dfdebb794 100644 --- a/star/test_suite/20M_pre_ms_to_core_collapse/inlist_remove_envelope +++ b/star/test_suite/20M_pre_ms_to_core_collapse/inlist_remove_envelope @@ -16,8 +16,6 @@ !relax_initial_mass_to_remove_H_env = .true. !lg_max_abs_mdot = -2 - !pgstar_flag = .false. - / ! end of star_job namelist &eos diff --git a/star/test_suite/20M_pre_ms_to_core_collapse/inlist_to_zams b/star/test_suite/20M_pre_ms_to_core_collapse/inlist_to_zams index d10dfafcf..b5d6d12aa 100644 --- a/star/test_suite/20M_pre_ms_to_core_collapse/inlist_to_zams +++ b/star/test_suite/20M_pre_ms_to_core_collapse/inlist_to_zams @@ -18,8 +18,6 @@ change_rotation_flag = .false. ! rotation off until near zams num_steps_to_relax_rotation = 50 - pgstar_flag = .false. - / ! end of star_job namelist diff --git a/star/test_suite/20M_pre_ms_to_core_collapse/standard_after_core_c_burn.mod b/star/test_suite/20M_pre_ms_to_core_collapse/standard_after_core_c_burn.mod index 27d843893..b7eb399f4 100644 --- a/star/test_suite/20M_pre_ms_to_core_collapse/standard_after_core_c_burn.mod +++ b/star/test_suite/20M_pre_ms_to_core_collapse/standard_after_core_c_burn.mod @@ -1,3 +1,3 @@ version https://git-lfs.github.com/spec/v1 -oid sha256:0438de9dfa9479211fd27de44844fb5d5e64a2ee279369d95b183034d3912f2f -size 785401 +oid sha256:27b3dbcdbeb520e582aa3aa0aa84505c8bfadc5392909415cc86fefe772c0d24 +size 771196 diff --git a/star/test_suite/20M_pre_ms_to_core_collapse/standard_after_core_he_burn.mod b/star/test_suite/20M_pre_ms_to_core_collapse/standard_after_core_he_burn.mod index b068bd47e..a4573cab3 100644 --- a/star/test_suite/20M_pre_ms_to_core_collapse/standard_after_core_he_burn.mod +++ b/star/test_suite/20M_pre_ms_to_core_collapse/standard_after_core_he_burn.mod @@ -1,3 +1,3 @@ version https://git-lfs.github.com/spec/v1 -oid sha256:cad819ff0177ba09b07f78a3a5788d4e04627701c186a96c18d3bfbf5925da59 -size 700871 +oid sha256:85949dee2e0eb16554fa788be125bb45510dcdfbf511dd0fb8a1dd2bc4d44f00 +size 697726 diff --git a/star/test_suite/20M_pre_ms_to_core_collapse/standard_late_pre_zams.mod b/star/test_suite/20M_pre_ms_to_core_collapse/standard_late_pre_zams.mod index e6246cfc8..16ff27646 100644 --- a/star/test_suite/20M_pre_ms_to_core_collapse/standard_late_pre_zams.mod +++ b/star/test_suite/20M_pre_ms_to_core_collapse/standard_late_pre_zams.mod @@ -1,3 +1,3 @@ version https://git-lfs.github.com/spec/v1 -oid sha256:e2da28aa6d45c618355a983ffab589b5ab5cf12b6b03a4c6de60fac8bb384631 -size 198408 +oid sha256:0048c37c1cc1f072de1566c8e512a20754e8e6a336838b4b2be5a537084ce762 +size 201118 diff --git a/star/test_suite/20M_pre_ms_to_core_collapse/standard_lgTmax.mod b/star/test_suite/20M_pre_ms_to_core_collapse/standard_lgTmax.mod index 543e1555d..e728df3c5 100644 --- a/star/test_suite/20M_pre_ms_to_core_collapse/standard_lgTmax.mod +++ b/star/test_suite/20M_pre_ms_to_core_collapse/standard_lgTmax.mod @@ -1,3 +1,3 @@ version https://git-lfs.github.com/spec/v1 -oid sha256:808de00ef2b333e36afbbf0a8dda7ebc146e78560f8faf79a8b12ce6288e4acf -size 1428461 +oid sha256:93614de8fd66541bec86d84ac636946c2354ae22b29ca40c47b9efaa35e0af33 +size 1494046 diff --git a/star/test_suite/20M_pre_ms_to_core_collapse/standard_removed_envelope.mod b/star/test_suite/20M_pre_ms_to_core_collapse/standard_removed_envelope.mod index e2e3b1d99..9c2fe229d 100644 --- a/star/test_suite/20M_pre_ms_to_core_collapse/standard_removed_envelope.mod +++ b/star/test_suite/20M_pre_ms_to_core_collapse/standard_removed_envelope.mod @@ -1,3 +1,3 @@ version https://git-lfs.github.com/spec/v1 -oid sha256:636aa0a7b04e6809ddd105569d27d629451f4a3cefb788b56e9908ccc0aece3c -size 700871 +oid sha256:8d26e81be912518fcb8bd8f7f5e4e27957123e4cec1eee96b06bc027851a6999 +size 697726 diff --git a/star/test_suite/20M_pre_ms_to_core_collapse/standard_zams.mod b/star/test_suite/20M_pre_ms_to_core_collapse/standard_zams.mod index 9a4ba47f3..c0d7a472f 100644 --- a/star/test_suite/20M_pre_ms_to_core_collapse/standard_zams.mod +++ b/star/test_suite/20M_pre_ms_to_core_collapse/standard_zams.mod @@ -1,3 +1,3 @@ version https://git-lfs.github.com/spec/v1 -oid sha256:133a5ee96544e548f4094310fe94ef2dbce7246c8bdbe51f9a6a7fb16ce119b0 -size 309031 +oid sha256:d3b979b176822035c37f0b8f93e2eecb8e38eba6310945da57d845ec27cc3149 +size 312206 diff --git a/star/test_suite/make_pre_ccsn_13bvn/inlist_massive_defaults b/star/test_suite/make_pre_ccsn_13bvn/inlist_massive_defaults index 3e7c184f9..03b2a27fe 100644 --- a/star/test_suite/make_pre_ccsn_13bvn/inlist_massive_defaults +++ b/star/test_suite/make_pre_ccsn_13bvn/inlist_massive_defaults @@ -39,6 +39,9 @@ &controls +! For test_suite testing +max_model_number = 6000 + ! Initial model controls initial_mass = 12.0 initial_z = 0.02 diff --git a/star/test_suite/make_pre_ccsn_13bvn/inlist_pgstar b/star/test_suite/make_pre_ccsn_13bvn/inlist_pgstar index 84b83cdcc..8bbcf1e46 100644 --- a/star/test_suite/make_pre_ccsn_13bvn/inlist_pgstar +++ b/star/test_suite/make_pre_ccsn_13bvn/inlist_pgstar @@ -1,6 +1,6 @@ &pgstar -pgstar_interval = 1 +pgstar_interval = 20 file_digits = 8 pgstar_show_age_in_years = .true. pgstar_show_age_in_seconds = .false. diff --git a/star/test_suite/make_pre_ccsn_13bvn/standard_after_si_burn.mod b/star/test_suite/make_pre_ccsn_13bvn/standard_after_si_burn.mod index 6b5987cdb..8291cda3a 100644 --- a/star/test_suite/make_pre_ccsn_13bvn/standard_after_si_burn.mod +++ b/star/test_suite/make_pre_ccsn_13bvn/standard_after_si_burn.mod @@ -1,3 +1,3 @@ version https://git-lfs.github.com/spec/v1 -oid sha256:f51b9324742131fd131f08c1144c7bcc40a95abfbb1f6df28f78f24182adf9fb -size 770391 +oid sha256:96e5c97b1277ec0fec955b0af21ac26d96278a55ba441a4c462cd9122cfafd1d +size 941828 diff --git a/star/test_suite/make_pre_ccsn_13bvn/standard_near_zams.mod b/star/test_suite/make_pre_ccsn_13bvn/standard_near_zams.mod index 3d1f0d59d..61c305a04 100644 --- a/star/test_suite/make_pre_ccsn_13bvn/standard_near_zams.mod +++ b/star/test_suite/make_pre_ccsn_13bvn/standard_near_zams.mod @@ -1,3 +1,3 @@ version https://git-lfs.github.com/spec/v1 -oid sha256:44a4791821d73651b594f0853a6365357f1c91e7954e1431aab4ffab8791ae13 -size 182983 +oid sha256:6452a46d109dc3f43f98e739a8c4f83a64e2b92e82fe5d482e5ecd6bc430b2f0 +size 182578 diff --git a/star/test_suite/make_pre_ccsn_13bvn/standard_ready.mod b/star/test_suite/make_pre_ccsn_13bvn/standard_ready.mod index cfbfb2530..2dd436841 100644 --- a/star/test_suite/make_pre_ccsn_13bvn/standard_ready.mod +++ b/star/test_suite/make_pre_ccsn_13bvn/standard_ready.mod @@ -1,3 +1,3 @@ version https://git-lfs.github.com/spec/v1 -oid sha256:605ccc1bff348eacfd7bbcef0c16789f86e980821798cb1c1e72cf28863bc797 -size 418755 +oid sha256:be5d7b8802fd5770c3617dcdd8af84d2b0ebb5477fc3299811daf9dc98a02e58 +size 340181 diff --git a/star/test_suite/make_pre_ccsn_13bvn/standard_ready_to_remove.mod b/star/test_suite/make_pre_ccsn_13bvn/standard_ready_to_remove.mod index de63bce4c..3bd285010 100644 --- a/star/test_suite/make_pre_ccsn_13bvn/standard_ready_to_remove.mod +++ b/star/test_suite/make_pre_ccsn_13bvn/standard_ready_to_remove.mod @@ -1,3 +1,3 @@ version https://git-lfs.github.com/spec/v1 -oid sha256:ef13d990bec377fee67829754112f3639565b86c0239104922c4b0af4399f925 -size 531139 +oid sha256:e4aa0b5a0b5854d80f6f7c8fbf6cee8d3e5fe8c20d16d8222ade6a97a8b0c01d +size 464418 diff --git a/star/test_suite/make_pre_ccsn_13bvn/standard_removed.mod b/star/test_suite/make_pre_ccsn_13bvn/standard_removed.mod index 9f02b7819..a13e13411 100644 --- a/star/test_suite/make_pre_ccsn_13bvn/standard_removed.mod +++ b/star/test_suite/make_pre_ccsn_13bvn/standard_removed.mod @@ -1,3 +1,3 @@ version https://git-lfs.github.com/spec/v1 -oid sha256:e4c3c9493370ea00bee2763ba0f905306ff199a1bd1ad711bfc2ef61bc9368a7 -size 461777 +oid sha256:25a58d9c96acfe1818142e18611f7acc982a1e7fe09d09d34f90169385345437 +size 392422 diff --git a/star_data/private/star_controls.inc b/star_data/private/star_controls.inc index 2d0c29d8d..b4e3b4174 100644 --- a/star_data/private/star_controls.inc +++ b/star_data/private/star_controls.inc @@ -335,6 +335,7 @@ logical :: do_premix_heating ! include eps_pre_mix in energy equation real(dp) :: overshoot_f(NUM_OVERSHOOT_PARAM_SETS) + real(dp) :: overshoot_f2(NUM_OVERSHOOT_PARAM_SETS) real(dp) :: overshoot_f0(NUM_OVERSHOOT_PARAM_SETS) real(dp) :: overshoot_D0(NUM_OVERSHOOT_PARAM_SETS) real(dp) :: overshoot_Delta0(NUM_OVERSHOOT_PARAM_SETS)