* add a new routine that computes dpde|rho and dpdrho|e

  at once.



git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEOS/EOS_Omni/EOS_Omni@23 8e189c6b-2ab8-4400-aa02-70a9cfce18b9

1 files changed, 4 insertions, 2 deletions

diff --git a/src/nuc_eos/findtemp.F90 b/src/nuc_eos/findtemp.F90
index d208a5f..646994d 100644
--- a/src/nuc_eos/findtemp.F90
+++ b/src/nuc_eos/findtemp.F90
@@ -16,10 +16,10 @@ subroutine findtemp(lr,lt0,y,epsin,keyerrt,rfeps)
   integer :: rl = 0
   integer itmax,i,keyerrt
   integer ii,jj,kk
-  
+
   keyerrt=0
 
-  tol=rfeps ! need to find energy to less than 1 in 10^-10
+  tol=rfeps ! need to find energy to less than 1 in 10^10
   itmax=20 ! use at most 20 iterations, then bomb
 
   lt=lt0
@@ -95,12 +95,14 @@ subroutine findtemp(lr,lt0,y,epsin,keyerrt,rfeps)
           keyerrt=0
           goto 12
        else
+#if 0
           ! total failure
           write(*,*) "EOS: Did not converge in findtemp!"
           write(*,*) "rl,logrho,logtemp0,ye,lt,eps,eps0,abs(eps-eps0)/eps0"
           write(*,"(i4,i4,1P10E19.10)") i,rl,lr,lt0,y,lt,eps,eps0,abs(eps-eps0)/eps0
           write(*,*) "Tried calling bisection... didn't help... :-/"
           write(*,*) "Bisection error: ",keyerrt
+#endif
        endif
     endif