Merge pull request #606 from MESAHub/fix_const_density_burn_solver

partial-Fix const density burn solver
MESAHub · Nov 22, 2023 · a563abf · a563abf
2 parents d8f68a1 + bf2e473
commit a563abf
Show file tree

Hide file tree

Showing 3 changed files with 15 additions and 4 deletions.
diff --git a/net/private/net_burn_const_density.f90 b/net/private/net_burn_const_density.f90
@@ -237,7 +237,8 @@ subroutine burn_const_density_1_zone( &
          stpmax = max_steps
 
          n% screening_mode = screening_mode
-
+         n% g => g
+
          if (dbg) write(*,2) 'call setup_net_info', iwork
          call setup_net_info(n) 
          if (dbg) write(*,*) 'done setup_net_info'

diff --git a/star/defaults/controls.defaults b/star/defaults/controls.defaults
@@ -9084,8 +9084,6 @@
 
       ! op_split_burn_min_T
       ! ~~~~~~~~~~~~~~~~~~~
-      ! op_split_burn_min_T_for_variable_T_solver
-      ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
       ! op_split_burn_eps
       ! ~~~~~~~~~~~~~~~~~
       ! op_split_burn_odescal
@@ -9096,7 +9094,6 @@
       ! ::
 
     op_split_burn_min_T = 2d9
-    op_split_burn_min_T_for_variable_T_solver = 1d99
     op_split_burn_eps = 1d-5
     op_split_burn_odescal = 1d-5
 

diff --git a/star/defaults/controls_dev.defaults b/star/defaults/controls_dev.defaults
@@ -161,3 +161,16 @@
     phase_separation_no_diffusion = .false.
 
 
+! split burn
+! ================
+
+
+      ! op_split_burn_min_T_for_variable_T_solver
+      ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+      ! use variable T, constant density solver for op_split_burn in cells with T >= this limit at start of step.
+      ! only adopted in cells with T >= op_split_burn_min_T at start of step.
+
+      ! ::  
+
+    op_split_burn_min_T_for_variable_T_solver = 1d99