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

[cott at tapir.caltech.edu]

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