&star_job

      mesa_dir = ''


      save_model_when_terminate = .true.
      !save_model_filename = 'end_h.mod'
      save_model_filename = 'end_oburn.mod'

      !load_saved_model = .true.
      !saved_model_name = 'end_h.mod'

      ! begin with a pre-main sequence model
      create_pre_main_sequence_model = .true.

      pre_ms_relax_num_steps = 200

      !## reaction network
      change_net = .true.
      new_net_name = 'basic.net' ! basic network; switches to bigger nets when
      ! certain conditions are met in order to reach oxygen depletion with the 
      ! 'auto_extend_net = .true.'
      ! new_net_name = 'o_burn_full.net' ! trunkated reaction network up to ca46
      ! new_net_name = 'approx21_plus_co56.net' ! use smaller network to make test for massive stars only

      !## opacity
      !kappa_file_prefix = 'gs98'
      !initial_zfracs = 3
      kappa_file_prefix = 'a09'
      kappa_lowT_prefix = 'lowT_fa05_a09p'
      initial_zfracs = 6
      kappa_CO_prefix = 'a09_co'

      !## reaction rates
      set_rates_preference = .true. ! for use by net + rates modules; default is false
      new_rates_preference = 2 ! latest jina includes most of the NACRE rates tbc
         ! 1 = NACRE rates -- this is the default
         ! 2 = jina reaclib rates -- to match jina where possible
 
      !jina_reaclib_filename = 'jina_reaclib_results_20171204default'
         !! latest version from JINA Reaclib Database
      !set_rate_c12ag = 'jina reaclib' ! jina_reaclib usually uses Kunz 
         !! one of 'NACRE', 'jina reaclib', 'Kunz', or 'CF88'
         !! note: original CF88 rate is actually multiplied by 1.7 as in Timmes' rates
      !set_rate_n14pg = 'jina reaclib' ! empty string (default) means ignore this control
         !! one of 'NACRE', 'jina reaclib', or 'CF88'
      !set_rate_3a = 'jina reaclib' ! empty string (default) means ignore this control
         !! one of 'NACRE', 'jina reaclib', or 'CF88


      !## modified history.data and profile.data
      history_columns_file = 'history_columns.list'
      profile_columns_file = 'profile_columns.list'

      ! change_v_flag = .true.
      ! new_v_flag = .true.

    ! set omega/omega_crit ratio
    new_omega_div_omega_crit = 0.6
    ! set_near_zams_omega_div_omega_crit_steps = 10 ! this option sets the new omega/omega_crit for x steps
    num_steps_to_relax_rotation = 100
    relax_omega_div_omega_crit = .true.
    near_zams_relax_omega_div_omega_crit = .true. ! this option relaxes the model until is has reached the 
      					     	     ! specified omega/omega_crit


/ ! end of star_job namelist
 
&controls

      ! steps_before_use_gold_tolerances = 1000
      !use_gold_tolerances = .false.

      warn_when_large_rel_run_E_err = 1.0d99
      extra_terminal_output_file = 'log'

      initial_mass = 75.0
      initial_z = 0.014d0
      !initial_y = 0.986d0
      !initial_he3 = 0d0

      ! atmosphere
      !which_atm_option = 'simple_photosphere'  
      Pextra_factor = 2 ! Pressure at optical depth tau is calculated as `P = tau*g/kap*(1 + Pextra)`

      ! convection
      mixing_length_alpha = 1.67d0
      !MLT_option = 'Henyey'

      !use_simple_es_for_kap= .true.
      use_Type2_opacities = .true.
      Zbase=0.014d0
      okay_to_reduce_gradT_excess = .true.
      !opacity_max = 0.5


      !## controls output
      photo_interval = 100 
      profile_interval = 50 
      history_interval = 1
      terminal_interval = 10
      max_num_profile_models = -1

      !## definition of core boundaries

      ! he core
      he_core_boundary_h1_fraction = 0.01
      ! c core
      c_core_boundary_he4_fraction = 0.01
      ! o core
      o_core_boundary_c12_fraction = 0.01
      ! si core
      si_core_boundary_o16_fraction = 0.01
      ! iron core
      fe_core_boundary_si28_fraction = 0.01

      ! wind mass loss
      ! 'Dutch' scheme combines several Dutch wind schemes: For O-stars the 'Vink' scheme and for
      ! wolf-rayet stars the one from Nugis & Lamers (2000)
      hot_wind_scheme = 'Dutch'
      cool_wind_AGB_scheme = 'Dutch'
      cool_wind_RGB_scheme = 'Dutch'
      Dutch_wind_lowT_scheme = 'de Jager' !'van Loon' gives a too strong mass loss! for Teff<1e4
      Dutch_scaling_factor = 0.49 ! as in GENEC (Maeder & Meynet (2001))

      !max_wind=1.e-99



      !## when to stop

      ! hydrogen depletion
      ! xa_central_lower_limit_species(1) = 'h1'
      ! xa_central_lower_limit(1) = 0.00001

      ! helium depletion
      ! xa_central_lower_limit_species(2) = 'he4'
      ! xa_central_lower_limit(2) = 0.001

      ! oxygen depletion
      xa_central_lower_limit_species(3) = 'o16'
      xa_central_lower_limit(3) = 0.00001
      
      ! silicon depletion
      ! xa_central_lower_limit_species(4) = 'si28'
      ! xa_central_lower_limit(4) = 0.01

      ! collapsing core
      ! fe_core_infall_limit = 1d8 ![cm/sec]

      fe56ec_fake_factor = 1d-4
      !rate_factor_for_fe56ec_to_cr56 does not exist in MESA 8118 anymore



      !## overshoot
      overshoot_f0_above_nonburn_core = 0.004
      overshoot_f_above_nonburn_core = 0.004
      overshoot_f0_above_nonburn_shell = 0.004
      overshoot_f_above_nonburn_shell = 0.004
      overshoot_f0_below_nonburn_shell = 0.004
      overshoot_f_below_nonburn_shell = 0.004     

      overshoot_f0_above_burn_h_core = 0.004
      overshoot_f_above_burn_h_core = 0.004
      overshoot_f0_above_burn_h_shell = 0.004
      overshoot_f_above_burn_h_shell = 0.004
      overshoot_f0_below_burn_h_shell = 0.004
      overshoot_f_below_burn_h_shell = 0.004

      overshoot_f0_above_burn_he_core = 0.004
      overshoot_f_above_burn_he_core = 0.004
      overshoot_f0_above_burn_he_shell = 0.004
      overshoot_f_above_burn_he_shell = 0.004
      overshoot_f0_below_burn_he_shell = 0.004
      overshoot_f_below_burn_he_shell = 0.004

      overshoot_f0_above_burn_z_core = 0.004
      overshoot_f_above_burn_z_core = 0.004
      overshoot_f0_above_burn_z_shell = 0.004
      overshoot_f_above_burn_z_shell = 0.004
      overshoot_f0_below_burn_z_shell = 0.004
      overshoot_f_below_burn_z_shell = 0.004

    ! rotation
    skip_rotation_in_convection_zones =.true.
    am_D_mix_factor = 0.03333333333333333d00 ! Heger et al. 2000

    ! AM lost through winds
    !actual_J_lost = adjust_J_fraction*mass_lost*j_rot_avg_surf + (1-adjust_J_fraction)*angular_momentum_removed
    do_adjust_J_lost = .true.
    adjust_J_fraction = 1d0

      ! rotationally induced instabilities (default is 0 for all)
      ! Each rotationally induced diffusion coefficient has a factor that lets you control it. Value of 1 gives 
      ! normal strength; value of 0 turns it off.
      ! These controls are used to determine D_mix, the diffusion coefficient for mixing material.
      ! Note that for regions with brunt N^2 < 0, we set Richardson number to 1, which is > Ri_critical and 
      ! therefore turns off DSI and SSI.
      ! dynamical shear instability
    D_DSI_factor = 0.0
      ! Solberg-Hoiland
    D_SH_factor  = 0.0
      ! secular shear
    D_SSI_factor = 1.0
      ! Eddington-Sweet circulation
    D_ES_factor  = 1.0
      ! Goldreich-Schubert-Fricke
    D_GSF_factor = 0.0
      ! Spruit-Tayler dynamo
    D_ST_factor  = 1.0
      ! Those factors are also used for angular momentum transport (diffusive in MESA) as long as 
      ! 'am_nu_*_factor' < 0

      ! if logT smaller than this -> mixing full on
    D_mix_rotation_max_logT_full_on = 9.4d0
      ! if logT larger than this -> mixing full_off
    D_mix_rotation_min_logT_full_off = 9.5d0

        ! smooth over grid lines
    	! default: 0
    !smooth_D_DSI = 0
    !smooth_D_SH  = 0
    !smooth_D_SSI = 0
    !smooth_D_ES  = 0
    !smooth_D_GSF = 0
    smooth_D_ST  = 2
    smooth_nu_ST = 2

        ! temporal smoothing. Set to 0 to turn it off
    angsmt_D_DSI = 0.0d0
    angsmt_D_SH  = 0.0d0
    angsmt_D_SSI = 0.0d0
    angsmt_D_ES  = 0.0d0
    angsmt_D_GSF = 0.0d0
    angsmt_D_ST  = 0.2d0
    angsmt_nu_ST = 0.2d0
    angsml       = 1d-3

    ! 'hooks' for artificial additional angular momentum transport
    !use_other_am_mixing = .false.
    !x_ctrl(4) = 1.0d5 ! additional viscosity; 1.0d5; 1.0d10
    !x_ctrl(5) = 0.0 ! additional mixing -- don't use for now


      !## resolution
      mesh_delta_coeff = 0.8 ! 0.9
      varcontrol_target = 0.5d-3 !1.0d-4 !0.5d-3 ! max. change between to steps
      min_timestep_limit = 1d-14

    omega_function_weight = 50

      ! until he dep
      ! dX_nuc_drop_limit = 0.1 ! limits max change in abundance mass fraction

!      ! beyond he burning
!      dX_nuc_drop_limit = 0.001
!      dX_nuc_drop_limit_at_high_T = 5.0d-3
!      dX_nuc_drop_min_X_limit = 1d-4

!    delta_lg_XC_cntr_limit = 0.01d0
!    delta_lg_XC_cntr_max   = 0.0d0
!    delta_lg_XC_cntr_min   = -6.0d0             ! default: -5
!    delta_lg_XC_cntr_hard_limit = 0.02d0


!?      delta_HR_limit = 0.0025d0


/ ! end of controls namelist

&pgstar
         
!pause = .true.
         
Profile_Panels1_xmax = 2.2 ! -101d0 ! only used if /= -101d0
Abundance_xmax = 2.2 ! -101 ! only used if /= -101d0

Profile_Panels1_xaxis_name = 'logR_cm'

!Grid2_win_flag = .true.

! file output
!Grid2_file_flag = .true.
   Grid2_file_dir = 'png'
   Grid2_file_prefix = 'grid_'
   Grid2_file_interval = 5 ! 1 ! output when mod(model_number,Grid2_file_interval)==0
   Grid2_file_width = -1 ! negative means use same value as for window
   Grid2_file_aspect_ratio = -1 ! negative means use same value as for window
      

/ ! end of pgstar namelist
