From 3c0323ec4a6d0cceefd1a2c1dcfd45144fea1f21 Mon Sep 17 00:00:00 2001 From: rhaas Date: Thu, 9 Aug 2012 06:10:15 +0000 Subject: Introduce a better treatment of the shocktube boundaries. * The 2d diagonal boundary conditions are correctly applied when using multiprocessors. * The ymin (ymax) surfaces have their points in the region indsum < minsum (indsum > maxsum) fixed after a first sync operation. A second sync guarantees the neighbour grid receives the correct value of the boundary on those regions. * Mirror symmetries applied at the faces intercepting the upper right corner assumed cubic local domain. However, rectangular local domains are common after domain decomposition. This commit fixes the indexes applying this symmetry for these cases. * Upper right corner/edge/surface condition should use sten instead of stenp1. Thanks for catching this Josh! * Attempt at a 3d diagonal boundary treatment (Josh's commit) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@135 ac85fae7-cede-4708-beff-ae01c7fa1c26 --- par/balsara1d.par | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) (limited to 'par') diff --git a/par/balsara1d.par b/par/balsara1d.par index 3d1980f..f688b16 100644 --- a/par/balsara1d.par +++ b/par/balsara1d.par @@ -81,7 +81,7 @@ Cactus::terminate="iteration" #Cactus::cctk_final_time = 0.02 #Cactus::cctk_final_time = 0.4 #cactus::cctk_itlast = 40 -cactus::cctk_itlast = 400 +cactus::cctk_itlast = 800 #cactus::cctk_itlast = 2 IO::out_dir = $parfile @@ -97,12 +97,12 @@ CarpetIOBasic::outInfo_every=1 ##CarpetIOASCII::out1D_vars = "HydroBase::rho HydroBase::press HydroBase::eps HydroBase::vel grhydro::psidc grhydro::dens grhydro::tau grhydro::scon HydroBase::Bvec" ##CarpetIOASCII::out1D_vars = "HydroBase::rho HydroBase::press HydroBase::eps HydroBase::vel grhydro::dens grhydro::tau grhydro::scon " -#CarpetIOHDF5::out_every = 1 -#CarpetIOHDF5::out_vars = "HydroBase::rho HydroBase::press HydroBase::eps HydroBase::vel HydroBase::w_lorentz grhydro::dens grhydro::divB grhydro::tau grhydro::scon HydroBase::Bvec" +CarpetIOHDF5::out_every = 10 +CarpetIOHDF5::out_vars = "HydroBase::rho HydroBase::press HydroBase::eps HydroBase::vel HydroBase::w_lorentz grhydro::dens grhydro::divB grhydro::tau grhydro::scon HydroBase::Bvec" ##CarpetIOHDF5::out_vars = "HydroBase::rho HydroBase::press HydroBase::eps HydroBase::vel HydroBase::w_lorentz grhydro::psidc grhydro::dens grhydro::divB grhydro::tau grhydro::scon HydroBase::Bvec" ##CarpetIOHDF5::out_vars = "HydroBase::rho HydroBase::press HydroBase::eps HydroBase::vel HydroBase::w_lorentz grhydro::dens grhydro::tau grhydro::scon " -CarpetIOHDF5::out2D_every = 1 +CarpetIOHDF5::out2D_every = 10 CarpetIOHDF5::out2D_xy = "yes" CarpetIOHDF5::out2D_xz = "no" CarpetIOHDF5::out2D_yz = "no" -- cgit v1.2.3