diff options
author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-09-15 16:47:44 +0000 |
---|---|---|
committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-09-15 16:47:44 +0000 |
commit | 374f96eabf2c562985209711c05b68624cd866e1 (patch) | |
tree | 8630463354ec9497dc9f405c45ae8caa09ba7264 | |
parent | b75ed5dcd756f0243be784abfb95ceb2c1095106 (diff) |
clean up variable definitions in Eigenproblem_Marquina
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@272 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
-rw-r--r-- | src/GRHydro_Con2Prim.F90 | 2 | ||||
-rw-r--r-- | src/GRHydro_Eigenproblem_Marquina.F90 | 36 |
2 files changed, 21 insertions, 17 deletions
diff --git a/src/GRHydro_Con2Prim.F90 b/src/GRHydro_Con2Prim.F90 index 220f92d..d3c9732 100644 --- a/src/GRHydro_Con2Prim.F90 +++ b/src/GRHydro_Con2Prim.F90 @@ -148,7 +148,7 @@ subroutine Conservative2Primitive(CCTK_ARGUMENTS) rho(i,j,k),eps(i,j,k),temperature(i,j,k),y_e(i,j,k),& press(i,j,k),keyerr,anyerr) else - ! use polytropic EOS + ! w_lorentz=1, so the expression for tau reduces to: keytemp = 0 call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,GRHydro_eos_rf_prec,n,& rho(i,j,k),eps(i,j,k),xtemp,xye,press(i,j,k),keyerr,anyerr) diff --git a/src/GRHydro_Eigenproblem_Marquina.F90 b/src/GRHydro_Eigenproblem_Marquina.F90 index 102e69e..970f7ac 100644 --- a/src/GRHydro_Eigenproblem_Marquina.F90 +++ b/src/GRHydro_Eigenproblem_Marquina.F90 @@ -91,15 +91,17 @@ subroutine eigenproblem_marquina(handle,rhor,velxr,velyr,& CCTK_REAL invulpfac, invulmfac, invurpfac, invurmfac ! begin EOS Omni vars - integer :: n = 1 - integer :: keytemp = 0 - integer :: anyerr = 0 - integer :: keyerr(1) = 0 - real*8 :: xpress = 0.0d0 - real*8 :: xeps = 0.0d0 - real*8 :: xtemp = 0.0d0 - real*8 :: xye = 0.0d0 + integer :: n + integer :: keytemp + integer :: anyerr + integer :: keyerr(1) + real*8 :: xpress + real*8 :: xeps + real*8 :: xtemp + real*8 :: xye ! end EOS Omni vars + n=1;keytemp=0;anyerr=0;keyerr(1)=0 + xpress=0.0d0;xeps=0.0d0;xtemp=0.0d0;xye=0.0d0 one = 1.0d0 two = 2.0d0 @@ -816,15 +818,17 @@ subroutine eigenproblem_marquina_hot(handle,rhor,velxr,velyr,& ! begin EOS Omni vars - integer :: n = 1 - integer :: keytemp = 0 - integer :: anyerr = 0 - integer :: keyerr(1) = 0 - real*8 :: xpress = 0.0d0 - real*8 :: xeps = 0.0d0 - real*8 :: xtemp = 0.0d0 - real*8 :: xye = 0.0d0 + integer :: n + integer :: keytemp + integer :: anyerr + integer :: keyerr(1) + real*8 :: xpress + real*8 :: xeps + real*8 :: xtemp + real*8 :: xye ! end EOS Omni vars + n=1;keytemp=0;anyerr=0;keyerr(1)=0 + xpress=0.0d0;xeps=0.0d0;xtemp=0.0d0;xye=0.0d0 one = 1.0d0 two = 2.0d0 |