Extract works fully with physical metric ... still some testing to do.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/Extract/trunk@62 5301f0c2-dbc4-4cee-b2f5-8d7afba4d129

70 files changed, 987 insertions, 26 deletions

diff --git a/src/ADMmass_integrand3D.F b/src/ADMmass_integrand3D.F
index a501ef1..566741a 100644
--- a/src/ADMmass_integrand3D.F
+++ b/src/ADMmass_integrand3D.F
@@ -89,7 +89,7 @@ c                 -------------------------------------
 
                  else 
 
-                    p = 1
+                    p = 1.0d0
 
                     dxx = gxx(i,j,k); dxy = gxy(i,j,k)
                     dxz = gxz(i,j,k); dyy = gyy(i,j,k)
diff --git a/src/D3_to_D2.F b/src/D3_to_D2.F
index 2948517..a7d22b2 100644
--- a/src/D3_to_D2.F
+++ b/src/D3_to_D2.F
@@ -288,7 +288,7 @@ c     ---------------------------------
 
          CALL momentum_integrand3D(origin,Dx,Dy,Dz,x,y,z,
      &        gxx,gxy,gxz,gyy,gyz,gzz,hxx,hxy,hxz,hyy,hyz,hzz,
-     &        Extract_temp3d,Psi,Psi_power)
+     &        Extract_temp3d,Psi,Psi_power,conformal_state)
          CALL CCTK_SyncGroup(ierror,cctkGH,"extract::temps")
          call CCTK_InterpGV (ierror, cctkGH, interp_handle, coord_system_handle,
      $                       npoints, 1, 1,
@@ -301,7 +301,7 @@ c     ---------------------------------
 
          CALL momentum_integrand3D(origin,Dx,Dy,Dz,x,y,z,
      &        gxx,gxy,gxz,gyy,gyz,gzz,hxx,hxy,hxz,hyy,hyz,hzz,
-     &        Extract_temp3d,Psi,Psi_power)
+     &        Extract_temp3d,Psi,Psi_power,conformal_state)
          CALL CCTK_SyncGroup(ierror,cctkGH,"extract::temps")
          call CCTK_InterpGV (ierror, cctkGH, interp_handle, coord_system_handle,
      $                       npoints, 1, 1,
@@ -315,7 +315,7 @@ c     ---------------------------------
 
          CALL momentum_integrand3D(origin,Dx,Dy,Dz,x,y,z,
      &        gxx,gxy,gxz,gyy,gyz,gzz,hxx,hxy,hxz,hyy,hyz,hzz,
-     &        Extract_temp3d,Psi,Psi_power)
+     &        Extract_temp3d,Psi,Psi_power,conformal_state)
          CALL CCTK_SyncGroup(ierror,cctkGH,"extract::temps")
          call CCTK_InterpGV (ierror, cctkGH, interp_handle, coord_system_handle,
      $                       npoints, 1, 1,
@@ -334,7 +334,7 @@ c     -----------------------------
 
          CALL spin_integrand3D(origin,x,y,z,
      &        gxx,gxy,gxz,gyy,gyz,gzz,hxx,hxy,hxz,hyy,hyz,hzz,
-     &        Extract_temp3d,Psi,Psi_power)
+     &        Extract_temp3d,Psi,Psi_power,conformal_state)
          CALL CCTK_SyncGroup(ierror,cctkGH,"extract::temps")
          call CCTK_InterpGV (ierror, cctkGH, interp_handle, coord_system_handle,
      $                       npoints, 1, 1,
@@ -347,7 +347,7 @@ c     -----------------------------
 
          CALL spin_integrand3D(origin,x,y,z,
      &        gxx,gxy,gxz,gyy,gyz,gzz,hxx,hxy,hxz,hyy,hyz,hzz,
-     &        Extract_temp3d,Psi,Psi_power)
+     &        Extract_temp3d,Psi,Psi_power,conformal_state)
          CALL CCTK_SyncGroup(ierror,cctkGH,"extract::temps")
          call CCTK_InterpGV (ierror, cctkGH, interp_handle, coord_system_handle,
      $                       npoints, 1, 1,
@@ -361,7 +361,7 @@ c     -----------------------------
 
          CALL spin_integrand3D(origin,x,y,z,
      &        gxx,gxy,gxz,gyy,gyz,gzz,hxx,hxy,hxz,hyy,hyz,hzz,
-     &        Extract_temp3d,Psi,Psi_power)
+     &        Extract_temp3d,Psi,Psi_power,conformal_state)
          CALL CCTK_SyncGroup(ierror,cctkGH,"extract::temps")
          call CCTK_InterpGV (ierror, cctkGH, interp_handle, coord_system_handle,
      $                       npoints, 1, 1,
diff --git a/src/momentum_integrand3D.F b/src/momentum_integrand3D.F
index 011fa68..cc481f9 100644
--- a/src/momentum_integrand3D.F
+++ b/src/momentum_integrand3D.F
@@ -5,7 +5,7 @@ c     ========================================================================
 
       SUBROUTINE momentum_integrand3D(origin,Dx,Dy,Dz,x,y,z,gxx,gxy,gxz,
      &     gyy,gyz,gzz,hxx,hxy,hxz,hyy,hyz,hzz,
-     &     momentum_int,Psi,Psi_power)
+     &     momentum_int,Psi,Psi_power,conformal_state)
 
 c     ------------------------------------------------------------------------
 c
@@ -18,7 +18,7 @@ c     -----------------------------------------------------------------------
 
 c     Input variables
       INTEGER,INTENT(IN) ::
-     &     Psi_power
+     &     Psi_power,conformal_state
       CCTK_REAL,INTENT(IN) :: 
      &     Dx,Dy,Dz,origin(3)
       CCTK_REAL,DIMENSION(:),INTENT(IN) :: 
@@ -81,11 +81,24 @@ c     Because other codes evolve Psi**4
 
 c                 Abbreviations for metric coefficients
 c                 -------------------------------------
-                  p = psi(i,j,k)**ip
+                  if (conformal_state > 0) then
+  
+                    p = psi(i,j,k)**ip
+
+                    dxx = p*gxx(i,j,k); dxy = p*gxy(i,j,k)
+                    dxz = p*gxz(i,j,k); dyy = p*gyy(i,j,k)
+                    dyz = p*gyz(i,j,k); dzz = p*gzz(i,j,k)
+
+                 else 
+
+                    p = 1.0d0
+
+                    dxx = gxx(i,j,k); dxy = gxy(i,j,k)
+                    dxz = gxz(i,j,k); dyy = gyy(i,j,k)
+                    dyz = gyz(i,j,k); dzz = gzz(i,j,k)
+
+                 end if
 
-                  dxx = p*gxx(i,j,k); dxy = p*gxy(i,j,k)
-                  dxz = p*gxz(i,j,k); dyy = p*gyy(i,j,k)
-                  dyz = p*gyz(i,j,k); dzz = p*gzz(i,j,k)
 
 c                 Determinant of 3-metric
 c                 -----------------------
diff --git a/src/momentum_integrand3D_int.F b/src/momentum_integrand3D_int.F
index e97174f..31a6f50 100644
--- a/src/momentum_integrand3D_int.F
+++ b/src/momentum_integrand3D_int.F
@@ -9,11 +9,11 @@ c     ------------------------------------------------------------------
 
           SUBROUTINE  momentum_integrand3D(origin,Dx,Dy,Dz,x,y,z,gxx,
      &        gxy,gxz,gyy,gyz,gzz,hxx,hxy,hxz,hyy,hyz,hzz,momentum_int,
-     &        Psi,Psi_power)
+     &        Psi,Psi_power,conformal_state)
 
           IMPLICIT NONE
           INTEGER,INTENT(IN) ::
-     &         Psi_power
+     &         Psi_power,conformal_state
           CCTK_REAL,INTENT(IN) :: 
      &        Dx,Dy,Dz,origin(3)
           CCTK_REAL,DIMENSION(:),INTENT(IN) :: 
diff --git a/src/spin_integrand3D.F b/src/spin_integrand3D.F
index cf617d8..7f93a0d 100644
--- a/src/spin_integrand3D.F
+++ b/src/spin_integrand3D.F
@@ -5,7 +5,7 @@ c     ========================================================================
 
       SUBROUTINE spin_integrand3D(origin,x,y,z,gxx,gxy,gxz,
      &     gyy,gyz,gzz,hxx,hxy,hxz,hyy,hyz,hzz,
-     &     spin_int,Psi,Psi_power)
+     &     spin_int,Psi,Psi_power,conformal_state)
 
 c     ------------------------------------------------------------------------
 c
@@ -19,7 +19,7 @@ c     -----------------------------------------------------------------------
 
 c     Input variables
       INTEGER,INTENT(IN) ::
-     &     Psi_power
+     &     Psi_power,conformal_state
       CCTK_REAL,INTENT(IN) :: 
      &     origin(3)
       CCTK_REAL,DIMENSION(:),INTENT(IN) :: 
@@ -81,11 +81,23 @@ c     ------------------------------------------------------------------------
 
 c                 Abbreviations for metric coefficients
 c                 -------------------------------------
-                  p = psi(i,j,k)**ip
+                  if (conformal_state > 0) then
+  
+                    p = psi(i,j,k)**ip
 
-                  dxx = p*gxx(i,j,k); dxy = p*gxy(i,j,k)
-                  dxz = p*gxz(i,j,k); dyy = p*gyy(i,j,k)
-                  dyz = p*gyz(i,j,k); dzz = p*gzz(i,j,k)
+                    dxx = p*gxx(i,j,k); dxy = p*gxy(i,j,k)
+                    dxz = p*gxz(i,j,k); dyy = p*gyy(i,j,k)
+                    dyz = p*gyz(i,j,k); dzz = p*gzz(i,j,k)
+
+                 else
+
+                    p = 1.0d0
+
+                    dxx = gxx(i,j,k); dxy = gxy(i,j,k)
+                    dxz = gxz(i,j,k); dyy = gyy(i,j,k)
+                    dyz = gyz(i,j,k); dzz = gzz(i,j,k)
+
+                 end if
 
 c                 Determinant of 3-metric
 c                 -----------------------
diff --git a/src/spin_integrand3D_int.F b/src/spin_integrand3D_int.F
index 1a31c78..4f9fc5e 100644
--- a/src/spin_integrand3D_int.F
+++ b/src/spin_integrand3D_int.F
@@ -9,11 +9,11 @@ c     ------------------------------------------------------------------
 
           SUBROUTINE  spin_integrand3D(origin,x,y,z,gxx,
      &        gxy,gxz,gyy,gyz,gzz,hxx,hxy,hxz,hyy,hyz,hzz,spin_int,
-     &        Psi,Psi_power)
+     &        Psi,Psi_power,conformal_state)
 
           IMPLICIT NONE
           INTEGER,INTENT(IN) ::
-     &         Psi_power
+     &         Psi_power,conformal_state
           CCTK_REAL,INTENT(IN) :: 
      &        origin(3)
           CCTK_REAL,DIMENSION(:),INTENT(IN) :: 
diff --git a/test/extract_physical.par b/test/extract_physical.par
new file mode 100644
index 0000000..1a20d93
--- /dev/null
+++ b/test/extract_physical.par
@@ -0,0 +1,54 @@
+!DESC "Extract waveforms from BH"
+
+
+ActiveThorns = "Boundary CartGrid3D Time IOUtil IOASCII PUGH PUGHSlab PUGHInterp Extract ADMBase StaticConformal IDAnalyticBH ADM ADMMacros CoordGauge"
+
+
+driver::global_nx = 15
+driver::global_ny = 15
+driver::global_nz = 15
+
+grid::type   = "bySpacing"
+grid::domain = "octant"
+grid::dxyz   = 0.3
+
+time::dtfac = 0.2  
+
+cactus::cctk_itlast = 5   
+IO::out_every       = 1      
+
+admbase::evolution_method ="adm"
+admbase::initial_data     = "schwarzschild" 
+admbase::metric_type = "physical"
+admbase::lapse_evolution_method          = "geodesic"
+idanalyticbh::mass	   = 1 
+
+IO::out_dir          = "extract_physical"
+io::out_fileinfo    = "none"
+ioascii::out1d_vars = "admbase::alp"
+
+### Extraction parameters
+
+extract::interpolation_order = 1
+extract::num_detectors = 2
+extract::itout         = 1
+
+extract::Nt        = 50
+extract::Np        = 50
+extract::origin_x  = 0.
+extract::origin_y  = 0.
+extract::origin_z  = 0.
+
+extract::all_modes  = "yes"
+extract::l_mode     = 4 
+extract::m_mode     = 0 
+
+extract::detector1 = 2.5 
+extract::detector2 = 3.0
+
+extract::doADMmass = "yes"
+extract::do_momentum = "yes"
+extract::do_spin = "yes"
+
+# do not use the new scheme for names of output files
+IO::new_filename_scheme = "no"
diff --git a/test/extract_physical/ADMmass_R1.tl b/test/extract_physical/ADMmass_R1.tl
new file mode 100644
index 0000000..768b4d2
--- /dev/null
+++ b/test/extract_physical/ADMmass_R1.tl
@@ -0,0 +1,7 @@
+"ADMmass_R1.tl 
+   0.00000000000000D+000   6.66832806766511D-002
+   6.00000000000000D-002   6.70419215988544D-002
+  0.120000000000000        6.81256065003433D-002
+  0.180000000000000        6.99578043210512D-002
+  0.240000000000000        7.25782235849593D-002
+  0.300000000000000        7.60436713612564D-002
diff --git a/test/extract_physical/ADMmass_R1.ul b/test/extract_physical/ADMmass_R1.ul
new file mode 100644
index 0000000..334834d
--- /dev/null
+++ b/test/extract_physical/ADMmass_R1.ul
@@ -0,0 +1,7 @@
+"ADMmass_R1.ul 
+   0.00000000000000D+000   6.66832806766511D-002
+   6.00012337233809D-002   6.70419215988544D-002
+  0.120002467446762        6.81256065003433D-002
+  0.180003701170143        6.99578043210512D-002
+  0.240004934893523        7.25782235849593D-002
+  0.300006168616904        7.60436713612564D-002
diff --git a/test/extract_physical/ADMmass_R2.tl b/test/extract_physical/ADMmass_R2.tl
new file mode 100644
index 0000000..c171640
--- /dev/null
+++ b/test/extract_physical/ADMmass_R2.tl
@@ -0,0 +1,7 @@
+"ADMmass_R2.tl 
+   0.00000000000000D+000   3.80317401216433D-002
+   6.00000000000000D-002   3.81730921853292D-002
+  0.120000000000000        3.85993978064040D-002
+  0.180000000000000        3.93174598328391D-002
+  0.240000000000000        4.03387987537457D-002
+  0.300000000000000        4.16799302676428D-002
diff --git a/test/extract_physical/ADMmass_R2.ul b/test/extract_physical/ADMmass_R2.ul
new file mode 100644
index 0000000..1a2a6ae
--- /dev/null
+++ b/test/extract_physical/ADMmass_R2.ul
@@ -0,0 +1,7 @@
+"ADMmass_R2.ul 
+   0.00000000000000D+000   3.80317401216433D-002
+   6.00012337233809D-002   3.81730921853292D-002
+  0.120002467446762        3.85993978064040D-002
+  0.180003701170143        3.93174598328391D-002
+  0.240004934893523        4.03387987537457D-002
+  0.300006168616904        4.16799302676428D-002
diff --git a/test/extract_physical/Qeven_R1_20.tl b/test/extract_physical/Qeven_R1_20.tl
new file mode 100644
index 0000000..88b98a3
--- /dev/null
+++ b/test/extract_physical/Qeven_R1_20.tl
@@ -0,0 +1,7 @@
+"Qeven_R1_20.tl
+   0.00000000000000D+000  -4.70935933465641D-006   0.00000000000000D+000
+   6.00000000000000D-002  -4.70863954630203D-006   0.00000000000000D+000
+  0.120000000000000       -4.70348115661468D-006   0.00000000000000D+000
+  0.180000000000000       -4.68454803026766D-006   0.00000000000000D+000
+  0.240000000000000       -4.63518565647178D-006   0.00000000000000D+000
+  0.300000000000000       -4.52991654628025D-006   0.00000000000000D+000
diff --git a/test/extract_physical/Qeven_R1_20.ul b/test/extract_physical/Qeven_R1_20.ul
new file mode 100644
index 0000000..e5a0f94
--- /dev/null
+++ b/test/extract_physical/Qeven_R1_20.ul
@@ -0,0 +1,7 @@
+"Qeven_R1_20.ul
+   0.00000000000000D+000  -4.70935933465641D-006   0.00000000000000D+000
+   6.00012337233809D-002  -4.70863954630203D-006   0.00000000000000D+000
+  0.120002467446762       -4.70348115661468D-006   0.00000000000000D+000
+  0.180003701170143       -4.68454803026766D-006   0.00000000000000D+000
+  0.240004934893523       -4.63518565647178D-006   0.00000000000000D+000
+  0.300006168616904       -4.52991654628025D-006   0.00000000000000D+000
diff --git a/test/extract_physical/Qeven_R1_22.tl b/test/extract_physical/Qeven_R1_22.tl
new file mode 100644
index 0000000..30e6203
--- /dev/null
+++ b/test/extract_physical/Qeven_R1_22.tl
@@ -0,0 +1,7 @@
+"Qeven_R1_22.tl
+   0.00000000000000D+000   7.63888023782705D-018   0.00000000000000D+000
+   6.00000000000000D-002  -1.79247954090654D-017   0.00000000000000D+000
+  0.120000000000000       -1.61683838617602D-017   0.00000000000000D+000
+  0.180000000000000        2.12098570696846D-017   0.00000000000000D+000
+  0.240000000000000        3.13472812608282D-017   0.00000000000000D+000
+  0.300000000000000        1.75721815945631D-017   0.00000000000000D+000
diff --git a/test/extract_physical/Qeven_R1_22.ul b/test/extract_physical/Qeven_R1_22.ul
new file mode 100644
index 0000000..96e1028
--- /dev/null
+++ b/test/extract_physical/Qeven_R1_22.ul
@@ -0,0 +1,7 @@
+"Qeven_R1_22.ul
+   0.00000000000000D+000   7.63888023782705D-018   0.00000000000000D+000
+   6.00012337233809D-002  -1.79247954090654D-017   0.00000000000000D+000
+  0.120002467446762       -1.61683838617602D-017   0.00000000000000D+000
+  0.180003701170143        2.12098570696846D-017   0.00000000000000D+000
+  0.240004934893523        3.13472812608282D-017   0.00000000000000D+000
+  0.300006168616904        1.75721815945631D-017   0.00000000000000D+000
diff --git a/test/extract_physical/Qeven_R1_40.tl b/test/extract_physical/Qeven_R1_40.tl
new file mode 100644
index 0000000..b0f20e3
--- /dev/null
+++ b/test/extract_physical/Qeven_R1_40.tl
@@ -0,0 +1,7 @@
+"Qeven_R1_40.tl
+   0.00000000000000D+000   1.44708663301951D-003   0.00000000000000D+000
+   6.00000000000000D-002   1.47950370238149D-003   0.00000000000000D+000
+  0.120000000000000        1.57809637147747D-003   0.00000000000000D+000
+  0.180000000000000        1.74691642027592D-003   0.00000000000000D+000
+  0.240000000000000        1.99279115345655D-003   0.00000000000000D+000
+  0.300000000000000        2.32538391111350D-003   0.00000000000000D+000
diff --git a/test/extract_physical/Qeven_R1_40.ul b/test/extract_physical/Qeven_R1_40.ul
new file mode 100644
index 0000000..7cea702
--- /dev/null
+++ b/test/extract_physical/Qeven_R1_40.ul
@@ -0,0 +1,7 @@
+"Qeven_R1_40.ul
+   0.00000000000000D+000   1.44708663301951D-003   0.00000000000000D+000
+   6.00012337233809D-002   1.47950370238149D-003   0.00000000000000D+000
+  0.120002467446762        1.57809637147747D-003   0.00000000000000D+000
+  0.180003701170143        1.74691642027592D-003   0.00000000000000D+000
+  0.240004934893523        1.99279115345655D-003   0.00000000000000D+000
+  0.300006168616904        2.32538391111350D-003   0.00000000000000D+000
diff --git a/test/extract_physical/Qeven_R1_42.tl b/test/extract_physical/Qeven_R1_42.tl
new file mode 100644
index 0000000..7da4bbf
--- /dev/null
+++ b/test/extract_physical/Qeven_R1_42.tl
@@ -0,0 +1,7 @@
+"Qeven_R1_42.tl
+   0.00000000000000D+000   9.01566387140683D-018   0.00000000000000D+000
+   6.00000000000000D-002   1.29624146133562D-017   0.00000000000000D+000
+  0.120000000000000       -4.64077338984613D-018   0.00000000000000D+000
+  0.180000000000000       -8.32151937915786D-018   0.00000000000000D+000
+  0.240000000000000       -1.06552891407775D-018   0.00000000000000D+000
+  0.300000000000000        4.15285974977847D-018   0.00000000000000D+000
diff --git a/test/extract_physical/Qeven_R1_42.ul b/test/extract_physical/Qeven_R1_42.ul
new file mode 100644
index 0000000..fba63bd
--- /dev/null
+++ b/test/extract_physical/Qeven_R1_42.ul
@@ -0,0 +1,7 @@
+"Qeven_R1_42.ul
+   0.00000000000000D+000   9.01566387140683D-018   0.00000000000000D+000
+   6.00012337233809D-002   1.29624146133562D-017   0.00000000000000D+000
+  0.120002467446762       -4.64077338984613D-018   0.00000000000000D+000
+  0.180003701170143       -8.32151937915786D-018   0.00000000000000D+000
+  0.240004934893523       -1.06552891407775D-018   0.00000000000000D+000
+  0.300006168616904        4.15285974977847D-018   0.00000000000000D+000
diff --git a/test/extract_physical/Qeven_R1_44.tl b/test/extract_physical/Qeven_R1_44.tl
new file mode 100644
index 0000000..e2a46b6
--- /dev/null
+++ b/test/extract_physical/Qeven_R1_44.tl
@@ -0,0 +1,7 @@
+"Qeven_R1_44.tl
+   0.00000000000000D+000   8.66573471858079D-004   0.00000000000000D+000
+   6.00000000000000D-002   8.85950122273846D-004   0.00000000000000D+000
+  0.120000000000000        9.44880281461343D-004   0.00000000000000D+000
+  0.180000000000000        1.04578083621531D-003   0.00000000000000D+000
+  0.240000000000000        1.19272400115100D-003   0.00000000000000D+000
+  0.300000000000000        1.39147298430074D-003   0.00000000000000D+000
diff --git a/test/extract_physical/Qeven_R1_44.ul b/test/extract_physical/Qeven_R1_44.ul
new file mode 100644
index 0000000..3134236
--- /dev/null
+++ b/test/extract_physical/Qeven_R1_44.ul
@@ -0,0 +1,7 @@
+"Qeven_R1_44.ul
+   0.00000000000000D+000   8.66573471858079D-004   0.00000000000000D+000
+   6.00012337233809D-002   8.85950122273846D-004   0.00000000000000D+000
+  0.120002467446762        9.44880281461343D-004   0.00000000000000D+000
+  0.180003701170143        1.04578083621531D-003   0.00000000000000D+000
+  0.240004934893523        1.19272400115100D-003   0.00000000000000D+000
+  0.300006168616904        1.39147298430074D-003   0.00000000000000D+000
diff --git a/test/extract_physical/Qeven_R2_20.tl b/test/extract_physical/Qeven_R2_20.tl
new file mode 100644
index 0000000..4341506
--- /dev/null
+++ b/test/extract_physical/Qeven_R2_20.tl
@@ -0,0 +1,7 @@
+"Qeven_R2_20.tl
+   0.00000000000000D+000   2.78006108622088D-006   0.00000000000000D+000
+   6.00000000000000D-002   2.78195307867046D-006   0.00000000000000D+000
+  0.120000000000000        2.78817516958580D-006   0.00000000000000D+000
+  0.180000000000000        2.80039311821707D-006   0.00000000000000D+000
+  0.240000000000000        2.82147161973111D-006   0.00000000000000D+000
+  0.300000000000000        2.85559882655036D-006   0.00000000000000D+000
diff --git a/test/extract_physical/Qeven_R2_20.ul b/test/extract_physical/Qeven_R2_20.ul
new file mode 100644
index 0000000..6d89d91
--- /dev/null
+++ b/test/extract_physical/Qeven_R2_20.ul
@@ -0,0 +1,7 @@
+"Qeven_R2_20.ul
+   0.00000000000000D+000   2.78006108622088D-006   0.00000000000000D+000
+   6.00012337233809D-002   2.78195307867046D-006   0.00000000000000D+000
+  0.120002467446762        2.78817516958580D-006   0.00000000000000D+000
+  0.180003701170143        2.80039311821707D-006   0.00000000000000D+000
+  0.240004934893523        2.82147161973111D-006   0.00000000000000D+000
+  0.300006168616904        2.85559882655036D-006   0.00000000000000D+000
diff --git a/test/extract_physical/Qeven_R2_22.tl b/test/extract_physical/Qeven_R2_22.tl
new file mode 100644
index 0000000..0afe5ae
--- /dev/null
+++ b/test/extract_physical/Qeven_R2_22.tl
@@ -0,0 +1,7 @@
+"Qeven_R2_22.tl
+   0.00000000000000D+000  -8.02668343797341D-018   0.00000000000000D+000
+   6.00000000000000D-002  -1.68584545514989D-017   0.00000000000000D+000
+  0.120000000000000       -1.40679415610302D-017   0.00000000000000D+000
+  0.180000000000000        1.34642410893566D-017   0.00000000000000D+000
+  0.240000000000000       -1.80637967904897D-017   0.00000000000000D+000
+  0.300000000000000       -6.38060173795871D-018   0.00000000000000D+000
diff --git a/test/extract_physical/Qeven_R2_22.ul b/test/extract_physical/Qeven_R2_22.ul
new file mode 100644
index 0000000..fb9f3ae
--- /dev/null
+++ b/test/extract_physical/Qeven_R2_22.ul
@@ -0,0 +1,7 @@
+"Qeven_R2_22.ul
+   0.00000000000000D+000  -8.02668343797341D-018   0.00000000000000D+000
+   6.00012337233809D-002  -1.68584545514989D-017   0.00000000000000D+000
+  0.120002467446762       -1.40679415610302D-017   0.00000000000000D+000
+  0.180003701170143        1.34642410893566D-017   0.00000000000000D+000
+  0.240004934893523       -1.80637967904897D-017   0.00000000000000D+000
+  0.300006168616904       -6.38060173795871D-018   0.00000000000000D+000
diff --git a/test/extract_physical/Qeven_R2_40.tl b/test/extract_physical/Qeven_R2_40.tl
new file mode 100644
index 0000000..fbca399
--- /dev/null
+++ b/test/extract_physical/Qeven_R2_40.tl
@@ -0,0 +1,7 @@
+"Qeven_R2_40.tl
+   0.00000000000000D+000   4.34645147888665D-004   0.00000000000000D+000
+   6.00000000000000D-002   4.42943704062136D-004   0.00000000000000D+000
+  0.120000000000000        4.68109800079682D-004   0.00000000000000D+000
+  0.180000000000000        5.10966356354637D-004   0.00000000000000D+000
+  0.240000000000000        5.72923598198710D-004   0.00000000000000D+000
+  0.300000000000000        6.56037215861202D-004   0.00000000000000D+000
diff --git a/test/extract_physical/Qeven_R2_40.ul b/test/extract_physical/Qeven_R2_40.ul
new file mode 100644
index 0000000..89c66c0
--- /dev/null
+++ b/test/extract_physical/Qeven_R2_40.ul
@@ -0,0 +1,7 @@
+"Qeven_R2_40.ul
+   0.00000000000000D+000   4.34645147888665D-004   0.00000000000000D+000
+   6.00012337233809D-002   4.42943704062136D-004   0.00000000000000D+000
+  0.120002467446762        4.68109800079682D-004   0.00000000000000D+000
+  0.180003701170143        5.10966356354637D-004   0.00000000000000D+000
+  0.240004934893523        5.72923598198710D-004   0.00000000000000D+000
+  0.300006168616904        6.56037215861202D-004   0.00000000000000D+000
diff --git a/test/extract_physical/Qeven_R2_42.tl b/test/extract_physical/Qeven_R2_42.tl
new file mode 100644
index 0000000..af77a81
--- /dev/null
+++ b/test/extract_physical/Qeven_R2_42.tl
@@ -0,0 +1,7 @@
+"Qeven_R2_42.tl
+   0.00000000000000D+000   1.00630187148561D-017   0.00000000000000D+000
+   6.00000000000000D-002   2.39466243453657D-017   0.00000000000000D+000
+  0.120000000000000        4.26486790280587D-018   0.00000000000000D+000
+  0.180000000000000        8.91660621376682D-018   0.00000000000000D+000
+  0.240000000000000       -3.14068235692405D-018   0.00000000000000D+000
+  0.300000000000000        1.15435752887110D-017   0.00000000000000D+000
diff --git a/test/extract_physical/Qeven_R2_42.ul b/test/extract_physical/Qeven_R2_42.ul
new file mode 100644
index 0000000..e7cf308
--- /dev/null
+++ b/test/extract_physical/Qeven_R2_42.ul
@@ -0,0 +1,7 @@
+"Qeven_R2_42.ul
+   0.00000000000000D+000   1.00630187148561D-017   0.00000000000000D+000
+   6.00012337233809D-002   2.39466243453657D-017   0.00000000000000D+000
+  0.120002467446762        4.26486790280587D-018   0.00000000000000D+000
+  0.180003701170143        8.91660621376682D-018   0.00000000000000D+000
+  0.240004934893523       -3.14068235692405D-018   0.00000000000000D+000
+  0.300006168616904        1.15435752887110D-017   0.00000000000000D+000
diff --git a/test/extract_physical/Qeven_R2_44.tl b/test/extract_physical/Qeven_R2_44.tl
new file mode 100644
index 0000000..d2adaae
--- /dev/null
+++ b/test/extract_physical/Qeven_R2_44.tl
@@ -0,0 +1,7 @@
+"Qeven_R2_44.tl
+   0.00000000000000D+000   2.55360577828519D-004   0.00000000000000D+000
+   6.00000000000000D-002   2.60294758509547D-004   0.00000000000000D+000
+  0.120000000000000        2.75258442065402D-004   0.00000000000000D+000
+  0.180000000000000        3.00741998318653D-004   0.00000000000000D+000
+  0.240000000000000        3.37585801584188D-004   0.00000000000000D+000
+  0.300000000000000        3.87014943801537D-004   0.00000000000000D+000
diff --git a/test/extract_physical/Qeven_R2_44.ul b/test/extract_physical/Qeven_R2_44.ul
new file mode 100644
index 0000000..69a0bfc
--- /dev/null
+++ b/test/extract_physical/Qeven_R2_44.ul
@@ -0,0 +1,7 @@
+"Qeven_R2_44.ul
+   0.00000000000000D+000   2.55360577828519D-004   0.00000000000000D+000
+   6.00012337233809D-002   2.60294758509547D-004   0.00000000000000D+000
+  0.120002467446762        2.75258442065402D-004   0.00000000000000D+000
+  0.180003701170143        3.00741998318653D-004   0.00000000000000D+000
+  0.240004934893523        3.37585801584188D-004   0.00000000000000D+000
+  0.300006168616904        3.87014943801537D-004   0.00000000000000D+000
diff --git a/test/extract_physical/alp.dl b/test/extract_physical/alp.dl
new file mode 100644
index 0000000..907c321
--- /dev/null
+++ b/test/extract_physical/alp.dl
@@ -0,0 +1,108 @@
+
+
+"Time = 0.0000000000000
+-0.2598076211353		1.0000000000000
+0.2598076211353		1.0000000000000
+0.7794228634060		1.0000000000000
+1.2990381056767		1.0000000000000
+1.8186533479473		1.0000000000000
+2.3382685902180		1.0000000000000
+2.8578838324886		1.0000000000000
+3.3774990747593		1.0000000000000
+3.8971143170300		1.0000000000000
+4.4167295593006		1.0000000000000
+4.9363448015713		1.0000000000000
+5.4559600438420		1.0000000000000
+5.9755752861126		1.0000000000000
+6.4951905283833		1.0000000000000
+7.0148057706540		1.0000000000000
+
+
+"Time = 0.0600000000000
+-0.2598076211353		1.0000000000000
+0.2598076211353		1.0000000000000
+0.7794228634060		1.0000000000000
+1.2990381056767		1.0000000000000
+1.8186533479473		1.0000000000000
+2.3382685902180		1.0000000000000
+2.8578838324886		1.0000000000000
+3.3774990747593		1.0000000000000
+3.8971143170300		1.0000000000000
+4.4167295593006		1.0000000000000
+4.9363448015713		1.0000000000000
+5.4559600438420		1.0000000000000
+5.9755752861126		1.0000000000000
+6.4951905283833		1.0000000000000
+7.0148057706540		1.0000000000000
+
+
+"Time = 0.1200000000000
+-0.2598076211353		1.0000000000000
+0.2598076211353		1.0000000000000
+0.7794228634060		1.0000000000000
+1.2990381056767		1.0000000000000
+1.8186533479473		1.0000000000000
+2.3382685902180		1.0000000000000
+2.8578838324886		1.0000000000000
+3.3774990747593		1.0000000000000
+3.8971143170300		1.0000000000000
+4.4167295593006		1.0000000000000
+4.9363448015713		1.0000000000000
+5.4559600438420		1.0000000000000
+5.9755752861126		1.0000000000000
+6.4951905283833		1.0000000000000
+7.0148057706540		1.0000000000000
+
+
+"Time = 0.1800000000000
+-0.2598076211353		1.0000000000000
+0.2598076211353		1.0000000000000
+0.7794228634060		1.0000000000000
+1.2990381056767		1.0000000000000
+1.8186533479473		1.0000000000000
+2.3382685902180		1.0000000000000
+2.8578838324886		1.0000000000000
+3.3774990747593		1.0000000000000
+3.8971143170300		1.0000000000000
+4.4167295593006		1.0000000000000
+4.9363448015713		1.0000000000000
+5.4559600438420		1.0000000000000
+5.9755752861126		1.0000000000000
+6.4951905283833		1.0000000000000
+7.0148057706540		1.0000000000000
+
+
+"Time = 0.2400000000000
+-0.2598076211353		1.0000000000000
+0.2598076211353		1.0000000000000
+0.7794228634060		1.0000000000000
+1.2990381056767		1.0000000000000
+1.8186533479473		1.0000000000000
+2.3382685902180		1.0000000000000
+2.8578838324886		1.0000000000000
+3.3774990747593		1.0000000000000
+3.8971143170300		1.0000000000000
+4.4167295593006		1.0000000000000
+4.9363448015713		1.0000000000000
+5.4559600438420		1.0000000000000
+5.9755752861126		1.0000000000000
+6.4951905283833		1.0000000000000
+7.0148057706540		1.0000000000000
+
+
+"Time = 0.3000000000000
+-0.2598076211353		1.0000000000000
+0.2598076211353		1.0000000000000
+0.7794228634060		1.0000000000000
+1.2990381056767		1.0000000000000
+1.8186533479473		1.0000000000000
+2.3382685902180		1.0000000000000
+2.8578838324886		1.0000000000000
+3.3774990747593		1.0000000000000
+3.8971143170300		1.0000000000000
+4.4167295593006		1.0000000000000
+4.9363448015713		1.0000000000000
+5.4559600438420		1.0000000000000
+5.9755752861126		1.0000000000000
+6.4951905283833		1.0000000000000
+7.0148057706540		1.0000000000000
diff --git a/test/extract_physical/alp.xl b/test/extract_physical/alp.xl
new file mode 100644
index 0000000..9a16fb2
--- /dev/null
+++ b/test/extract_physical/alp.xl
@@ -0,0 +1,108 @@
+
+
+"Time = 0.0000000000000
+-0.1500000000000		1.0000000000000
+0.1500000000000		1.0000000000000
+0.4500000000000		1.0000000000000
+0.7500000000000		1.0000000000000
+1.0500000000000		1.0000000000000
+1.3500000000000		1.0000000000000
+1.6500000000000		1.0000000000000
+1.9500000000000		1.0000000000000
+2.2500000000000		1.0000000000000
+2.5500000000000		1.0000000000000
+2.8500000000000		1.0000000000000
+3.1500000000000		1.0000000000000
+3.4500000000000		1.0000000000000
+3.7500000000000		1.0000000000000
+4.0500000000000		1.0000000000000
+
+
+"Time = 0.0600000000000
+-0.1500000000000		1.0000000000000
+0.1500000000000		1.0000000000000
+0.4500000000000		1.0000000000000
+0.7500000000000		1.0000000000000
+1.0500000000000		1.0000000000000
+1.3500000000000		1.0000000000000
+1.6500000000000		1.0000000000000
+1.9500000000000		1.0000000000000
+2.2500000000000		1.0000000000000
+2.5500000000000		1.0000000000000
+2.8500000000000		1.0000000000000
+3.1500000000000		1.0000000000000
+3.4500000000000		1.0000000000000
+3.7500000000000		1.0000000000000
+4.0500000000000		1.0000000000000
+
+
+"Time = 0.1200000000000
+-0.1500000000000		1.0000000000000
+0.1500000000000		1.0000000000000
+0.4500000000000		1.0000000000000
+0.7500000000000		1.0000000000000
+1.0500000000000		1.0000000000000
+1.3500000000000		1.0000000000000
+1.6500000000000		1.0000000000000
+1.9500000000000		1.0000000000000
+2.2500000000000		1.0000000000000
+2.5500000000000		1.0000000000000
+2.8500000000000		1.0000000000000
+3.1500000000000		1.0000000000000
+3.4500000000000		1.0000000000000
+3.7500000000000		1.0000000000000
+4.0500000000000		1.0000000000000
+
+
+"Time = 0.1800000000000
+-0.1500000000000		1.0000000000000
+0.1500000000000		1.0000000000000
+0.4500000000000		1.0000000000000
+0.7500000000000		1.0000000000000
+1.0500000000000		1.0000000000000
+1.3500000000000		1.0000000000000
+1.6500000000000		1.0000000000000
+1.9500000000000		1.0000000000000
+2.2500000000000		1.0000000000000
+2.5500000000000		1.0000000000000
+2.8500000000000		1.0000000000000
+3.1500000000000		1.0000000000000
+3.4500000000000		1.0000000000000
+3.7500000000000		1.0000000000000
+4.0500000000000		1.0000000000000
+
+
+"Time = 0.2400000000000
+-0.1500000000000		1.0000000000000
+0.1500000000000		1.0000000000000
+0.4500000000000		1.0000000000000
+0.7500000000000		1.0000000000000
+1.0500000000000		1.0000000000000
+1.3500000000000		1.0000000000000
+1.6500000000000		1.0000000000000
+1.9500000000000		1.0000000000000
+2.2500000000000		1.0000000000000
+2.5500000000000		1.0000000000000
+2.8500000000000		1.0000000000000
+3.1500000000000		1.0000000000000
+3.4500000000000		1.0000000000000
+3.7500000000000		1.0000000000000
+4.0500000000000		1.0000000000000
+
+
+"Time = 0.3000000000000
+-0.1500000000000		1.0000000000000
+0.1500000000000		1.0000000000000
+0.4500000000000		1.0000000000000
+0.7500000000000		1.0000000000000
+1.0500000000000		1.0000000000000
+1.3500000000000		1.0000000000000
+1.6500000000000		1.0000000000000
+1.9500000000000		1.0000000000000
+2.2500000000000		1.0000000000000
+2.5500000000000		1.0000000000000
+2.8500000000000		1.0000000000000
+3.1500000000000		1.0000000000000
+3.4500000000000		1.0000000000000
+3.7500000000000		1.0000000000000
+4.0500000000000		1.0000000000000
diff --git a/test/extract_physical/alp.yl b/test/extract_physical/alp.yl
new file mode 100644
index 0000000..9a16fb2
--- /dev/null
+++ b/test/extract_physical/alp.yl
@@ -0,0 +1,108 @@
+
+
+"Time = 0.0000000000000
+-0.1500000000000		1.0000000000000
+0.1500000000000		1.0000000000000
+0.4500000000000		1.0000000000000
+0.7500000000000		1.0000000000000
+1.0500000000000		1.0000000000000
+1.3500000000000		1.0000000000000
+1.6500000000000		1.0000000000000
+1.9500000000000		1.0000000000000
+2.2500000000000		1.0000000000000
+2.5500000000000		1.0000000000000
+2.8500000000000		1.0000000000000
+3.1500000000000		1.0000000000000
+3.4500000000000		1.0000000000000
+3.7500000000000		1.0000000000000
+4.0500000000000		1.0000000000000
+
+
+"Time = 0.0600000000000
+-0.1500000000000		1.0000000000000
+0.1500000000000		1.0000000000000
+0.4500000000000		1.0000000000000
+0.7500000000000		1.0000000000000
+1.0500000000000		1.0000000000000
+1.3500000000000		1.0000000000000
+1.6500000000000		1.0000000000000
+1.9500000000000		1.0000000000000
+2.2500000000000		1.0000000000000
+2.5500000000000		1.0000000000000
+2.8500000000000		1.0000000000000
+3.1500000000000		1.0000000000000
+3.4500000000000		1.0000000000000
+3.7500000000000		1.0000000000000
+4.0500000000000		1.0000000000000
+
+
+"Time = 0.1200000000000
+-0.1500000000000		1.0000000000000
+0.1500000000000		1.0000000000000
+0.4500000000000		1.0000000000000
+0.7500000000000		1.0000000000000
+1.0500000000000		1.0000000000000
+1.3500000000000		1.0000000000000
+1.6500000000000		1.0000000000000
+1.9500000000000		1.0000000000000
+2.2500000000000		1.0000000000000
+2.5500000000000		1.0000000000000
+2.8500000000000		1.0000000000000
+3.1500000000000		1.0000000000000
+3.4500000000000		1.0000000000000
+3.7500000000000		1.0000000000000
+4.0500000000000		1.0000000000000
+
+
+"Time = 0.1800000000000
+-0.1500000000000		1.0000000000000
+0.1500000000000		1.0000000000000
+0.4500000000000		1.0000000000000
+0.7500000000000		1.0000000000000
+1.0500000000000		1.0000000000000
+1.3500000000000		1.0000000000000
+1.6500000000000		1.0000000000000
+1.9500000000000		1.0000000000000
+2.2500000000000		1.0000000000000
+2.5500000000000		1.0000000000000
+2.8500000000000		1.0000000000000
+3.1500000000000		1.0000000000000
+3.4500000000000		1.0000000000000
+3.7500000000000		1.0000000000000
+4.0500000000000		1.0000000000000
+
+
+"Time = 0.2400000000000
+-0.1500000000000		1.0000000000000
+0.1500000000000		1.0000000000000
+0.4500000000000		1.0000000000000
+0.7500000000000		1.0000000000000
+1.0500000000000		1.0000000000000
+1.3500000000000		1.0000000000000
+1.6500000000000		1.0000000000000
+1.9500000000000		1.0000000000000
+2.2500000000000		1.0000000000000
+2.5500000000000		1.0000000000000
+2.8500000000000		1.0000000000000
+3.1500000000000		1.0000000000000
+3.4500000000000		1.0000000000000
+3.7500000000000		1.0000000000000
+4.0500000000000		1.0000000000000
+
+
+"Time = 0.3000000000000
+-0.1500000000000		1.0000000000000
+0.1500000000000		1.0000000000000
+0.4500000000000		1.0000000000000
+0.7500000000000		1.0000000000000
+1.0500000000000		1.0000000000000
+1.3500000000000		1.0000000000000
+1.6500000000000		1.0000000000000
+1.9500000000000		1.0000000000000
+2.2500000000000		1.0000000000000
+2.5500000000000		1.0000000000000
+2.8500000000000		1.0000000000000
+3.1500000000000		1.0000000000000
+3.4500000000000		1.0000000000000
+3.7500000000000		1.0000000000000
+4.0500000000000		1.0000000000000
diff --git a/test/extract_physical/alp.zl b/test/extract_physical/alp.zl
new file mode 100644
index 0000000..9a16fb2
--- /dev/null
+++ b/test/extract_physical/alp.zl
@@ -0,0 +1,108 @@
+
+
+"Time = 0.0000000000000
+-0.1500000000000		1.0000000000000
+0.1500000000000		1.0000000000000
+0.4500000000000		1.0000000000000
+0.7500000000000		1.0000000000000
+1.0500000000000		1.0000000000000
+1.3500000000000		1.0000000000000
+1.6500000000000		1.0000000000000
+1.9500000000000		1.0000000000000
+2.2500000000000		1.0000000000000
+2.5500000000000		1.0000000000000
+2.8500000000000		1.0000000000000
+3.1500000000000		1.0000000000000
+3.4500000000000		1.0000000000000
+3.7500000000000		1.0000000000000
+4.0500000000000		1.0000000000000
+
+
+"Time = 0.0600000000000
+-0.1500000000000		1.0000000000000
+0.1500000000000		1.0000000000000
+0.4500000000000		1.0000000000000
+0.7500000000000		1.0000000000000
+1.0500000000000		1.0000000000000
+1.3500000000000		1.0000000000000
+1.6500000000000		1.0000000000000
+1.9500000000000		1.0000000000000
+2.2500000000000		1.0000000000000
+2.5500000000000		1.0000000000000
+2.8500000000000		1.0000000000000
+3.1500000000000		1.0000000000000
+3.4500000000000		1.0000000000000
+3.7500000000000		1.0000000000000
+4.0500000000000		1.0000000000000
+
+
+"Time = 0.1200000000000
+-0.1500000000000		1.0000000000000
+0.1500000000000		1.0000000000000
+0.4500000000000		1.0000000000000
+0.7500000000000		1.0000000000000
+1.0500000000000		1.0000000000000
+1.3500000000000		1.0000000000000
+1.6500000000000		1.0000000000000
+1.9500000000000		1.0000000000000
+2.2500000000000		1.0000000000000
+2.5500000000000		1.0000000000000
+2.8500000000000		1.0000000000000
+3.1500000000000		1.0000000000000
+3.4500000000000		1.0000000000000
+3.7500000000000		1.0000000000000
+4.0500000000000		1.0000000000000
+
+
+"Time = 0.1800000000000
+-0.1500000000000		1.0000000000000
+0.1500000000000		1.0000000000000
+0.4500000000000		1.0000000000000
+0.7500000000000		1.0000000000000
+1.0500000000000		1.0000000000000
+1.3500000000000		1.0000000000000
+1.6500000000000		1.0000000000000
+1.9500000000000		1.0000000000000
+2.2500000000000		1.0000000000000
+2.5500000000000		1.0000000000000
+2.8500000000000		1.0000000000000
+3.1500000000000		1.0000000000000
+3.4500000000000		1.0000000000000
+3.7500000000000		1.0000000000000
+4.0500000000000		1.0000000000000
+
+
+"Time = 0.2400000000000
+-0.1500000000000		1.0000000000000
+0.1500000000000		1.0000000000000
+0.4500000000000		1.0000000000000
+0.7500000000000		1.0000000000000
+1.0500000000000		1.0000000000000
+1.3500000000000		1.0000000000000
+1.6500000000000		1.0000000000000
+1.9500000000000		1.0000000000000
+2.2500000000000		1.0000000000000
+2.5500000000000		1.0000000000000
+2.8500000000000		1.0000000000000
+3.1500000000000		1.0000000000000
+3.4500000000000		1.0000000000000
+3.7500000000000		1.0000000000000
+4.0500000000000		1.0000000000000
+
+
+"Time = 0.3000000000000
+-0.1500000000000		1.0000000000000
+0.1500000000000		1.0000000000000
+0.4500000000000		1.0000000000000
+0.7500000000000		1.0000000000000
+1.0500000000000		1.0000000000000
+1.3500000000000		1.0000000000000
+1.6500000000000		1.0000000000000
+1.9500000000000		1.0000000000000
+2.2500000000000		1.0000000000000
+2.5500000000000		1.0000000000000
+2.8500000000000		1.0000000000000
+3.1500000000000		1.0000000000000
+3.4500000000000		1.0000000000000
+3.7500000000000		1.0000000000000
+4.0500000000000		1.0000000000000
diff --git a/test/extract_physical/extract_physical.par b/test/extract_physical/extract_physical.par
new file mode 100644
index 0000000..1a20d93
--- /dev/null
+++ b/test/extract_physical/extract_physical.par
@@ -0,0 +1,54 @@
+!DESC "Extract waveforms from BH"
+
+
+ActiveThorns = "Boundary CartGrid3D Time IOUtil IOASCII PUGH PUGHSlab PUGHInterp Extract ADMBase StaticConformal IDAnalyticBH ADM ADMMacros CoordGauge"
+
+
+driver::global_nx = 15
+driver::global_ny = 15
+driver::global_nz = 15
+
+grid::type   = "bySpacing"
+grid::domain = "octant"
+grid::dxyz   = 0.3
+
+time::dtfac = 0.2  
+
+cactus::cctk_itlast = 5   
+IO::out_every       = 1      
+
+admbase::evolution_method ="adm"
+admbase::initial_data     = "schwarzschild" 
+admbase::metric_type = "physical"
+admbase::lapse_evolution_method          = "geodesic"
+idanalyticbh::mass	   = 1 
+
+IO::out_dir          = "extract_physical"
+io::out_fileinfo    = "none"
+ioascii::out1d_vars = "admbase::alp"
+
+### Extraction parameters
+
+extract::interpolation_order = 1
+extract::num_detectors = 2
+extract::itout         = 1
+
+extract::Nt        = 50
+extract::Np        = 50
+extract::origin_x  = 0.
+extract::origin_y  = 0.
+extract::origin_z  = 0.
+
+extract::all_modes  = "yes"
+extract::l_mode     = 4 
+extract::m_mode     = 0 
+
+extract::detector1 = 2.5 
+extract::detector2 = 3.0
+
+extract::doADMmass = "yes"
+extract::do_momentum = "yes"
+extract::do_spin = "yes"
+
+# do not use the new scheme for names of output files
+IO::new_filename_scheme = "no"
diff --git a/test/extract_physical/mass_R1.tl b/test/extract_physical/mass_R1.tl
new file mode 100644
index 0000000..ddbe79a
--- /dev/null
+++ b/test/extract_physical/mass_R1.tl
@@ -0,0 +1,7 @@
+"mass_R1.tl    
+   0.00000000000000D+000   6.61218534601497D-002
+   6.00000000000000D-002   6.64629086221119D-002
+  0.120000000000000        6.74926863635501D-002
+  0.180000000000000        6.92311174898341D-002
+  0.240000000000000        7.17117113214540D-002
+  0.300000000000000        7.49819375473559D-002
diff --git a/test/extract_physical/mass_R1.ul b/test/extract_physical/mass_R1.ul
new file mode 100644
index 0000000..7b1f6c5
--- /dev/null
+++ b/test/extract_physical/mass_R1.ul
@@ -0,0 +1,7 @@
+"mass_R1.ul    
+   0.00000000000000D+000   6.61218534601497D-002
+   6.00012337233809D-002   6.64629086221119D-002
+  0.120002467446762        6.74926863635501D-002
+  0.180003701170143        6.92311174898341D-002
+  0.240004934893523        7.17117113214540D-002
+  0.300006168616904        7.49819375473559D-002
diff --git a/test/extract_physical/mass_R2.tl b/test/extract_physical/mass_R2.tl
new file mode 100644
index 0000000..9669c18
--- /dev/null
+++ b/test/extract_physical/mass_R2.tl
@@ -0,0 +1,7 @@
+"mass_R2.tl    
+   0.00000000000000D+000   3.79970161063922D-002
+   6.00000000000000D-002   3.81353009437551D-002
+  0.120000000000000        3.85522585175652D-002
+  0.180000000000000        3.92542438588434D-002
+  0.240000000000000        4.02520014140427D-002
+  0.300000000000000        4.15608959438484D-002
diff --git a/test/extract_physical/mass_R2.ul b/test/extract_physical/mass_R2.ul
new file mode 100644
index 0000000..0cf82bd
--- /dev/null
+++ b/test/extract_physical/mass_R2.ul
@@ -0,0 +1,7 @@
+"mass_R2.ul    
+   0.00000000000000D+000   3.79970161063922D-002
+   6.00012337233809D-002   3.81353009437551D-002
+  0.120002467446762        3.85522585175652D-002
+  0.180003701170143        3.92542438588434D-002
+  0.240004934893523        4.02520014140427D-002
+  0.300006168616904        4.15608959438484D-002
diff --git a/test/extract_physical/momentum_x_R1.tl b/test/extract_physical/momentum_x_R1.tl
new file mode 100644
index 0000000..95b9275
--- /dev/null
+++ b/test/extract_physical/momentum_x_R1.tl
@@ -0,0 +1,7 @@
+"momentum_x_R1.
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00000000000000D-002   1.52177558310929D-003
+  0.120000000000000        3.07348618677577D-003
+  0.180000000000000        4.68584017545205D-003
+  0.240000000000000        6.39109977136692D-003
+  0.300000000000000        8.22388170610372D-003
diff --git a/test/extract_physical/momentum_x_R1.ul b/test/extract_physical/momentum_x_R1.ul
new file mode 100644
index 0000000..f6c00e3
--- /dev/null
+++ b/test/extract_physical/momentum_x_R1.ul
@@ -0,0 +1,7 @@
+"momentum_x_R1.
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00012337233809D-002   1.52177558310929D-003
+  0.120002467446762        3.07348618677577D-003
+  0.180003701170143        4.68584017545205D-003
+  0.240004934893523        6.39109977136692D-003
+  0.300006168616904        8.22388170610372D-003
diff --git a/test/extract_physical/momentum_x_R2.tl b/test/extract_physical/momentum_x_R2.tl
new file mode 100644
index 0000000..d640208
--- /dev/null
+++ b/test/extract_physical/momentum_x_R2.tl
@@ -0,0 +1,7 @@
+"momentum_x_R2.
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00000000000000D-002   7.40186023914552D-004
+  0.120000000000000        1.49082992352751D-003
+  0.180000000000000        2.26261288068105D-003
+  0.240000000000000        3.06666767190871D-003
+  0.300000000000000        3.91481790514123D-003
diff --git a/test/extract_physical/momentum_x_R2.ul b/test/extract_physical/momentum_x_R2.ul
new file mode 100644
index 0000000..ffb0d99
--- /dev/null
+++ b/test/extract_physical/momentum_x_R2.ul
@@ -0,0 +1,7 @@
+"momentum_x_R2.
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00012337233809D-002   7.40186023914552D-004
+  0.120002467446762        1.49082992352751D-003
+  0.180003701170143        2.26261288068105D-003
+  0.240004934893523        3.06666767190871D-003
+  0.300006168616904        3.91481790514123D-003
diff --git a/test/extract_physical/momentum_y_R1.tl b/test/extract_physical/momentum_y_R1.tl
new file mode 100644
index 0000000..c954e16
--- /dev/null
+++ b/test/extract_physical/momentum_y_R1.tl
@@ -0,0 +1,7 @@
+"momentum_y_R1.
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00000000000000D-002   1.52177558310929D-003
+  0.120000000000000        3.07348618677577D-003
+  0.180000000000000        4.68584017545204D-003
+  0.240000000000000        6.39109977136692D-003
+  0.300000000000000        8.22388170610371D-003
diff --git a/test/extract_physical/momentum_y_R1.ul b/test/extract_physical/momentum_y_R1.ul
new file mode 100644
index 0000000..065f926
--- /dev/null
+++ b/test/extract_physical/momentum_y_R1.ul
@@ -0,0 +1,7 @@
+"momentum_y_R1.
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00012337233809D-002   1.52177558310929D-003
+  0.120002467446762        3.07348618677577D-003
+  0.180003701170143        4.68584017545204D-003
+  0.240004934893523        6.39109977136692D-003
+  0.300006168616904        8.22388170610371D-003
diff --git a/test/extract_physical/momentum_y_R2.tl b/test/extract_physical/momentum_y_R2.tl
new file mode 100644
index 0000000..2ae683b
--- /dev/null
+++ b/test/extract_physical/momentum_y_R2.tl
@@ -0,0 +1,7 @@
+"momentum_y_R2.
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00000000000000D-002   7.40186023914552D-004
+  0.120000000000000        1.49082992352751D-003
+  0.180000000000000        2.26261288068105D-003
+  0.240000000000000        3.06666767190871D-003
+  0.300000000000000        3.91481790514123D-003
diff --git a/test/extract_physical/momentum_y_R2.ul b/test/extract_physical/momentum_y_R2.ul
new file mode 100644
index 0000000..4f8411c
--- /dev/null
+++ b/test/extract_physical/momentum_y_R2.ul
@@ -0,0 +1,7 @@
+"momentum_y_R2.
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00012337233809D-002   7.40186023914552D-004
+  0.120002467446762        1.49082992352751D-003
+  0.180003701170143        2.26261288068105D-003
+  0.240004934893523        3.06666767190871D-003
+  0.300006168616904        3.91481790514123D-003
diff --git a/test/extract_physical/momentum_z_R1.tl b/test/extract_physical/momentum_z_R1.tl
new file mode 100644
index 0000000..f7618f7
--- /dev/null
+++ b/test/extract_physical/momentum_z_R1.tl
@@ -0,0 +1,7 @@
+"momentum_z_R1.
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00000000000000D-002   1.52199082739009D-003
+  0.120000000000000        3.07392128045055D-003
+  0.180000000000000        4.68650455527377D-003
+  0.240000000000000        6.39200816662428D-003
+  0.300000000000000        8.22505484590828D-003
diff --git a/test/extract_physical/momentum_z_R1.ul b/test/extract_physical/momentum_z_R1.ul
new file mode 100644
index 0000000..ecce202
--- /dev/null
+++ b/test/extract_physical/momentum_z_R1.ul
@@ -0,0 +1,7 @@
+"momentum_z_R1.
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00012337233809D-002   1.52199082739009D-003
+  0.120002467446762        3.07392128045055D-003
+  0.180003701170143        4.68650455527377D-003
+  0.240004934893523        6.39200816662428D-003
+  0.300006168616904        8.22505484590828D-003
diff --git a/test/extract_physical/momentum_z_R2.tl b/test/extract_physical/momentum_z_R2.tl
new file mode 100644
index 0000000..0d39b2c
--- /dev/null
+++ b/test/extract_physical/momentum_z_R2.tl
@@ -0,0 +1,7 @@
+"momentum_z_R2.
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00000000000000D-002   7.40220781384088D-004
+  0.120000000000000        1.49090041767084D-003
+  0.180000000000000        2.26272109611529D-003
+  0.240000000000000        3.06681665079707D-003
+  0.300000000000000        3.91501182491871D-003
diff --git a/test/extract_physical/momentum_z_R2.ul b/test/extract_physical/momentum_z_R2.ul
new file mode 100644
index 0000000..31284d2
--- /dev/null
+++ b/test/extract_physical/momentum_z_R2.ul
@@ -0,0 +1,7 @@
+"momentum_z_R2.
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00012337233809D-002   7.40220781384088D-004
+  0.120002467446762        1.49090041767084D-003
+  0.180003701170143        2.26272109611529D-003
+  0.240004934893523        3.06681665079707D-003
+  0.300006168616904        3.91501182491871D-003
diff --git a/test/extract_physical/rsch_R1.tl b/test/extract_physical/rsch_R1.tl
new file mode 100644
index 0000000..b61efbd
--- /dev/null
+++ b/test/extract_physical/rsch_R1.tl
@@ -0,0 +1,7 @@
+"rsch_R1.tl    
+   0.00000000000000D+000   2.51624637785126     
+   6.00000000000000D-002   2.51628379065722     
+  0.120000000000000        2.51639673790065     
+  0.180000000000000        2.51658736254029     
+  0.240000000000000        2.51685929058176     
+  0.300000000000000        2.51721771239714     
diff --git a/test/extract_physical/rsch_R1.ul b/test/extract_physical/rsch_R1.ul
new file mode 100644
index 0000000..1bc63b8
--- /dev/null
+++ b/test/extract_physical/rsch_R1.ul
@@ -0,0 +1,7 @@
+"rsch_R1.ul    
+   0.00000000000000D+000   2.51624637785126     
+   6.00012337233809D-002   2.51628379065722     
+  0.120002467446762        2.51639673790065     
+  0.180003701170143        2.51658736254029     
+  0.240004934893523        2.51685929058176     
+  0.300006168616904        2.51721771239714     
diff --git a/test/extract_physical/rsch_R2.tl b/test/extract_physical/rsch_R2.tl
new file mode 100644
index 0000000..2c453f1
--- /dev/null
+++ b/test/extract_physical/rsch_R2.tl
@@ -0,0 +1,7 @@
+"rsch_R2.tl    
+   0.00000000000000D+000   3.00939932577571     
+   6.00000000000000D-002   3.00941438151010     
+  0.120000000000000        3.00945975854980     
+  0.180000000000000        3.00953609075923     
+  0.240000000000000        3.00964444914157     
+  0.300000000000000        3.00978636400156     
diff --git a/test/extract_physical/rsch_R2.ul b/test/extract_physical/rsch_R2.ul
new file mode 100644
index 0000000..371caf2
--- /dev/null
+++ b/test/extract_physical/rsch_R2.ul
@@ -0,0 +1,7 @@
+"rsch_R2.ul    
+   0.00000000000000D+000   3.00939932577571     
+   6.00012337233809D-002   3.00941438151010     
+  0.120002467446762        3.00945975854980     
+  0.180003701170143        3.00953609075923     
+  0.240004934893523        3.00964444914157     
+  0.300006168616904        3.00978636400156     
diff --git a/test/extract_physical/spin_x_R1.tl b/test/extract_physical/spin_x_R1.tl
new file mode 100644
index 0000000..3d05443
--- /dev/null
+++ b/test/extract_physical/spin_x_R1.tl
@@ -0,0 +1,7 @@
+"spin_x_R1.tl  
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00000000000000D-002   1.07995317097359D-008
+  0.120000000000000        2.15685453004565D-008
+  0.180000000000000        3.22985012382555D-008
+  0.240000000000000        4.30287011908392D-008
+  0.300000000000000        5.38795017757992D-008
diff --git a/test/extract_physical/spin_x_R1.ul b/test/extract_physical/spin_x_R1.ul
new file mode 100644
index 0000000..55c4f71
--- /dev/null
+++ b/test/extract_physical/spin_x_R1.ul
@@ -0,0 +1,7 @@
+"spin_x_R1.ul  
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00000000000000D-002   1.07995317097359D-008
+  0.120000000000000        2.15685453004565D-008
+  0.180000000000000        3.22985012382555D-008
+  0.240000000000000        4.30287011908392D-008
+  0.300000000000000        5.38795017757992D-008
diff --git a/test/extract_physical/spin_x_R2.tl b/test/extract_physical/spin_x_R2.tl
new file mode 100644
index 0000000..86967a8
--- /dev/null
+++ b/test/extract_physical/spin_x_R2.tl
@@ -0,0 +1,7 @@
+"spin_x_R2.tl  
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00000000000000D-002  -4.28444012444931D-009
+  0.120000000000000       -8.47345323446115D-009
+  0.180000000000000       -1.24708358206800D-008
+  0.240000000000000       -1.61790562144519D-008
+  0.300000000000000       -1.94992310935978D-008
diff --git a/test/extract_physical/spin_x_R2.ul b/test/extract_physical/spin_x_R2.ul
new file mode 100644
index 0000000..bc06b9c
--- /dev/null
+++ b/test/extract_physical/spin_x_R2.ul
@@ -0,0 +1,7 @@
+"spin_x_R2.ul  
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00000000000000D-002  -4.28444012444931D-009
+  0.120000000000000       -8.47345323446115D-009
+  0.180000000000000       -1.24708358206800D-008
+  0.240000000000000       -1.61790562144519D-008
+  0.300000000000000       -1.94992310935978D-008
diff --git a/test/extract_physical/spin_y_R1.tl b/test/extract_physical/spin_y_R1.tl
new file mode 100644
index 0000000..4485f11
--- /dev/null
+++ b/test/extract_physical/spin_y_R1.tl
@@ -0,0 +1,7 @@
+"spin_y_R1.tl  
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00000000000000D-002  -1.07995317097865D-008
+  0.120000000000000       -2.15685453003211D-008
+  0.180000000000000       -3.22985012372215D-008
+  0.240000000000000       -4.30287011887272D-008
+  0.300000000000000       -5.38795017726367D-008
diff --git a/test/extract_physical/spin_y_R1.ul b/test/extract_physical/spin_y_R1.ul
new file mode 100644
index 0000000..5fd108f
--- /dev/null
+++ b/test/extract_physical/spin_y_R1.ul
@@ -0,0 +1,7 @@
+"spin_y_R1.ul  
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00000000000000D-002  -1.07995317097865D-008
+  0.120000000000000       -2.15685453003211D-008
+  0.180000000000000       -3.22985012372215D-008
+  0.240000000000000       -4.30287011887272D-008
+  0.300000000000000       -5.38795017726367D-008
diff --git a/test/extract_physical/spin_y_R2.tl b/test/extract_physical/spin_y_R2.tl
new file mode 100644
index 0000000..d870f39
--- /dev/null
+++ b/test/extract_physical/spin_y_R2.tl
@@ -0,0 +1,7 @@
+"spin_y_R2.tl  
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00000000000000D-002   4.28444012441639D-009
+  0.120000000000000        8.47345323442142D-009
+  0.180000000000000        1.24708358206241D-008
+  0.240000000000000        1.61790562143525D-008
+  0.300000000000000        1.94992310935669D-008
diff --git a/test/extract_physical/spin_y_R2.ul b/test/extract_physical/spin_y_R2.ul
new file mode 100644
index 0000000..83a2716
--- /dev/null
+++ b/test/extract_physical/spin_y_R2.ul
@@ -0,0 +1,7 @@
+"spin_y_R2.ul  
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00000000000000D-002   4.28444012441639D-009
+  0.120000000000000        8.47345323442142D-009
+  0.180000000000000        1.24708358206241D-008
+  0.240000000000000        1.61790562143525D-008
+  0.300000000000000        1.94992310935669D-008
diff --git a/test/extract_physical/spin_z_R1.tl b/test/extract_physical/spin_z_R1.tl
new file mode 100644
index 0000000..b4e19e5
--- /dev/null
+++ b/test/extract_physical/spin_z_R1.tl
@@ -0,0 +1,7 @@
+"spin_z_R1.tl  
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00000000000000D-002   4.76755153562323D-020
+  0.120000000000000       -1.32461207639645D-019
+  0.180000000000000       -1.00712928011351D-018
+  0.240000000000000       -2.18426444150420D-018
+  0.300000000000000       -3.09606680651214D-018
diff --git a/test/extract_physical/spin_z_R1.ul b/test/extract_physical/spin_z_R1.ul
new file mode 100644
index 0000000..ecb6fb0
--- /dev/null
+++ b/test/extract_physical/spin_z_R1.ul
@@ -0,0 +1,7 @@
+"spin_z_R1.ul  
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00000000000000D-002   4.76755153562323D-020
+  0.120000000000000       -1.32461207639645D-019
+  0.180000000000000       -1.00712928011351D-018
+  0.240000000000000       -2.18426444150420D-018
+  0.300000000000000       -3.09606680651214D-018
diff --git a/test/extract_physical/spin_z_R2.tl b/test/extract_physical/spin_z_R2.tl
new file mode 100644
index 0000000..a2b1688
--- /dev/null
+++ b/test/extract_physical/spin_z_R2.tl
@@ -0,0 +1,7 @@
+"spin_z_R2.tl  
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00000000000000D-002   3.41875628705859D-020
+  0.120000000000000        4.70946098859468D-020
+  0.180000000000000        7.85374632740441D-020
+  0.240000000000000        4.60707880322157D-020
+  0.300000000000000        5.26509184760700D-020
diff --git a/test/extract_physical/spin_z_R2.ul b/test/extract_physical/spin_z_R2.ul
new file mode 100644
index 0000000..cf0648a
--- /dev/null
+++ b/test/extract_physical/spin_z_R2.ul
@@ -0,0 +1,7 @@
+"spin_z_R2.ul  
+   0.00000000000000D+000   0.00000000000000D+000
+   6.00000000000000D-002   3.41875628705859D-020
+  0.120000000000000        4.70946098859468D-020
+  0.180000000000000        7.85374632740441D-020
+  0.240000000000000        4.60707880322157D-020
+  0.300000000000000        5.26509184760700D-020
diff --git a/test/test_extract.par b/test/test_extract.par
index 1463283..abcfb21 100644
--- a/test/test_extract.par
+++ b/test/test_extract.par
@@ -46,5 +46,9 @@ extract::m_mode     = 0
 extract::detector1 = 2.5 
 extract::detector2 = 3.0
 
+extract::doADMmass = "yes"
+extract::do_momentum = "yes"
+extract::do_spin = "yes"
+
 # do not use the new scheme for names of output files
 IO::new_filename_scheme = "no"
\ No newline at end of file
diff --git a/test/test_extract/test_extract.par b/test/test_extract/test_extract.par
index abcfb21..1a20d93 100644
--- a/test/test_extract/test_extract.par
+++ b/test/test_extract/test_extract.par
@@ -19,11 +19,11 @@ IO::out_every       = 1
 
 admbase::evolution_method ="adm"
 admbase::initial_data     = "schwarzschild" 
-admbase::metric_type = "static conformal"
+admbase::metric_type = "physical"
 admbase::lapse_evolution_method          = "geodesic"
 idanalyticbh::mass	   = 1 
 
-IO::out_dir          = "test_extract"
+IO::out_dir          = "extract_physical"
 io::out_fileinfo    = "none"
 ioascii::out1d_vars = "admbase::alp"
 
@@ -51,4 +51,4 @@ extract::do_momentum = "yes"
 extract::do_spin = "yes"
 
 # do not use the new scheme for names of output files
-IO::new_filename_scheme = "no"
\ No newline at end of file
+IO::new_filename_scheme = "no"