create envelope of about 10^-11 Msun pure Fe56
for 1.4 Msun neutron star with radius 10 km.
L from core = 1e33 erg/sec




edit inlist to specify inlist_ns1.
run until it saves ns1.mod and stops.

then edit inlist to specify inlist_ns2.
run until it saves ns2.mod and stops.

keep going until you've done all of the inlist_ns's.



1 change net and tau factor

      change_net = .true.      
      new_net_name = 'cno_extras_plus_fe56.net'

      set_tau_factor = .true. ! change tau_factor without reconverging.
      set_to_this_tau_factor = 75

2 remove big part of inner mass

      remove_inner_fraction_q = 0.9

3 remove more mass

      remove_inner_fraction_q = 0.99

4 remove more and convert all to fe56
      
      remove_inner_fraction_q = 0.99
      
      set_to_xa_for_accretion = .true.
      set_nzlo = 1
      set_nzhi = 99999

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

5 remove more mass
      remove_inner_fraction_q = 0.999
      
6 remove more mass
      remove_inner_fraction_q = 0.999
      
7 change M_center
      relax_M_center = .true. 
      new_mass = 1.4
      dlgm_per_step = 1d-3
      relax_M_center_dt = 3.1558149984d-1 ! secyer*1d-8
      
8 change R_center
      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

9 change L_center
      relax_L_center = .true.
      new_L_center = 1e34 ! in ergs/second
      dlgL_per_step = 1d-2
      relax_L_center_dt = 3.1558149984d-1 ! secyer*1d-8


