[Draft] [debug] debug 20Msun crash

net/private/net_approx21.f90

@@ -1683,7 +1683,17 @@ subroutine approx21_eps_info( &
                xx = xx + a1*a2
             end do
             eps_total_q = -m3(avo,clight,clight) * xx
+            ! XXX error here??? eps_total is really negative?
             eps_total = eps_total_q
+!            if (n% zone >= 1816 .and. n% zone <= 1819) then
+!               xx = 0.0_qp
+!               do i=1,species(plus_co56)
+!                  a1 = dydt(i)
+!                  a2 = mion(i)
+!                  xx = xx + a1*a2
+!                  write(*,*) n% zone, eps_total, i, a1, a2, xx
+!               end do
+!            end if
 
             fe56ec_fake_factor = eval_fe56ec_fake_factor( &
                n% g% fe56ec_fake_factor, n% g% min_T_for_fe56ec_fake_factor, n% temp)

net/private/net_eval.f90

@@ -384,6 +384,8 @@ subroutine eval_net_approx21_procs(n,just_dxdt, ierr)
 
          if (ierr /= 0) return
          n% eps_nuc = n% eps_total - n% eps_neu_total
+         ! XXX test:
+         !n% eps_nuc = n% eps_total 
 
          do i=1, num_isos
             n% dxdt(i) = chem_isos% Z_plus_N(g% chem_id(i)) * n% dydt1(i)

star/dev_cases_test_solver/test_20M_near_cc_approx21/.gitattributes

@@ -0,0 +1,2 @@
+plots/* filter=lfs diff=lfs merge=lfs -text
+photos/* filter=lfs diff=lfs merge=lfs -text

star/dev_cases_test_solver/test_20M_near_cc_approx21/README.rst

@@ -0,0 +1,45 @@
+.. _20M_pre_ms_to_core_collapse:
+
+***************************
+20M_pre_ms_to_core_collapse
+***************************
+
+This test suite evolves a solar metalicity 20 |MSun| model from the pre-ms to core collapse. 
+For bit for bit convergence, we recomended to run by using the ./run_all script instead of restarting from models,
+see https://github.com/MESAHub/mesa/issues/610.
+
+This test_suite has been tested up to 80 solar masses, up to solar metallicity, with mass loss, and produces reasonable HR-tracks.
+Note that for higher masses at solar metallicity, some combination of Pextra_factor, mass-loss, and/or superadiabatic convection reduction (e.g. mlt++)
+might be necessary to stabilize the surface and avoid numerical issues. See the 80Msun_zams_to_cc test_suite as an example.
+
+For production science we recommend adopting tighter mesh and timestep controls, such as those suggested in the comments of inlist_common.
+
+Physical checks
+===============
+
+None
+
+Inlists
+=======
+
+This test case has seven parts.
+
+* Part 1 (``inlist_make_late_pre_zams``) creates a 20 |Msun|, Z=1.42*10^-2 metallicity, pre-main sequence model and evolves it for 100 years.
+
+* Part 2 (``inlist_to_zams``) evolves the model to the zero age main sequence.
+
+* Part 3 (``inlist_to_end_core_he_burn``) takes the model to core helium depletion.
+
+* Part 4 (``inlist_remove_envelope``) removes the remianing hydrogen envelope. (optional)
+
+* Part 5 (``inlist_to_end_core_c_burn``) takes the model to core carbon depletion.
+
+* Part 6 (``inlist_to_lgTmax``) evolves the model until the core temperature reaches log T =9.60 (approximately silicon-shell burning)
+
+* Part 7 (``inlist_to_cc``) evolves until core collapse.
+
+
+
+
+Last-Updated: 18Dec2023 by EbF
+

star/dev_cases_test_solver/test_20M_near_cc_approx21/ck

@@ -0,0 +1,7 @@
+#!/bin/bash
+
+# this provides the definition of check_one
+# check_one
+source "${MESA_DIR}/star/test_suite/test_suite_helpers"
+
+check_one

star/dev_cases_test_solver/test_20M_near_cc_approx21/clean

@@ -0,0 +1,4 @@
+#!/bin/bash
+
+cd make
+make clean