User: cott Date: 2010/12/04 03:29 PM

Modified: /trunk/src/ GRHydro_Prim2Con.F90

Log: * correct Prim2Con for Y_e

File Changes:

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

File [modified]: GRHydro_Prim2Con.F90 Delta lines: +12 -12 =================================================================== --- trunk/src/GRHydro_Prim2Con.F90 2010-12-04 07:36:15 UTC (rev 187) +++ trunk/src/GRHydro_Prim2Con.F90 2010-12-04 21:29:56 UTC (rev 188) @@ -319,18 +319,6 @@ 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,7 +361,19 @@ !$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