[mesa-users] Low metalicity star
Liliana Novais
linovais at gmail.com
Sat Oct 31 11:08:46 EDT 2015
Hello,
I am making a model of a low metallicity star and I am trying to save
the plots to analyze after it runs but it gives me back an error:
save LOGS/profile1.data for model 1
%PGPLOT, Too many active plots.
PGSTAR: failed to open /xwin
please check that X Window System is running
and that you have set the environment variables
DISPLAY :0
PGPLOT_DEV /xwin
and also that you have set PGPLOT_DIR and PGPLOT_FONT
Grid1 ierr -1
init_pgstar ierr -1
onScreen_Plots ierr -1
update_pgstar_plots ierr -1
after_step_loop ierr -1
The inlist is attached... I don't know what is going wrong...
Thanks,
Liliana
-------------- next part --------------
! this is the master inlist that MESA reads when it starts.
! This file tells MESA to go look elsewhere for its configuration
! info. This makes changing between different inlists easier, by
! allowing you to easily change the name of the file that gets read.
&star_job
create_pre_main_sequence_model = .true.
pre_ms_relax_num_steps = 10
set_initial_age = .true.
initial_age = 0 ! in years
set_initial_model_number = .true.
initial_model_number = 0
change_initial_net = .true.
new_net_name = 'o18_and_ne22.net'
new_omega_div_omega_crit = 0.0
set_near_zams_omega_div_omega_crit_steps = 50
pgstar_flag = .true.
/ !end of star_job namelist
&controls
! starting specifications
initial_mass = 0.83 ! in Msun units
initial_z = 0.004
initial_y = 0.28
max_num_profile_models = -1
RGB_wind_scheme = 'Reimers'
Reimers_wind_eta = 0.5
AGB_wind_scheme = 'Blocker'
Blocker_wind_eta = 0.1
profile_interval = 100
history_interval = 1
terminal_cnt = 10000000
write_header_frequency = 1000000
photostep = 500000
max_age = 9d10
varcontrol_target = 1d-4
mesh_delta_coeff = 1.0
dX_nuc_drop_min_X_limit = 1d-4
mixing_length_alpha = 2.0
MLT_option = 'Cox'
am_D_mix_factor = 0.033333333d0
D_DSI_factor = 1
D_SH_factor = 1
D_SSI_factor = 1
D_ES_factor = 1
D_GSF_factor = 1
D_ST_factor = 1
/ ! end of controls namelist
&pgstar
read_extra_pgstar_inlist1 = .true.
extra_pgstar_inlist1_name = 'inlist_pgstar'
Grid1_win_flag = .true.
Grid1_win_width = 8
Grid1_win_aspect_ratio = 0.75 ! aspect_ratio = height/width
show_TRho_Profile_eos_regions = .true.
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 = .true.
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 = .true.
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
-------------- next part --------------
&pgstar
! MESA uses PGPLOT for live plotting and gives the user a tremendous
! amount of control of the presentation of the information.
! show HR diagram
! this plots the history of L,Teff over many timesteps
HR_win_flag = .true.
! set static plot bounds
HR_logT_min = 3.5
HR_logT_max = 4.6
HR_logL_min = 2.0
HR_logL_max = 6.0
! set window size (aspect_ratio = height/width)
HR_win_width = 6
HR_win_aspect_ratio = 1.0
! show temperature/density profile
! this plots the internal structure at single timestep
TRho_Profile_win_flag = .true.
! add legend explaining colors
show_TRho_Profile_legend = .true.
! display numerical info about the star
show_TRho_Profile_text_info = .true.
! set window size (aspect_ratio = height/width)
TRho_Profile_win_width = 8
TRho_Profile_win_aspect_ratio = 0.75
/ ! end of pgstar namelist
-------------- next part --------------
!This in an Inlist to evolve a ow metalicity solar type star called HD 10700
!First atempt on doing it
&star_job
show_log_description_at_start = .false.
!Save the model
save_model_when_terminate = .true.
save_model_filename = 'HD10700.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 = 6 ! AGSS09_zfracs = 6
initial_h1 = 6.9977251770969628D-01
initial_h2 = 0
initial_he3 = 2.7978169353764870D-05
initial_he4 = 2.7975371536829496D-01
!initial_z = 6.5D-03
! initial_y = 0.27978169353764870; FeH = 0.10524293223273457
kappa_file_prefix = 'OP_a09_nans_removed_by_hand'
kappa_lowT_prefix = 'lowT_fa05_a09p' ! 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.
save_pgstar_files_when_terminate=.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 = 9.5d9 ! 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 = .true.
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 = .true.
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 = .true.
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