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

6 Jul 2013

User: rhaas
Date: 2013/07/06 01:10 PM
Modified:
 /trunk/
  schedule.ccl
 /trunk/src/
  GRHydro_CalcBcom.F90
Log:
 GRHydro: Make "CalcBcom" compatible with multipatch by enforcing scheduling after local->global trafo of primitives.
From: Christian Reisswig reisswig@tapir.caltech.edu
File Changes:
Directory: /trunk/src/
======================
File [modified]: GRHydro_CalcBcom.F90
Delta lines: +0 -4
===================================================================

--- trunk/src/GRHydro_CalcBcom.F90	2013-07-06 18:10:28 UTC (rev 545)
+++ trunk/src/GRHydro_CalcBcom.F90	2013-07-06 18:10:32 UTC (rev 546)
@@ -40,10 +40,6 @@
      call CCTK_WARN(0,"No storage for GRHydro::bcom")
    end if
-   if (GRHydro_UseGeneralCoordinates(cctkGH).ne.0) then
-     call CCTK_WARN(0,"MP not yet supported in GRHydro_CalcBcom")
-   end if
-
    !$OMP PARALLEL DO PRIVATE(i,j,k,velxlow, velylow, velzlow,&
    !$OMP Bvecxlow,Bvecylow,Bveczlow, bdotv,dum1,dum2,b2,bxlow,bylow,bzlow)
Directory: /trunk/
==================
File [modified]: schedule.ccl
Delta lines: +31 -1
===================================================================
--- trunk/schedule.ccl	2013-07-06 18:10:28 UTC (rev 545)
+++ trunk/schedule.ccl	2013-07-06 18:10:32 UTC (rev 546)
@@ -1363,6 +1363,36 @@
}  # sync_conserved_only
+
+#SCHEDULE GRHydro_TestSync in HydroBase_PostStep AFTER HydroBase_Boundaries BEFORE HydroBase_Con2Prim
+#{
+#  LANG: C
+#  OPTIONS: LEVEL
+#} "Sync test variables"
+
+#SCHEDULE GROUP Bla in HydroBase_PostStep AFTER HydroBase_Boundaries AFTER HydroBase_Con2Prim
+#{
+#} "bla"
+
+#SCHEDULE GRHydro_TestSync in Bla
+#{
+#  LANG: C
+#  OPTIONS: LEVEL
+#} "Sync test variables"
+
+#SCHEDULE GRHydro_TestSync in CCTK_POSTRESTRICT AFTER Refluxing_CorrectState BEFORE MoL_PostStep
+#{
+#  LANG: C
+#  OPTIONS: LEVEL
+#} "Sync test variables"
+
+#SCHEDULE GRHydro_TestSync in CCTK_POSTREGRID AFTER Refluxing_Reset
+#{
+#  LANG: C
+#  OPTIONS: LEVEL
+#} "Sync test variables"
+
+
 ############################################################
 ### Compute first differences of rho for mesh refinement ###
 ############################################################
@@ -1509,7 +1539,7 @@
 # When Tmunu is not calculated, then the following 
 # routine is used to evaluate bcom, if requested.
   if (calculate_bcom){
-    schedule GRHydro_CalcBcom IN HydroBase_PostStep
+    schedule GRHydro_CalcBcom IN HydroBase_PostStep AFTER GRHydroTransformPrimToGlobalBasis
     {
       LANG: Fortran
     } "Compute comoving magnetic field, pressure, etc..."

    