Add routines to copy the three-metric into other variables, so that

more than one timelevels can be accessed.


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/CalcK/trunk@4 f47d718b-0e4f-0410-8445-f2f96c8ccefb

3 files changed, 141 insertions, 2 deletions

diff --git a/par/calck.par b/par/calck.par
index fc2446a..bb3934f 100644
--- a/par/calck.par
+++ b/par/calck.par
@@ -22,13 +22,15 @@ IO::filereader_ID_files = "
 	gxx_3d gxy_3d gxz_3d gyy_3d gyz_3d gzz_3d
 	alp_3d
 	betax_3d betay_3d betaz_3d
+        gxx_prev_3d gxy_prev_3d gxz_prev_3d gyy_prev_3d gyz_prev_3d gzz_prev_3d
+        gxx_next_3d gxy_next_3d gxz_next_3d gyy_next_3d gyz_next_3d gzz_next_3d
 "
 IO::filereader_ID_vars = "
 	ADMBase::lapse[cctk_iteration=101]
 	ADMBase::shift[cctk_iteration=101]
 	ADMBase::metric[cctk_iteration=101]
-	CalcK::metric_prev[cctk_iteration=100]
-	CalcK::metric_next[cctk_iteration=102]
+	CalcK::metric_prev
+	CalcK::metric_next
 "
 
 ActiveThorns = "IOHDF5 IOHDF5Util"
@@ -44,6 +46,7 @@ ActiveThorns = "ADMBase ADMCoupling ADMMacros CoordGauge SpaceMask StaticConform
 ADMBase::initial_shift = "zero"
 
 ActiveThorns = "CalcK"
+CalcK::calc_extcurv = yes
 
 ActiveThorns = "IOBasic"
 IOBasic::outInfo_every      = 1
diff --git a/par/copy_to_next.par b/par/copy_to_next.par
new file mode 100644
index 0000000..ec68881
--- /dev/null
+++ b/par/copy_to_next.par
@@ -0,0 +1,68 @@
+# $Header$
+
+Cactus::cctk_run_title     = "Copy to next"
+Cactus::cctk_full_warnings = yes
+Cactus::cctk_timer_output  = full
+
+Cactus::cctk_itlast = 0
+
+ActiveThorns = "always"
+
+ActiveThorns = "PUGH PUGHInterp PUGHReduce PUGHSlab"
+driver::global_nx = 80
+driver::global_ny = 80
+driver::global_nz = 40
+
+ActiveThorns = "NaNChecker"
+
+ActiveThorns = "IOUtil"
+IO::out_dir = $parfile
+
+IO::filereader_ID_files = "
+	gxx_3d gxy_3d gxz_3d gyy_3d gyz_3d gzz_3d
+"
+IO::filereader_ID_vars = "
+	ADMBase::metric[cctk_iteration=102]
+"
+
+ActiveThorns = "IOHDF5 IOHDF5Util"
+
+ActiveThorns = "Boundary CartGrid3D CoordBase"
+CartGrid3D::type         = byspacing
+CartGrid3D::domain       = bitant
+CartGrid3D::dxyz         = 0.08
+CartGrid3D::avoid_origin = no
+
+ActiveThorns = "ADMBase ADMCoupling ADMMacros CoordGauge SpaceMask StaticConformal"
+
+ADMBase::initial_shift = "zero"
+
+ActiveThorns = "CalcK"
+CalcK::copy_to_next = yes
+
+ActiveThorns = "IOBasic"
+IOBasic::outInfo_every      = 1
+IOBasic::outInfo_reductions = "minimum maximum"
+IOBasic::outInfo_vars       = "
+        ADMBase::gxx
+        CalcK::gxx_next
+"
+
+IOBasic::outScalar_every      = 1
+IOBasic::outScalar_reductions = "minimum maximum mean norm1 norm2 norm_inf"
+IOBasic::outScalar_style      = gnuplot
+IOBasic::outScalar_vars       = "
+        ADMBase::metric
+"
+
+ActiveThorns = "IOASCII"
+IOASCII::out1D_every = 1
+IOASCII::out1D_style = "gnuplot f(x)"
+IOASCII::out1D_vars  = "
+        ADMBase::metric
+"
+
+IOHDF5::out_every = 1
+IOHDF5::out_vars  = "
+        CalcK::metric_next
+"
diff --git a/par/copy_to_prev.par b/par/copy_to_prev.par
new file mode 100644
index 0000000..85113f0
--- /dev/null
+++ b/par/copy_to_prev.par
@@ -0,0 +1,68 @@
+# $Header$
+
+Cactus::cctk_run_title     = "Copy to prev"
+Cactus::cctk_full_warnings = yes
+Cactus::cctk_timer_output  = full
+
+Cactus::cctk_itlast = 0
+
+ActiveThorns = "always"
+
+ActiveThorns = "PUGH PUGHInterp PUGHReduce PUGHSlab"
+driver::global_nx = 80
+driver::global_ny = 80
+driver::global_nz = 40
+
+ActiveThorns = "NaNChecker"
+
+ActiveThorns = "IOUtil"
+IO::out_dir = $parfile
+
+IO::filereader_ID_files = "
+	gxx_3d gxy_3d gxz_3d gyy_3d gyz_3d gzz_3d
+"
+IO::filereader_ID_vars = "
+	ADMBase::metric[cctk_iteration=100]
+"
+
+ActiveThorns = "IOHDF5 IOHDF5Util"
+
+ActiveThorns = "Boundary CartGrid3D CoordBase"
+CartGrid3D::type         = byspacing
+CartGrid3D::domain       = bitant
+CartGrid3D::dxyz         = 0.08
+CartGrid3D::avoid_origin = no
+
+ActiveThorns = "ADMBase ADMCoupling ADMMacros CoordGauge SpaceMask StaticConformal"
+
+ADMBase::initial_shift = "zero"
+
+ActiveThorns = "CalcK"
+CalcK::copy_to_prev = yes
+
+ActiveThorns = "IOBasic"
+IOBasic::outInfo_every      = 1
+IOBasic::outInfo_reductions = "minimum maximum"
+IOBasic::outInfo_vars       = "
+        ADMBase::gxx
+        CalcK::gxx_prev
+"
+
+IOBasic::outScalar_every      = 1
+IOBasic::outScalar_reductions = "minimum maximum mean norm1 norm2 norm_inf"
+IOBasic::outScalar_style      = gnuplot
+IOBasic::outScalar_vars       = "
+        ADMBase::metric
+"
+
+ActiveThorns = "IOASCII"
+IOASCII::out1D_every = 1
+IOASCII::out1D_style = "gnuplot f(x)"
+IOASCII::out1D_vars  = "
+        ADMBase::metric
+"
+
+IOHDF5::out_every = 1
+IOHDF5::out_vars  = "
+        CalcK::metric_prev
+"