oldmod_mod.f

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

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

Variables on this page (12)

hcompp, hdo, hlo, hlsuro, hpo, hro, hss, hto, json, lco, lczo, tefflo

Purpose

OLDMOD_MOD holds the model at the start of the step (the "old" model): its structure, composition, effective temperature, and surface luminosity, shared by 17 routines. It replaces COMMON/OLDMOD/ (MHP via Claude, 10/14/26); JSON comes from DIMENSIONS_MOD.

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

HPO, HTO, HRO, HLO, HDO (JSON each): $\log P$, $\log T$, $\log r$, $L$ and $\log\rho$ of the old model. Saved by HPOINT (when it remaps the model onto a new grid, interpolating them with OSPLIN), by the main program ("save old model for PTIME") and by STARIN. Read by PTIME (the time step from the structure changes), MIDMOD (interpolation within the step), LIRATE88 (the old rates) and LIBURN, LIBURN2.
HCOMPP(15,JSON): the old composition, saved by the main program, HPOINT, GETW and BURSMIX; read by the burning routines (KEMCOM, EQBURN, DBURN, DBURNM, LIRATE88), MIDMOD, CHECKC and CRRECT.
HSS(JSON): the old mass coordinates, against which HPOINT interpolates (and the new grid while HPOINT builds it).
LCO, LCZO (JSON each): the old convective flags, saved by HPOINT (LCZO also by STARIN); LCO is read by MIDMOD, LCZO by CRRECT, which passes it to GETROT.
TEFFLO, HLSURO: the old $\log T_\mathrm{eff}$ ("JVS 04/14") and surface luminosity, saved by the main program, HPOINT (before it remaps the model) and STARIN; read by ENTIME, which compares them with the new model.

Design

A blanket USE OLDMOD_MOD. HPOINT also uses HPO, HTO and HDO as scratch arrays ("note HDO is being used as a dummy array here"), so between HPOINT and the next save they do not hold the old model. At rp_1, MIX and SECULR also declared the block without using it; they no longer do.

Definitions

moduleOLDMOD_MOD

lines 19–43 · GitHub

Uses modules DIMENSIONS_MOD

Used by (17 units)

(file level), BURSMIX, CHECKC, CRRECT, DBURN, DBURNM, ENTIME, EQBURN, GETW, HPOINT, KEMCOM, LIBURN, LIBURN2, LIRATE88, MIDMOD, PTIME, STARIN

Change log

The changes since rp_1

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

Changed since a8a5950. The point count MO is replaced by HLSURO, the old model's surface luminosity, which is what ENTIME needs: HLO(MO) was not the surface once HPOINT had changed the number of points. Since CRRECT passes LCZO to GETROT (crrect.f changed, not this file), LCZO is read, which settles former question 1.

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C OLDMOD_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/OLDMOD/ into this module. 8C 9C Imported with a blanket USE, the convention for the conversion. 10C All 17 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 OLDMOD_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 :: HPO(JSON) 29 REAL*8 :: HTO(JSON) 30 REAL*8 :: HRO(JSON) 31 REAL*8 :: HLO(JSON) 32 REAL*8 :: HDO(JSON) 33 REAL*8 :: HCOMPP(15,JSON) 34 REAL*8 :: HSS(JSON) 35 LOGICAL*4 :: LCO(JSON) 36 LOGICAL*4 :: LCZO(JSON) 37 REAL*8 :: TEFFLO 38C 10/29/26 MHP via Claude 39C code changed to replace the point count MO with the previous model's 40C surface luminosity HLSURO, which is what ENTIME needs. 41 REAL*8 :: HLSURO 42 43 END MODULE OLDMOD_MOD

1"""Python version of oldmod_mod.f (YREC 7be4133): the model at the start of 2the step, formerly COMMON/OLDMOD/. 3 4Translates: oldmod_mod.f#OLDMOD_MOD (7be4133) 5 6Routines that USE OLDMOD_MOD import this module and write 7``oldmod_mod.<name>``. Saved by main, hpoint and starin (hpoint also uses 8hpo, hto, hdo as scratch arrays). 9""" 10 11import numpy as np 12 13from pyrec import dimensions_mod 14 15_JSON = dimensions_mod.json 16 17# The old model, its composition and flags (Fortran lines 28–41). 18hpo = np.zeros(_JSON + 1) # log P 19hto = np.zeros(_JSON + 1) # log T 20hro = np.zeros(_JSON + 1) # log r 21hlo = np.zeros(_JSON + 1) # L 22hdo = np.zeros(_JSON + 1) # log rho 23hcompp = np.zeros((16, _JSON + 1)) # composition, hcompp[k, j] 24hss = np.zeros(_JSON + 1) # mass coordinates 25lco = np.zeros(_JSON + 1, dtype=bool) # convective flags (read by midmod) 26lczo = np.zeros(_JSON + 1, dtype=bool) # read by crrect (for getrot) 27tefflo = 0.0 # log Teff 28hlsuro = 0.0 # surface L (for entime)