[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 272)

[roland.haas at physics.gatech.edu]

User: rhaas
Date: 2011/09/15 11:47 AM

Modified:
 /trunk/src/
  GRHydro_Con2Prim.F90, GRHydro_Eigenproblem_Marquina.F90

Log:
 clean up variable definitions in Eigenproblem_Marquina

File Changes:

Directory: /trunk/src/
======================

File [modified]: GRHydro_Con2Prim.F90
Delta lines: +1 -1
===================================================================
--- trunk/src/GRHydro_Con2Prim.F90	2011-09-15 16:47:00 UTC (rev 271)
+++ trunk/src/GRHydro_Con2Prim.F90	2011-09-15 16:47:44 UTC (rev 272)
@@ -148,7 +148,7 @@
                    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)

File [modified]: GRHydro_Eigenproblem_Marquina.F90
Delta lines: +20 -16
===================================================================
--- trunk/src/GRHydro_Eigenproblem_Marquina.F90	2011-09-15 16:47:00 UTC (rev 271)
+++ trunk/src/GRHydro_Eigenproblem_Marquina.F90	2011-09-15 16:47:44 UTC (rev 272)
@@ -91,15 +91,17 @@
   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 @@
 
 
 ! 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