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@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!)