alatm01_mod.f

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

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

Variables on this page (14)

alphas, fehs, gls, loldnextgen, nalpha, nfeh, nga, ngl, nta, nteff, p100ls, pls, t100ls, teffls

Purpose

ALATM01_MOD holds the Allard model-atmosphere table, shared by the four Allard routines. It replaces COMMON/ALATM01/ (MHP via Claude, 10/14/26); the sizes NGA = 25 and NTA = 250 are private, since each user declares its own.

What its variables hold, who sets them, who reads them. All are set by ALFILEIN, which reads the table file:

TEFFLs(NTA), GLs(NGA): the grid in $\log T_\mathrm{eff}$ and $\log g$, with their sizes nTEFF, nGL.
FeHs(NGA), ALPHAs(NGA), with nFeH, nALPHA: the metallicities and $\alpha$-enhancements found in the file, used only within ALFILEIN to choose the table.
PLs, P100Ls, T100Ls (NTA × NGA): $\log P$ where $T = T_\mathrm{eff}$, and $\log P$ and $\log T$ at $\tau = 100$ (for the option LALTPTau100).
LOldNextGen: the file's format, set and not read elsewhere.

Read by ALTABINIT (the table's limits), ALPRINT and ALSURFP (the interpolation).

Definitions

moduleALATM01_MOD

lines 23–44 · GitHub

Used by (4 units)

ALFILEIN, ALPRINT, ALSURFP, ALTABINIT

Change log

The changes since rp_1

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

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C ALATM01_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/ALATM01/ into this module. 8C 9C Imported with a blanket USE, the convention for the conversion. 10C All 4 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 ALATM01_MOD 24 25 IMPLICIT NONE 26 SAVE 27 28 INTEGER*4, PARAMETER, PRIVATE :: NGA = 25 29 INTEGER*4, PARAMETER, PRIVATE :: NTA = 250 30 31 REAL*8 :: TEFFLs(NTA) 32 REAL*8 :: GLs(NGA) 33 REAL*8 :: FeHs(NGA) 34 REAL*8 :: ALPHAs(NGA) 35 REAL*8 :: PLs(NTA,NGA) 36 REAL*8 :: P100Ls(NTA,NGA) 37 REAL*8 :: T100Ls(NTA,NGA) 38 LOGICAL*4 :: LOldNextGen 39 INTEGER*4 :: nTEFF 40 INTEGER*4 :: nGL 41 INTEGER*4 :: nFeH 42 INTEGER*4 :: nALPHA 43 44 END MODULE ALATM01_MOD

1"""Python version of alatm01_mod.f (YREC 7be4133): the Allard model-atmosphere table, formerly COMMON/ALATM01/. 2 3Translates: alatm01_mod.f#ALATM01_MOD (7be4133) 4 5Module-level names, as in the Fortran module, in lower case. Routines that 6USE ALATM01_MOD import this module and write ``alatm01_mod.<name>``. 7""" 8 9import numpy as np 10 11# Sizes, PRIVATE in the Fortran (Fortran lines 28–29). 12_NGA = 25 13_NTA = 250 14 15# All set by alfilein (Fortran lines 31–42). 16teffls = np.zeros(_NTA + 1) # log Teff grid 17gls = np.zeros(_NGA + 1) # log g grid 18fehs = np.zeros(_NGA + 1) # [Fe/H] values in the file 19alphas = np.zeros(_NGA + 1) # [alpha/Fe] values in the file 20pls = np.zeros((_NTA + 1, _NGA + 1)) # log P where T = Teff 21p100ls = np.zeros((_NTA + 1, _NGA + 1)) # log P at tau = 100 22t100ls = np.zeros((_NTA + 1, _NGA + 1)) # log T at tau = 100 23loldnextgen = False # the file's format 24nteff = 0 25ngl = 0 26nfeh = 0 27nalpha = 0