diff options
| author | cott <cott@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2010-12-04 21:29:56 +0000 |
|---|---|---|
| committer | cott <cott@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2010-12-04 21:29:56 +0000 |
| commit | 9b147a8c397b5bb7a0c1e8c1537f870474c71a74 (patch) | |
| tree | 13f7a74dd1aa83f5b72bab69fb7825a605bf5ff3 /src/GRHydro_Prim2Con.F90 | |
| parent | 658e7648d888ae72d7b52a297e3bc11b3bcf6c55 (diff) | |
* correct Prim2Con for Y_e
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@188 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_Prim2Con.F90')
| -rw-r--r-- | src/GRHydro_Prim2Con.F90 | 24 |
1 files changed, 12 insertions, 12 deletions
diff --git a/src/GRHydro_Prim2Con.F90 b/src/GRHydro_Prim2Con.F90 index aef6f3a..dffeba1 100644 --- a/src/GRHydro_Prim2Con.F90 +++ b/src/GRHydro_Prim2Con.F90 @@ -319,18 +319,6 @@ subroutine Primitive2ConservativeCells(CCTK_ARGUMENTS) CCTK_REAL :: det character(len=512) :: warnline - if(evolve_Y_e.ne.0) then - !$OMP PARALLEL DO PRIVATE(i, j) - do k = GRHydro_stencil,cctk_lsh(3)-GRHydro_stencil+1 - do j = GRHydro_stencil,cctk_lsh(2)-GRHydro_stencil+1 - do i = GRHydro_stencil,cctk_lsh(1)-GRHydro_stencil+1 - Y_e_con(i,j,k) = Y_e(i,j,k) * dens(i,j,k) - enddo - enddo - enddo - !$OMP END PARALLEL DO - endif - if(evolve_temper.ne.1) then !$OMP PARALLEL DO PRIVATE(i, j) do k = GRHydro_stencil,cctk_lsh(3)-GRHydro_stencil+1 @@ -373,6 +361,18 @@ subroutine Primitive2ConservativeCells(CCTK_ARGUMENTS) !$OMP END PARALLEL DO endif + if(evolve_Y_e.ne.0) then + !$OMP PARALLEL DO PRIVATE(i, j) + do k = GRHydro_stencil,cctk_lsh(3)-GRHydro_stencil+1 + do j = GRHydro_stencil,cctk_lsh(2)-GRHydro_stencil+1 + do i = GRHydro_stencil,cctk_lsh(1)-GRHydro_stencil+1 + Y_e_con(i,j,k) = Y_e(i,j,k) * dens(i,j,k) + enddo + enddo + enddo + !$OMP END PARALLEL DO + endif + end subroutine Primitive2ConservativeCells |