atmstruct_mod.f

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

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

Variables on this page (13)

abetas, adels, afxions, agam1, aqdt, atmocp, atmod, atmoo, atmop, atmor, atmot, json, numatm

Purpose

ATMSTRUCT_MOD holds the structure of the atmosphere, point by point, for the stitched output of the whole star. It replaces COMMON/ATMSTRUCT/ (MHP via Claude, 10/14/26); JSON comes from DIMENSIONS_MOD.

What its variables hold, who sets them, who reads them. All are set by ENVINT during the atmosphere integration, and read only by STITCH, which joins them to the envelope and the interior: NUMATM, the number of points; at each point $\log P$, $\log T$, $\log\rho$ (ATMOP, ATMOT, ATMOD) and ATMOR, not $\log r$ but the step DELTR in cm (ENVINT), which STITCH adds to the envelope's top radius, $10^{\tt ENVR} + {\tt ATMOR}$, the three gradients ADELS(3,·), ABETAS ($\beta$), AGAM1 ($\Gamma_1$), AQDT, the ionization fractions AFXIONS(3,·), the opacity ATMOO and $c_P$ ATMOCP.

Definitions

moduleATMSTRUCT_MOD

lines 19–41 · GitHub

Uses modules DIMENSIONS_MOD

Used by (2 units)

ENVINT, STITCH

Change log

The changes since rp_1

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

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C ATMSTRUCT_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/ATMSTRUCT/ 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. 19 MODULE ATMSTRUCT_MOD 20C 10/20/26 MHP via Claude 21C code changed to take JSON from DIMENSIONS_MOD instead of declaring 22C a PRIVATE copy here. 23 USE DIMENSIONS_MOD 24 25 IMPLICIT NONE 26 SAVE 27 28 REAL*8 :: ATMOP(JSON) 29 REAL*8 :: ATMOT(JSON) 30 REAL*8 :: ATMOD(JSON) 31 REAL*8 :: ATMOR(JSON) 32 REAL*8 :: ADELS(3,JSON) 33 REAL*8 :: ABETAS(JSON) 34 REAL*8 :: AGAM1(JSON) 35 REAL*8 :: AQDT(JSON) 36 REAL*8 :: AFXIONS(3,JSON) 37 REAL*8 :: ATMOO(JSON) 38 REAL*8 :: ATMOCP(JSON) 39 INTEGER*4 :: NUMATM 40 41 END MODULE ATMSTRUCT_MOD

1"""Python version of atmstruct_mod.f (YREC 7be4133): the atmosphere's 2structure, point by point, for the stitched output, formerly COMMON/ATMSTRUCT/. 3 4Translates: atmstruct_mod.f#ATMSTRUCT_MOD (7be4133) 5 6Module-level names, as in the Fortran module, in lower case. Routines that 7USE ATMSTRUCT_MOD import this module and write ``atmstruct_mod.<name>``. 8""" 9 10import numpy as np 11 12from pyrec import dimensions_mod 13 14_JSON = dimensions_mod.json 15 16# All set by envint, read by stitch; index 0 unused (Fortran lines 28–39). 17atmop = np.zeros(_JSON + 1) # log P 18atmot = np.zeros(_JSON + 1) # log T 19atmod = np.zeros(_JSON + 1) # log rho 20atmor = np.zeros(_JSON + 1) # log r 21adels = np.zeros((4, _JSON + 1)) # the three gradients 22abetas = np.zeros(_JSON + 1) # beta 23agam1 = np.zeros(_JSON + 1) # Gamma_1 24aqdt = np.zeros(_JSON + 1) 25afxions = np.zeros((4, _JSON + 1)) # ionization fractions 26atmoo = np.zeros(_JSON + 1) # opacity 27atmocp = np.zeros(_JSON + 1) # c_P 28numatm = 0 # the number of points