atmprt_mod.f

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

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

Variables on this page (6)

ad, afxion, ao, ap, at, taul

Purpose

ATMPRT_MOD holds the state at the last point of the atmosphere integration (or the fitting point of a tabulated atmosphere), for the envelope integration and the output. It replaces COMMON/ATMPRT/ (MHP via Claude, 10/14/26), where the optical depth was named TAU; it is now TAUL.

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

TAUL ($\log\tau$), AP ($\log P$), AT ($\log T$): set by QATM at each point of the atmosphere integration; AP and AT also by the table routines SURFP and KCSURFP, and AT by ALSURFP.
AD ($\log\rho$), AO (the opacity), AFXION(3) (the ionization fractions): set by QATM only when the model is stored or written for pulsation (LSAVE, LPUMOD).
Read by ENVINT, which writes them with the stored model, and by MASSLOSS (AP, AT) and the table routines.

Definitions

moduleATMPRT_MOD

lines 19–31 · GitHub

Used by (6 units)

ALSURFP, ENVINT, KCSURFP, MASSLOSS, QATM, SURFP

Change log

The changes since rp_1

Changed since rp_1. New file: COMMON/ATMPRT/ as a module, TAU renamed TAUL.

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C ATMPRT_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/ATMPRT/ into this module. 8C 9C Imported with a blanket USE, the convention for the conversion. 10C All 6 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. 19 MODULE ATMPRT_MOD 20 21 IMPLICIT NONE 22 SAVE 23 24 REAL*8 :: TAUL 25 REAL*8 :: AP 26 REAL*8 :: AT 27 REAL*8 :: AD 28 REAL*8 :: AO 29 REAL*8 :: AFXION(3) 30 31 END MODULE ATMPRT_MOD

1"""Python version of atmprt_mod.f (YREC 7be4133): the state at the last point 2of the atmosphere, formerly COMMON/ATMPRT/. 3 4Translates: atmprt_mod.f#ATMPRT_MOD (7be4133) 5 6Module-level names, as in the Fortran module, in lower case. Routines that 7USE ATMPRT_MOD import this module and write ``atmprt_mod.<name>``. 8""" 9 10import numpy as np 11 12# Fortran lines 24–29. 13taul = 0.0 # log tau (qatm; formerly TAU) 14ap = 0.0 # log P (qatm, surfp, kcsurfp) 15at = 0.0 # log T (qatm, surfp, kcsurfp, alsurfp) 16ad = 0.0 # log rho (qatm, when the model is stored) 17ao = 0.0 # opacity (qatm, likewise) 18afxion = np.zeros(4) # ionization fractions, index 0 unused (qatm)