[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 371)
rhaas at tapir.caltech.edu
rhaas at tapir.caltech.edu
Fri Jun 22 15:37:45 CDT 2012
User: rhaas
Date: 2012/06/22 03:37 PM
Modified:
/trunk/src/
GRHydro_ParamCheck.F90
Log:
Revert "Revert "GRHydro: do not access whichpsidcspeed in ParamCheck unless evolve_MHD""
Ahem. I had reverted the wrong commit. This one was blameless.
This reverts commit 0b5a561872841b50ea32da2b665a04b2984a8c16.
File Changes:
Directory: /trunk/src/
======================
File [modified]: GRHydro_ParamCheck.F90
Delta lines: +5 -1
===================================================================
--- trunk/src/GRHydro_ParamCheck.F90 2012-06-22 02:29:21 UTC (rev 370)
+++ trunk/src/GRHydro_ParamCheck.F90 2012-06-22 20:37:44 UTC (rev 371)
@@ -126,7 +126,7 @@
end if
- if(clean_divergence.ne.0) then
+ if(evolve_MHD.ne.0 .and. clean_divergence.ne.0) then
if (CCTK_EQUALS(psidcspeed,"char speed")) then
whichpsidcspeed = 0
else if (CCTK_EQUALS(psidcspeed,"light speed")) then
@@ -138,6 +138,10 @@
end if
end if
+ if(clean_divergence.eq.0 .and. evolve_MHD.ne.0) then
+ call CCTK_WARN(1, "You activated divergence cleaning but do not evolve a magnetic field. At best, nothing should happen, be preapered for a segfault otherwise.")
+ end if
+
if (Tmunu_damping_radius_min .gt. Tmunu_damping_radius_max) then
call CCTK_WARN(0, "Minimum damping radius is greater than maximum damping radius!")
end if
More information about the Commits
mailing list