#
# Sperical particle distribution: initial data parameter file
# Usage: ./sedov_3d_generator <this_file>.par to generate initial data
#
  initial_data_prefix = "wd_initial_n32"
  
  # geometry: spherical domain with radius R = 1
  domain_type = 1   # 0:box, 1:sphere
  density_profile = "from file"
  input_density_file = "density_profile_nr8000.dat"
  sphere_radius = 4.79283e+08

  # icosahedra lattice with small perturbations
  lattice_nx =  32  # particle lattice dimension
  lattice_type = 3  # 0:rectangular, 1:hcp, 2:fcc, 3:icosahedral
  lattice_perturbation_amplitude = 0.09 # in units of sm. length
  
  # equation of state type and parameters
  eos_type = "polytropic"
  #poly_gamma = 1.66667   # polytropic index
  poly_gamma = 0.99
  
  # density and pressure for relaxation stage
  rho_initial = 5.2e+6
  #pressure_initial = 1.56085e+23
  pressure_initial = 1e+28

  # since we only need spherical distribution of particles,
  # set Sedov energy to zero
  sedov_blast_energy = 0.0 
  
  # use a good kernel
  sph_kernel = "Wendland C6"
  sph_eta = 1.5