# Temperature difference between hot and cold model
dtemp                   real              3.81e8_rt                 y

# Theta at which the blending of the hot and cold model begins
theta_half_max          real              1.745e-2_rt               y

# Controls theta width at which the blending stops.
# Blending stops at theta = theta_half_max + theta_half_width
theta_half_width        real              4.9e-3_rt                 y

# Highest temperature achieved during the tanh transition.
T_hi                    real              5.e8_rt                   y

# Isothermal temperature of the base region, i.e. neutron star temp
T_star                  real              1.e8_rt                   y

# Temperature of the isothermal fluff region on top of the atmosphere
T_lo                    real              5.e7_rt                   y

# Density of the isothermal fluff region on top of the atmosphere
low_density_cutoff      real              1.e-4_rt                  y

# Density of the base layer, at problo(0) + H_star.
# This is used as the starting point for HSE integration.
dens_base               real              2.e6_rt                   y

# height above which we start tanh transition to blend T_star and T_hi
# This is relative to problo(0).
H_star                  real              500.e0_rt                 y

# Characteristic height of the tanh transition blending layer
# between T_hi and T_star. Transition stops at H_star + 3 atm_delta
atm_delta               real              25.e0_rt                  y

# Initialization that defines the composition of the accretion layer
fuel1_name              string            "helium-4"                y

fuel2_name              string            ""                        y

fuel3_name              string            ""                        y

fuel4_name              string            ""                        y

ash1_name               string            "iron-56"                 y

ash2_name               string            ""                        y

ash3_name               string            ""                        y

fuel1_frac              real              1.0_rt                    y

fuel2_frac              real              0.0_rt                    y

fuel3_frac              real              0.0_rt                    y

fuel4_frac              real              0.0_rt                    y

ash1_frac               real              1.0_rt                    y

ash2_frac               real              0.0_rt                    y

ash3_frac               real              0.0_rt                    y

smallx                  real              1.e-10_rt                 y

### Parameters for controlling grid refinement ###

# do we dynamically refine based on density and temp? or based on height?
tag_by_thermodynamics   integer           1                         y

# If tag_by_thermodynamic=0, then Tag everything below this height.
refine_height           real              3600                      y

# used for tagging if tag_by_thermodynamic = 1.
# Tag when density and temperature are both above the cutoff
cutoff_density          real              500.e0_rt                 y
cutoff_temp             real              1.e8_rt                   y

## Allow different tagging levels for tag_by_thermodynamic=1 ##

# Only cutoff_density and cutoff_temp need to be satisfied for this level.
max_base_tagging_level  integer           1                         y

# Allow for additional refinement level to better resolve initial
# burning phase. Aside from cutoff_temp and cutoff_dens,
# also require the region to be within a certain cylindrical-r distance,
# this is controlled by refine_distance and helium-4 mass fraction to be
# greater than refine_X_min.
max_hse_tagging_level   integer           2                         y

# Cylindrical-r distance below which are refined.
# See description for max_hse_tagging_level
refine_distance         real              1.e30_rt                  y

# cutoff helium-4 mass fraction to allow for additional refinement.
X_min                   real              1.e-4_rt                  y
