more tweaks and new controls for mixing near surface

Author: Debraheem

star/defaults/controls.defaults

@@ -4135,6 +4135,27 @@
     D_mix_rotation_min_logT_full_off = 9.5d0
 
 
+      ! D_mix_rotation_min_tau_full_off
+      ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+      ! Drop rotational components of ``D_mix`` for locations where tau <= this.
+      ! For numerical stability, turn off rotational part of ``D_mix`` at very low tau.
+
+      ! ::
+      
+    D_mix_rotation_min_tau_full_off = 0d0
+    
+          ! D_mix_rotation_max_tau_full_on
+      ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+      ! Use rotational components of ``D_mix`` for locations where tau >= this.
+      ! For numerical stability, turn off rotational part of ``D_mix`` at very low tau.
+
+      ! ::
+      
+    D_mix_rotation_min_tau_full_on = 0d0
+
+
       ! D_omega_max_replacement_fraction
       ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
       ! D_omega_growth_rate

star/dev_cases_test_TDC/dev_TDC_through_ppisn/he_dep.mod

@@ -14,11 +14,12 @@
    Grid2_win_flag = .true.
    !Grid2_win_flag = .false.
 
-   Grid2_win_width = 13 !15
+   Grid2_win_width = 21 !15
    Grid2_win_aspect_ratio = 0.65 ! aspect_ratio = height/width
 
    Grid2_plot_name(4) = 'Mixing'
-   
+   Mixing_show_rotation_details = .true.
+
    Grid2_num_cols = 7 ! divide plotting region into this many equal width cols
    Grid2_num_rows = 8 ! divide plotting region into this many equal height rows
 
@@ -202,7 +203,7 @@
    Grid2_file_dir = 'png'
    Grid2_file_prefix = 'grid_'
    Grid2_file_interval = 10 ! 1 ! output when mod(model_number,Grid2_file_interval)==0
-   Grid2_file_width = -1 ! negative means use same value as for window
+   Grid2_file_width = 27 ! negative means use same value as for window
    Grid2_file_aspect_ratio = -1 ! negative means use same value as for window
 
    !show_TRho_Profile_eos_regions = .true.

star/dev_cases_test_TDC/dev_TDC_through_ppisn/inlist_pgstar

@@ -32,7 +32,7 @@
    ! this are options for the relaxation of the model in-between pulses
    timescale_for_relax_entropy = 1d-15
    max_dt_for_relax_entropy = 5d-17
-   num_timescales_for_relax_entropy = 30
+   num_timescales_for_relax_entropy = 30 ! probably need 50 or more, takes a while.
    max_steps_to_relax_entropy = 10000
    relax_omega_max_yrs_dt = 1d-5
 
@@ -64,6 +64,12 @@
 
 &controls
 
+! remove limits on rotational mixing in the core
+    D_mix_rotation_max_logT_full_on = 1d11
+    D_mix_rotation_max_logT_full_off = 1d12
+
+
+
 ! atmosphere
     !  Pextra_factor = 1 
     !  ! need this extra pressure to stabilize the atmosphere during core He burning
@@ -117,24 +123,24 @@
 
    !relax_use_gold_tolerances = .false.
    use_gold_tolerances = .true.
-   use_gold2_tolerances = .false.
-   maxT_for_gold_tolerances = 3.7d9 ! this is to ease core-collapse (~logT=9.7)
+   use_gold2_tolerances = .true.
+   maxT_for_gold_tolerances = 1d99!3.7d9 ! this is to ease core-collapse (~logT=9.7)
    gold_iter_for_resid_tol2 = 10
    gold_iter_for_resid_tol3 = 10
    gold_tol_residual_norm3 = 1d-6
    gold_tol_max_residual3 = 1d-3
    gold_solver_iters_timestep_limit = 20
    solver_iters_timestep_limit = 50
-   ignore_too_large_correction = .true.
+   ignore_too_large_correction = .false.
    scale_max_correction = 0.1d0
    !corr_coeff_limit = 0.2d0
    ignore_min_corr_coeff_for_scale_max_correction = .true.
    ignore_species_in_max_correction = .true.
 
       ! Can be helpful to decrease op_split_burn_min_T 
       ! to 3d9 or even 1d9 for large networks to avoid crashes
-      op_split_burn = .false.
-      op_split_burn_min_T = 3.0d9 
+      op_split_burn = .true.
+      op_split_burn_min_T = 1d9 
       burn_steps_limit = 150
       burn_steps_hard_limit = 250
       op_split_burn_eps = 1d-5
@@ -238,7 +244,7 @@
    ! mesh controls
    okay_to_remesh = .true.
    max_dq = 1d-3 ! RECOMMENDED 1d-3
-   mesh_delta_coeff = 0.8d0 ! RECOMMENDED 0.8d0
+   mesh_delta_coeff = 1d0 ! RECOMMENDED 0.8d0
    !! this one is turned on in run_star_extras
    !use_split_merge_amr = .true.
    split_merge_amr_log_zoning = .true.

star/dev_cases_test_TDC/dev_TDC_through_ppisn/inlist_ppisn

@@ -15,16 +15,13 @@
    save_model_filename = 'he_dep.mod'
    required_termination_code_string = 'xa_central_lower_limit'
 
-   read_extra_star_job_inlist(1) = .true.
-   extra_star_job_inlist_name(1) = 'inlist_ppisn'
+
 
 / !end of star_job namelist
 
 
 &controls
 
-   read_extra_controls_inlist(1) = .true.
-   extra_controls_inlist_name(1)= 'inlist_ppisn'
 
   ! stop when the center mass fraction of h1 drops below this limit
     xa_central_lower_limit_species(1) = 'he4'

star/dev_cases_test_TDC/dev_TDC_through_ppisn/inlist_to_he_dep

@@ -1,10 +1,13 @@
 
 &star_job
 
+      read_extra_star_job_inlist(1) = .true.
+      extra_star_job_inlist_name(1) = 'inlist_ppisn'
 
       read_extra_star_job_inlist(2) = .true.
       extra_star_job_inlist_name(2) = 'inlist_to_he_dep'
 
+
 / ! end of star_job namelist
 
 
@@ -28,7 +31,10 @@
 &controls
 
       read_extra_controls_inlist(1) = .true.
-      extra_controls_inlist_name(1)= 'inlist_to_he_dep'
+      extra_controls_inlist_name(1)= 'inlist_ppisn'
+
+      read_extra_controls_inlist(2) = .true.
+      extra_controls_inlist_name(2)= 'inlist_to_he_dep'
 
 / ! end of controls namelist
 

star/dev_cases_test_TDC/dev_TDC_through_ppisn/inlist_to_he_dep_header

@@ -10,12 +10,12 @@ source "${MESA_DIR}/star/test_suite/test_suite_helpers"
 date "+DATE: %Y-%m-%d%nTIME: %H:%M:%S"
 
 # check if can skip building starting model
-#if [ -n "$MESA_SKIP_OPTIONAL" ]; then
-#   cp standard_he_dep.mod he_dep.mod
-#else
-#   do_one inlist_to_he_dep_header he_dep.mod
-#   cp he_dep.mod standard_he_dep.mod
-#fi
+if [ -n "$MESA_SKIP_OPTIONAL" ]; then
+   cp standard_he_dep.mod he_dep.mod
+else
+   do_one inlist_to_he_dep_header he_dep.mod
+   cp he_dep.mod standard_he_dep.mod
+fi
 
 do_one inlist_pulses_header final.mod
 

star/dev_cases_test_TDC/dev_TDC_through_ppisn/rn

@@ -866,7 +866,7 @@ integer function extras_start_step(id)
          if (.not. in_inlist_pulses .and. .not. s% lxtra(lx_he_zams)) then
             lburn_div_lsurf  = abs(s% L_nuc_burn_total*Lsun/s% L(1))
             if (lburn_div_lsurf > 0d0) then
-               if((abs(log10(lburn_div_lsurf))) < 0.01 .and. &
+               if((abs(log10(lburn_div_lsurf))) < 0.1 .and. &
                   (s% star_age > 1d3 .or. s% center_he4 < 0.98d0)) then
                   s% use_other_before_struct_burn_mix = .false.
                   call star_relax_uniform_omega(id, 1, s% job% new_omega_div_omega_crit,&
@@ -1272,9 +1272,13 @@ subroutine my_before_struct_burn_mix(id, dt, res)
          if(s% u_flag) then
             call star_read_controls(id, 'inlist_hydro_on', ierr)
             if (s% xtra(x_time_start_pulse) > 0d0) then
-               s% v_drag = s% u(1)
-               s% v_drag_factor = 1d0
-               write(*,*) 'using drag'
+               !s% v_drag = 2d4 !s% u(1)
+               !s% v_drag_factor = 1d0
+               !s% v_drag = 2d9
+               !s% v_drag_factor = 1d0
+               s% D_mix_rotation_min_tau_full_off = 1d0
+               s% D_mix_rotation_min_tau_full_on = 10d0
+               write(*,*) 'u_flag is on, pulse has begun'
                if (.not. s% lxtra(lx_using_bb_bcs)) then
                   s% max_timestep = max_dt_during_pulse
                else
@@ -1289,13 +1293,17 @@ subroutine my_before_struct_burn_mix(id, dt, res)
                end do
                if (k > 1) then
                  s% max_q_for_convection_with_hydro_on = s% q(k)
+                 !if (s% tau(k)> 10d0) then
+                 !s% D_mix_rotation_min_tau_full_off = s% tau(k)
+                 !end if
                end if
             else
+               write (*,*) 'u_flag is on, Pulse has not begun.'
                s% max_timestep = max_dt_before_pulse
                !!this prevents surface from going mad right when hydro is turned on
                !right after a relax
-               s% v_drag = 0d0
-               s% v_drag_factor = 1d0
+               !s% v_drag = 0d0
+               !s% v_drag_factor = 1d0
                !s% max_timestep = 1d99
                s% max_q_for_convection_with_hydro_on = 0.999d0
             end if
@@ -1390,7 +1398,7 @@ subroutine my_before_struct_burn_mix(id, dt, res)
             end if
             if (safe_log10(abs(power_photo)) > max_Lphoto_for_lgLnuc_limit) then
                s% delta_lgL_power_photo_limit = delta_lgLnuc_limit
-               s% delta_lgL_power_photo_hard_limit = 2d0*delta_lgLnuc_limit
+               s% delta_lgL_power_photo_hard_limit = 10d0*delta_lgLnuc_limit
             else
                s% delta_lgL_power_photo_limit = -1d0
                s% delta_lgL_power_photo_hard_limit = -1d0