eeeos06_mod.f

34 lines · 1 definitions · on GitHub (7be4133)

Python version: eeeos06_mod.py, 20 lines, on the right (switch to Python above the code).

Variables on this page (5)

epl, mx, nr, nt, xx

Purpose

EEEOS06_MOD holds workspace of the OPAL 2006 equation-of-state interpolation: the table values of one quantity being interpolated, and the $X$ values of the tables. It replaces COMMON/EEEOS06/ (MHP via Claude, 10/20/26). The table dimensions MX = 5, NT = 197, NR = 169 are PRIVATE copies, as in AEOS06_MOD.

What its variables hold, who sets them, who reads them.

EPL(MX,NT,NR): for each $X$ index, the table values of the quantity being interpolated (XZ(MIDX,IV,IT,IR) of AEOS06_MOD), copied by ESAC06 for each variable IV and read there in the interpolation.
XX(MX): the $X$ values of the five tables, copied from XA by SETEOS06; read by ESAC06 and by READCOEOS06, which forms the inverse spacings DFSX from them.
T6RINTEOS06 also uses the module but reads neither member.

Design

A blanket USE EEEOS06_MOD. XX duplicates XA.

Definitions

moduleEEEOS06_MOD

lines 22–34 · GitHub

Used by (4 units)

ESAC06, READCOEOS06, SETEOS06, T6RINTEOS06

Change log

The changes since rp_1

Changed since rp_1. New file: COMMON/EEEOS06/ as a module.

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C EEEOS06_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/20/26 MHP via Claude 7C code changed to convert COMMON/EEEOS06/ into this module, completing 8C stages 3 and 4 (KNOWN_ISSUES.md issue 37). 9C 10C Imported with a blanket USE, the convention for the conversion. 11C All 4 declarations of this block were identical, so the USE 12C imports exactly the names the declaration already made global in 13C each scoping unit and cannot change meaning. 14C 15C Both members are REAL*8, as every declaring unit already gave them. 16C SAVE is explicit rather than relying on module variables being 17C implicitly saved, which holds only from Fortran 2008 while ifort 18C Classic predates it. That matches COMMON's semantics. 19C 20C The dimension parameters are PRIVATE: they are OPAL06 table 21C dimensions, which each importing unit also declares (issue 10). 22 MODULE EEEOS06_MOD 23 24 IMPLICIT NONE 25 SAVE 26 27 INTEGER*4, PARAMETER, PRIVATE :: MX = 5 28 INTEGER*4, PARAMETER, PRIVATE :: NR = 169 29 INTEGER*4, PARAMETER, PRIVATE :: NT = 197 30 31 REAL*8 :: EPL(MX,NT,NR) 32 REAL*8 :: XX(MX) 33 34 END MODULE EEEOS06_MOD

1"""Python version of eeeos06_mod.f (YREC 7be4133): workspace of the OPAL 22006 EOS interpolation, formerly COMMON/EEEOS06/. 3 4Translates: eeeos06_mod.f#EEEOS06_MOD (7be4133) 5 6Routines that USE EEEOS06_MOD import this module and write 7``eeeos06_mod.<name>``. Arrays keep Fortran's index order, with index 0 8unused. 9""" 10 11import numpy as np 12 13# Table dimensions (PRIVATE in the Fortran) (Fortran lines 27–29). 14_MX = 5 15_NR = 169 16_NT = 197 17 18# Fortran lines 31–32. 19epl = np.zeros((_MX + 1, _NT + 1, _NR + 1)) # one quantity's tables (esac06) 20xx = np.zeros(_MX + 1) # X of the tables, from xa (seteos06)