[Commits] [svn:einsteintoolkit] GRHydro/branches/hot_and_MHD_temp_dev/src/ (Rev. 187)

cott at tapir.caltech.edu cott at tapir.caltech.edu
Sat Dec 4 01:36:16 CST 2010


User: cott
Date: 2010/12/04 01:36 AM

Modified:
 /branches/hot_and_MHD_temp_dev/src/
  GRHydro_Prim2Con.F90

Log:
 * fix small bug with Y_e in Prim2ConCells

File Changes:

Directory: /branches/hot_and_MHD_temp_dev/src/
==============================================

File [modified]: GRHydro_Prim2Con.F90
Delta lines: +13 -12
===================================================================
--- branches/hot_and_MHD_temp_dev/src/GRHydro_Prim2Con.F90	2010-11-26 20:02:43 UTC (rev 186)
+++ branches/hot_and_MHD_temp_dev/src/GRHydro_Prim2Con.F90	2010-12-04 07:36:15 UTC (rev 187)
@@ -49,6 +49,7 @@
   CCTK_REAL :: xtemp
 
   if(evolve_temper.ne.1) then
+     
      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
@@ -319,18 +320,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 +362,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
 
 



More information about the Commits mailing list