atmos2_mod.f

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

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

Variables on this page (6)

atmgl, atmpl, atmtl, atmz, ng, nt

Purpose

ATMOS2_MOD holds the surface-pressure table read when KTTAU is 3 or 4 (used only with 3; see below): $\log P$ on a grid of 57 values of $\log T_\mathrm{eff}$ and 11 of $\log g$. It replaces COMMON/ATMOS2/ (MHP via Claude, 10/14/26); the sizes NG, NT are private.

What its variables hold, who sets them, who reads them. SETUPS reads them all from the file FATM: ATMZ, the table's metallicity; ATMTL(NT), the $\log T_\mathrm{eff}$ grid; ATMGL(NG), the $\log g$ grid; ATMPL(NT,NG), $\log P$. SURFP interpolates in them. ATMZ is read by no routine; SETUPS also stores in it the metallicity of the KTTAU = 5 table (ATMOS2C_MOD).

Why KTTAU = 4 loads it. With KTTAU = 4 (Allard atmospheres, ALFILEIN, ALSURFP) SETUPS reads this Kurucz table too, but no routine then uses it: above the Allard table's top $T_\mathrm{eff}$ (AL_TEFFLmax, ALATM05_MOD), or off its gravities, SURFBC and ENVINT switch to the gray atmosphere (KTTAU = 0), not to Kurucz, although ALSURFP's header names Kurucz (KTTAU = 3) as a fallback the caller could choose.

Definitions

moduleATMOS2_MOD

lines 23–36 · GitHub

Used by (2 units)

SETUPS, SURFP

Change log

The changes since rp_1

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

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C ATMOS2_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/ATMOS2/ into this module. 8C 9C Imported with a blanket USE, the convention for the conversion. 10C All 2 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 ATMOS2_MOD 24 25 IMPLICIT NONE 26 SAVE 27 28 INTEGER*4, PARAMETER, PRIVATE :: NG = 11 29 INTEGER*4, PARAMETER, PRIVATE :: NT = 57 30 31 REAL*8 :: ATMPL(NT,NG) 32 REAL*8 :: ATMTL(NT) 33 REAL*8 :: ATMGL(NG) 34 REAL*8 :: ATMZ 35 36 END MODULE ATMOS2_MOD

1"""Python version of atmos2_mod.f (YREC 7be4133): the surface-pressure table 2for KTTAU = 3 or 4, formerly COMMON/ATMOS2/. 3 4Translates: atmos2_mod.f#ATMOS2_MOD (7be4133) 5 6Module-level names, as in the Fortran module, in lower case. Routines that 7USE ATMOS2_MOD import this module and write ``atmos2_mod.<name>``. 8""" 9 10import numpy as np 11 12# Fortran lines 28–29. 13_NG = 11 # PRIVATE in the Fortran 14_NT = 57 15 16# All read by setups from the file FATM (Fortran lines 31–34). 17atmpl = np.zeros((_NT + 1, _NG + 1)) # log P 18atmtl = np.zeros(_NT + 1) # log Teff grid 19atmgl = np.zeros(_NG + 1) # log g grid 20atmz = 0.0 # the table's Z (read by no routine)