################################################################
# headon NS (conformal initial data)
#################################################################

ActiveThorns = "Mahc_Evolve Mahc_Init_Data conf_adm Time einstein pugh pughslab
  pughreduce cartgrid3d ioascii ioutil iobasic Mahc_Boundary ADM_Constraints "

einstein::evolution_system = "Conf_ADM"     
einstein::initial_data     = "boost_multiple_star"
einstein::einstein_register_slicing = "no"
einstein::slicing         = "ch_oneplog"
einstein::use_conformal = "no"
einstein::use_conformal_derivs = "no"
einstein::shift = "zero"

conf_adm::conf_adm_method = "icn"
Mahc_Evolve::mahc_spacetime_icn = "yes"

Mahc_Evolve::mahc_evolution = "yes"
Mahc_Evolve::hydro_order = "fluxlim"
Mahc_Evolve::mahc_stencil_size = 2
Mahc_Evolve::eos_gamma = 2.0
Mahc_Evolve::eos_k = 80.0
Mahc_Evolve::bound = "static"
Mahc_Evolve::centden = 2.43e-3

Mahc_Init_Data::boost_multi_star_num = 2
Mahc_Init_Data::bms_centden_1 = 2.43e-3
Mahc_Init_Data::bms_centden_2 = 2.43e-3
Mahc_Init_Data::bms_xpos_1 = 0.0
Mahc_Init_Data::bms_ypos_1 = 0.0
Mahc_Init_Data::bms_zpos_1 = -9.48
Mahc_Init_Data::bms_xpos_2 = 0.0
Mahc_Init_Data::bms_ypos_2 = 0.0
Mahc_Init_Data::bms_zpos_2 = 9.48
Mahc_Init_Data::bms_boost_x_1 = 0.0
Mahc_Init_Data::bms_boost_y_1 = 0.0
Mahc_Init_Data::bms_boost_z_1 = 0.2
Mahc_Init_Data::bms_boost_x_2 = 0.0
Mahc_Init_Data::bms_boost_y_2 = 0.0
Mahc_Init_Data::bms_boost_z_2 = -0.2

Mahc_Boundary::mahc_bound = "robin"
Mahc_Boundary::mahc_bnd_stencil_size = 2


driver::global_nx = 30
driver::global_ny = 30
driver::global_nz = 30
driver::ghost_size = 2

time::dtfac = 0.3
grid::type        = "byspacing"
grid::dxyz        = 2.0
grid::domain      = "full"

cactus::cctk_itlast          =  10

IO::out_fileinfo = "none"
IO::outdir = "test_headonns_conformal"
IOBasic::outInfo_every  =  1        
IOBasic::outInfo_vars       =  "einstein::gxx gen_rel_perfect_fluid::rho "

IOASCII::out1D_every =  5
IOASCII::out1D_vars  = "einstein::gxx einstein::kxx einstein::alp
   einstein::gzz einstein::gyy einstein::gxy einstein::gxz einstein::gyz
   einstein::kzz einstein::kyy einstein::kxy einstein::kxz einstein::kyz
                        gen_rel_perfect_fluid::dens
   gen_rel_perfect_fluid::tau gen_rel_perfect_fluid::sx
   gen_rel_perfect_fluid::sy  gen_rel_perfect_fluid::sz
   gen_rel_perfect_fluid::Mahc_prim_vars 
   conf_adm::ch_axx conf_adm::ch_axy conf_adm::ch_axz
   conf_adm::ch_ayy conf_adm::ch_ayz conf_adm::ch_azz
   conf_adm::ch_gxx conf_adm::ch_gxy conf_adm::ch_gxz
   conf_adm::ch_gyy conf_adm::ch_gyz conf_adm::ch_gzz
   conf_adm::ch_gammax conf_adm::ch_gammay conf_adm::ch_gammaz
   conf_adm::ch_trk conf_adm::ch_psi 
   adm_constraints::ADMconstraints "

IOBasic::outScalar_every = 1
IOBasic::outScalar_vars = "einstein::gzz conf_adm::ch_gzz
   einstein::kzz conf_adm::ch_azz 
   gen_rel_perfect_fluid::rho "