From e6f4a06735c8a0d38251cd984147a06abea2df16 Mon Sep 17 00:00:00 2001 From: knarf Date: Thu, 16 Dec 2010 21:02:36 +0000 Subject: remove some usage of uninitialized memory. Tidy up code a bit git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEOS/EOS_Omni/trunk@31 8e189c6b-2ab8-4400-aa02-70a9cfce18b9 --- src/nuc_eos/nuc_eos.F90 | 28 ++++++++++++---------------- 1 file changed, 12 insertions(+), 16 deletions(-) (limited to 'src') diff --git a/src/nuc_eos/nuc_eos.F90 b/src/nuc_eos/nuc_eos.F90 index f537930..b7c9312 100644 --- a/src/nuc_eos/nuc_eos.F90 +++ b/src/nuc_eos/nuc_eos.F90 @@ -233,13 +233,13 @@ subroutine nuc_eos_short(xrho,xtemp,xye,xenr,xprs,xent,xcs2,xdedt,& lr = log10(xrho) lt = log10(xtemp) y = xye - xeps = xenr + energy_shift - leps = log10(max(xeps,1.0d0)) keyerr = 0 if(keytemp.eq.0) then !need to find temperature based on xeps + xeps = xenr + energy_shift + leps = log10(max(xeps,1.0d0)) call findtemp(lr,lt,y,leps,keyerrt,rfeps) if(keyerrt.ne.0) then keyerr = keyerrt @@ -349,34 +349,30 @@ subroutine nuc_eos_press_eps(xrho,xtemp,xye,xenr,xprs,& endif endif - lr = log10(xrho) - lt = log10(xtemp) - y = xye - xeps = xenr + energy_shift - leps = log10(max(xeps,1.0d0)) - keyerr = 0 + if(keytemp.gt.1) then + stop "eos_nuc_press does not support keytemp other than 0 and 1" + endif + + lr = log10(xrho) + lt = log10(xtemp) + y = xye if(keytemp.eq.0) then !need to find temperature based on xeps + xeps = xenr + energy_shift + leps = log10(max(xeps,1.0d0)) call findtemp(lr,lt,y,leps,keyerrt,rfeps) if(keyerrt.ne.0) then keyerr = keyerrt return endif xtemp = 10.0d0**lt - - elseif(keytemp.gt.1) then - stop "eos_nuc_press does not support keytemp other than 0 and 1" endif - ! have temperature, proceed: call findall_press_eps(lr,lt,y,ff) xprs = 10.0d0**ff(1) - - if(keytemp.eq.1) then - xenr = 10.0d0**ff(2) - energy_shift - endif + xenr = 10.0d0**ff(2) - energy_shift end subroutine nuc_eos_press_eps -- cgit v1.2.3