&star_job

    mesa_dir = '' 

    save_model_when_terminate = .true.
    save_model_filename = '3M_FeHm030_ZACHeB.mod'

    history_columns_file = 'history_columns.list'
    profile_columns_file = 'profile_columns.list'

    create_pre_main_sequence_model = .true.
    pre_ms_T_c = 4d5
    
    save_model_when_terminate = .true.
    save_model_filename = 'M3_FeHm030_ZACHeB.mod'

    change_net = .true. 
    new_net_name = 'sagb_NeNa_MgAl.net'

    set_rates_preference = .true.
    new_rates_preference = 2

    set_uniform_initial_xa_from_file = .true.
    file_for_uniform_xa = 'init_comp_fehm0.3_y27_ap4_a09.data'

    kappa_file_prefix = 'a09'
    kappa_CO_prefix = 'a09_co'
    kappa_lowT_prefix = 'lowT_fa05_a09p'

    pgstar_flag = .false.

/ !end of star_job namelist


&controls

    use_other_D_mix = .false.

    use_Ledoux_criterion = .false.

    thermohaline_option = 'Kippenhahn'
    thermohaline_coeff = 667

    initial_mass = 3.0 

    initial_y = 0.27
    initial_z = 0.0066129064

    cubic_interpolation_in_Z = .true.
    use_Type2_opacities = .true.
    Zbase = 0.0066129064 

    max_num_profile_models = 1000
    history_interval = 1
    profile_interval = 100
    photo_interval = 500

    do_element_diffusion = .false.
    diffusion_dt_limit = 3.15d11 

    HB_limit = 0.95

    MLT_option = 'Henyey'

    mixing_length_alpha = 2.0

    atm_T_tau_relation = 'Krishna_Swamy'

    mesh_delta_coeff = 1.0

    min_thermohaline_gap = 0.01

    min_thermohaline_dropout = +1
    max_dropout_gradL_sub_grada = 1d-3

    overshoot_f = 0.015
    overshoot_f0 = 0.001
    overshoot_scheme = 'exponential'
    overshoot_zone_type = 'any'
    overshoot_zone_loc = 'any'
    overshoot_bdy_loc = 'any'

/ ! end of controls namelist