llot95e_mod.f

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

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

Variables on this page (4)

cappae, numd, numt, xop95

Purpose

LLOT95E_MOD holds the OPAL95 opacity table interpolated to the surface composition (one $X$ and $Z$), and that $X$, shared by 3 routines. It replaces COMMON/LLOT95E/ (MHP via Claude, 10/14/26), with its member XTAB renamed XOP95 (a name collision resolved before the conversion). The table dimensions NUMT = 70, NUMD = 19 are PRIVATE.

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

CAPPAE(NUMT,NUMD): $\log\kappa$ on the ($\log T$, $\log R$) grid at $X$ = XOP95 and $Z$ = ZTAB (LLOT95_MOD), interpolated in $X$ (four points) by OP95XTAB, which SETUPOPAC and SURFOPAC call with the envelope's $X$. Read by OP952D (interpolation in $T$ and $R$).
XOP95: the $X$ of that table, set by OP95XTAB; read by GETOPAL95, which uses the 2D interpolation when the local $X$ and $Z$ match the table's to $10^{-4}$.

Design

A blanket USE LLOT95E_MOD. At rp_1, main.f declared the block with its own NUMTT (also 70); that declaration is gone.

Definitions

moduleLLOT95E_MOD

lines 23–34 · GitHub

Used by (3 units)

GETOPAL95, OP952D, OP95XTAB

Change log

The changes since rp_1

Changed since rp_1. New file: COMMON/LLOT95E/ as a module, XTAB renamed XOP95.

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C LLOT95E_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/LLOT95E/ into this module. 8C 9C Imported with a blanket USE, the convention for the conversion. 10C All 3 files declared this block identically, so the USE 11C imports exactly the names the declaration already made global in 12C each scoping unit and cannot change meaning. 13C 14C The types restate what the implicit rules supplied: L is 15C LOGICAL*4, I/J/K/M/N are INTEGER*4, everything else REAL*8. 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: the importing files 21C declare their own PARAMETER statements for these, and an 22C exported copy would collide with them. 23 MODULE LLOT95E_MOD 24 25 IMPLICIT NONE 26 SAVE 27 28 INTEGER*4, PARAMETER, PRIVATE :: NUMD = 19 29 INTEGER*4, PARAMETER, PRIVATE :: NUMT = 70 30 31 REAL*8 :: CAPPAE(NUMT,NUMD) 32 REAL*8 :: XOP95 33 34 END MODULE LLOT95E_MOD

1"""Python version of llot95e_mod.f (YREC 7be4133): the OPAL95 opacity table 2at the surface X and Z, formerly COMMON/LLOT95E/ (XTAB now XOP95). 3 4Translates: llot95e_mod.f#LLOT95E_MOD (7be4133) 5 6Routines that USE LLOT95E_MOD import this module and write 7``llot95e_mod.<name>``. Index 0 unused. 8""" 9 10import numpy as np 11 12# Table dimensions (PRIVATE in the Fortran) (Fortran lines 28–29). 13_NUMD = 19 14_NUMT = 70 15 16# The table at X = xop95 (Fortran lines 31–32). 17cappae = np.zeros((_NUMT + 1, _NUMD + 1)) # log kappa at X = xop95, Z = ztab (op95xtab) 18xop95 = 0.0 # the X of cappae (op95xtab)