[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 341)
rhaas at tapir.caltech.edu
rhaas at tapir.caltech.edu
Tue Jun 5 15:28:24 CDT 2012
User: rhaas
Date: 2012/06/05 03:28 PM
Modified:
/trunk/src/
GRHydro_Boundaries.F90
Log:
GRHydro: sync those primitives which are needed for initial Con2Prim guesses
Patch by Christian Reisswig.
File Changes:
Directory: /trunk/src/
======================
File [modified]: GRHydro_Boundaries.F90
Delta lines: +33 -30
===================================================================
--- trunk/src/GRHydro_Boundaries.F90 2012-06-05 20:25:44 UTC (rev 340)
+++ trunk/src/GRHydro_Boundaries.F90 2012-06-05 20:28:24 UTC (rev 341)
@@ -378,15 +378,18 @@
!!$Flat boundaries if required
+! The commented out code are those primitives which do not require and explicit
+! initial guess for Con2Prim. The guesses are either not needed or computed from other quantities!
+
if (CCTK_EQUALS(bound,"flat")) then
+ !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")
+ !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, &
@@ -405,34 +408,34 @@
endif
endif
- if(evolve_tracer.ne.0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "GRHydro::GRHydro_tracers", "Flat")
- endif
+ !if(evolve_tracer.ne.0) then
+ ! ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ ! "GRHydro::GRHydro_tracers", "Flat")
+ !endif
- if(evolve_y_e.ne.0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::Y_e", "Flat")
- endif
+ !if(evolve_y_e.ne.0) then
+ ! ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ ! "HydroBase::Y_e", "Flat")
+ !endif
if(evolve_temper.ne.0) then
ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
"HydroBase::temperature", "Flat")
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::entropy", "Flat")
+ !ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ ! "HydroBase::entropy", "Flat")
endif
endif
if (CCTK_EQUALS(bound,"none")) then
+ !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")
+ !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, &
@@ -453,23 +456,23 @@
endif
- if(evolve_tracer.ne.0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "GRHydro::GRHydro_tracers", "None")
-
- endif
+ !if(evolve_tracer.ne.0) then
+ ! ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ ! "GRHydro::GRHydro_tracers", "None")
+ !
+ !endif
- if(evolve_y_e.ne.0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::Y_e", "None")
-
- endif
+ !if(evolve_y_e.ne.0) then
+ ! ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ ! "HydroBase::Y_e", "None")
+ !
+ !endif
if(evolve_temper.ne.0) then
ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
"HydroBase::temperature", "None")
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::entropy", "None")
+ !ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ ! "HydroBase::entropy", "None")
endif
end if
More information about the Commits
mailing list