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

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2012/06/20 11:19 AM

Modified:
 /trunk/src/
  GRHydro_ParamCheck.F90

Log:
 GRHydro: do not access whichpsidcspeed in ParamCheck unless evolve_MHD

File Changes:

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

File [modified]: GRHydro_ParamCheck.F90
Delta lines: +5 -1
===================================================================
--- trunk/src/GRHydro_ParamCheck.F90	2012-06-20 16:18:04 UTC (rev 359)
+++ trunk/src/GRHydro_ParamCheck.F90	2012-06-20 16:19:10 UTC (rev 360)
@@ -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