diff options
| -rw-r--r-- | schedule.ccl | 4 | ||||
| -rw-r--r-- | src/GRHydro_Boundaries.F90 | 20 |
2 files changed, 10 insertions, 14 deletions
diff --git a/schedule.ccl b/schedule.ccl index 1068ed9..cc5ce5a 100644 --- a/schedule.ccl +++ b/schedule.ccl @@ -1150,8 +1150,8 @@ else { 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 diff --git a/src/GRHydro_Boundaries.F90 b/src/GRHydro_Boundaries.F90 index d32dd76..c04994c 100644 --- a/src/GRHydro_Boundaries.F90 +++ b/src/GRHydro_Boundaries.F90 @@ -513,14 +513,12 @@ subroutine GRHydro_SelectPrimitiveBoundaries(CCTK_ARGUMENTS) 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 @@ -561,14 +559,12 @@ subroutine GRHydro_SelectPrimitiveBoundaries(CCTK_ARGUMENTS) 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 |