atmos2c_mod.f

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

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

Variables on this page (5)

atmglc, atmplc, atmtlc, ngc, ntc

Purpose

ATMOS2C_MOD holds the surface-pressure table used when KTTAU is 5: $\log P$ on a grid of 76 values of $\log T_\mathrm{eff}$ and 11 of $\log g$. It replaces COMMON/ATMOS2C/ (MHP via Claude, 10/14/26); the sizes NGC, NTC are private.

What its variables hold, who sets them, who reads them. SETUPS reads ATMTLC(NTC) ($\log T_\mathrm{eff}$), ATMGLC(NGC) ($\log g$) and ATMPLC(NTC,NGC) ($\log P$) from the file FATM (the table's metallicity goes to ATMZ in ATMOS2_MOD); KCSURFP interpolates in them.

Definitions

moduleATMOS2C_MOD

lines 23–35 · GitHub

Used by (2 units)

KCSURFP, SETUPS

Change log

The changes since rp_1

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

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C ATMOS2C_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/ATMOS2C/ 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 ATMOS2C_MOD 24 25 IMPLICIT NONE 26 SAVE 27 28 INTEGER*4, PARAMETER, PRIVATE :: NGC = 11 29 INTEGER*4, PARAMETER, PRIVATE :: NTC = 76 30 31 REAL*8 :: ATMPLC(NTC,NGC) 32 REAL*8 :: ATMTLC(NTC) 33 REAL*8 :: ATMGLC(NGC) 34 35 END MODULE ATMOS2C_MOD

1"""Python version of atmos2c_mod.f (YREC 7be4133): the surface-pressure table for KTTAU = 5, formerly COMMON/ATMOS2C/. 2 3Translates: atmos2c_mod.f#ATMOS2C_MOD (7be4133) 4 5Module-level names, as in the Fortran module, in lower case. Routines that 6USE ATMOS2C_MOD import this module and write ``atmos2c_mod.<name>``. 7""" 8 9import numpy as np 10 11# Fortran lines 28–29. 12_NGC = 11 # PRIVATE in the Fortran 13_NTC = 76 14 15# All read by setups from the file FATM (Fortran lines 31–33). 16atmplc = np.zeros((_NTC + 1, _NGC + 1)) # log P 17atmtlc = np.zeros(_NTC + 1) # log Teff grid 18atmglc = np.zeros(_NGC + 1) # log g grid