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

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2012/10/25 03:35 PM

Modified:
 /trunk/
  schedule.ccl
 /trunk/src/
  GRHydro_Boundaries.F90

Log:
 GRHydro: undo revision 426 since it causes test failures since poltytropes use
 pressure as an initial guess

File Changes:

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

File [modified]: GRHydro_Boundaries.F90
Delta lines: +8 -12
===================================================================
--- trunk/src/GRHydro_Boundaries.F90	2012-10-25 16:37:02 UTC (rev 427)
+++ trunk/src/GRHydro_Boundaries.F90	2012-10-25 20:35:02 UTC (rev 428)
@@ -513,15 +513,13 @@
 
   if (CCTK_EQUALS(bound,"flat")) then
     
-    ! Lorentz factor is not needed since it is computed in Con2Prim which runs after any sync!
-    !ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-    !      "HydroBase::w_lorentz", "Flat")
     ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+          "HydroBase::w_lorentz", "Flat")
+    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
           "HydroBase::rho", "Flat")
-    ! Pressure is not needed since it is computed in Con2Prim which runs after any sync!
-    !ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-    !      "HydroBase::press", "Flat")
     ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+          "HydroBase::press", "Flat")
+    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
           "HydroBase::eps", "Flat")
     if (general_coordinates .ne. 0) then
        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
@@ -561,15 +559,13 @@
   endif
 
   if (CCTK_EQUALS(bound,"none")) then      
-    ! Lorentz factor is not needed since it is computed in Con2Prim which runs after any sync!
-    !ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-    !      "HydroBase::w_lorentz", "None")
     ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+          "HydroBase::w_lorentz", "None")
+    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
           "HydroBase::rho", "None")
-    ! Pressure is not needed since it is computed in Con2Prim which runs after any sync!
-    !ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-    !      "HydroBase::press", "None")
     ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+          "HydroBase::press", "None")
+    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
           "HydroBase::eps", "None")
     if (general_coordinates .ne. 0) then
        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &

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

File [modified]: schedule.ccl
Delta lines: +2 -2
===================================================================
--- trunk/schedule.ccl	2012-10-25 16:37:02 UTC (rev 427)
+++ trunk/schedule.ccl	2012-10-25 20:35:02 UTC (rev 428)
@@ -1150,8 +1150,8 @@
       {
 	 LANG: Fortran
 	 OPTIONS: LEVEL
-#         SYNC: HydroBase::w_lorentz   # don't need that since computed in Con2Prim which runs in postregrid!
-#         SYNC: HydroBase::press       # don't need that since computed in Con2Prim which runs in postregrid!
+         SYNC: HydroBase::w_lorentz
+         SYNC: HydroBase::press
 #         SYNC: HydroBase::vel
 #         SYNC: HydroBase::Bvec
          SYNC: hydrobase::entropy