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

[rhaas at tapir.caltech.edu]

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 at 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)