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

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2014/04/15 02:49 PM

Modified:
 /trunk/src/
  GRHydro_CalcUpdate.F90

Log:
 GRHydro: ignore transport_constraints unless evolve_mhd is set

File Changes:

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

File [modified]: GRHydro_CalcUpdate.F90
Delta lines: +3 -2
===================================================================
--- trunk/src/GRHydro_CalcUpdate.F90	2014-04-15 19:49:40 UTC (rev 629)
+++ trunk/src/GRHydro_CalcUpdate.F90	2014-04-15 19:49:43 UTC (rev 630)
@@ -242,7 +242,7 @@
   
   else if (CCTK_EQUALS(method_type, "Flux split FD")) then
 
-    if (transport_constraints .ne. 0) then
+    if (evolve_mhd.ne.0 .and. transport_constraints .ne. 0) then
       call CCTK_WARN(0, "Not supported")
     end if
 
@@ -296,7 +296,8 @@
       
   end if
 
-  if (transport_constraints.ne.0 .and. flux_direction.eq.1) then ! HACK: x direction is last
+  if (evolve_mhd.ne.0 .and. transport_constraints.ne.0 .and. &
+      flux_direction.eq.1) then ! HACK: x direction is last
     ! FIXME: I think one could wrap all of this into a single do loop and remove the
     !        Evec storage
     !$OMP PARALLEL DO PRIVATE(i,j,k)