entrop_mod.f

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

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

Variables on this page (5)

hll0, hpp0, hrr0, htt0, json

Purpose

ENTROP_MOD holds the rates of change of the structure variables in the last model, used for a better first guess of the next one (LNEWS), shared by the main program and HPOINT. It replaces COMMON/ENTROP/ (MHP via Claude, 10/14/26); JSON comes from DIMENSIONS_MOD.

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

HTT0, HPP0 (JSON each): the changes of $\log T$ and $\log P$ over the last step, divided by the step DELTS; HLL0: the relative change of $L$, $2(L - L_\mathrm{old})/(L + L_\mathrm{old})$, per unit time; HRR0: the change of $\log r$ per unit time.
The main program zeroes them at the start of a run and when LNEWS is off or the step is not positive, sets them after each converged model ("7/91 store changes in the structure"), and reads them to estimate the next model's structure: $\Delta T = $ HTT0*DELTS, and so on.
HPOINT interpolates them onto the new mass grid when points are added or removed ("7/91 add entropy term interpolation").

Design

A blanket USE ENTROP_MOD. The comments call them "entropy terms," but they hold time derivatives of the structure.

Definitions

moduleENTROP_MOD

lines 19–33 · GitHub

Uses modules DIMENSIONS_MOD

Used by (2 units)

(file level), HPOINT

Change log

The changes since rp_1

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

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C ENTROP_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/ENTROP/ 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 ENTROP_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 :: HTT0(JSON) 29 REAL*8 :: HPP0(JSON) 30 REAL*8 :: HLL0(JSON) 31 REAL*8 :: HRR0(JSON) 32 33 END MODULE ENTROP_MOD

1"""Python version of entrop_mod.f (YREC 7be4133): the rates of change of 2the structure in the last model, formerly COMMON/ENTROP/. 3 4Translates: entrop_mod.f#ENTROP_MOD (7be4133) 5 6Routines that USE ENTROP_MOD import this module and write 7``entrop_mod.<name>``. Set by main, remapped by hpoint. 8""" 9 10import numpy as np 11 12from pyrec import dimensions_mod 13 14_JSON = dimensions_mod.json 15 16# The rates of change of the structure (Fortran lines 28–31). 17htt0 = np.zeros(_JSON + 1) # d log T / dt in the last step 18hpp0 = np.zeros(_JSON + 1) # d log P / dt 19hll0 = np.zeros(_JSON + 1) # relative dL / dt 20hrr0 = np.zeros(_JSON + 1) # d log r / dt