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)