diff options
| author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-01-14 14:23:37 +0000 |
|---|---|---|
| committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-01-14 14:23:37 +0000 |
| commit | 5199c61af4721444a26c4e7c1c30d3b8fdedf54c (patch) | |
| tree | 9428f573d53e23f6b2970cbe664ba96cfb2eb2ec /src/GRHydro_ParamCheck.F90 | |
| parent | 48b053d41b2aa9ce720d443d288eef535fb9651d (diff) | |
GRHydro: Add basic vector potential support
Basic cell-centered, algebraic gauge vector potential
method with place-holders for lorenz gauge.
Initial Avec constrained to poloidal at the moment.
From: Tanja Bode <tanja.bode@physics.gatech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@456 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_ParamCheck.F90')
| -rw-r--r-- | src/GRHydro_ParamCheck.F90 | 10 |
1 files changed, 9 insertions, 1 deletions
diff --git a/src/GRHydro_ParamCheck.F90 b/src/GRHydro_ParamCheck.F90 index c738f79..1ab2a40 100644 --- a/src/GRHydro_ParamCheck.F90 +++ b/src/GRHydro_ParamCheck.F90 @@ -88,12 +88,20 @@ subroutine GRHydro_ParamCheck(CCTK_ARGUMENTS) call CCTK_PARAMWARN("GRHydro::bound = 'static' is no longer supported, use 'none' instead"); end if - if (CCTK_EQUALS(Bvec_evolution_method,"GRHydro")) then + if (CCTK_EQUALS(Bvec_evolution_method,"GRHydro").or.CCTK_EQUALS(Bvec_evolution_method,"GRHydro_Avec")) then evolve_MHD = 1 else evolve_MHD = 0 endif + if ( CCTK_EQUALS(Bvec_evolution_method,"GRHydro_Avec") ) then + if ( CCTK_EQUALS(Avec_gauge,"Lorenz") ) then + evolve_Lorenz_gge = 1 + else + evolve_Lorenz_gge = 0 + endif + endif + if (CCTK_EQUALS(Y_e_evolution_method,"GRHydro")) then evolve_Y_e = 1 else |