[Commits] [svn:einsteintoolkit] DistortedBHIVP/trunk/ (Rev. 43)

[knarf at cct.lsu.edu]

User: knarf
Date: 2010/05/24 01:02 PM

Modified:
 /trunk/
  param.ccl
 /trunk/src/
  DistortedBHIVP.F, Stab3d.F

Log:
 - add parameter 'tolerance' to be able to control this without recompile
 - correct intent of several function arguments

File Changes:

Directory: /trunk/src/
======================

File [modified]: DistortedBHIVP.F
Delta lines: +1 -1
===================================================================
--- trunk/src/DistortedBHIVP.F	2010-05-24 18:02:08 UTC (rev 42)
+++ trunk/src/DistortedBHIVP.F	2010-05-24 18:02:56 UTC (rev 43)
@@ -68,7 +68,7 @@
       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
 

File [modified]: Stab3d.F
Delta lines: +13 -11
===================================================================
--- trunk/src/Stab3d.F	2010-05-24 18:02:08 UTC (rev 42)
+++ trunk/src/Stab3d.F	2010-05-24 18:02:56 UTC (rev 43)
@@ -9,7 +9,7 @@
       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
       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 @@
       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
 *

Directory: /trunk/
==================

File [modified]: param.ccl
Delta lines: +4 -1
===================================================================
--- trunk/param.ccl	2010-05-24 18:02:08 UTC (rev 42)
+++ trunk/param.ccl	2010-05-24 18:02:56 UTC (rev 43)
@@ -73,4 +73,7 @@
   1:3 :: "Choose between first, second, and third-order"
 } 1
 
-
+REAL tolerance "Tolerance"
+{
+ 0:* :: "Tolerance"
+} 1.e-9