
! inlist_co_nova


&star_job

      show_log_description_at_start = .false.
      
      load_saved_model = .true.
      saved_model_name = '1.30Msun_Tc_3e7_init.mod' 

      save_model_when_terminate = .true.
      save_model_filename = 'final.mod'
      
      write_profile_when_terminate = .true.
      filename_for_profile_when_terminate = 'final_profile.data'

      change_initial_net = .true.      
      new_net_name = 'cno_extras_o18_to_mg26_plus_fe56.net'
      auto_extend_net = .false.

      change_v_flag = .true.
      new_v_flag = .true.

      set_initial_age = .true.
      initial_age = 0 ! in years

      set_initial_model_number = .true.
      initial_model_number = 0

      set_tau_factor = .true.
      set_to_this_tau_factor = 30
      
!      set_initial_dt = .true.
!      years_for_initial_dt = 1d-4

       !pgstar_flag = .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 = 20
         max_number_retries = 100
         max_model_number = 2000

      initial_mass = 1.30
      initial_z = 0.02d0
      
      use_Type2_opacities = .true.
      Zbase = 0.02d0

      mass_change_full_on_dt = 1d-6 ! (seconds)
      mass_change_full_off_dt = 1d-6 ! (seconds)

      ! NOTE: for super eddington wind,
      ! we use Ledd averaged over mass to optical depth tau = 100
      super_eddington_wind_eta = 1
         ! parameter for mass loss driven by super Eddington luminosity        
      super_eddington_wind_Ledd_factor = 1
         ! multiply Ledd by this factor when computing super Eddington wind
         ! e.g., if this is 2, then only get wind when L > 2*Ledd


      mixing_length_alpha = 1.9
      MLT_option = 'Henyey'

      use_Ledoux_criterion = .true.
      alpha_semiconvection = 1d-3
      thermohaline_coeff = 1d3
      
         

      logQ_limit = 100
      gamma_center_limit = 1000 ! FOR THIS TEST
            
      star_mass_max_limit = 1.301d0
   
      mass_change = 1d-9

      photostep = 50
      profile_interval = 50
      history_interval = 5
      terminal_cnt = 5
      write_header_frequency = 10
      
      min_timestep_limit = 1d-12

                  
      RGB_wind_scheme = ''
      AGB_wind_scheme = ''

      xa_function_species(1) = 'h1'  ! name of nuclide as defined in chem_def
      xa_function_weight(1) = 20
      xa_function_param(1) = 1d-6
      xa_function_species(2) = 'he4'  ! name of nuclide as defined in chem_def
      xa_function_weight(2) = 20
      xa_function_param(2) = 1d-2

!      delta_lgL_H_limit = 0.01 ! limit for magnitude of change in lgL_H
      delta_lgL_H_limit = 0.1 ! limit for magnitude of change in lgL_H
      lgL_H_burn_min = -0.5 ! ignore changes in lgL_H if value is less than this

!      delta_lgL_He_limit = 0.025 ! limit for magnitude of change in lgL_He
      delta_lgL_He_limit = 0.25 ! limit for magnitude of change in lgL_He
      lgL_He_burn_min = 2.5 ! ignore changes in lgL_He if value is less than this
      
      delta_lgRho_limit = 1 ! limit for magnitude of max change in log10 density at any point
      delta_lgRho_hard_limit = -1

      delta_lgT_limit = 0.5
      delta_lgL_limit = 0.05
      delta_lgTeff_limit = 0.10 



      accrete_same_as_surface = .false. 
      accretion_h1 = 0.7000 ! Estimate Nomoto 2007
      accretion_h2 = 0
      accretion_he3 = 0.292910D-04
      accretion_he4 = 0.279971 ! Estimated Nomoto 2007
      accretion_zfracs = 4 ! Lodders 03

      which_atm_option = 'grey_and_kap'
      Pextra_factor = -1

      varcontrol_target = 3d-4

      mesh_delta_coeff = 1.5
      max_allowed_nz = 20000
      
      tol_max_correction = 3d-2

      ! 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
      
      !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

/ ! end of controls namelist



&pgstar
         Grid8_win_flag = .true.
         Grid8_win_width = 7
         Summary_Burn_xaxis_name = 'logxq' 
         Summary_Burn_xaxis_reversed = .true.
         Summary_Burn_xmin = -14 ! -101d0 ! only used if /= -101d0
         Summary_Burn_xmax = -3 ! -101d0 ! only used if /= -101d0
         Abundance_xaxis_name = 'logxq' 
         Abundance_xaxis_reversed = .true.
         Abundance_xmin = -14 ! -101d0 ! only used if /= -101d0
         Abundance_xmax = -3 ! -101d0 ! only used if /= -101d0
/ ! end of pgstar namelist

