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

[rhaas at tapir.caltech.edu]

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 at 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