From 3ce67a02f07dea1eed9b5a4406acf798edb5fd53 Mon Sep 17 00:00:00 2001 From: Ian Hinder Date: Thu, 30 Sep 2010 01:01:33 +0100 Subject: Examples: Add hydro parameter files --- Examples/euler_shock.par | 141 +++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 141 insertions(+) create mode 100644 Examples/euler_shock.par (limited to 'Examples/euler_shock.par') diff --git a/Examples/euler_shock.par b/Examples/euler_shock.par new file mode 100644 index 0000000..a452db2 --- /dev/null +++ b/Examples/euler_shock.par @@ -0,0 +1,141 @@ +#!/usr/bin/perl -W + +ActiveThorns = " +Boundary +Carpet +CarpetIOASCII +CarpetIOBasic +CarpetIOScalar +CarpetLib +CarpetReduce +CarpetSlab +CartGrid3d +CoordBase +Dissipation +GenericFD +IOUtil +LoopControl +MoL +NanChecker +Periodic +Euler +Slab +SpaceMask +SphericalSurface +SymBase +Time +" + +############################################################# +# Grid +############################################################# + +CoordBase::domainsize = minmax + +CoordBase::xmin = -1 +CoordBase::ymin = 0 +CoordBase::zmin = 0 + +CoordBase::xmax = 2 +CoordBase::ymax = 0.1 +CoordBase::zmax = 0.1 + +CoordBase::dx = 0.04 +CoordBase::dy = 0.1 +CoordBase::dz = 0.1 + +CoordBase::boundary_size_x_lower = 2 +CoordBase::boundary_size_y_lower = 2 +CoordBase::boundary_size_z_lower = 2 +CoordBase::boundary_shiftout_x_lower = 1 +CoordBase::boundary_shiftout_y_lower = 1 +CoordBase::boundary_shiftout_z_lower = 1 + +CoordBase::boundary_size_x_upper = 2 +CoordBase::boundary_size_y_upper = 2 +CoordBase::boundary_size_z_upper = 2 +CoordBase::boundary_shiftout_x_upper = 0 +CoordBase::boundary_shiftout_y_upper = 0 +CoordBase::boundary_shiftout_z_upper = 0 + +CartGrid3D::type = "coordbase" +CartGrid3D::domain = "full" +CartGrid3D::avoid_origin = "no" + +Periodic::periodic = "yes" + +############################################################# +# Carpet +############################################################# + +Carpet::ghost_size = 2 +Carpet::domain_from_coordbase = "yes" +Carpet::max_refinement_levels = 1 +#Carpet::init_each_timelevel = "yes" +Carpet::num_integrator_substeps = 4 + +############################################################# +# Time integration +############################################################# + +Cactus::terminate = "time" +Cactus::cctk_final_time = 0.2 + +Time::dtfac = 0.25 +MethodOfLines::ode_method = "RK4" +MethodOfLines::MoL_Intermediate_Steps = 4 +MethodOfLines::MoL_Num_Scratch_Levels = 1 +MethodOfLines::MoL_NaN_Check = "yes" + +# MethodOfLines::ode_method = "generic" +# MethodOfLines::generic_type = "RK" +# MethodOfLines::MoL_Intermediate_Steps = 1 + +Dissipation::vars = "" +Dissipation::order = 1 +Dissipation::epsdis = 0.5 + +Euler::rhoL0 = 1 +Euler::rhoR0 = 0.125 +Euler::vL0 = 0 +Euler::vR0 = 0 +Euler::pL0 = 1 +Euler::pR0 = 0.1 +Euler::alpha = 1 +Euler::gamma = 1.4 +Euler::initial_data = "shock" + +############################################################# +# Boundary conditions +############################################################# + +Euler::Den_group_bound = "none" +Euler::S_group_bound = "none" +Euler::En_group_bound = "none" + +Carpet::poison_new_timelevels = yes +Carpet::poison_value = 100 +CarpetLib::poison_new_memory = yes +CarpetLib::poison_value = 200 +Carpet::check_for_poison = no + +############################################################# +# Output +############################################################# + +IO::out_dir = $parfile +IO::out_fileinfo = "all" + +CarpetIOBasic::outInfo_every = 1 +CarpetIOBasic::outInfo_vars = "Euler::rho Euler::DenF Euler::v1 Euler::p Euler::Den Euler::S1 Euler::En" + +CarpetIOScalar::outScalar_every = 1 +CarpetIOScalar::outScalar_vars = "" +CarpetIOScalar::outScalar_reductions = "minimum maximum norm2" + +IOASCII::out1D_every = 1 +IOASCII::out1D_x = "yes" +IOASCII::out1D_vars = "Euler::rho Euler::v1 Euler::p Euler::Den Euler::S1 Euler::En" + +CarpetIOASCII::out_precision = 19 +CarpetIOASCII::out3D_ghosts = "yes" -- cgit v1.2.3