#
# Overwrite quantities before relaxation step
# Usage: ./sedov_3d_generator <this_file>.par
# WARNING: overwrites initial_data_prefix !!
#
  modify_initial_data = yes
  initial_data_prefix = "wd_modified_n32"
  initial_iteration = 1000

  # geometry: spherical domain with radius R = 1
  domain_type = 1   # 0:box, 1:sphere
  sphere_radius = 4.79283e+08

  # equation of state type and parameters
  eos_type = "polytropic"
  poly_gamma = 1.66667   # polytropic index

  # reset thermodynamical quantities
  density_profile = "from file"
  input_density_file = "density_profile_nr8000.dat"
  external_force_type = "spherical density support"
  rho_initial = 5.2e+6
  pressure_initial = 1.56085e+23

  sedov_blast_energy = 0.0
  sedov_blast_radius = 0.1 # whatever

  # good kernel
  sph_kernel = "Wendland C6"
  sph_eta = 1.5