diff options
author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-01-14 14:23:38 +0000 |
---|---|---|
committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-01-14 14:23:38 +0000 |
commit | 3d793e03c6cd0a1e452b0b4d9f1a8bff38e5569e (patch) | |
tree | 1f6bbe0664867181cabd402ae3e883fb1f197e53 /src/GRHydro_Tmunu.F90 | |
parent | 5199c61af4721444a26c4e7c1c30d3b8fdedf54c (diff) |
GRHydro: Improved Tmunu damping.
* faster falloff to make sure to damp to atmosphere at
Tmunu_damping_radius_max
From: Christian Reisswig <reisswig@scriwalker.(none)>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@457 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_Tmunu.F90')
-rw-r--r-- | src/GRHydro_Tmunu.F90 | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/src/GRHydro_Tmunu.F90 b/src/GRHydro_Tmunu.F90 index bc7724c..6a2eee3 100644 --- a/src/GRHydro_Tmunu.F90 +++ b/src/GRHydro_Tmunu.F90 @@ -122,7 +122,7 @@ if ((Tmunu_damping_radius_min .gt. 0) .and. (r(i,j,k) .gt. Tmunu_damping_radius_min)) then ! 0.5 * (1.0 - tanh(4.0*(x-x0)/sigma0)) if (r(i,j,k) .lt. Tmunu_damping_radius_max) then - dampfac = 0.5 * (1.0 - tanh(4.0*(r(i,j,k)-0.5*(Tmunu_damping_radius_max+Tmunu_damping_radius_min))/(Tmunu_damping_radius_max-Tmunu_damping_radius_min))) + dampfac = 0.5d0 * (1.0d0 - tanh((8.0d0*r(i,j,k)-4.0d0*(Tmunu_damping_radius_max+Tmunu_damping_radius_min))/(Tmunu_damping_radius_max-Tmunu_damping_radius_min))) else dampfac = 0.0 continue ! no need to add anything to Tmunu at the current point (it's zero anyway!) |