[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 457)

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2013/01/14 08:23 AM

Modified:
 /trunk/src/
  GRHydro_Tmunu.F90

Log:
 GRHydro: Improved Tmunu damping.
 
 * faster falloff to make sure to damp to atmosphere at
   Tmunu_damping_radius_max
 
 From: Christian Reisswig <reisswig at scriwalker.(none)>

File Changes:

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

File [modified]: GRHydro_Tmunu.F90
Delta lines: +1 -1
===================================================================
--- trunk/src/GRHydro_Tmunu.F90	2013-01-14 14:23:37 UTC (rev 456)
+++ trunk/src/GRHydro_Tmunu.F90	2013-01-14 14:23:38 UTC (rev 457)
@@ -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!)