diff options
author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-01-14 14:23:33 +0000 |
---|---|---|
committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-01-14 14:23:33 +0000 |
commit | 48b053d41b2aa9ce720d443d288eef535fb9651d (patch) | |
tree | bb5274de4dfb94de758bafc0ee77e6b5f499fb06 /src/GRHydro_ParamCheck.F90 | |
parent | d6e45e79d92f3d2adf3eb4ebcc7641db07119f30 (diff) |
GRHydro: Implement entropy evolution (separated from temp evolution)
From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@455 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_ParamCheck.F90')
-rw-r--r-- | src/GRHydro_ParamCheck.F90 | 15 |
1 files changed, 15 insertions, 0 deletions
diff --git a/src/GRHydro_ParamCheck.F90 b/src/GRHydro_ParamCheck.F90 index 085fac1..c738f79 100644 --- a/src/GRHydro_ParamCheck.F90 +++ b/src/GRHydro_ParamCheck.F90 @@ -106,6 +106,12 @@ subroutine GRHydro_ParamCheck(CCTK_ARGUMENTS) evolve_temper = 0 endif + if (CCTK_EQUALS(entropy_evolution_method,"GRHydro")) then + evolve_entropy = 1 + else + evolve_entropy = 0 + endif + call CCTK_IsImplementationActive(coordinates_is_active, "Coordinates") ! this test is somewhat overzealous but we cannot access ! coordinates::general_coordinates yet since it is only set in BaseGrid @@ -153,5 +159,14 @@ subroutine GRHydro_ParamCheck(CCTK_ARGUMENTS) call CCTK_WARN(0, "Minimum damping radius is greater than maximum damping radius!") end if + if(evolve_entropy.ne.0) then + if(evolve_MHD.eq.0) then + call CCTK_PARAMWARN("Entropy evolution not properly implemented yet for unmagnetized fluids!") + end if + if (CCTK_EQUALS(recon_method,"ppm")) then + call CCTK_PARAMWARN("Entropy evolution not implemented yet with PPM reconstruction!") + end if + end if + end subroutine GRHydro_ParamCheck |