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