diff options
Diffstat (limited to 'src/nuc_eos/findtemp.F90')
-rw-r--r-- | src/nuc_eos/findtemp.F90 | 6 |
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 |