
! inlist_wd_ignition


&star_job
      show_log_description_at_start = .false.
      
      load_saved_model = .true.
      saved_model_name = 'wd_ignite.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 = 'co_burn.net'
      
      set_rate_3a = 'FL87' ! Fushiki and Lamb, Apj, 317, 368-388, 1987
      set_rate_1212 = 'G05' ! Gasques, et al. Phys Review C, 72, 025806, 2005

      show_log_description_at_start = .false. 

      ! shut off PC since get bad convergence for high gamma
      set_eos_PC_parameters = .true.
      mass_fraction_limit_for_PC = 1d-2 ! skip species if abundance < this
      logRho1_PC_limit = 20 ! 3.7d0 ! okay for pure PC for logRho > this
      logRho2_PC_limit = 20 ! 2.8d0 ! don't use PC for logRho < this (>= 2.8 or so because of PPT)
      log_Gamma_all_HELM = 1.60206d0 ! HELM for log_Gamma <= this
      log_Gamma_all_PC = 1.90309d0 ! PC for log_Gamma >= this

      change_v_flag = .true.
      new_v_flag = .true.

      set_Z_all_HELM = .true. ! if false, just use the eos default
      Z_all_HELM = 1.000001d0


      ! move the surface down to larger optical depth
      set_tau_factor = .true.
      set_to_this_tau_factor = 300 


      !pgstar_flag = .true.
      
/ ! end of star_job namelist


&controls
      
         
      num_trace_history_values = 1
      trace_history_value_name(1) = 'split_mixing_choice'

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

      max_model_number = 2000

      initial_mass =  0.45
      initial_z = 0.02d0
      
      use_Type2_opacities = .true.
      Zbase = 0.02d0
      
      power_nuc_burn_upper_limit = 1d8 
      
      log_center_density_limit = 15
      
      varcontrol_target = 1d-3
      mesh_delta_coeff = 0.5

      logQ_limit = 1d6
      gamma_center_limit = 1d6

      photostep = 50
      profile_interval = 100
      history_interval = 10
      terminal_cnt = 10
      write_header_frequency = 10
      
      min_timestep_limit = 1d-4 ! seconds

      delta_lgL_nuc_limit = 0.05 ! limit for magnitude of change in lgL_nuc

      T_mix_limit = 1d4
      
      which_atm_option = 'grey_and_kap'

      accrete_same_as_surface = .false. 
      accrete_given_mass_fractions = .true. 
      num_accretion_species = 2
      accretion_species_id(1) = 'c12'
      accretion_species_xa(1) = 0.25
      accretion_species_id(2) = 'o16'
      accretion_species_xa(2) = 0.75
      
      mass_change = 1d-9 ! rate of accretion (Msun/year)


      ! 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_matrix_input = .true. ! matrix before factor
      !trace_newton_bcyclic_matrix_output = .true. ! matrix after factor
      
      trace_newton_bcyclic_steplo = 905 ! 1st model number to trace
      trace_newton_bcyclic_stephi = 905 ! last model number to trace
      
      trace_newton_bcyclic_iterlo = 5 ! 1st newton iter to trace
      trace_newton_bcyclic_iterhi = 5 ! 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 = 1852

/ ! end of controls namelist



&pgstar

!pause = .true.

         TRho_Profile_win_flag = .true.
         TRho_Profile_xmin = -2.5
         TRho_Profile_xmax = 11.0
         TRho_Profile_ymin = 5.1
         TRho_Profile_ymax = 10.0        
         
         !History_Panels1_win_flag = .true.
         History_Panels1_yaxis_name(3) = 'log_R' 
         History_Panels1_other_yaxis_name(3) = 'v_div_csound_surf' 

/ ! end of pgstar namelist
