User: bmundim Date: 2011/08/04 12:05 PM

Modified: /trunk/ param.ccl /trunk/src/ GRHydro_SourceM.F90

Log: RIT GRMHD dev:

Change legacy divergence cleaning parameters ch_dc and cp_dc, introduced for the special relativistic case, to kap_dc, the one introduced in the covariant formulation (Liebling et al.).

File Changes:

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

File [modified]: GRHydro_SourceM.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_SourceM.F90 2011-08-03 18:00:47 UTC (rev 262) +++ trunk/src/GRHydro_SourceM.F90 2011-08-04 17:05:23 UTC (rev 263) @@ -422,7 +422,7 @@ taurhs(i,j,k) = alp(i,j,k)*sqrtdet*tau_source

if(clean_divergence.ne.0) then - psidcrhs(i,j,k) = -1.d0 * ( ch_dc*ch_dc/cp_dc/cp_dc*alp(i,j,k) + & + psidcrhs(i,j,k) = -1.d0 * (kap_dc*alp(i,j,k) + & dx_betax + dy_betay + dz_betaz ) * psidc(i,j,k) + & Bconsx(i,j,k) * ( dx_alp - half*alp(i,j,k) * & ( uxx*dx_gxx + uyy*dx_gyy + uzz*dx_gzz + 2.d0*uxy*dx_gxy + &

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

File [modified]: param.ccl Delta lines: +2 -7 =================================================================== --- trunk/param.ccl 2011-08-03 18:00:47 UTC (rev 262) +++ trunk/param.ccl 2011-08-04 17:05:23 UTC (rev 263) @@ -506,16 +506,11 @@ { } "no"

-CCTK_REAL ch_dc "The c_h parameter for divergence cleaning" +CCTK_REAL kap_dc "The kap parameter for divergence cleaning" { - 0:* :: "Any value, but one is preferred" + 0:* :: "Any non-negative value, but 1/144 to 1 seems preferred" } 1.0

-CCTK_REAL cp_dc "The c_p parameter for divergence cleaning" -{ - 0:* :: "Any value, but one to 12 is preferred" -} 1.0 - KEYWORD psidcspeed "Which speed to set for psidc" { "char speed" :: "Based on the characteristic speeds"