alot06_mod.f

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

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

Variables on this page (8)

capalex06, id, it, ix, numd, numt, xe, ze

Purpose

ALOT06_MOD holds the table, for the current composition, of the low-temperature opacities of Ferguson et al. (2005), called ALEX06 in the code (READALEX06), with the composition and the starting indices of the search. It replaces COMMON/ALOT06/ (MHP via Claude, 10/14/26); the sizes NUMD = 19 and NUMT = 85 are private. Five members carry initial values, moved from READALEX06's DATA statement (at rp_1); the module's header says they came from PARMIN's DATA statements and calls them namelist variables, but they are in no namelist (a comment copied from the other conversions).

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

CAPALEX06(NUMT,NUMD): $\log\kappa$ on the $(T, R)$ grid, interpolated to XE, ZE by ALEX06TAB (READALEX06 first fills it with the first composition's table); read by GETALEX06.
XE, ZE (initially 0): the composition of the current table; GETALEX06 resets them, and calls ALEX06TAB, when the point's $X$ or $Z$ differs by more than $10^{-8}$.
IX (initially 4): the first of the $X$ tables used, set by ALEX06TAB; IT (43) and ID (10): the starting indices in $T$ and $R$ of the search in GETALEX06, which updates them, so each search starts from the last point.

Definitions

moduleALOT06_MOD

lines 29–44 · GitHub

Used by (3 units)

ALEX06TAB, GETALEX06, READALEX06

Change log

The changes since rp_1

Changed since rp_1. New file: COMMON/ALOT06/ as a module, with the initial values of XE, ZE, IX, IT, ID moved from READALEX06.

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C ALOT06_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/ALOT06/ 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 imports 11C exactly the names the declaration already made global in each 12C 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 18C ifort Classic predates it. That matches COMMON's semantics. 19C 20C 5 of the 6 members carry an initializer, migrated 21C from parmin.f's DATA statements because DATA on a USE-associated 22C name is rejected. The literal text is copied verbatim: these 23C are namelist variables whose DATA supplied a default the 24C namelist read could override, and a module initializer is 25C applied once at load time, before PARMIN runs, which is the 26C same thing. Rewriting a literal such as 0.9828 as 0.9828D0 27C would change its stored value, so that is deliberately not 28C done here; see KNOWN_ISSUES.md issue 23. 29 MODULE ALOT06_MOD 30 31 IMPLICIT NONE 32 SAVE 33 34 INTEGER*4, PARAMETER, PRIVATE :: NUMD = 19 35 INTEGER*4, PARAMETER, PRIVATE :: NUMT = 85 36 37 REAL*8 :: CAPALEX06(NUMT,NUMD) 38 REAL*8 :: XE = 0.0D0 39 REAL*8 :: ZE = 0.0D0 40 INTEGER*4 :: IX = 4 41 INTEGER*4 :: IT = 43 42 INTEGER*4 :: ID = 10 43 44 END MODULE ALOT06_MOD

1"""Python version of alot06_mod.f (YREC 7be4133): the low-temperature opacity 2table for the current composition, formerly COMMON/ALOT06/. 3 4Translates: alot06_mod.f#ALOT06_MOD (7be4133) 5 6Module-level names, as in the Fortran module, in lower case. Routines that 7USE ALOT06_MOD import this module and write ``alot06_mod.<name>``. 8""" 9 10import numpy as np 11 12_NUMD = 19 # PRIVATE in the Fortran (lines 34–35) 13_NUMT = 85 14 15# The table for the current composition (Fortran line 37). 16capalex06 = np.zeros((_NUMT + 1, _NUMD + 1)) # log kappa (alex06tab) 17 18# Initializers, from readalex06's DATA (Fortran lines 38–42). 19xe = 0.0 # composition of the current table 20ze = 0.0 21ix = 4 # first X table used (alex06tab) 22it = 43 # starting indices of getalex06's search 23id = 10