From 09ce7a16b9ede9ff720607cf87d2ade17f217717 Mon Sep 17 00:00:00 2001
From: knarf <knarf@850bcc8b-0e4f-0410-8c26-8d28fbf1eda9>
Date: Mon, 20 Jun 2005 10:39:59 +0000
Subject: -new parameter ah_max_epsdis -bugfix for ah not at (0,0,0)

git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/Dissipation/trunk@13 850bcc8b-0e4f-0410-8c26-8d28fbf1eda9
---
 src/setup_epsdis.c | 9 +++++++--
 1 file changed, 7 insertions(+), 2 deletions(-)

(limited to 'src')

diff --git a/src/setup_epsdis.c b/src/setup_epsdis.c
index 5acb649..6239e6c 100644
--- a/src/setup_epsdis.c
+++ b/src/setup_epsdis.c
@@ -154,7 +154,8 @@ setup_epsdis (CCTK_ARGUMENTS)
         if (verbose) {
           CCTK_VInfo(CCTK_THORNSTRING,
                      "Invalid Surface: s=%d, sf_va=%d, surf_no=%d, ah_no=%d",
-                     s,sf_valid[surface_number[s]],surface_number[s],horizon_number[s]);
+                     s,sf_valid[surface_number[s]],
+                     surface_number[s],horizon_number[s]);
         }
         continue;
       }
@@ -217,9 +218,13 @@ setup_epsdis (CCTK_ARGUMENTS)
       assert(!ierr);
 
       for (i=0;i<npts;i++) {
-        radp=sqrt(xa[i]*xa[i]+ya[i]*ya[i]+za[i]*za[i]);
+        radp=sqrt((xa[i]-odx)*(xa[i]-odx)+(ya[i]-ody)*(ya[i]-ody)+
+                  (za[i]-odz)*(za[i]-odz));
         if (radp<=rads[i]) {
           epsdisA[inds[i]]=epsdis+ ah_slope*(rads[i]-radp);
+          if (epsdisA[inds[i]] > ah_max_epsdis) {
+            epsdisA[inds[i]] = ah_max_epsdis;
+          }
         }
       }
 
-- 
cgit v1.2.3