envprt_mod.f

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

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

Variables on this page (10)

ebeta, ed, edel, efxion, eo, ep, er, es, et, evel

Purpose

ENVPRT_MOD holds the state at the current point of the envelope integration, passed from QENV (the right-hand side) to ENVINT (the integrator), which stores it. It replaces COMMON/ENVPRT/ (MHP via Claude, 10/14/26). All members are REAL*8.

What its variables hold, who sets them, who reads them. QENV sets them at each call ("07/02 always store the basic structure variables"):

EP, ET, ED, ER: $\log P$, $\log T$, $\log\rho$, $\log r$; ES: $\log(m/M_*)$ (SL - STOTAL).
EO: the opacity; EBETA: $\beta$; EVEL: the convective velocity; EDEL(3): the radiative and adiabatic gradients and the actual one (DELR, DELA, DEL); EFXION(3): the ionization fractions.
ENVINT reads them to print the point and to store the envelope's run in ENVSTRUCT_MOD.

Design

A blanket USE ENVPRT_MOD. QENV is called by the ODE solver through a fixed argument list, so its extra outputs go through the module. At rp_1, GETNEWENV and STARIN also declared the block without using it.

Definitions

moduleENVPRT_MOD

lines 19–35 · GitHub

Used by (2 units)

ENVINT, QENV

Change log

The changes since rp_1

Changed since rp_1. New file: COMMON/ENVPRT/ as a module (two unused declarations dropped).

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C ENVPRT_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/ENVPRT/ 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 ENVPRT_MOD 20 21 IMPLICIT NONE 22 SAVE 23 24 REAL*8 :: EP 25 REAL*8 :: ET 26 REAL*8 :: ER 27 REAL*8 :: ES 28 REAL*8 :: ED 29 REAL*8 :: EO 30 REAL*8 :: EBETA 31 REAL*8 :: EDEL(3) 32 REAL*8 :: EFXION(3) 33 REAL*8 :: EVEL 34 35 END MODULE ENVPRT_MOD

1"""Python version of envprt_mod.f (YREC 7be4133): the state at the current 2point of the envelope integration, formerly COMMON/ENVPRT/. 3 4Translates: envprt_mod.f#ENVPRT_MOD (7be4133) 5 6Routines that USE ENVPRT_MOD import this module and write 7``envprt_mod.<name>``. Set by qenv, read by envint. 8""" 9 10import numpy as np 11 12# Fortran lines 24–33. 13ep = 0.0 # log P 14et = 0.0 # log T 15er = 0.0 # log r 16es = 0.0 # log(m/M) 17ed = 0.0 # log rho 18eo = 0.0 # opacity 19ebeta = 0.0 # beta 20edel = np.zeros(4) # radiative, adiabatic, actual gradient; index 0 unused 21efxion = np.zeros(4) # ionization fractions; index 0 unused 22evel = 0.0 # convective velocity