User: rhaas Date: 2013/01/14 08:23 AM

Modified: /trunk/src/ GRHydro_CalcBcom.F90

Log: GRHydro: use distinct dummys for calc_vlow_blow

also check actual timelevels information for storage of bcom

From: Roland Haas roland.haas@physics.gatech.edu

File Changes:

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

File [modified]: GRHydro_CalcBcom.F90 Delta lines: +7 -5 =================================================================== --- trunk/src/GRHydro_CalcBcom.F90 2013-01-14 14:23:47 UTC (rev 462) +++ trunk/src/GRHydro_CalcBcom.F90 2013-01-14 14:23:48 UTC (rev 463) @@ -30,12 +30,14 @@

CCTK_REAL :: velxlow, velylow, velzlow CCTK_REAL :: Bvecxlow,Bvecylow,Bveczlow - CCTK_REAL :: bdotv,b2,bxlow,bylow,bzlow,dum + CCTK_REAL :: bdotv,b2,bxlow,bylow,bzlow,dum1,dum2 CCTK_INT :: i,j,k CCTK_INT :: GRHydro_UseGeneralCoordinates + integer :: timelevels

- if(stress_energy_state.ne.0) then - return + call CCTK_ActiveTimeLevels(timelevels, cctkGH, "GRHydro::bcom") + if(timelevels.eq.0) then + call CCTK_WARN(0,"No storage for GRHydro::bcom") end if

if (GRHydro_UseGeneralCoordinates(cctkGH).ne.0) then @@ -43,7 +45,7 @@ end if

!$OMP PARALLEL DO PRIVATE(i,j,k,velxlow, velylow, velzlow,& - !$OMP Bvecxlow,Bvecylow,Bveczlow, bdotv,dum,b2,bxlow,bylow,bzlow) + !$OMP Bvecxlow,Bvecylow,Bveczlow, bdotv,dum1,dum2,b2,bxlow,bylow,bzlow)

do k = 1, cctk_lsh(3) do j = 1, cctk_lsh(2) @@ -53,7 +55,7 @@ gyy(i,j,k),gyz(i,j,k),gzz(i,j,k), & velx(i,j,k),vely(i,j,k),velz(i,j,k),Bvecx(i,j,k),Bvecy(i,j,k),Bvecz(i,j,k), & velxlow,velylow,velzlow,Bvecxlow,Bvecylow,Bveczlow, & - bdotv,b2,dum,dum,bxlow,bylow,bzlow) + bdotv,b2,dum1,dum2,bxlow,bylow,bzlow)

bcom_sq(i,j,k) = b2 bcom0(i,j,k) = w_lorentz(i,j,k)*bdotv/alp(i,j,k)