[ci optional] avoid hard floors on HELM, fall back to ideal instead

Author: evbauer

eos/private/helm.f90

@@ -257,6 +257,13 @@ subroutine helmeos2aux( &
          return
        end if
 
+       ! error out for very low pgas,egas,sgas (previously just set hard floor at this value)
+       if(pgas < 1d-20 .or. egas < 1d-20 .or. sgas < 1d-20) then
+          ierr = 1
+          if (dbg) write(*,*) 'failed in helm, sums too small'
+          return
+       end if
+
 !..compute the derivative quantities (cv, gamma1 ...etc)
        include 'helm_gammas.dek'
        if (ierr /= 0) then

eos/private/helm_sum_totals.dek

@@ -1,4 +1,4 @@
-       pgas    = max(pion + pele + ppos + pcoul, 1d-20)
+       pgas    = pion + pele + ppos + pcoul
 
        dpgasdd  = dpiondd  + dpepdd  + dpcouldd 
        dpgasdt  = dpiondt  + dpepdt  + dpcouldt
@@ -15,7 +15,7 @@
        dpgasdaz = dpiondaz + dpepdaz + dpcouldaz 
        dpgasdzz = dpiondzz + dpepdzz + dpcouldzz 
 
-       egas    = max(eion + eele + epos + ecoul, 1d-20)
+       egas    = eion + eele + epos + ecoul
 
        if (.false.) then
        write(*,'(a6,e15.5)') 'egas', egas
@@ -41,7 +41,7 @@
        degasdaz = deiondaz + deepdaz + decouldaz 
        degasdzz = deiondzz + deepdzz + decouldzz 
 
-       sgas    = max(sion + sele + spos + scoul, 1d-20)
+       sgas    = sion + sele + spos + scoul
 
        dsgasdd  = dsiondd  + dsepdd  + dscouldd
        dsgasdt  = dsiondt  + dsepdt  + dscouldt