- add parameter 'tolerance' to be able to control this without recompile

- correct intent of several function arguments


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/DistortedBHIVP/trunk@43 971fb155-194f-0410-9daf-e2eca44e59f5

3 files changed, 18 insertions, 13 deletions

diff --git a/param.ccl b/param.ccl
index 225bfba..db18e45 100644
--- a/param.ccl
+++ b/param.ccl
@@ -73,4 +73,7 @@ INT interpolation_order "Order for interpolation" STEERABLE = ALWAYS
   1:3 :: "Choose between first, second, and third-order"
 } 1
 
-
+REAL tolerance "Tolerance"
+{
+ 0:* :: "Tolerance"
+} 1.e-9
diff --git a/src/DistortedBHIVP.F b/src/DistortedBHIVP.F
index 50139ed..d8bbdca 100644
--- a/src/DistortedBHIVP.F
+++ b/src/DistortedBHIVP.F
@@ -68,7 +68,7 @@ c@@   */
       CCTK_POINTER, dimension(10) :: in_arrays, out_arrays
       CCTK_INT,     dimension(10), parameter :: type_codes = CCTK_VARIABLE_REAL
 
-      dbh_eps = 1.0d-9
+      dbh_eps = tolerance
 
 c Check if we should create and store conformal factor stuff
 
diff --git a/src/Stab3d.F b/src/Stab3d.F
index aad61b6..6215ac1 100644
--- a/src/Stab3d.F
+++ b/src/Stab3d.F
@@ -9,7 +9,7 @@ c
       integer,intent(in) :: im,jm,km
       real*8,intent(inout) :: cc(im,jm,km),cn(im,jm,km),cs(im,jm,km),
      $     ce(im,jm,km),cw(im,jm,km),ct(im,jm,km),cb(im,jm,km)
-      real*8,intent(out) :: eps
+      real*8,intent(in) :: eps
       real*8,intent(out) :: rmax
       real*8 :: u(im,jm,km),rhs(im,jm,km)
 c     Local variable
@@ -88,11 +88,12 @@ c
 c
       integer,intent(in) :: im,jm,km
       integer,intent(in) :: ncyc
-      real*8,intent(out) :: cc(im*jm*km),cn(im*jm*km)
-      real*8,intent(out) :: cs(im*jm*km),ce(im*jm*km)
-      real*8,intent(out) :: cw(im*jm*km),ct(im*jm*km)
-      real*8,intent(out) :: cb(im*jm*km)
-      real*8,intent(out) :: tol,rnorm
+      real*8,intent(in) :: cc(im*jm*km),cn(im*jm*km)
+      real*8,intent(in) :: cs(im*jm*km),ce(im*jm*km)
+      real*8,intent(in) :: cw(im*jm*km),ct(im*jm*km)
+      real*8,intent(in) :: cb(im*jm*km)
+      real*8,intent(in)  :: tol
+      real*8,intent(out) :: rnorm
       integer,intent(out) :: ier
       real*8 x(im*jm*km),r(im*jm*km)
 c     Local variables
@@ -257,11 +258,12 @@ c
       implicit none
 c
       integer,intent(in) :: im,jm,km
-      real*8,intent(out) :: cc(im*jm*km),cn(im*jm*km)
-      real*8,intent(out) :: cs(im*jm*km),ce(im*jm*km)
-      real*8,intent(out) :: cw(im*jm*km),ct(im*jm*km)
-      real*8,intent(out) :: cb(im*jm*km)
-      real*8 :: x(im*jm*km), y(im*jm*km)
+      real*8,intent(in) :: cc(im*jm*km),cn(im*jm*km)
+      real*8,intent(in) :: cs(im*jm*km),ce(im*jm*km)
+      real*8,intent(in) :: cw(im*jm*km),ct(im*jm*km)
+      real*8,intent(in) :: cb(im*jm*km)
+      real*8,intent(in) :: x(im*jm*km)
+      real*8,intent(out) :: y(im*jm*km)
 c     Local variables
       integer :: i, j, k
 *