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

roland.haas at physics.gatech.edu roland.haas at physics.gatech.edu
Sun May 13 14:12:48 CDT 2012


User: rhaas
Date: 2012/05/13 02:12 PM

Modified:
 /trunk/src/
  GRHydro_Boundaries.F90

Log:
 GRHydro: Fixed minor bugs in registration in Boundaries
 
 Some variables were stil set to "Flat" boundaries even though "None" is
 supposed to be set.
 
 Patch by Christian Reisswig.

File Changes:

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

File [modified]: GRHydro_Boundaries.F90
Delta lines: +9 -3
===================================================================
--- trunk/src/GRHydro_Boundaries.F90	2012-05-13 19:12:00 UTC (rev 324)
+++ trunk/src/GRHydro_Boundaries.F90	2012-05-13 19:12:48 UTC (rev 325)
@@ -260,8 +260,14 @@
          "HydroBase::press", "None")
     ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
          "HydroBase::eps", "None")
-    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-         "HydroBase::vel", "None")
+    if (general_coordinates .ne. 0) then
+       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "GRHydro::lvel", "None")
+    else
+       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "HydroBase::vel", "None")
+    endif
+    
     if(evolve_mhd.ne.0) then
        if (general_coordinates .ne. 0) then
           ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
@@ -299,7 +305,7 @@
             "HydroBase::entropy", "None")
     endif
 
-  end if
+  endif
 
   if (CCTK_EQUALS(bound,"scalar")) then
     call CCTK_WARN(0, "Until somebody uses this I see no reason to support it")



More information about the Commits mailing list