Re: [mesa-users] astero...

[Frank T. <fx...@ma...>]

thanks. 

see warricks's message:
> you've added the astero stuff to &controls in the usual MESA inlist, rather than to &astero_search_controls.


fxt



> On May 16, 2017, at 7:04 AM, Bruno Lustosa <bru...@gm...> wrote:
> 
> mesa-r8845
> 
> 
> ##########################################################################
> 
> &star_job
> 
>       pause_before_terminate = .true.
>       create_pre_main_sequence_model = .true.
>       pre_ms_relax_num_steps = 50
>       pre_ms_T_c = 9e5 
> 
>       
>      
>       !change_initial_net = .true.      
>       !new_net_name = 'agb.net'
>       !set_abundance = .true.
>       set_initial_abundance = .true.
>       chem_name = 'li7'
>       new_frac = 7.77594697D-09
>       !new_frac = 1.06927581475e-08
>       !set_abundance_nzlo = -1
>       !set_abundance_nzhi = -1
> 
> 
>       change_net=.true.
>       eos_file_prefix = 'mesa'
>       new_net_name = 'pp_extras.net'
>       kap_cache_dir = '/home/bruno/mesa-r8845/data/kap_data/cache'
>       kappa_file_prefix = 'OP_gs98'
>       dump_missing_metals_into_heaviest = .false.  
>       set_uniform_initial_composition = .true.
>       initial_h1 = 0.6997725177
>       !initial y = 0.2797816935
>       initial_he3 = 0.0000279782
>       initial_he4 = 0.2797537154
>       initial_h2 = 0
>       
>       set_rates_preference = .true.
>       new_rates_preference = 1 
>      
>       
>       new_rotation_flag = .true.
>       change_rotation_flag = .false.
>       change_initial_rotation_flag = .true.
>       
> 
>       initial_zfracs = 7
>       new_surface_rotation_v = 2 ! solar (km sec^1)
>       set_near_zams_surface_rotation_v_steps = 10
>          
> 
> 
>       pgstar_flag = .true.
> 
> 
> / ! end of star_job namelist
> 
> &controls
>       !max_age = 12000000000   
>        stop_at_TP = .true.
>       ! check for retries and backups as part of test_suite
>       ! you can/should delete this for use outside of test_suite
>          max_number_backups = 200
>          max_number_retries = 500
>          
>          
>       max_model_number = 6000 !958 max
> 
>       
>       initial_mass = 1.08
>       initial_z = 0.0204457888 
>       !initial_z = 0.02 
>       initial_y = 0.2797816935
>       !use_Type2_opacities = .true.
>       Zbase = 0.0204457888 
> 
>       am_nu_visc_factor = 0
>       am_D_mix_factor = 0.0333333333333333d0
>       D_DSI_factor = 0
>       D_SH_factor = 1
>       D_SSI_factor = 1
>       D_ES_factor = 1
>       D_GSF_factor = 1
>       D_ST_factor = 1
>       
>       varcontrol_target = 1d-3
>       mesh_delta_coeff = 1.5
> 
>       
>       photo_interval = 10
>       profile_interval = 50
>       history_interval = 1
>       terminal_interval = 10
>       write_header_frequency = 10
> 
>       write_pulse_info_with_profile = .true. !.true. ! if true, write pulse info file when write profile
>       pulse_info_format = 'FGONG'
> 
>       log_L_lower_limit = -1
>       mixing_length_alpha = 1.9107649497
>       MLT_option = 'Cox'
> 		! options are:
> 		!    'Cox'        MLT as developed in Cox & Giuli 1968, Chapter 14.
> 		!    'ML1'        Bohm-Vitense 1958
> 		!    'ML2'        Bohm and Cassinelli 1971
> 		!    'Mihalas'    Mihalas 1978, Kurucz 1979
> 		!    'Henyey'     Henyey, Vardya, and Bodenheimer 1965
> 		! 'Cox' option assumes optically thick material.
> 		! the other options are various ways of extending to include optically thin material. 
> 
>       !which_atm_option = 'simple_photosphere'
>       smooth_convective_bdy = .true.                  
>       convective_bdy_weight = 1
>       D_visc_factor = 0  
>      
>       cool_wind_full_on_T = 9d6
>       cool_wind_full_off_T = 1d7 
>       cool_wind_RGB_scheme = 'Reimers'
>       cool_wind_AGB_scheme = 'Blocker'
>       RGB_to_AGB_wind_switch = 1d-4
>       Reimers_scaling_factor = 0.5d0  
>       Blocker_scaling_factor = 0.5d0  
> 
> 
>       ! FOR DEBUGGING
> 
>       !report_hydro_solver_progress = .true. ! set true to see info about newton iterations
>       !report_ierr = .true. ! if true, produce terminal output when have some internal error
>       !hydro_show_correction_info = .true.
>       
>       !max_years_for_timestep = 3.67628942044319d-05
> 
>       !report_why_dt_limits = .true.
>       !report_all_dt_limits = .true.
>       
>       !show_mesh_changes = .true.
>       !mesh_dump_call_number = 5189
>       !okay_to_remesh = .false.
>       
>       trace_evolve = .true.
>             
> 
>       ! hydro debugging
>       !hydro_check_everything = .true.
>       !hydro_inspectB_flag = .true.
>       
>       !hydro_numerical_jacobian = .true.
>       !hydro_save_numjac_plot_data = .true.
>       !small_mtx_decsol = 'lapack'
>       !large_mtx_decsol = 'lapack'
>       !hydro_dump_call_number = 195
> 
>       !trace_newton_bcyclic_solve_input = .true. ! input is "B" j k iter B(j,k)
>       !trace_newton_bcyclic_solve_output = .true. ! output is "X" j k iter X(j,k)
>       
>       !trace_newton_bcyclic_steplo = 1 ! 1st model number to trace
>       !trace_newton_bcyclic_stephi = 1 ! last model number to trace
>       
>       !trace_newton_bcyclic_iterlo = 2 ! 1st newton iter to trace
>       !trace_newton_bcyclic_iterhi = 2 ! last newton iter to trace
>       
>       !trace_newton_bcyclic_nzlo = 1 ! 1st cell to trace
>       !trace_newton_bcyclic_nzhi = 10000 ! last cell to trace; if < 0, then use nz as nzhi
>       
>       !trace_newton_bcyclic_jlo = 1 ! 1st var to trace
>       !trace_newton_bcyclic_jhi = 100 ! last var to trace; if < 0, then use nvar as jhi
>       
>       !trace_k = 0
>       
>       ! misc
> 	    add_atmosphere_to_pulse_info = .true.
>    	add_center_point_to_pulse_info = .true. ! add point for r=0
>    	stop_at_TP = .true.  
>      
> 
> 
>       !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
>       ! seismic data (frequencies in microHz)
> 
>       !nu_max = 1010
>          ! nu_max is needed when 
>             ! chi2_seismo_nu_max_fraction > 0 or correction_factor > 0 (see below)
>       !nu_max_sigma = -1
>          ! nu_max_sigma is needed when chi2_seismo_nu_max_fraction > 0
> 
>       !delta_nu = 56.28
>       !delta_nu_sigma = 1.0
>          
>          !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
> / ! end of controls namelist
> 
> 
> 
> &pgstar
>             
>         Power_win_flag = .false.
>         Power_win_width = 9
> 	HR_win_flag = .true.
> 	HR_win_width = 10
> 	HR_win_aspect_ratio = 1
>         logg_logT_win_flag = .false. 
> 	logg_logT_win_width = 8
>         logg_logT_win_aspect_ratio = 1.0
> 	Abundance_win_flag = .false.
> 	Abundance_win_width = 6
>         Abundance_win_aspect_ratio = 1
>         HR_file_flag = .true.
>         HR_file_dir = 'png'
>         HR_file_prefix = 'HR_'
> 	HR_title = 'DHR - Bruno'
> 
>  ! Mode Propagation window
>       
>          Mode_Prop_win_flag = .false.
>       
>          Mode_Prop_win_width = 10
>          Mode_Prop_win_aspect_ratio = 0.62 ! aspect_ratio = height/width
>          
>          Mode_Prop_nu_max_obs = 1000
>       
>          Mode_Prop_xaxis_name = 'r_div_R' 
>          
>          ! power xaxis limits -- to override system default selections
>          Mode_Prop_xmin = -101 ! only used if > -100
>          Mode_Prop_xmax = 0.3 !-101 ! only used if > -100
>          
>          ! power yaxis limits -- to override system default selections
>          Mode_Prop_ymin = 2 ! -101 ! only used if > -100
>          Mode_Prop_ymax = 5 ! -101 ! only used if > -100
>          
>          ! file output
>          Mode_Prop_file_flag = .false.
>          Mode_Prop_file_dir = 'pgstar_out'
>          Mode_Prop_file_prefix = 'mode_prop'
>          Mode_Prop_file_interval = 5 ! output when mod(model_number,Mode_Prop_file_interval)==0
>          Mode_Prop_file_width = -1 ! negative means use same value as for window
>          Mode_Prop_file_aspect_ratio = -1 ! negative means use same value as for window   
>      
>    
> / ! end of pgstar namelist
> <inlist_2.0>