diff options
| author | bmundim <bmundim@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-04-28 21:33:20 +0000 |
|---|---|---|
| committer | bmundim <bmundim@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-04-28 21:33:20 +0000 |
| commit | 80d7ee804818a46db34f102c5f08b913544a34c2 (patch) | |
| tree | a1b3781cf1bf8f22c4cf5cdb92e0bc41bc87aeef /src/GRHydro_HLLEM.F90 | |
| parent | 67ebaf5166748ff5aeed70883df73163047bad6d (diff) | |
JOsh's fixes.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@246 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_HLLEM.F90')
| -rw-r--r-- | src/GRHydro_HLLEM.F90 | 11 |
1 files changed, 2 insertions, 9 deletions
diff --git a/src/GRHydro_HLLEM.F90 b/src/GRHydro_HLLEM.F90 index 8cb2e49..e841450 100644 --- a/src/GRHydro_HLLEM.F90 +++ b/src/GRHydro_HLLEM.F90 @@ -111,8 +111,6 @@ subroutine GRHydro_HLLEM(CCTK_ARGUMENTS) cons_m(7) = Bconsyminus(i+xoffset,j+yoffset,k+zoffset) cons_m(8) = Bconszminus(i+xoffset,j+yoffset,k+zoffset) -!!$ if(i.eq.75.and.j.eq.5.and.k.eq.5)write(6,*)'HLLEM0:',cons_m(6),cons_p(6),cons_m(7),cons_p(7),cons_m(8),cons_p(8) - prim_p(1) = rhoplus(i,j,k) prim_p(2) = velxplus(i,j,k) prim_p(3) = velyplus(i,j,k) @@ -353,6 +351,7 @@ subroutine GRHydro_HLLEM(CCTK_ARGUMENTS) fminus(7),fminus(8),fminus(6),& vytm,vztm,vxtm,pressstarm,bylowm,bzlowm,bxlowm,ab0m,wm, & avg_det,avg_alp,avg_beta) + if(clean_divergence.ne.0) then fminus(6)=fminus(6) + avg_alp*sdet*uxyh*psidcm - cons_m(7)*avg_betax fminus(7)=fminus(7) + avg_alp*sdet*uyyh*psidcm - cons_m(7)*avg_betay @@ -386,6 +385,7 @@ subroutine GRHydro_HLLEM(CCTK_ARGUMENTS) fminus(8),fminus(6),fminus(7), & vztm,vxtm,vytm,pressstarm,bzlowm,bxlowm,bylowm,ab0m,wm, & avg_det,avg_alp,avg_beta) + if(clean_divergence.ne.0) then fminus(6)=fminus(6) + avg_alp*sdet*uxzh*psidcm - cons_m(8)*avg_betax fminus(7)=fminus(7) + avg_alp*sdet*uyzh*psidcm - cons_m(8)*avg_betay @@ -414,8 +414,6 @@ subroutine GRHydro_HLLEM(CCTK_ARGUMENTS) charpm = charmax - charmin -!!$ if(i.eq.75.and.j.eq.5.and.k.eq.5)write(6,*)'HLLEM:',fplus(6),fminus(6),fplus(8),fminus(8),charmax,charmin,cons_p(6),cons_m(6) - !!$ Calculate flux by standard formula do m = 1,8 @@ -427,9 +425,6 @@ subroutine GRHydro_HLLEM(CCTK_ARGUMENTS) end do -!!$ if(i.eq.75.and.j.eq.5.and.k.eq.5)write(6,*)'HLLEM2:',fplus(6),fminus(6),fplus(8),fminus(8),charmax,charmin,qdiff(6),qdiff(8),f1(6),f1(8) - - if(clean_divergence.ne.0) then psidcdiff = psidcm - psidcp @@ -456,8 +451,6 @@ subroutine GRHydro_HLLEM(CCTK_ARGUMENTS) Bconsyflux(i, j, k) = f1(7) Bconszflux(i, j, k) = f1(8) -!!$ if(i.eq.75.and.j.eq.5.and.k.eq.5)write(6,*)'HLLEM2:',fplus(6),fminus(6),fplus(8),fminus(8),f1(6),f1(8) - if(clean_divergence.ne.0) then psidcflux(i,j,k) = psidcf endif |