[Commits] [svn:einsteintoolkit] GRHydro_InitData/trunk/ (Rev. 215)

rhaas at tapir.caltech.edu rhaas at tapir.caltech.edu
Sat Jul 6 13:12:16 CDT 2013


User: rhaas
Date: 2013/07/06 01:12 PM

Modified:
 /trunk/
  schedule.ccl
 /trunk/src/
  GRHydro_PoloidalMagFieldM.F90

Log:
 GRHydro_InitData: Reschedule to make things work with multipatch. Spit out a compile-time warning since
 poloidal ID B-field routine is not yet adapted to MP.
 
 From: Christian Reisswig <reisswig at tapir.caltech.edu>

File Changes:

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

File [modified]: GRHydro_PoloidalMagFieldM.F90
Delta lines: +9 -1
===================================================================
--- trunk/src/GRHydro_PoloidalMagFieldM.F90	2013-05-27 22:32:37 UTC (rev 214)
+++ trunk/src/GRHydro_PoloidalMagFieldM.F90	2013-07-06 18:12:15 UTC (rev 215)
@@ -85,7 +85,8 @@
 
   write(*,*)'GRHydro_InitData: Setting up initial poloidal magnetic field'
 
-  
+! Initialize to zero
+  Bvec = 0.0d0
 
   do i=2,nx-1
    do j=2,ny-1
@@ -106,6 +107,7 @@
 !!     write(*,*)'delPcut = ', delPcut
 !!     write(*,*)'maxP_Pcut = ', maxP_Pcut 
 
+#warning "This algorithm does only work on Cartesian grids!!"
      rho_dx = 0.5d0*(rho(i+1,j,k)-rho(i-1,j,k))/dx 
      rho_dy = 0.5d0*(rho(i,j+1,k)-rho(i,j-1,k))/dy 
      rho_dz = 0.5d0*(rho(i,j,k+1)-rho(i,j,k-1))/dz 
@@ -183,6 +185,12 @@
   if (clean_divergence .ne. 0) then
     psidcrhs =0.0
   endif
+  
+  !Bvec = 0
+  !lBvec = 0
+  !Bcons = 0
+  !Avec = 0
+  
   return
   
 end subroutine GRHydro_PoloidalMagFieldM

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

File [modified]: schedule.ccl
Delta lines: +4 -2
===================================================================
--- trunk/schedule.ccl	2013-05-27 22:32:37 UTC (rev 214)
+++ trunk/schedule.ccl	2013-07-06 18:12:15 UTC (rev 215)
@@ -270,11 +270,13 @@
 if (CCTK_EQUALS(initial_Bvec, "poloidalmagfield") || CCTK_EQUALS(initial_Avec, "poloidalmagfield"))
 {
 #  SCHEDULE GRHydro_PoloidalMagFieldM AT CCTK_INITIAL AFTER IN HydroBase_Initial AFTER rnsid_init AFTER TOV_Initial_Data after CCCC_StarMapper_InitialData
-  SCHEDULE GRHydro_PoloidalMagFieldM AT CCTK_POSTINITIAL 
+#  SCHEDULE GRHydro_PoloidalMagFieldM AT CCTK_POSTINITIAL 
+  SCHEDULE GRHydro_PoloidalMagFieldM AT CCTK_INITIAL AFTER HydroBase_Initial BEFORE GRHydroTransformPrimToLocalBasis
   {
     LANGUAGE: Fortran
   } "Set up a poloidal magnetic field. It expects the other fluid variables already to be set, as for example in the TOV solution"
-  SCHEDULE group HydroBase_Boundaries IN HydroBase_Initial AFTER GRHydro_PoloidalMagFieldM 
+#  SCHEDULE group HydroBase_Boundaries IN HydroBase_Initial AFTER GRHydro_PoloidalMagFieldM 
+  SCHEDULE group HydroBase_Boundaries IN CCTK_INITIAL AFTER GRHydro_PoloidalMagFieldM BEFORE GRHydroTransformPrimToLocalBasis
   {
   } "Call boundary conditions after magnetic field initial data setup"
 }



More information about the Commits mailing list