Skip to content

fix low density/Pressure SCVH EOS tables #997

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
Debraheem wants to merge 4 commits into
base: main
Choose a base branch
from
Open

fix low density/Pressure SCVH EOS tables #997

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
b282ad2
fix low density/Pressure SCVH EOS tables with smooth extrapolation be…
Debraheem May 14, 2026
aba2d31
update test output
Debraheem May 14, 2026
51b48f0
update test additional test output
Debraheem May 14, 2026
11c6859
blend to ideal, update eos boundaries, fix small eos bugs
Debraheem May 28, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
8 changes: 4 additions & 4 deletions atm/test/test_output
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -85,17 +85,17 @@

test_irradiated

test_grey_irradiated: kap 4.3649748642304882D-03
test_grey_irradiated: kap 4.3649748858080083D-03
M/M_jupiter 1.4999999999999998D+00
R/R_jupiter 1.5743735095633800D+00
R/Rsun 1.6178684913857289D-01
kap_v 4.0000000000000001D-03
P 1.0000000000000000D+06

tau 2.9101126947637423D+00
tau 2.9101127094055270D+00
Teff 5.7759999999999891D+03
T_eq 1.0000000000000000D+03
T_int 9.0000000000000000D+02
T 1.2929349604416070D+03
logT 3.1115766787614829D+00
T 1.2929349615774468D+03
logT 3.1115766791430097D+00

20 changes: 20 additions & 0 deletions docs/source/changelog.rst
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -39,6 +39,26 @@ The parameter ``report_max_infall_inside_fe_core`` was ignored in versions r25.1

``fe_core_infall_limit`` now obeys ``when_to_stop_rtol`` and ``when_to_stop_atol`` again (broken since r11532).

Fixed a bug in the FreeEOS table builder where MESA fallback rows could write
``lnPgas``, ``lnE``, and ``lnS`` values into table columns that are read as
base-10 logarithms. This could create discontinuities and noisy EOS partial
derivatives.

At the low-density OPAL/SCVH edge, MESA now blends the OPAL/SCVH component to
the ideal EOS over :math:`-10 \leq \log_{10}\rho \leq -8` in the low-T
OPAL/SCVH region where MESA was previously falling back to HELM.

Fixed the OPAL/SCVH blend between ``Z = 0.035`` and ``Z = 0.04``. This blend
was not working properly, so the EOS was not actually blending in Z over that
interval. For more details, see
:ref:`the known bugs entry <freeeos_fallback_log_units_bug>` and
`gh-997 <https://github.com/MESAHub/mesa/pull/997>`_.

.. figure:: changelog_plots/eos_regions_x070_z002_compare.png
:alt: EOS regions in mesa-r26.4.1 and this release for X = 0.7, Z = 0.02

EOS regions for a solar-like composition in mesa-r26.4.1 and this release.


.. note:: Before releasing a new version of MESA, move `Changes in main` to a new section below with the version number as the title, and add a new `Changes in main` section at the top of the file (see ``changelog_template.rst``).

Expand Down
Binary file added docs/source/changelog_plots/eos_regions_x070_z002_compare.png
View file
Open in desktop
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified docs/source/eos/eos_regions_x000_z100.png
View file
Open in desktop
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified docs/source/eos/eos_regions_x070_z002.png
View file
Open in desktop
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
27 changes: 27 additions & 0 deletions docs/source/known_bugs.rst
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -13,6 +13,33 @@ issue, but it may not be complete.
r26.4.1
=======

.. _freeeos_fallback_log_units_bug:

EOS: FreeEOS fallback rows and OPAL/SCVH blends
-----------------------------------------------

Fixed an issue in the FreeEOS table builder where fallback rows marked by
``MESA = 1`` could write the MESA EOS values ``lnPgas``, ``lnE``, and ``lnS``
into table columns that are read as base-10 logarithms. This could create
discontinuities in the FreeEOS support table and noisy EOS partial derivatives,
including near the low-density ``logT = 3.0 - 3.1`` FreeEOS boundary.

This also exposed a low-density OPAL/SCVH edge where MESA was falling back to
HELM. The fix updates this region so OPAL/SCVH blends to the ideal EOS at low
density instead.

Another OPAL/SCVH blend issue affected the high-Z transition between
``Z = 0.035`` and ``Z = 0.04``. In practice these regions were not being
blended in Z: OPAL/SCVH remained fully active below ``Z = 0.04`` wherever it
could be used, then switched discontinuously to the HELM/ideal layer at
``Z = 0.04``. This happened because the Z ramp multiplied the existing fallback
fraction, which was zero in pure OPAL/SCVH regions. The fix fades OPAL/SCVH
into the HELM/ideal layer across this Z interval.

This is fixed after ``r26.4.1``. See
`gh-921 <https://github.com/MESAHub/mesa/issues/921>`_. Users who need the
patch directly can refer to `gh-997 <https://github.com/MESAHub/mesa/pull/997>`_.

.. _report_max_infall_inside_fe_core_bug:

Controls: ``report_max_infall_inside_fe_core`` is ignored
Expand Down
4 changes: 2 additions & 2 deletions eos/defaults/eos.defaults
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -101,8 +101,8 @@
logRho_max_FreeEOS_hi = 10.1d0
logRho_max_FreeEOS_lo = 10.0d0

logT_min_FreeEOS_hi = 3.1d0
logT_min_FreeEOS_lo = 3.0d0
logT_min_FreeEOS_hi = 3.12d0
logT_min_FreeEOS_lo = 3.02d0
logT_max_FreeEOS_hi = 8.2d0 ! FreeEOS tabulations go up to here
logT_max_FreeEOS_lo = 8.1d0

Expand Down
18 changes: 9 additions & 9 deletions eos/eosDT_builder/Makefile
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3,13 +3,13 @@ include ../../make/defaults-module.mk
# Build

MODULE_NAME := eos-dt-builder
SRCS := src/azbar.f \
src/create_eos_files.f \
SRCS := src/azbar.f90 \
src/create_eos_files.f90 \
src/eos5_xtrin_h-he.f \
src/helm_opal_scvh_driver.f \
src/helm_opal_scvh_driver.f90 \
src/opal_core.f \
src/opal_scvh_driver.f \
src/scvh_core.f
src/opal_scvh_driver.f90 \
src/scvh_core.f90
SRCS_GENERATED := \
private/helm_alloc.f90 \
private/helm_polynomials.f90 \
Expand All @@ -32,14 +32,14 @@ include ../../make/Makefile
# Generated sources
# This just copies over private implementations from eos

$(BUILD_DIR_MODULE)/private/helm_alloc.f90: ../private/helm_alloc.f90 | $(BUILD_DIR_MODULE)/private/
$(BUILD_DIR_MODULE)/private/helm_alloc.f90: ../private/helm_alloc.f90 | $(BUILD_DIR_MODULE)/private/.
@cp $^ $@

$(BUILD_DIR_MODULE)/private/helm_polynomials.f90: ../private/helm_polynomials.f90 | $(BUILD_DIR_MODULE)/private/
$(BUILD_DIR_MODULE)/private/helm_polynomials.f90: ../private/helm_polynomials.f90 | $(BUILD_DIR_MODULE)/private/.
@cp $^ $@

$(BUILD_DIR_MODULE)/private/gauss_fermi.f90: ../private/gauss_fermi.f90 | $(BUILD_DIR_MODULE)/private/
$(BUILD_DIR_MODULE)/private/gauss_fermi.f90: ../private/gauss_fermi.f90 | $(BUILD_DIR_MODULE)/private/.
@cp $^ $@

$(BUILD_DIR_MODULE)/private/helm.f90: ../private/helm.f90 | $(BUILD_DIR_MODULE)/private/
$(BUILD_DIR_MODULE)/private/helm.f90: ../private/helm.f90 | $(BUILD_DIR_MODULE)/private/.
@cp $^ $@
7 changes: 4 additions & 3 deletions eos/eosDT_builder/src/create_eos_files.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -32,7 +32,7 @@ program create_eosDT_files
real(dp) :: logQ_min = -10.09d0, logQ_max = 5.69d0, logRho_min = -14d0
real(dp) :: logT_max

integer, parameter :: version_number = 51 ! update this to force rebuilding of caches
integer, parameter :: version_number = 52 ! update this to force rebuilding of caches
! update min_version in eosDT_load_tables to force rebuild of data files

integer :: ix, io_unit, ios, info, irad
Expand Down Expand Up @@ -97,7 +97,8 @@ subroutine Make_EoS_Files(Z_in, X_in, include_radiation)
opal_only = .false., scvh_only = .false., search_for_SCVH = .true.

!opal_only = .true.
scvh_only = .true.
opal_scvh_only = .true.
scvh_only = .false.

dir = 'data' ! where to put the new data files
io_unit = 40
Expand Down Expand Up @@ -178,7 +179,7 @@ subroutine Make_EoS_Files(Z_in, X_in, include_radiation)
call do_stop('failed in helm_opal_scvh')
end if

write (io_unit, '(f4.2,3(f10.5),7(1pe13.5),1(0pf9.5),4(0pf10.5),2(0pf11.5))') &
write (io_unit, '(f4.2,3(1x,f10.5),7(1x,1pe13.5),1x,0pf9.5,1x,f12.5,3(1x,1pe13.5),2(1x,0pf12.5))') &
Copy link
Copy Markdown
Member

@evbauer evbauer May 18, 2026

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Why did the write format change?

Copy link
Copy Markdown
Member Author

@Debraheem Debraheem May 27, 2026

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Some values overflowed the old fixed-width f10.5 columns when I regenerated the eosDT tables as part of this PR. I still need to check why this overflow happened. It might be symptoms of another bug i might have introduced in this pr, so i'm looking closer.

logT, logPgas, logE, logS, chiRho, chiT, Cp, Cv, dE_dRho, dS_dT, dS_dRho, mu, &
log10(max(1d-99, free_e)), gamma1, gamma3, grad_ad, eta, HELM_fraction
end do
Expand Down
2 changes: 1 addition & 1 deletion eos/eosDT_builder/src/eos_regions_code.dek
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
if (logT >= logT1 .or. logT <= logT7 .or. logRho >= logRho1) then
if (logT >= logT1 .or. logT < logT7 .or. logRho >= logRho1) then
iregion = pure_helm
else if (logRho - 2*logT + 12 < logQmin) then
iregion = pure_helm
Expand Down
13 changes: 7 additions & 6 deletions eos/eosDT_builder/src/eos_regions_defs.dek
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -15,10 +15,11 @@
double precision, parameter :: logT2_no_rad_default = 7.7d0
! problems with HELM blend if logT2_no_rad_default is < logT1_no_rad_default

double precision, parameter :: logRho5 = -9.0d0
double precision, parameter :: logRho6 = -9.99d0
double precision, parameter :: logRho7 = -13.5d0
double precision, parameter :: logRho8 = -13.9d0
! keep the generated table out of the warm low-density FreeEOS side
double precision, parameter :: logRho5 = -14.299d0
double precision, parameter :: logRho6 = -14.3d0
double precision, parameter :: logRho7 = -14.90d0
double precision, parameter :: logRho8 = -14.99d0
double precision, parameter :: logT4 = 3.6d0
double precision, parameter :: logT5 = 3.5d0
! 2.10 is lowest in SCVH; make mesa table use all of SCVH range at lowT end
Expand All @@ -33,7 +34,7 @@
double precision, parameter :: logRho4 = logQ2 + 2*logT6 - 12

double precision, parameter :: logT_HELM = 3d0 ! smallest logT for HELM
double precision, parameter :: logQmin = -8.2d0
double precision, parameter :: logQmin = -8.05d0
double precision, parameter :: logQmax = logQ1

double precision :: alfa, beta, c_dx, c_dy, logRho_lo, logRho_hi, logT1, logT2
double precision :: alfa, beta, c_dx, c_dy, logRho_lo, logRho_hi, logT1, logT2
86 changes: 71 additions & 15 deletions eos/eosDT_builder/src/helm_opal_scvh_driver.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -151,7 +151,7 @@ subroutine helm_opal_scvh(helm_only, opal_scvh_only, opal_only, scvh_only, &
use helm
use opal_scvh_driver
use utils_lib, only: is_bad_num
use const_def, only: crad, ln10
use const_def, only: crad, kerg, amu

!..logT and logRho are log10's of temp and den
implicit none
Expand Down Expand Up @@ -371,14 +371,18 @@ subroutine helm_opal_scvh(helm_only, opal_scvh_only, opal_only, scvh_only, &

subroutine check_results
include 'formats'

! Only reject non-finite or unphysical OPAL/SCVH support points.
! Falling back to HELM for finite pressure-ionization structure
! causes noisy mu features in regenerated tables.
if (is_bad_num(logPgas_opalscvh) .or. is_bad_num(logE_opalscvh) .or. is_bad_num(logS_opalscvh) .or. &
is_bad_num(chiRho_opalscvh) .or. is_bad_num(chiT_opalscvh) .or. is_bad_num(Cp_opalscvh) .or. &
is_bad_num(Cv_opalscvh) .or. is_bad_num(dE_dRho_opalscvh) .or. is_bad_num(dS_dT_opalscvh) .or. &
is_bad_num(dS_dRho_opalscvh) .or. is_bad_num(mu_opalscvh) .or. is_bad_num(logNe_opalscvh) .or. &
is_bad_num(gamma1_opalscvh) .or. is_bad_num(gamma3_opalscvh) .or. is_bad_num(grad_ad_opalscvh) .or. &
is_bad_num(eta_opalscvh) .or. Cv_opalscvh <= 0d0 .or. logPgas_opalscvh <= -50d0 .or. &
logE_opalscvh <= -50d0 .or. logS_opalscvh <= -50d0 .or. logPgas_opalscvh <= -50d0 .or. &
gamma1_opalscvh <= 1d0 .or. gamma1_opalscvh >= 2d0 .or. grad_ad_opalscvh >= 50d0) then
gamma1_opalscvh <= 0d0) then
ierr = -1
!write(*,*) 'bail to HELM for logT logRho logQ', logT, logRho, logRho - 2*logT + 12
if (.false.) then
Expand Down Expand Up @@ -414,19 +418,14 @@ subroutine get_helmeos(include_radiation)
call helmeos2(temp, logT, den, logRho, abar, zbar, 1d6, 1d3, &
helm_res, clip_to_table_boundaries, .true., include_elec_pos, off_table, ierr)
if (ierr /= 0) then
write (*, *) 'failed in helmeos2'
write (*, 1) 'temp', temp
write (*, 1) 'logT', logT
write (*, 1) 'den', den
write (*, 1) 'logRho', logRho
write (*, 1) 'Z', Z
write (*, 1) 'X', X
write (*, 1) 'abar', abar
write (*, 1) 'zbar', zbar
write (*, *) 'clip_to_table_boundaries', clip_to_table_boundaries
write (*, *) 'include_radiation', include_radiation
write (*, *) 'stop in get_helmeos'
stop 1
call helmeos2(temp, logT, den, logRho, abar, zbar, 1d6, 1d3, &
helm_res, clip_to_table_boundaries, .true., .false., off_table, ierr)
end if
if (ierr /= 0) then
! Last resort for rectangular eosDT support points outside HELM coverage.
call get_ideal_helmeos(include_radiation)
ierr = 0
return
end if

logPgas_helm = log10(helm_res(h_pgas))
Expand Down Expand Up @@ -460,6 +459,63 @@ subroutine get_helmeos(include_radiation)

end subroutine get_helmeos

! Used only when HELM cannot fill an eosDT support point.
subroutine get_ideal_helmeos(include_radiation)
logical, intent(in) :: include_radiation
double precision :: mu_ideal, Rgas, Pgas, Prad, P, egas, erad, ener, &
sgas, srad, entr, dpressdd, dpressdt, dedt, dedd, dsdt, dsdd, &
gamma3_minus1

mu_ideal = abar/(1d0 + zbar)
Rgas = kerg/(mu_ideal*amu)

Pgas = den*Rgas*temp
egas = 1.5d0*Rgas*temp
sgas = Rgas*(1.5d0*log(temp) - log(den) + 50d0)

Prad = 0d0
erad = 0d0
srad = 0d0
if (include_radiation) then
Prad = crad*temp**4/3d0
erad = 3d0*Prad/den
srad = 4d0*Prad/(den*temp)
end if

P = Pgas + Prad
ener = egas + erad
entr = sgas + srad

dpressdd = Rgas*temp
dpressdt = den*Rgas + 4d0*Prad/temp

dedt = 1.5d0*Rgas + 4d0*erad/temp
dedd = -erad/den

dsdt = 1.5d0*Rgas/temp + 3d0*srad/temp
dsdd = -(Rgas + srad)/den

logPgas_helm = log10(Pgas)
logE_helm = log10(ener)
logS_helm = log10(entr)

chiRho_helm = dpressdd*den/P
chiT_helm = dpressdt*temp/P
Cv_helm = dedt
Cp_helm = Cv_helm + P*chiT_helm**2/(den*temp*chiRho_helm)
gamma3_minus1 = dpressdt/(den*dedt)
gamma3_helm = 1d0 + gamma3_minus1
gamma1_helm = chiRho_helm + chiT_helm*gamma3_minus1
grad_ad_helm = gamma3_minus1/gamma1_helm
dE_dRho_helm = dedd
dS_dRho_helm = dsdd
dS_dT_helm = dsdt
logNe_helm = log10(max(1d-99, den*avo*zbar/abar))
mu_helm = mu_ideal
eta_helm = -20d0

end subroutine get_ideal_helmeos

end subroutine helm_opal_scvh

end module helm_opal_scvh_driver
12 changes: 8 additions & 4 deletions eos/eosDT_builder/src/scvh_core.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -402,6 +402,11 @@ subroutine interp_vals_bicub(only_densities, search_for_SCVH, den_h, ener_h, ent
max_logP = logPs(kP)
end if

if (logP < logP_min) then
info = -1
return
end if

den_h = 10.0**interp_value(num_h_pts, denlog_h_si); if (info /= 0) return
dddpress_ct_h = d_dlogP
dddt_cp_h = d_dlogT
Expand Down Expand Up @@ -654,7 +659,6 @@ subroutine interpolate_scvh(include_radiation, search_for_SCVH, logT, logRho, T,
si0t, si1t, si2t, si0mt, si1mt, si2mt, si0p, si1p, si2p, si0mp, si1mp, si2mp, logRho_test, &
logP_guess, epslogP, epslogRho, x1, x3, y1, y3, dfdlogP, &
z, fi_h(36), fi_he(36), w0t, w1t, w2t, w0mt, w1mt, w2mt, w0p, w1p, w2p, w0mp, w1mp, w2mp

include 'formats'
info = 0
ipar => ipar_array
Expand Down Expand Up @@ -936,10 +940,10 @@ subroutine interpolate_scvh(include_radiation, search_for_SCVH, logT, logRho, T,
chiT = dpressdt*T/P

gamma3 = 1 + dpressdt/(Rho*dedt)
grad_ad = dtdpress_cs_hhe/(T*inv_P)
gamma1 = (gamma3 - 1)/grad_ad ! C&G 9.88 & 9.89

Cp = Cv + P*chiT**2/(Rho*T*chiRho) ! C&G 9.86
! Use total derivatives here since radiation may have been added above.
gamma1 = chiRho + chiT*(gamma3 - 1) ! C&G 9.88
grad_ad = (gamma3 - 1)/gamma1 ! C&G 9.89

xnhp = max(0d0, min(1d0, 1 - (xnh2 + xnh)))
xnhepp = max(0d0, min(1d0, 1 - (xnhep + xnhe)))
Expand Down
4 changes: 2 additions & 2 deletions eos/eosDT_data.tar.xz
View file
Open in desktop
Git LFS file not shown
2 changes: 1 addition & 1 deletion eos/eosFreeEOS_builder/inlist_co
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -4,7 +4,7 @@
mass_list = 'mass_frac_CO.txt'

!MESA/eos table version number
table_version = 51
table_version = 52

!FreeEOS option suite, see FreeEOS README for details
eos_version = 1
Expand Down
2 changes: 1 addition & 1 deletion eos/eosFreeEOS_builder/inlist_solar
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -4,7 +4,7 @@
mass_list = 'mass_frac_solar.txt'

!MESA/eos table version number
table_version = 51
table_version = 52

!FreeEOS option suite, see FreeEOS README for details
eos_version = 1
Expand Down
9 changes: 7 additions & 2 deletions eos/eosFreeEOS_builder/src/free_eos_table.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -248,7 +248,7 @@ subroutine read_namelist

! set defaults
mass_list = 'mass_frac.txt'
table_version = 51
table_version = 52
eos_version = 1

!set T range
Expand Down Expand Up @@ -564,7 +564,12 @@ subroutine mesa_eos_eval( logRho0, logT, mass_Frac, eos_result)
end if

eos_result(1:num_eos_basic_results) = res
eos_result(i_lnRho) = logRho

! FreeEOS data files store these logarithms in base 10.
eos_result(i_lnPgas) = res(i_lnPgas)/ln10
eos_result(i_lnE) = res(i_lnE)/ln10
eos_result(i_lnS) = res(i_lnS)/ln10
eos_result(i_lnRho) = log10Rho
Copy link
Copy Markdown
Member

@evbauer evbauer May 18, 2026

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Great catch! I think you're right about this. Do you know if it has an impact on anything?

In practice, I think we cut out most of the region where FreeEOS was failing and falling back to MESA with the "logQ_cut"-related boundary controls. Anyway, I've always been a bit uncomfortable with the idea that our "FreeEOS" tables that were meant to replace the older EOS might end up falling back to the older EOS without being very clear about that. So I like the fix, but I don't really like the existence of the thing that we're fixing here...

Copy link
Copy Markdown
Member Author

@Debraheem Debraheem May 27, 2026

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I've definitely seen models crash from this, in fact this issue from jaime comes to mind : #822

Copy link
Copy Markdown
Member Author

@Debraheem Debraheem May 27, 2026

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

after asking an ai to check this further:

It affects every row marked MESA = 1 in the FreeEOS tables, since those rows were written through the mesa_eos_eval fallback path. It only showed up at the low-T boundary because that was the place where runtime interpolation of the FreeEOS data was actually reaching those MESA = 1 rows.

total numeric rows: 14,399,221
MESA fallback rows: 2,618,641

examples pulled from the tables comparing the old and new free eos tables:

Here are a few matched v51/v52 fallback rows from mesa-FreeEOS_02z70x.data. In all of these rows, the v51 values are larger than the v52 values by ln(10), which is the sign that v51 wrote natural-log values into fields that are read as base-10 logs.

logT T [K] logRho rho [g/cm^3] v51 logP/logE/logS v52 logP/logE/logS
2.10 1.26e2 -17.89 1.29e-18 -10.0294 / 22.6406 / 20.9587 -4.3557 / 9.8327 / 9.1022
2.50 3.16e2 -9.20 6.31e-10 1.9463 / 23.8743 / 20.7341 0.8453 / 10.3685 / 9.0047
3.00 1.00e3 -0.37 4.27e-1 23.4272 / 24.6846 / 19.4259 10.1743 / 10.7204 / 8.4366
3.50 3.16e3 0.18 1.51e0 29.2003 / 28.4440 / 19.7273 12.6815 / 12.3531 / 8.5675
4.00 1.00e4 1.09 1.23e1 33.3487 / 30.8600 / 19.6948 14.4831 / 13.4023 / 8.5533
5.00 1.00e5 1.65 4.47e1 35.7939 / 32.1696 / 20.1509 15.5451 / 13.9711 / 8.7514

|

eos_result(i_dpe) = 0._dp
eos_result(i_dsp) = 0._dp
eos_result(i_dse) = 0._dp
Expand Down
Loading
Loading