[mesa-users] solar mode with AGSS09

Sun Aug 30 13:51:08 EDT 2015

Hi,

I am trying to change the solar model from GS98 to AGSS09 I found that I 
had to change from:

|initial_zfracs = 3 to ||initial_zfracs = 6 |

But then I have a question on the changes that I need to do more. Is it 
possible to see if my inlist is ok???? thanks, Liliana
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150830/1cd43c30/attachment.html>
-------------- 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 = 6 ! AGSS09_zfracs = 6
      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_AGSS09'
      kappa_lowT_prefix = 'lowT_fa05_AGSS09' ! 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 = 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 = .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