add dev test cases for radial pulsations

Author: Debraheem

star/dev_cases_test_RSP2/dev_rsp2_Cepheid_6M/inlist_rsp2_Cepheid

@@ -51,7 +51,14 @@
    kap_CO_prefix = 'a09_co'
 / ! end of kap namelist
 
-&controls  
+&controls
+   num_trace_history_values = 2
+   trace_history_value_name(1) = 'log_rel_run_E_err'
+   trace_history_value_name(2) = 'rel_E_err'
+   warn_when_large_rel_run_E_err = 1d-2
+
+!convergence_ignore_equL_residuals = .true.
+
    !save_controls_namelist = .true.
    !controls_namelist_name = 'rsp2_controls.txt'
 
@@ -67,8 +74,16 @@
 
 ! surf BCs
    use_zero_Pgas_outer_BC = .true.
-   RSP2_use_L_eqn_at_surface = .true.
-
+   !RSP2_use_L_eqn_at_surface = .true.
+   !use_compression_outer_BC = .true.
+
+RSP2_Lsurf_factor = 0.5d0
+ RSP2_use_L_eqn_at_surface = .true.
+ RSP2_assume_HSE = .true.
+ RSP2_use_RSP_eqn_for_Y_face = .true.
+ RSP2_use_mass_interp_face_values = .false.
+ RSP2_use_Stellingwerf_Lr = .true.
+ 
 ! timesteps
 !  RSP_max_dt_times_min_dr_div_cs = 2d0
    dt_div_min_dr_div_cs_limit = 2d0
@@ -123,29 +138,33 @@
 
    energy_eqn_option = 'dedt'
 
-   RSP2_assume_HSE = .true.
    include_w_in_correction_limits = .true.
    include_v_in_correction_limits = .true.
    include_L_in_correction_limits = .true.
 
+   RSP2_w_fix_if_neg = 5d-5
+RSP2_w_min_for_damping = 1d2
+
    RSP2_alfap =     0.0d0         ! turbulent pressure; Pt ~ alfap     
    RSP2_alfat =     0.0d0         ! turbulent flux; Lt ~ RSP2_alfat
    RSP2_alfam =     0.25d0        ! turbulent viscosity; Eq and Uq ~ RSP2_alfam
    RSP2_alfar =     0.0d0         ! radiative dissipation of turbulent energy; Dr ~ RSP2_alfar   
    RSP2_alfad =     1d0
 
-   RSP2_num_outermost_cells_forced_nonturbulent = 1
-   RSP2_num_innermost_cells_forced_nonturbulent = 5
+  ! RSP2_num_outermost_cells_forced_nonturbulent = 1
+  ! RSP2_num_innermost_cells_forced_nonturbulent = 5
    alt_scale_height_flag = .false.
 
-   tol_max_correction = 5d-7
+   !tol_max_correction = 5d-7
 
 ! output
-   x_integer_ctrl(1) = -1 ! which period to check
+   x_integer_ctrl(1) = 0 ! which period to check
    x_ctrl(1) = 20 ! expected period (in days) 
 
+ x_logical_ctrl(5) = .false.
+
    x_logical_ctrl(7) = .true. ! doing pulses
-   x_integer_ctrl(7) = 0 ! 3 ! which period to check (<= 0 means don't check any)
+   x_integer_ctrl(7) = -1 ! 3 ! which period to check (<= 0 means don't check any)
    x_ctrl(7) = 20 ! expected period (in days) 
 !   x_ctrl(8) = 1.0 ! min_deltaR_for_periods (Rsun)
    x_ctrl(9) = 0.1 ! KE_growth_avg_abs_frac_new
@@ -163,7 +182,7 @@
       photo_interval = 2000
       profile_interval = 2000
       history_interval = 40
-      terminal_interval = 2000
+      terminal_interval = 100 !2000
 
       !photo_interval = 100
       !profile_interval = 400
@@ -182,7 +201,7 @@
 
       !max_number_retries = 0
 
-      !report_solver_progress = .true. ! set true to see info about solver iterations
+      report_solver_progress = .true. ! set true to see info about solver iterations
       !report_ierr = .true. ! if true, produce terminal output when have some internal error
       !stop_for_bad_nums = .true.
 
@@ -272,7 +291,7 @@
    History_Panels2_xaxis_name='day'
    History_Panels2_max_width = 40 ! only used if > 0.  causes xmin to move with xmax.
 
-   !Grid2_file_flag = .true.
+   Grid2_file_flag = .true.
    file_digits = 7
    Grid2_file_dir = 'png'
    Grid2_file_prefix = 'grid2'

star/dev_cases_test_TDC/dev_TDC_Cepheid_6M/README.rst

@@ -0,0 +1,21 @@
+.. _ppisn:
+
+*****
+ppisn
+*****
+
+This test case evolves a very massive helium star from the He-ZAMS
+up to the ocurrence of a pulsational pair-instability event (see |Marchant2019|).
+
+.. |Marchant2019| replace:: `Marchant et al. 2019 <https://ui.adsabs.harvard.edu/abs/2019ApJ...882...36M/abstract>`__
+
+Initialization of the model
+===========================
+The initial mass of the helium star is set in ``inlist_extra``
+
+.. literalinclude:: ../../../star/test_suite/ppisn/inlist_extra
+
+In this case we use a :math:`72 M_\odot`
+
+Last-Updated: 2019-11-12 (mesa r12413) by Pablo Marchant
+

star/dev_cases_test_TDC/dev_TDC_Cepheid_6M/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_TDC/dev_TDC_Cepheid_6M/clean

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

star/dev_cases_test_TDC/dev_TDC_Cepheid_6M/gyre.in

@@ -0,0 +1,33 @@
+&model
+  add_center = .FALSE.
+/
+
+&mode
+   l = 0
+/
+
+&osc
+   nonadiabatic = .TRUE.
+  inner_bound = 'ZERO_R'
+/
+
+&rot
+/
+
+&num
+   diff_scheme = 'MAGNUS_GL2' 
+/
+
+&scan
+   grid_type = 'LINEAR'
+   freq_min = 0.5
+   freq_max = 5.0
+   freq_min_units = 'ACOUSTIC_DELTA'
+   freq_max_units = 'ACOUSTIC_DELTA'
+   n_freq = 50
+/
+
+&grid
+   w_osc = 10
+   w_exp = 2
+/