User: rhaas Date: 2013/01/11 09:04 AM

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

Log: GRHydro_InitData: Add basic vector potential support

Initial Avec constrained to poloidal at the moment.

From: Tanja Bode tanja.bode@physics.gatech.edu

File Changes:

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

File [modified]: GRHydro_PoloidalMagFieldM.F90 Delta lines: +17 -1 =================================================================== --- trunk/src/GRHydro_PoloidalMagFieldM.F90 2013-01-11 15:04:11 UTC (rev 194) +++ trunk/src/GRHydro_PoloidalMagFieldM.F90 2013-01-11 15:04:13 UTC (rev 195) @@ -29,6 +29,9 @@ #define Bconsx(i,j,k) Bcons(i,j,k,1) #define Bconsy(i,j,k) Bcons(i,j,k,2) #define Bconsz(i,j,k) Bcons(i,j,k,3) +#define Avecx(i,j,k) Avec(i,j,k,1) +#define Avecy(i,j,k) Avec(i,j,k,2) +#define Avecz(i,j,k) Avec(i,j,k,3)

/*@@ @routine GRHydro_PoloidalMagFieldM @@ -56,7 +59,7 @@ DECLARE_CCTK_PARAMETERS DECLARE_CCTK_FUNCTIONS

- CCTK_INT :: i, j, k, nx, ny, nz + CCTK_INT :: i, j, k, nx, ny, nz, set_Avec CCTK_REAL :: det CCTK_REAL :: sdet CCTK_REAL :: dx,dy,dz @@ -75,8 +78,15 @@ dy = CCTK_DELTA_SPACE(2) dz = CCTK_DELTA_SPACE(3)

+ set_Avec = 0 + if ( CCTK_EQUALS(Bvec_evolution_method,"GRHydro_Avec") ) then + set_Avec = 1 + end if + write(*,*)'GRHydro_InitData: Setting up initial poloidal magnetic field'

+ + do i=2,nx-1 do j=2,ny-1 do k=2,nz-1 @@ -123,6 +133,12 @@ Bvecy(i,j,k) = - Ax_dz/sdet Bvecz(i,j,k) = (Ax_dy-Ay_dx)/sdet

+ if ( set_Avec.gt.0 ) then + Avecx(i,j,k) = Ax + Avecy(i,j,k) = Ay + Avecz(i,j,k) = Az + end if + !Bvecx(i,j,k) = 0.0d0 !Bvecy(i,j,k) = 0.0d0 !Bvecz(i,j,k) = 0.00000001/sdet

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

File [modified]: param.ccl Delta lines: +5 -0 =================================================================== --- trunk/param.ccl 2013-01-11 15:04:11 UTC (rev 194) +++ trunk/param.ccl 2013-01-11 15:04:13 UTC (rev 195) @@ -35,6 +35,11 @@ "magnetized Bondi" :: "radial magnetic field appropriate for Bondi test" }

+EXTENDS KEYWORD initial_Avec +{ + "poloidalmagfield" :: "Poloidal Magnetic Field" +} + EXTENDS KEYWORD initial_entropy { "magnetized Bondi" :: "Initial entropy for a radial magnetic field appropriate for Bondi test"

File [modified]: schedule.ccl Delta lines: +1 -1 =================================================================== --- trunk/schedule.ccl 2013-01-11 15:04:11 UTC (rev 194) +++ trunk/schedule.ccl 2013-01-11 15:04:13 UTC (rev 195) @@ -250,7 +250,7 @@ } "setup GRHydro vars for the magnetized Bondi solution" }

-if (CCTK_EQUALS(initial_Bvec, "poloidalmagfield")) +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