#
# Sperical particle distribution: initial data parameter file
# Usage: ./sedov_[123]d_generator <this_file>.par to generate initial data
#
  initial_data_prefix = "mesa_initial"
  
  # geometry: spherical domain with uniform density
  domain_type = 1   # 0:box, 1:sphere
  
  # pick the radius of the uniform sphere in such a way that total mass
  # is the same mass as the target density distrubution; in this example 
  # the target density profile is "mesa" profile (flat with a smooth rim):
  sphere_radius = 0.922984 # 3D
  # sphere_radius = 0.918559 # 2D
  # sphere_radius = 0.916667 # 1D

  # icosahedra lattice with small perturbations
  lattice_nx =  24  # particle lattice dimension
  lattice_type = 2  # 0:rectangular, 1:hcp, 2:fcc, 3:icosahedral
  lattice_perturbation_amplitude = 0.1 # in units of sm. length
  
  # equation of state type and parameters
  eos_type = "polytropic"
  poly_gamma = 1.4         # polytropic index
  
  # density and pressure for relaxation stage
  rho_initial = 1.0
  pressure_initial = 0.01

  # 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.6