envstruct_mod.f

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

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

Variables on this page (19)

ebetas, edels, efxions, egam1, envd, envl, envo, envp, envr, envs, envt, envx, envz, eqcp, eqdt, evels, json, lcenv, numenv

Purpose

ENVSTRUCT_MOD holds the run of the envelope's structure, point by point, as stored by the envelope integration, shared by 7 routines. It replaces COMMON/ENVSTRUCT/ (MHP via Claude, 10/14/26); JSON comes from DIMENSIONS_MOD.

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

NUMENV (INTEGER*4): the number of stored points.
At each point (JSON each): ENVP, ENVT, ENVD, ENVR ($\log P$, $\log T$, $\log\rho$, $\log r$), ENVS ($\log(m/M_*)$), ENVX, ENVZ ($X$, $Z$), ENVO (opacity), ENVL (luminosity), LCENV (convective, EVEL > 0), EVELS (convective velocity), EDELS(3,JSON) (the three gradients), EBETAS, EGAM1 ($\Gamma_1$), EQCP, EQDT (thermodynamic derivatives), EFXIONS(3,JSON) (ionization fractions).
ENVINT fills them from ENVPRT_MOD, starting at the surface, and then reverses their order, so that point 1 is the bottom; it passes them to TAUCAL.
GETNEWENV and STARIN shift the structure variables when the fitting point moves (points above the new interior mass kept, offset by the jump in each variable), and use them to add envelope points to the interior model.
Read also by CALCAD, GETTAU, TAUINT (convective velocity beyond the interior) and STITCH (which writes them out).

Design

A blanket USE ENVSTRUCT_MOD. At rp_1, tauintnew.f also declared the block; it no longer does.

Definitions

moduleENVSTRUCT_MOD

lines 19–47 · GitHub

Uses modules DIMENSIONS_MOD

Used by (7 units)

CALCAD, ENVINT, GETNEWENV, GETTAU, STARIN, STITCH, TAUINT

Change log

The changes since rp_1

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

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C ENVSTRUCT_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/ENVSTRUCT/ into this module. 8C 9C Imported with a blanket USE, the convention for the conversion. 10C All 7 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 ENVSTRUCT_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 :: ENVP(JSON) 29 REAL*8 :: ENVT(JSON) 30 REAL*8 :: ENVS(JSON) 31 REAL*8 :: ENVD(JSON) 32 REAL*8 :: ENVR(JSON) 33 REAL*8 :: ENVX(JSON) 34 REAL*8 :: ENVZ(JSON) 35 LOGICAL*4 :: LCENV(JSON) 36 REAL*8 :: EDELS(3,JSON) 37 REAL*8 :: EVELS(JSON) 38 REAL*8 :: EBETAS(JSON) 39 REAL*8 :: EGAM1(JSON) 40 REAL*8 :: EQCP(JSON) 41 REAL*8 :: EFXIONS(3,JSON) 42 REAL*8 :: ENVO(JSON) 43 REAL*8 :: ENVL(JSON) 44 REAL*8 :: EQDT(JSON) 45 INTEGER*4 :: NUMENV 46 47 END MODULE ENVSTRUCT_MOD

1"""Python version of envstruct_mod.f (YREC 7be4133): the envelope's 2structure, point by point, formerly COMMON/ENVSTRUCT/. 3 4Translates: envstruct_mod.f#ENVSTRUCT_MOD (7be4133) 5 6Routines that USE ENVSTRUCT_MOD import this module and write 7``envstruct_mod.<name>``. Filled by envint (point 1 at the bottom after it 8reverses them); shifted by getnewenv and starin. 9""" 10 11import numpy as np 12 13from pyrec import dimensions_mod 14 15_JSON = dimensions_mod.json 16 17# Fortran lines 28–45. 18envp = np.zeros(_JSON + 1) # log P 19envt = np.zeros(_JSON + 1) # log T 20envs = np.zeros(_JSON + 1) # log(m/M) 21envd = np.zeros(_JSON + 1) # log rho 22envr = np.zeros(_JSON + 1) # log r 23envx = np.zeros(_JSON + 1) # X 24envz = np.zeros(_JSON + 1) # Z 25lcenv = np.zeros(_JSON + 1, dtype=bool) # convective point 26edels = np.zeros((4, _JSON + 1)) # the three gradients, edels[k, j] 27evels = np.zeros(_JSON + 1) # convective velocity 28ebetas = np.zeros(_JSON + 1) # beta 29egam1 = np.zeros(_JSON + 1) # Gamma_1 30eqcp = np.zeros(_JSON + 1) # thermodynamic derivative QQCP 31efxions = np.zeros((4, _JSON + 1)) # ionization fractions, efxions[k, j] 32envo = np.zeros(_JSON + 1) # opacity 33envl = np.zeros(_JSON + 1) # luminosity 34eqdt = np.zeros(_JSON + 1) # thermodynamic derivative QQDT 35numenv = 0 # number of stored points