! inlist_ns


&star_job
      show_log_description_at_start = .false.
      
      !show_net_reactions_info = .true.
      
      !create_pre_main_sequence_model = .true.
      !profile_starting_model = .true.

      !save_model_number = 11
      !save_model_filename = '4.0_z2m2_pre_ms.mod'
      !save_model_filename = 'wd.mod'
      !save_model_filename = 'ns_1.4M.mod'
      !save_model_filename = 'ns_1.4M_h_cold.mod'
      save_model_filename = 'ns_1.4M_c.mod'

      load_saved_model = .true.
      !saved_model_name = '../../data/star_data/white_dwarf_models/1.000_Tc_7.6.mod'
      saved_model_name = 'ns_1.4M_c.mod'


      change_net = .true.      
      !new_net_name = 'cno_extras_plus_fe56.net'
      new_net_name = 'approx21.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


      log_columns_file = 'log_columns.list' 
      profile_columns_file = 'profile_columns.list' 


      ! 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

      !remove_center_at_cell_k = 657

      !relax_M_center = .true. 
      new_mass = 1.4
      dlgm_per_step = 1d-3
      relax_M_center_dt = 3.1558149984d-1 ! secyer*1d-8
       
      !relax_R_center = .true.
      new_R_center = 1d6 ! in cm
      dlgR_per_step = 1d-2
      relax_R_center_dt = 3.1558149984d-1 ! secyer*1d-8
         
      !relax_L_center = .true.
      new_L_center = 1e35 ! in ergs/second
      dlgL_per_step = 1d-2
      relax_L_center_dt = 3.1558149984d-1 ! secyer*1d-8
         ! mdot = 1.7d-9 goes with 1e35 L_center
      
      
      !report_mass_not_fe56 = .true.
      !report_cell_for_xm = 1.9892d23
      
      !set_to_xa_for_accretion = .true.
      set_nzlo = 1
      set_nzhi = 99999
      

      set_to_this_tau_factor = 300 ! tau_factor of 1 puts outer cell at photosphere.
      set_tau_factor = .true. ! change tau_factor without reconverging.

      kappa_file_prefix = 'gs98'
      

      !pgstar_flag = .true.



/ ! end of star_job namelist

&controls


      use_lnE_for_eps_grav = .false.


      D_norm_kappa = 1
      D_norm_alpha = 2 
      D_norm_atol = 1d-3
      D_norm_rtol = 1d-3


      
      initial_mass = 1
      initial_z = 0.02d0
      
      photostep = 100
      profile_interval = 100
      log_cnt = 100
      terminal_cnt = 10
      write_header_frequency = 10
      
      min_timestep_limit = 1d-9
      gamma_center_limit = 1d99
      
      varcontrol_target = 1d-3

      power_c_burn_upper_limit = 1d3
      !when_to_stop_rtol = 1d-4
      !when_to_stop_atol = 1d-4
      
      merge_if_dlnR_too_small = .true.
      mesh_min_dlnR = 1d-9
      
      !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

      report_why_dt_limits = .false.
      report_all_dt_limits = .false.
      
      !show_mesh_changes = .true.
      !mesh_dump_call_number = 131
      
      !max_model_number = 21100
      !max_model_number = 11000
      !max_model_number = 18000 ! << increase HR xmax

      !trace_evolve = .true.
      
      ! hydro debugging
       !hydro_numerical_jacobian = .true.
      !hydro_check_everything = .true.
      !hydro_show_rcond = .true.
      !hydro_inspectB_flag = .true.
      !hydro_save_numjac_plot_data = .true.
      !hydro_dump_call_number = 2
      !okay_to_remesh = .false.

      mesh_delta_coeff = 0.5
      
      which_atm_option = 'grey_and_kap'
      
      varcontrol_target = 1d-3 
      
      RGB_wind_scheme = ''
      AGB_wind_scheme = ''
      RGB_to_AGB_wind_switch = 1d-4
      Reimers_wind_eta = 0.5d0  
      Blocker_wind_eta = 0.1d0 
      
      screening_mode = 'extended'
      
      use_GR_factors = .true.
      
      !super_eddington_wind_eta = 1
         
      supersonic_wind_eta = 0
         mass_change_full_on_dt = 1d-4 ! (seconds)
         mass_change_full_off_dt = 1d-6 ! (seconds)

      
      accrete_same_as_surface = .false. 
      accrete_given_mass_fractions = .true. 
      num_accretion_species = 1
      accretion_species_xa(1) = 1
      !accretion_species_id(1) = 'he4'
      accretion_species_id(1) = 'c12'
      !accretion_species_id(1) = 'fe56'

      !accrete_given_mass_fractions = .false. 
      accretion_h1 = 0.72 ! mass fraction
      accretion_h2 = 0 ! if no h2 in current net, then this is automatically added to h1
      accretion_he3 = 0
      accretion_he4 = 0.26
      accretion_zfracs = 3 ! one of the identifiers for different Z fractions from chem_def


      use_CO_enhanced_opacities = .false.

      net_pp_rate_factor = 0
      net_cno_rate_factor = 0

      mass_change = 1.7d-9  ! rate of accretion (Msun/year)
         ! mdot = 1.7d-11 goes with 1e33 L_center
         ! mdot = 1.7d-13 goes with 1e31 L_center
         ! mdot = 5.1d-13 goes with 3e31 L_center
      
      !max_years_for_timestep = 1d-2

      ! when reach target mass, turn these off
      	!max_years_for_timestep = 1d-8

      max_abar_for_burning = 55



/ ! end of controls namelist






&pgstar
         
      ! top level controls

         !pause = .true. 
            ! if true, the code waits for user to enter a RETURN on the command line
      
      ! main window
      
         MAIN_win_flag = .true.
         
         show_TRho_Profile_with_main = .true.
      
         show_HR_TRho_with_main = .false. ! if false, use aspect ratio > 1
         MAIN_win_width = 7.5
         MAIN_win_aspect_ratio = 1.2 ! aspect_ratio = height/width
      
         xaxis_by = 'by_log_column_depth' ! select xaxis for main window

         show_main_win_text_info = .false.
                     
         ! xaxis limits -- to override system default selections
         main_xmin = 2 !-18 !0.4
         main_xmax = 10.2
         
         logxq_cutoff = -9.1 ! min value when using logxm for xaxis
         !main_logL_max = 7
         !main_logL_min = 0
         
         ! control for plot showing abundances
         log_mass_frac_ymax = 0.7 ! making this > 0 helps readability
         log_mass_frac_ymin = -4.01
         num_abundance_line_labels = 5
      
         log_abund_vary = 0.1 ! don't show species if it would vary by < this in plot
      
         ! file output
         !MAIN_file_flag = .true.
         MAIN_file_dir = 'png'
         MAIN_file_prefix = 'main'
         MAIN_file_cnt = 2 ! output when mod(model_number,main_file_cnt)==0
               
               
      ! TRho Profile window -- current model in T-Rho plane
      
         TRho_Profile_win_flag = .false.
         
         TRho_switch_to_Column_Depth = .true.

         TRho_Profile_win_width = 7.5
         TRho_Profile_win_aspect_ratio = 0.618 ! aspect_ratio = height/width
            
         show_TRho_Profile_text_info = .false.
         show_TRho_Profile_legend = .true.
         show_HR_TRho_with_TRho_Profile = .false.
         show_TRho_Profile_burn_labels = .false.
         TRho_Profile_show_logQ_limit = .true.
         !show_TRho_Profile_cross_hair = .true.
         show_TRho_Profile_mass_locs = .false.

         TRho_Profile_legend_disp1 = -2.0
         TRho_Profile_legend_del_disp = -1.3

         show_TRho_Profile_kap_regions = .false.
         show_TRho_Profile_eos_regions = .false.
         show_TRho_Profile_degeneracy_line = .true.
         show_TRho_Profile_Pgas_Prad_line = .true.
         show_TRho_Profile_burn_lines = .true.

         !show_TRho_Profile_annotation1 = .true.
         !show_TRho_Profile_annotation2 = .true.
         !show_TRho_Profile_annotation3 = .true.
      
         ! axis limits
         TRho_Profile_xmin = -3
         TRho_Profile_xmax = 11
         TRho_Profile_ymin = 6
         TRho_Profile_ymax = 10

         TRho_Profile_xmin = 2
         TRho_Profile_xmax = 10.5
         TRho_Profile_ymin = 7
         TRho_Profile_ymax = 9.4
         
         ! file output
         !TRho_Profile_file_flag = .true.
         TRho_Profile_file_dir = 'png'
         TRho_Profile_file_prefix = 'trho_profile'
         TRho_Profile_file_cnt = 1 ! output when mod(model_number,TRho_Profile_file_cnt)==0
      
      
      ! Profile window
      
         !Profile_win_flag = .true.

         Profile_win_width = 7.5
         Profile_win_aspect_ratio = 0.618 ! aspect_ratio = height/width
            
         show_Profile_legend = .false.
         Profile_legend_coord = 0.65
         Profile_legend_fjust = 0.0
         Profile_legend_disp1 = -7.5
         Profile_legend_del_disp = -1.5
         Profile_show_mass_locs = .false.
         Profile_show_decorated_line = .false.
         
         show_HR_TRho_with_Profile = .false.
         show_Profile_text_info = .false.
         Profile_text_info_xfac = 0.78 ! controls x location
         Profile_text_info_dxfac = 0.02 ! controls x spacing to value from text
         Profile_text_info_yfac = 0.95 ! controls y location of 1st line
         Profile_text_info_dyfac = -0.04 ! controls line spacing
         show_Profile_cross_hair = .false.

         show_Profile_annotation1 = .false.
         show_Profile_annotation2 = .false.
         show_Profile_annotation3 = .false.
         
         ! axis choices
            ! to get a list of axis identifiers, edit your star_job controls
               ! to set show_profile_column_numbers = .true.
            ! you can use any identifier for x or y
         Profile_xaxis_name = 'log_column_depth'
         Profile_xaxis_reversed = .false.
         Profile_xmin = 2 ! only used if > -100
         Profile_xmax = 10 !-101 ! only used if > -100
         
         Profile_yaxis_name = 'logT' !'grada'
         Profile_yaxis_reversed = .false.
         Profile_ymin = 6.8 ! only used if > -100
         Profile_ymax = 9.8 !8.6 ! only used if > -100    
         Profile_dymin = 0.025    
         
         Profile_other_yaxis_name = 'mu'
         Profile_other_yaxis_reversed = .false.
         Profile_other_ymin = 1.25 ! only used if > -100
         Profile_other_ymax = 1.9 ! only used if > -100        
         
         ! file output
         !Profile_file_flag = .true.
         Profile_file_dir = 'png_profile'
         Profile_file_prefix = 'profile'
         Profile_file_cnt = 2 ! output when mod(model_number,Profile_file_cnt)==0
      
      
      ! Convection window -- history of convection and more.
      
         CONV_win_flag = .false.
      
         CONV_win_width = 9
         CONV_win_aspect_ratio = 0.62 ! aspect_ratio = height/width
      
         ! axis limits
         CONV_xmax = -1 ! step number.  negative means use default.
         CONV_xmin = -1 ! step number.  negative means use default.
         CONV_max_width = 0 ! only used if > 0.  causes xmin to move with xmax.
         CONV_mmax = 1.00007  !-1 ! (Msun units) negative means use start initial mass
         CONV_mmin = 1  !0 ! (Msun units)
         CONV_lgLmax = -101 ! only used if > -100; (L in Lsun units)
         CONV_lgLmin = -2 ! only used if > -100; (L in Lsun units)

         CONV_show_burn = .true.
         CONV_show_mixing = .true.
         CONV_show_log_radius = .true.
         CONV_show_luminosities = .true.
         CONV_show_mass_boundaries = .true.

         show_CONV_annotation1 = .false.
         show_CONV_annotation2 = .false.
         show_CONV_annotation3 = .false.
         
         ! file output
         CONV_file_flag = .false.
         CONV_file_dir = 'pgstar_out'
         CONV_file_prefix = 'conv'
         CONV_file_cnt = 5 ! output when mod(model_number,CONV_file_cnt)==0
         CONV_file_width = -1 ! negative means use same value as for window
         CONV_file_aspect_ratio = -1 ! negative means use same value as for window
      
      
         
         
      ! Surface History window
      
         Surf_Hist_win_flag = .false.
      
         Surf_Hist_win_width = 6.5
         Surf_Hist_win_aspect_ratio = 1.2 ! aspect_ratio = height/width
      
         ! axis limits
         Surf_Hist_xmax = -1 ! step number.  negative means use default.
         Surf_Hist_xmin = -1 ! step number.  negative means use default.
         Surf_Hist_max_width = 0 ! only used if > 0.  causes xmin to move with xmax.
         
         ! abundance limits
         Surf_Hist_mass_frac_min = 0
         Surf_Hist_mass_frac_max = 1.1 ! > 0 helps legibility

         show_Surf_Hist_annotation1 = .false.
         show_Surf_Hist_annotation2 = .false.
         show_Surf_Hist_annotation3 = .false.
         
         ! file output
         Surf_Hist_file_flag = .false.
         Surf_Hist_file_dir = 'pgstar_out'
         Surf_Hist_file_prefix = 'conv'
         Surf_Hist_file_cnt = 5 ! output when mod(model_number,Surf_Hist_file_cnt)==0
         Surf_Hist_file_width = -1 ! negative means use same value as for window
         Surf_Hist_file_aspect_ratio = -1 ! negative means use same value as for window
      
         
         
      ! Power window
      
         Power_win_flag = .false.
      
         Power_win_width = 8.5
         Power_win_aspect_ratio = 0.618 ! aspect_ratio = height/width
      
         Power_xaxis_by = 6 ! select xaxis
            ! 1 = by_mass
            ! 2 = by_grid
            ! 3 = by_radius
            ! 4 = by_logR
            ! 5 = by_logP
            ! 6 = by_logxm
         
         ! power xaxis limits -- to override system default selections
         Power_xmin = -101 ! only used if > -100
         Power_xmax = -101 ! only used if > -100
         
         ! power yaxis limits -- to override system default selections
         Power_ymin = -101 ! only used if > -100
         Power_ymax = -101 ! only used if > -100

         show_Power_cross_hair = .false.
         
         ! file output
         Power_file_flag = .false.
         Power_file_dir = 'pgstar_out'
         Power_file_prefix = 'power'
         Power_file_cnt = 5 ! output when mod(model_number,Power_file_cnt)==0
         Power_file_width = -1 ! negative means use same value as for window
         Power_file_aspect_ratio = -1 ! negative means use same value as for window
         
         
      ! Abundance window
      
         Abundance_win_flag = .false.
      
         Abundance_win_width = 8.5
         Abundance_win_aspect_ratio = 0.618 ! aspect_ratio = height/width
      
         Abundance_xaxis_by = 6 ! select xaxis
            ! 1 = by_mass
            ! 2 = by_grid
            ! 3 = by_radius
            ! 4 = by_logR
            ! 5 = by_logP
            ! 6 = by_logxm
         
         ! power xaxis limits -- to override system default selections
         Abundance_xmin = -101 ! only used if > -100
         Abundance_xmax = -101 ! only used if > -100
         
         ! power yaxis limits -- to override system default selections
         Abundance_log_mass_frac_min = 1 ! only used if < 0
         Abundance_log_mass_frac_max = 1 ! only used if < 0

         show_Abundance_cross_hair = .false.
         
         ! file output
         Abundance_file_flag = .false.
         Abundance_file_dir = 'pgstar_out'
         Abundance_file_prefix = 'abund'
         Abundance_file_cnt = 5 ! output when mod(model_number,Abundance_file_cnt)==0
         Abundance_file_width = -1 ! negative means use same value as for window
         Abundance_file_aspect_ratio = -1 ! negative means use same value as for window
         
         


/ ! end of pgstar namelist
