User: rhaas Date: 2013/07/06 01:10 PM

Modified: /trunk/src/ GRHydro_CalcBcom.F90

Log: GRHydro: Don't compute Bcom twice.

From: Christian Reisswig reisswig@tapir.caltech.edu

File Changes:

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

File [modified]: GRHydro_CalcBcom.F90 Delta lines: +5 -0 =================================================================== --- trunk/src/GRHydro_CalcBcom.F90 2013-07-06 18:10:32 UTC (rev 546) +++ trunk/src/GRHydro_CalcBcom.F90 2013-07-06 18:10:37 UTC (rev 547) @@ -40,6 +40,11 @@ call CCTK_WARN(0,"No storage for GRHydro::bcom") end if

+! Only compute if Tmunu is not computed! (Bcom is also computed with Tmunu) + if (stress_energy_state .ne. 0) then + return + endif + !$OMP PARALLEL DO PRIVATE(i,j,k,velxlow, velylow, velzlow,& !$OMP Bvecxlow,Bvecylow,Bveczlow, bdotv,dum1,dum2,b2,bxlow,bylow,bzlow)