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

[bcmsma at astro.rit.edu]

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"