aeos06_mod.f

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

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

Variables on this page (16)

dfs, dfsr, dfsx, esk, esk2, mf, midx, mv, mx, nr, nt, rho, t6a, t6list, xa, xz

Purpose

AEOS06_MOD holds the OPAL 2006 equation-of-state tables and their grids, shared by 5 routines. It replaces COMMON/AEOS06/ (MHP via Claude, 10/20/26), whose member M is now named MIDX; the table dimensions (MX = 5, MV = 10, NR = 169, NT = 197) are private to the module, since each user declares its own.

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

XZ(MX,MV,NT,NR): the tables, for each of the five $X$ values, room for ten quantities on the $(T_6, \rho)$ grid, of which the file fills nine (readcoeos06.f, lines 124–126 ↗; the rest stays at $10^{35}$); read in by READCOEOS06, read by ESAC06 and RHOOFP06.
T6LIST(NR,NT), RHO(NR), T6A(NT): the grids, from READCOEOS06; read by ESAC06, RHOOFP06, EQBOUND06 and T6RINTEOS06.
DFSX(MX), DFS(NT), DFSR(NR): inverse grid spacings in $X$, $T_6$ and $\rho$, set by READCOEOS06; read by ESAC06 and T6RINTEOS06.
ESK(NT,NR), ESK2(NT,NR): values interpolated in $X$, set by ESAC06 and read by it and T6RINTEOS06.
MIDX: a loop index over the $X$ tables, in READCOEOS06 (as it reads each table) and ESAC06; MF: the $X$ index, set in ESAC06.
XA(MX): the five $X$ values, 0, 0.2, 0.4, 0.6, 0.8 (a DATA statement moved from ESAC06); read by ESAC06 and RHOOFP06.

Design

A blanket USE AEOS06_MOD; the table alone holds $5 \times 10 \times 197 \times 169$ values (about 13 MB).

Definitions

moduleAEOS06_MOD

lines 27–52 · GitHub

Used by (6 units)

EQBOUND06, ESAC06, READCOEOS06, RHOOFP06, SETEOS06, T6RINTEOS06

Change log

The changes since rp_1

Changed since rp_1. New file: COMMON/AEOS06/ as a module, M renamed MIDX, XA's values moved from ESAC06.

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C AEOS06_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/20/26 MHP via Claude 7C code changed to convert COMMON/AEOS06/ 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 6 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 The types restate what every declaring unit already gave them 16C explicitly: MIDX and MF are INTEGER*4, everything else REAL*8. 17C SAVE is explicit rather than relying on module variables being 18C implicitly saved, which holds only from Fortran 2008 while ifort 19C Classic predates it. That matches COMMON's semantics. 20C 21C The dimension parameters are PRIVATE: they are OPAL06 table 22C dimensions, which each importing unit also declares (issue 10). 23C 24C XA's initializer moved here from esac06.f, because DATA on a 25C USE-associated name is rejected. The literal text is copied 26C verbatim, so the stored values are unchanged. 27 MODULE AEOS06_MOD 28 29 IMPLICIT NONE 30 SAVE 31 32 INTEGER*4, PARAMETER, PRIVATE :: MX = 5 33 INTEGER*4, PARAMETER, PRIVATE :: MV = 10 34 INTEGER*4, PARAMETER, PRIVATE :: NR = 169 35 INTEGER*4, PARAMETER, PRIVATE :: NT = 197 36 37 REAL*8 :: XZ(MX,MV,NT,NR) 38 REAL*8 :: T6LIST(NR,NT) 39 REAL*8 :: RHO(NR) 40 REAL*8 :: T6A(NT) 41 REAL*8 :: ESK(NT,NR) 42 REAL*8 :: ESK2(NT,NR) 43 REAL*8 :: DFSX(MX) 44 REAL*8 :: DFS(NT) 45 REAL*8 :: DFSR(NR) 46 INTEGER*4 :: MIDX 47 INTEGER*4 :: MF 48 REAL*8 :: XA(MX) 49 50 DATA XA/0.0,0.2,0.4,0.6,0.8/ 51 52 END MODULE AEOS06_MOD

1"""Python version of aeos06_mod.f (YREC 7be4133): the OPAL 2006 2equation-of-state tables and grids, formerly COMMON/AEOS06/. 3 4Translates: aeos06_mod.f#AEOS06_MOD (7be4133) 5 6Module-level names, as in the Fortran module, in lower case. Routines that 7USE AEOS06_MOD import this module and write ``aeos06_mod.<name>``. 8""" 9 10import numpy as np 11 12# Table dimensions (PRIVATE in the Fortran) (Fortran lines 32–35). 13_MX = 5 14_MV = 10 15_NR = 169 16_NT = 197 17 18# Fortran lines 37–47. 19xz = np.zeros((_MX + 1, _MV + 1, _NT + 1, _NR + 1)) # the tables (readcoeos06) 20t6list = np.zeros((_NR + 1, _NT + 1)) 21rho = np.zeros(_NR + 1) # density grid (readcoeos06) 22t6a = np.zeros(_NT + 1) # T6 grid (readcoeos06) 23esk = np.zeros((_NT + 1, _NR + 1)) # values interpolated in X (esac06) 24esk2 = np.zeros((_NT + 1, _NR + 1)) 25dfsx = np.zeros(_MX + 1) # inverse spacings in X, T6, rho 26dfs = np.zeros(_NT + 1) 27dfsr = np.zeros(_NR + 1) 28midx = 0 # formerly M (readcoeos06) 29mf = 0 # the X index (esac06) 30# The X values of the tables, index 0 unused (DATA, from esac06; the 31# Fortran's single-precision literals, as stored in REAL*8) (Fortran lines 32# 48–50). 33xa = np.array([0.0, 0.0, np.float32(0.2), np.float32(0.4), 34 np.float32(0.6), np.float32(0.8)], dtype=float)