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@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" }