Pgplot residuals

star/Makefile

@@ -202,6 +202,8 @@ SRCS += \
    private/pgstar_profile_panels.f90 \
    private/pgstar_grid.f90 \
    private/pgstar_network.f90 \
+   private/pgstar_equation_residuals.f90 \
+   private/pgstar_kipp_residuals.f90 \
    other/sample_pgstar_plot.f90 \
    private/pgstar_full.f90
 EXTERNAL_DEPENDS_ON := pgplot
@@ -237,11 +239,13 @@ MODULES += pgstar_abundance.mod \
            pgstar_decorator.mod \
            pgstar_dpg_dnu.mod \
            pgstar_dynamo.mod \
+		   pgstar_equation_residuals.mod \
            pgstar_grid.mod \
            pgstar_history_panels.mod \
            pgstar_hist_track.mod \
            pgstar_hr.mod \
            pgstar_kipp.mod \
+           pgstar_kipp_residuals.mod \
            pgstar_logg_logt.mod \
            pgstar_logg_teff.mod \
            pgstar_logl_r.mod \

star/defaults/pgstar.defaults

@@ -7814,6 +7814,66 @@
     Grid9_file_width = -1
     Grid9_file_aspect_ratio = -1
 
+
+      ! Max_eq_resid
+      ! ------------
+
+      ! ::
+
+   Max_eq_resid_win_flag = .false.
+   Max_eq_resid_win_width = 7
+   Max_eq_resid_win_aspect_ratio = 5d-1
+
+   Max_eq_resid_xleft = 0.15
+   Max_eq_resid_xright = 0.85
+   Max_eq_resid_ybot = 0.15
+   Max_eq_resid_ytop = 0.85
+   Max_eq_resid_txt_scale = 1
+   Max_eq_resid_title = 'Max equation residuals'
+   Max_eq_resid_max_width = -1
+
+      ! **File output**
+
+      ! ::
+
+   Max_eq_resid_file_flag = .false.
+   Max_eq_resid_file_dir = 'png'
+   Max_eq_resid_file_prefix = 'max_eq_resid_'
+   Max_eq_resid_file_width = -1
+   Max_eq_resid_file_aspect_ratio = -1
+
+
+      ! Kipp_residuals
+      ! ------------
+
+      ! ::
+
+   Kipp_residuals_win_flag = .false.
+   Kipp_residuals_win_width = 9
+   Kipp_residuals_win_aspect_ratio = 6d-1
+
+      ! ::
+
+   Kipp_residuals_xleft = 0.15
+   Kipp_residuals_xright = 0.85
+   Kipp_residuals_ybot = 0.15
+   Kipp_residuals_ytop = 0.85
+   Kipp_residuals_txt_scale = 1
+   Kipp_residuals_title = 'Kippenhahn residuals'
+   Kipp_residuals_max_width = -1
+   Kipp_residuals_yaxis_name = 'mass'
+
+      ! **File output**
+
+      ! ::
+
+   Kipp_residuals_file_flag = .false.
+   Kipp_residuals_file_dir = 'png'
+   Kipp_residuals_file_prefix = 'kipp_residuals_'
+   Kipp_residuals_file_width = -1
+   Kipp_residuals_file_aspect_ratio = -1
+
+
 ! Annotation
 ! ==========
 

star/dev_cases_massive_stars/15M_binary_stripped_envelope/.gitignore

@@ -0,0 +1,9 @@
+*.mod
+*.data
+*.dat
+png/*
+LOGS*/*
+photos/*
+*cache*
+build/*
+output.txt

star/dev_cases_massive_stars/15M_binary_stripped_envelope/Makefile

@@ -0,0 +1,21 @@
+ifeq ($(MESA_DIR),)
+  $(error MESA_DIR environment variable is not set)
+endif
+
+include $(MESA_DIR)/make/defaults-module.mk
+
+MODULE_NAME := binary
+SRCS := \
+   src/binary_run.f90 \
+   src/run_star_extras.f90 \
+   src/run_binary_extras.f90 \
+   $(MESA_DIR)/binary/job/run_binary.f90
+
+include $(MAKE_DIR)/work.mk
+
+run: $(OBJ_OUT)
+	@date "+DATE: %Y-%m-%d%nTIME: %H:%M:%S"
+	@$(OBJ_OUT)
+	@date "+DATE: %Y-%m-%d%nTIME: %H:%M:%S"
+
+.PHONY: run

star/dev_cases_massive_stars/15M_binary_stripped_envelope/README.md

@@ -0,0 +1,47 @@
+#  MESA setup for Stripped Envelope Supernova Progenitor
+
+Developed for the [TARDIS connector workshop 2025](https://kasukabedefenceforce.github.io/tardis-con/2025) at CCA/Flatiron.
+
+The setup template is in <./template_binary>. It computes the
+evolution of a 15M<sub>☉</sub> donor star + 12M<sub>☉</sub> point-mass
+companion with an initial period of 20 days through case B (post
+donor-main-sequence) fully conservative mass transfer until carbon
+core depletion and core-collapse.
+
+Originally developed for [Renzo et al.
+2021](https://ui.adsabs.harvard.edu/abs/2021ApJ...923..277R/abstract),
+[2023](https://ui.adsabs.harvard.edu/abs/2023ApJ...942L..32R/abstract),
+ported to MESA `r24.03.1` or later by
+[@neevshah101](https://github.com/neevshah101). The post-carbon
+depletion evolution (see below) is based on a template inlist provided
+by [@Debraheem](https://github.com/Debraheem). These have been
+iterated upon by trial-and-error until viable models run.
+
+Output files generated by using this setup are saved in the
+run<sub>outputs</sub> dir for each directory.
+
+
+# Note on evolution beyond C depletion
+
+At carbon depletion ($X_{center}(^{12}C)_{} <0.01$, defined in
+`extras_finish_step` in `run_star_extras.f90`) we read on-the-fly a
+`star_job` and `controls` from an inlist whose path is **hard coded**.
+This inlist lives in
+[inlist<sub>to</sub><sub>cc</sub>](./template_binary/inlist_to_cc).
+and it helps to reach the onset of core collapse.
+
+As usual, most computing time is spent post-carbon depletion. This
+model uses `use_drag = .true.` post-carbon depletion to successfully go
+to core-collapse without spurious envelope velocities. It also uses
+`op_split = .true.` for T&ge;4&times;10<sup>9</sup> K, and changes overshooting parameters
+for late burning (this has been found to be necessary through a
+non-systematic but extensive experimentation).
+
+
+# TODO
+
+-   test it reaches onset of core-collapse with MESA `r24.08.1`
+    -   [X] 22 iso network
+    -   [ ] large network (work in progress in separate folder/branches)
+-   [ ] check CPUh cost: &sim; 3h on 12 cores
+

star/dev_cases_massive_stars/15M_binary_stripped_envelope/clean

@@ -0,0 +1,3 @@
+#!/usr/bin/env bash
+
+make clean