diff options
-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 |