Re: [mesa-users] astero...

SourceForge Headquarters 225 Broadway Suite 1600 San Diego, CA 92101 +1 (858) 422-6466

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