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

14 Jan 2013

User: rhaas
Date: 2013/01/14 08:23 AM
Modified:
 /trunk/
  schedule.ccl
 /trunk/src/
  GRHydro_Boundaries.F90, GRHydro_Con2PrimAM.F90
Log:
 GRHydro: Starting to push vector potential changes
* incl boundary calls for Avec
 * disable call to BvecfromAvec from within loop, this one is a schedule
   function call now.
From: Philipp Moesta pmoesta@tapir.caltech.edu
File Changes:
Directory: /trunk/src/
======================
File [modified]: GRHydro_Boundaries.F90
Delta lines: +120 -134
===================================================================

--- trunk/src/GRHydro_Boundaries.F90	2013-01-14 14:23:38 UTC (rev 457)
+++ trunk/src/GRHydro_Boundaries.F90	2013-01-14 14:23:40 UTC (rev 458)
@@ -85,11 +85,6 @@
       call SetCartSymGN(ierr, cctkGH, sym, "HydroBase::entropy")
     !endif
   endif
-
-  if(evolve_entropy.ne.0) then
-      call SetCartSymGN(ierr, cctkGH, sym, "HydroBase::entropy")
-      call SetCartSymGN(ierr, cctkGH, sym, "GRHydro::entropycons")
-  endif
call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::atmosphere_mask_real")
@@ -99,9 +94,9 @@
!if (sync_conserved_only .eq. 0) then
       if (general_coordinates .ne. 0) then
-        call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::lvel[0]")
+	 call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::lvel[0]")
       else
-        call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::vel[0]")
+	 call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::vel[0]")
       endif
   !endif
   call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::scon[0]")
@@ -118,15 +113,15 @@
!if (sync_conserved_only .eq. 0) then
       if (general_coordinates .ne. 0) then
-         call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::lvel[1]")
+	 call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::lvel[1]")
       else
-         call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::vel[1]")
+	 call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::vel[1]")
       endif
   !endif
   call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::scon[1]")
   if(evolve_mhd.ne.0) then
      !if (sync_conserved_only .eq. 0) then
-     call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::Bvec[1]")
+	 call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::Bvec[1]")
      !endif
      call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::Bcons[1]")
   endif
@@ -137,9 +132,9 @@
!if (sync_conserved_only .eq. 0) then
       if (general_coordinates .ne. 0) then
-        call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::lvel[2]")
+	 call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::lvel[2]")
       else
-        call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::vel[2]")
+	 call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::vel[2]")
       endif
   !endif
   call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::scon[2]")
@@ -206,31 +201,31 @@
          "GRHydro::scon", "Flat")
     if (sync_conserved_only .eq. 0) then
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-           "HydroBase::w_lorentz", "Flat")
+	    "HydroBase::w_lorentz", "Flat")
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-           "HydroBase::rho", "Flat")
+	    "HydroBase::rho", "Flat")
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-           "HydroBase::press", "Flat")
+	    "HydroBase::press", "Flat")
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-           "HydroBase::eps", "Flat")
+	    "HydroBase::eps", "Flat")
       if (general_coordinates .ne. 0) then
-        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "GRHydro::lvel", "Flat")
+	 ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+	       "GRHydro::lvel", "Flat")
       else
-        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "HydroBase::vel", "Flat")
+	 ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+	       "HydroBase::vel", "Flat")
       endif
     endif
     if(evolve_mhd.ne.0) then
-      if (sync_conserved_only .eq. 0) then
-        if (general_coordinates .ne. 0) then
-          ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-               "GRHydro::lBvec", "Flat")
-        else
-          ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-               "HydroBase::Bvec", "Flat")
-        endif
-      endif
+       if (sync_conserved_only .eq. 0) then
+	 if (general_coordinates .ne. 0) then
+	    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+		  "GRHydro::lBvec", "Flat")
+	 else
+	    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+		  "HydroBase::Bvec", "Flat")
+	 endif
+       endif
        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
             "GRHydro::Bcons", "Flat")
        if(clean_divergence.ne.0) then
@@ -239,10 +234,22 @@
        endif
     endif
+    if (CCTK_Equals(Bvec_evolution_method, "GRHydro_Avec")) then 
+      ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "HydroBase::Avec", "Flat")
+      if (general_coordinates .ne. 0) then
+        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+               "GRHydro::lBvec", "Flat")
+      else
+        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+               "HydroBase::Bvec", "Flat")
+      endif
+    endif 
+
     if(evolve_tracer.ne.0) then 
       if (sync_conserved_only .eq. 0) then
-        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "GRHydro::GRHydro_tracers", "Flat")
+	 ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+	       "GRHydro::GRHydro_tracers", "Flat")
       endif
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
             "GRHydro::GRHydro_cons_tracers", "Flat")
@@ -250,8 +257,8 @@
if(evolve_y_e.ne.0) then
       if (sync_conserved_only .eq. 0) then
-        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "HydroBase::Y_e", "Flat")
+	 ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+	       "HydroBase::Y_e", "Flat")
       endif
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
             "GRHydro::Y_e_con", "Flat")
@@ -259,21 +266,13 @@
if(evolve_temper.ne.0) then
       if (sync_conserved_only .eq. 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::temperature", "Flat")
+	 ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+	       "HydroBase::entropy", "Flat")
       endif
     endif
-    if(evolve_entropy.ne.0) then
-      if (sync_conserved_only .eq. 0) then
-        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "HydroBase::entropy", "Flat")
-      endif
-      ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-           "GRHydro::entropycons", "Flat")
-    endif
endif
@@ -287,44 +286,56 @@
if (sync_conserved_only .eq. 0) then
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-           "HydroBase::w_lorentz", "None")
+	    "HydroBase::w_lorentz", "None")
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-           "HydroBase::rho", "None")
+	    "HydroBase::rho", "None")
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-           "HydroBase::press", "None")
+	    "HydroBase::press", "None")
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-           "HydroBase::eps", "None")
+	    "HydroBase::eps", "None")
       if (general_coordinates .ne. 0) then
-        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "GRHydro::lvel", "None")
+	 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")
+	 ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+	       "HydroBase::vel", "None")
       endif
     endif
if(evolve_mhd.ne.0) then
-      if (sync_conserved_only .eq. 0) then
-        if (general_coordinates .ne. 0) then
+       if (sync_conserved_only .eq. 0) then
+	 if (general_coordinates .ne. 0) then
+	    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+		  "GRHydro::lBvec", "None")
+	 else
+	    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+		  "HydroBase::Bvec", "None")
+	 endif
+       endif
+       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+            "GRHydro::Bcons", "None")
+       if(clean_divergence.ne.0) then
           ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-               "GRHydro::lBvec", "None")
-        else
-          ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-               "HydroBase::Bvec", "None")
-        endif
-      endif
+               "GRHydro::psidc", "None")
+       endif
+    endif
+
+    if (CCTK_Equals(Bvec_evolution_method, "GRHydro_Avec")) then 
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-           "GRHydro::Bcons", "None")
-      if(clean_divergence.ne.0) then
-         ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-              "GRHydro::psidc", "None")
+             "HydroBase::Avec", "None")
+      if (general_coordinates .ne. 0) then
+        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+               "GRHydro::lBvec", "Flat")
+      else
+        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+               "HydroBase::Bvec", "Flat")
       endif
-    endif
+    endif
if(evolve_tracer.ne.0) then 
       if (sync_conserved_only .eq. 0) then 
-        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "GRHydro::GRHydro_tracers", "None")
+	 ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+	       "GRHydro::GRHydro_tracers", "None")
       endif
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
             "GRHydro::GRHydro_cons_tracers", "None")
@@ -332,8 +343,8 @@
if(evolve_y_e.ne.0) then
       if (sync_conserved_only .eq. 0) then
-        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "HydroBase::Y_e", "None")
+	 ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+	       "HydroBase::Y_e", "None")
       endif
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
             "GRHydro::Y_e_con", "None")
@@ -341,22 +352,13 @@
if(evolve_temper.ne.0) then
       if (sync_conserved_only .eq. 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::temperature", "None")
+	 ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+	       "HydroBase::entropy", "None")
       endif
     endif
-    if(evolve_entropy.ne.0) then
-      if (sync_conserved_only .eq. 0) then
-        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "HydroBase::entropy", "None")
-      endif
-      ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-           "GRHydro::entropycons", "None")
-    endif
-
   endif
if (CCTK_EQUALS(bound,"scalar")) then
@@ -415,13 +417,13 @@
              "HydroBase::vel", "Flat")
     endif
     if(evolve_mhd.ne.0) then
-      if (general_coordinates .ne. 0) then
-        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "GRHydro::lBvec", "Flat")
-      else
-        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "HydroBase::Bvec", "Flat")
-      endif
+	 if (general_coordinates .ne. 0) then
+	    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+		  "GRHydro::lBvec", "Flat")
+	 else
+	    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+		  "HydroBase::Bvec", "Flat")
+	 endif
     endif
!if(evolve_tracer.ne.0) then 
@@ -441,10 +443,6 @@
        !      "HydroBase::entropy", "Flat")
     endif
-    if(evolve_entropy.ne.0) then
-       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "HydroBase::entropy", "Flat")
-    endif
endif
@@ -466,13 +464,14 @@
     endif
if(evolve_mhd.ne.0) then
-      if (general_coordinates .ne. 0) then
-        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "GRHydro::lBvec", "None")
-      else
-        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "HydroBase::Bvec", "None")
-      endif
+	 if (general_coordinates .ne. 0) then
+	    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+		  "GRHydro::lBvec", "None")
+	 else
+	    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+		  "HydroBase::Bvec", "None")
+	 endif
+       
     endif
!if(evolve_tracer.ne.0) then 
@@ -494,11 +493,6 @@
        !      "HydroBase::entropy", "None")
     endif
-    if(evolve_entropy.ne.0) then
-       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "HydroBase::entropy", "None")
-    endif
-
   end if
if (CCTK_EQUALS(bound,"scalar")) then
@@ -558,25 +552,24 @@
        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
              "HydroBase::vel", "Flat")
     endif
-
     if(evolve_mhd.ne.0) then
-      if (general_coordinates .ne. 0) then
-        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "GRHydro::lBvec", "Flat")
-      else
-        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "HydroBase::Bvec", "Flat")
-      endif
+	 if (general_coordinates .ne. 0) then
+	    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+		  "GRHydro::lBvec", "Flat")
+	 else
+	    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+		  "HydroBase::Bvec", "Flat")
+	 endif
     endif
if(evolve_tracer.ne.0) then 
-      ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-           "GRHydro::GRHydro_tracers", "Flat")
+    	 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")
+    	 ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+    	       "HydroBase::Y_e", "Flat")
     endif
if(evolve_temper.ne.0) then
@@ -586,10 +579,6 @@
              "HydroBase::entropy", "Flat")
     endif
-    if(evolve_entropy.ne.0) then
-       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "HydroBase::entropy", "Flat")
-    endif
endif
@@ -611,18 +600,20 @@
     endif
if(evolve_mhd.ne.0) then
-      if (general_coordinates .ne. 0) then
-        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "GRHydro::lBvec", "None")
-      else
-        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "HydroBase::Bvec", "None")
-      endif
+	 if (general_coordinates .ne. 0) then
+	    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+		  "GRHydro::lBvec", "None")
+	 else
+	    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+		  "HydroBase::Bvec", "None")
+	 endif
+       
     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
@@ -638,11 +629,6 @@
              "HydroBase::entropy", "None")
     endif
-    if(evolve_entropy.ne.0) then
-       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "HydroBase::entropy", "None")
-    endif
-
   end if
if (CCTK_EQUALS(bound,"scalar")) then
File [modified]: GRHydro_Con2PrimAM.F90
Delta lines: +0 -3
===================================================================
--- trunk/src/GRHydro_Con2PrimAM.F90	2013-01-14 14:23:38 UTC (rev 457)
+++ trunk/src/GRHydro_Con2PrimAM.F90	2013-01-14 14:23:40 UTC (rev 458)
@@ -102,9 +102,6 @@
   !as conformally flat only for con2prim inversion purposes.
   posdef = .true.
-  !! Get Bvec from Avec before rest, so we can reuse pointwise MHD recovery
-  call GRHydro_BvecfromAvec(CCTK_PASS_FTOF)
-
   if (GRHydro_UseGeneralCoordinates(cctkGH).ne.0) then
     g11 => gaa
     g12 => gab
Directory: /trunk/
==================
File [modified]: schedule.ccl
Delta lines: +11 -6
===================================================================
--- trunk/schedule.ccl	2013-01-14 14:23:38 UTC (rev 457)
+++ trunk/schedule.ccl	2013-01-14 14:23:40 UTC (rev 458)
@@ -1066,21 +1066,26 @@
     {
       LANG: Fortran
     } "Convert initial data given in primive variables to conserved variables - MHD version"
-  } 
-  else if(CCTK_Equals(Bvec_evolution_method,"GRHydro_Avec")) 
+  } else if(CCTK_Equals(Bvec_evolution_method,"GRHydro_Avec")) 
   {
-    schedule Conservative2PrimitiveAM IN HydroBase_Con2Prim AS Con2Prim IF GRHydro::execute_MoL_PostStep
+    schedule GRHydro_BvecfromAvec IN HydroBase_Con2Prim
     {
       LANG: Fortran
+    } "Populate Bvec from Avec"
+#   SCHEDULE group HydroBase_Boundaries AT CCTK_POSTINITIAL AFTER GRHydro_PoloidalMagFieldM AFTER GRHydro_BvecfromAvec
+    SCHEDULE group HydroBase_Boundaries IN HydroBase_Con2Prim AFTER GRHydro_BvecfromAvec
+    {
+    } "Call boundary conditions after magnetic field initial data setup"
+    schedule Conservative2PrimitiveAM IN HydroBase_Con2Prim AS Con2Prim AFTER HydroBase_Boundaries IF GRHydro::execute_MoL_PostStep
+    {
+      LANG: Fortran
     } "Convert back to primitive variables (general) - MHD with Avec version"
schedule Primitive2ConservativeCellsAM IN HydroBase_Prim2ConInitial
     {
       LANG: Fortran
     } "Convert initial data given in primive variables to conserved variables - MHD with Avec version"
-  }
-  else
-  {
+  } else {
     schedule Conservative2Primitive IN HydroBase_Con2Prim AS Con2Prim IF GRHydro::execute_MoL_PostStep
     {
       LANG: Fortran

    