Three component phase separation for ultramassive white dwarfs, part 2

star/defaults/controls.defaults

@@ -7547,6 +7547,9 @@
       !   phase diagram of `Blouin & Daligault (2021a) <https://ui.adsabs.harvard.edu/abs/2021PhRvE.103d3204B>`_.
       ! + ``'ONe'`` : oxygen-neon phase separation using the two-component
       !   phase diagram of `Blouin & Daligault (2021b) <https://ui.adsabs.harvard.edu/abs/2021ApJ...919...87B>`_.
+      ! + ``'3c'`` : three-component phase separation option dynamically selects dominant three components out
+      !   of C/O/Ne/Na/Mg and applies three-component PS procedure based on
+      !   `Castro-Tapia & Cumming (2025) <https://ui.adsabs.harvard.edu/abs/2025ApJ...991...64C>`_.
 
       ! ::
 
@@ -7563,6 +7566,17 @@
 
     do_phase_separation_heating = .true.
 
+      ! smooth_phase_separation_heating
+      ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+      ! if true, redistribute all heating associated with phase
+      ! separation evenly throughout the inner half of the white dwarf
+      ! (prevents unphysical localized spikes at small time steps)
+
+      ! ::
+
+    smooth_phase_separation_heating = .true.
+
       ! phase_separation_mixing_use_brunt
       ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 

star/private/ctrls_io.f90

@@ -323,6 +323,7 @@ module ctrls_io
     do_phase_separation, &
     phase_separation_option, &
     do_phase_separation_heating, &
+    smooth_phase_separation_heating, &
     phase_separation_mixing_use_brunt, &
     phase_separation_no_diffusion, &
 
@@ -1812,6 +1813,7 @@ subroutine store_controls(s, ierr)
  s% do_phase_separation = do_phase_separation
  s% phase_separation_option = phase_separation_option
  s% do_phase_separation_heating = do_phase_separation_heating
+ s% smooth_phase_separation_heating = smooth_phase_separation_heating
  s% phase_separation_mixing_use_brunt = phase_separation_mixing_use_brunt
  s% phase_separation_no_diffusion = phase_separation_no_diffusion
 
@@ -3472,6 +3474,7 @@ subroutine set_controls_for_writing(s, ierr)
  do_phase_separation = s% do_phase_separation
  phase_separation_option = s% phase_separation_option
  do_phase_separation_heating = s% do_phase_separation_heating
+ smooth_phase_separation_heating = s% smooth_phase_separation_heating
  phase_separation_mixing_use_brunt = s% phase_separation_mixing_use_brunt
  phase_separation_no_diffusion = s% phase_separation_no_diffusion