[mesa-users] Standard sollar model
Liliana Novais
linovais at gmail.com
Tue Aug 18 07:08:27 EDT 2015
Hello,
Yesterday I was offline.
Here it goes the inlist.
Thanks,
Liliana
Às 08:13 de 17-08-2015, Aaron Dotter escreveu:
> Hi Liliana,
>
> Please include your inlist so that we can see what you are doing.
> What do you want the initial values of CNO to be?
>
> Aaron
>
> On Monday, August 17, 2015, Liliana Novais <linovais at gmail.com
> <mailto:linovais at gmail.com>> wrote:
>
> Hi i run the program but the initial quantities of c n and o aee
> zero, i need it to have the "real" quantities.
> Thanks
>
> Liliana Novais
>
> Em 17/08/2015 03:24, "Josiah Schwab" <jwschwab at berkeley.edu
> <javascript:_e(%7B%7D,'cvml','jwschwab at berkeley.edu');>> escreveu:
>
> Hello Liliana,
>
> On 16 August 2015 at 13:10 PDT, Liliana Novais wrote:
>
> > I am new at Mesa, finally I understand how it works and I
> started to
> > run the examples that came with it. The example of the solar
> model has
> > the cno cycle starting at zero, that it isn’t true… Where
> can I find
> > the files for the ssm??
>
> Unfortunately, your message does not contain enough
> information for us
> to be able to provide a useful response.
>
> The example_solar_model test case, available in
>
> $MESA_DIR/star/test_suite/example_solar_model
>
> uses the 'pp_and_cno_extras.net
> <http://pp_and_cno_extras.net>', which has both pp chain and
> cno cycle
> reactions.
>
> Could you provide more information about what you've done, what
> difficulties you've encountered, and what you're hoping to
> accomplish?
>
> Best,
> Josiah
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150818/8308b098/attachment.html>
-------------- next part --------------
&star_job
read_extra_star_job_inlist1 = .true.
extra_star_job_inlist1_name = '../inlist_test_suite'
read_extra_star_job_inlist2 = .true.
extra_star_job_inlist2_name = 'inlist_solar_model'
/ ! end of star_job namelist
&controls
read_extra_controls_inlist1 = .true.
extra_controls_inlist1_name = 'inlist_solar_model'
/ ! end of controls namelist
&pgstar
read_extra_pgstar_inlist1 = .true.
extra_pgstar_inlist1_name = 'inlist_solar_model'
/ ! end of pgstar namelist
-------------- next part --------------
&star_job
show_log_description_at_start = .false.
load_saved_model = .true.
saved_model_name = 'solar_calibration_input.mod'
save_model_when_terminate = .true.
save_model_filename = 'solar.mod'
write_profile_when_terminate = .true.
filename_for_profile_when_terminate = 'final_profile.data'
change_net = .true. ! switch nuclear reaction network
new_net_name = 'pp_and_cno_extras.net'
set_uniform_initial_composition = .true.
initial_zfracs = 3 ! GS98_zfracs = 3
initial_h1 = 6.9977251770969628D-01
initial_h2 = 0
initial_he3 = 2.7978169353764870D-05
initial_he4 = 2.7975371536829496D-01
! initial_z = 2.0445788752654973D-02
! initial_y = 0.27978169353764870; FeH = 0.10524293223273457
kappa_file_prefix = 'OP_gs98'
kappa_lowT_prefix = 'lowT_fa05_gs98' ! for lower temperatures.
! THIS IS A TEMPORARY PATCH TO MAKE THE TEST_CASE WORK USING OLD SOLAR CALIBRATION RESULTS
kappa_blend_logT_upper_bdy = 4.1d0 ! new default is better. 3.88d0
kappa_blend_logT_lower_bdy = 4.0d0 ! new default is better. 3.80d0
set_rate_c12ag = 'Kunz'
set_rate_n14pg = 'jina reaclib'
change_lnPgas_flag = .true.
new_lnPgas_flag = .true.
set_initial_model_number = .true.
initial_model_number = 0
!pgstar_flag = .true.
/ ! end of star_job namelist
&controls
! check for retries and backups as part of test_suite
! you can/should delete this for use outside of test_suite
max_number_backups = 0
max_number_retries = 0
max_model_number = 1500
! timesteps
max_years_for_timestep = 1d7
varcontrol_target = 1d-3
max_timestep_factor = 2d0
delta_lgT_cntr_limit = 0.1
delta_lgRho_cntr_limit = 0.5
max_age = 4.61d9 ! 4.57d9 + 4d7 for pre-ms
num_adjusted_dt_steps_before_max_age = 200
dt_years_for_steps_before_max_age = 1d6
mixing_length_alpha = 1.9107649497310391D+00
overshoot_f_above_nonburn_core = 2.9091163079927995D-03
overshoot_f_above_nonburn_shell = 2.9091163079927995D-03
overshoot_f_below_nonburn_shell = 2.9091163079927995D-03
overshoot_f_above_burn_h_core = 2.9091163079927995D-03
overshoot_f_above_burn_h_shell = 2.9091163079927995D-03
overshoot_f_below_burn_h_shell = 2.9091163079927995D-03
overshoot_f0_above_nonburn_core = 2.9091163079927995D-03
overshoot_f0_above_nonburn_shell = 2.9091163079927995D-03
overshoot_f0_below_nonburn_shell = 2.9091163079927995D-03
overshoot_f0_above_burn_h_core = 2.9091163079927995D-03
overshoot_f0_above_burn_h_shell = 2.9091163079927995D-03
overshoot_f0_below_burn_h_shell = 2.9091163079927995D-03
! controls for output
photostep = 50
profile_interval = 100
history_interval = 10
terminal_cnt = 10
write_header_frequency = 10
! atmosphere
which_atm_option = 'photosphere_tables'
! atomic diffusion
do_element_diffusion = .true. ! determines whether or not we do diffusion
diffusion_dt_limit = 7d11 ! no element diffusion if dt < this limit (in seconds)
diffusion_T_full_on = 1d3
diffusion_T_full_off = 1d3
diffusion_calculates_ionization = .true.
diffusion_num_classes = 4 ! number of classes of species for diffusion calculations
diffusion_class_representative(1) = 'h1'
diffusion_class_representative(2) = 'he4'
diffusion_class_representative(3) = 'o16'
diffusion_class_representative(4) = 'fe56'
! in ascending order. species goes into 1st class with A_max >= species A
diffusion_class_A_max(1) = 2
diffusion_class_A_max(2) = 4
diffusion_class_A_max(3) = 16
diffusion_class_A_max(4) = 10000
diffusion_use_isolve = .true.
diffusion_rtol_for_isolve = 1d-4
diffusion_atol_for_isolve = 1d-5
diffusion_maxsteps_for_isolve = 1000
diffusion_isolve_solver = 'ros2_solver'
! mesh adjustment
mesh_delta_coeff = 0.8
P_function_weight = 25
T_function1_weight = 75
xtra_coef_czb_full_on = 1
xtra_coef_czb_full_off = 1
!xtra_coef_czb_full_on_if_center_he4_below_this = 1
!xtra_coef_czb_full_off_if_center_he4_above_this = 1
xtra_coef_a_l_nb_czb = 0.15 ! above lower nonburn convective boundary
xtra_dist_a_l_nb_czb = 20 ! above lower nonburn convective boundary
xtra_coef_b_l_nb_czb = 0.15 ! below lower nonburn convective boundary
xtra_dist_b_l_nb_czb = 20 ! below lower nonburn convective boundary
xa_function_species(1) = 'he4' ! name of nuclide as defined in chem_def
xa_function_weight(1) = 130
xa_function_param(1) = 1d-2
xa_function_species(2) = 'he3' ! name of nuclide as defined in chem_def
xa_function_weight(2) = 60
xa_function_param(2) = 1d-5
! opacity
use_Type2_opacities = .false.
cubic_interpolation_in_Z = .true.
! 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
/ ! end of controls namelist
&pgstar
Grid1_win_flag = .true.
Grid1_win_width = 8
Grid1_win_aspect_ratio = 0.75 ! aspect_ratio = height/width
show_TRho_Profile_eos_regions = .false.
TRho_Profile_xmin = -5.1
TRho_Profile_xmax = 10.2
TRho_Profile_ymin = 5.1
TRho_Profile_ymax = 10.2
!Grid1_file_flag = .true.
Grid1_file_dir = 'png'
Grid1_file_prefix = 'Grid1_'
Grid1_file_cnt = 5 ! output when mod(model_number,Grid1_file_cnt)==0
Grid1_file_width = -1 ! (inches) negative means use same value as for window
Grid1_file_aspect_ratio = -1 ! negative means use same value as for window
Abundance_win_flag = .true.
Abundance_win_width = 6
Abundance_win_aspect_ratio = 1 ! aspect_ratio = height/width
Abundance_xaxis_name = 'mass'
Abundance_xaxis_reversed = .false.
Abundance_xmin = 0.00 ! -101d0 ! only used if /= -101d0
Abundance_xmax = -101 ! only used if /= -101d0
Abundance_log_mass_frac_min = -6 ! only used if < 0
!Abundance_file_flag = .true.
Abundance_file_dir = 'png'
Abundance_file_prefix = 'abund_'
Abundance_file_cnt = 5 ! output when mod(model_number,Abundance_file_cnt)==0
Abundance_file_width = -1 ! (inches) negative means use same value as for window
Abundance_file_aspect_ratio = -1 ! negative means use same value as for window
Profile_Panels1_win_flag = .true.
Profile_Panels1_win_width = 6
Profile_Panels1_win_aspect_ratio = 1 ! aspect_ratio = height/width
Profile_Panels1_title = 'Profile Panels1'
Profile_Panels1_num_panels = 3
Profile_Panels1_yaxis_name(1) = 'logT'
Profile_Panels1_other_yaxis_name(1) = 'entropy'
Profile_Panels1_yaxis_name(2) = 'logRho'
Profile_Panels1_other_yaxis_name(2) = 'opacity'
Profile_Panels1_yaxis_name(3) = 'zone'
Profile_Panels1_ymin(3) = -101 ! only used if /= -101d0
Profile_Panels1_ymax(3) = -101 ! only used if /= -101d0
Profile_Panels1_ymargin(3) = 0.1
Profile_Panels1_other_yaxis_name(3) = 'grada'
Profile_Panels1_other_ymin(3) = -101 ! only used if /= -101d0
Profile_Panels1_other_ymax(3) = -101d0 ! only used if /= -101d0
Profile_Panels1_other_ymargin(:) = 0.1
Profile_Panels1_xaxis_name = 'mass'
Profile_Panels1_xaxis_reversed = .false.
Profile_Panels1_xmin = -101d0 ! only used if /= -101d0
Profile_Panels1_xmax = -101d0 ! only used if /= -101d0
Profile_Panels1_xmargin = 0d0
! file output
Profile_Panels1_file_flag = .false.
Profile_Panels1_file_dir = 'png'
Profile_Panels1_file_prefix = 'profile_panels1_'
Profile_Panels1_file_cnt = 5 ! output when mod(model_number,Profile_Panels1_file_cnt)==0
Profile_Panels1_file_width = -1 ! (inches) negative means use same value as for window
Profile_Panels1_file_aspect_ratio = -1 ! negative means use same value as for window
History_Panels1_win_flag = .true.
History_Panels1_win_width = 7
History_Panels1_win_aspect_ratio = 1 ! aspect_ratio = height/width
History_Panels1_title = 'History_Panels1'
History_Panels1_num_panels = 3
History_Panels1_yaxis_name(1) = 'log_center_T'
History_Panels1_yaxis_reversed(1) = .false.
History_Panels1_ymin(1) = -101d0 ! only used if /= -101d0
History_Panels1_ymax(1) = -101d0 ! only used if /= -101d0
History_Panels1_dymin(1) = -1
History_Panels1_other_yaxis_name(1) = 'log_center_Rho'
History_Panels1_other_yaxis_reversed(1) = .false.
History_Panels1_other_ymin(1) = -101d0 ! only used if /= -101d0
History_Panels1_other_ymax(1) = -101d0 ! only used if /= -101d0
History_Panels1_other_dymin(1) = -1
History_Panels1_yaxis_name(2) = 'log_L'
History_Panels1_yaxis_reversed(2) = .false.
History_Panels1_ymin(2) = -101d0 ! only used if /= -101d0
History_Panels1_ymax(2) = -101d0 ! only used if /= -101d0
History_Panels1_dymin(2) = 0.01
History_Panels1_other_yaxis_name(2) = 'log_Teff'
History_Panels1_other_yaxis_reversed(2) = .false.
History_Panels1_other_ymin(2) = -101d0 ! only used if /= -101d0
History_Panels1_other_ymax(2) = -101d0 ! only used if /= -101d0
History_Panels1_other_dymin(2) = 0.01
History_Panels1_yaxis_name(3) = 'log_dt'
History_Panels1_yaxis_reversed(3) = .false.
History_Panels1_ymin(3) = -101d0 ! only used if /= -101d0
History_Panels1_ymax(3) = -101d0 ! only used if /= -101d0
History_Panels1_dymin(3) = -1
History_Panels1_other_yaxis_name(3) = 'num_zones'
History_Panels1_other_yaxis_reversed(3) = .false.
History_Panels1_other_ymin(3) = -101d0 ! only used if /= -101d0
History_Panels1_other_ymax(3) = -101d0 ! only used if /= -101d0
History_Panels1_other_dymin(3) = -1
! file output
History_Panels1_file_flag = .false.
History_Panels1_file_dir = 'png'
History_Panels1_file_prefix = 'History_Panels1_'
History_Panels1_file_cnt = 5 ! output when mod(model_number,History_Panels1_file_cnt)==0
History_Panels1_file_width = -1 ! (inches) negative means use same value as for window
History_Panels1_file_aspect_ratio = -1 ! negative means use same value as for window
/ ! end of pgstar namelist
More information about the Mesa-users
mailing list