diff options
author | knarf <knarf@850bcc8b-0e4f-0410-8c26-8d28fbf1eda9> | 2005-06-20 10:39:59 +0000 |
---|---|---|
committer | knarf <knarf@850bcc8b-0e4f-0410-8c26-8d28fbf1eda9> | 2005-06-20 10:39:59 +0000 |
commit | 09ce7a16b9ede9ff720607cf87d2ade17f217717 (patch) | |
tree | ebc43b564f836a33140c04739ac77a4c6087c913 /src | |
parent | 92c17f4c2d24acdd48915f47c99290d5238c61a5 (diff) |
-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
Diffstat (limited to 'src')
-rw-r--r-- | src/setup_epsdis.c | 9 |
1 files changed, 7 insertions, 2 deletions
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; + } } } |