intvar_mod.f

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

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

Variables on this page (12)

alm, delami, delmi, dm, epsilm, hgm, hs3, json, pm, qdtmi, rm, tm

Purpose

INTVAR_MOD holds the structure variables at the interfaces between shells (the midpoints in radius), for the rotation routines, shared by SETUPV, VCIRC, GETFC and SECULR. It replaces COMMON/INTVAR/ (MHP via Claude, 10/14/26); JSON comes from DIMENSIONS_MOD. All members are arrays of JSON.

What its variables hold, who sets them, who reads them. All are set by SETUPV:

RM: the radius of the interface between shells I-1 and I, the mean of the two shells' radii.
At each interface, by four-point Lagrangian interpolation in radius (the first and last interfaces from the first and last four points): PM, TM, DM ($P$, $T$, $\rho$, unlogged), ALM ($L$, cgs), HS3 (the mass), HGM (gravity), DELMI, DELAMI (the radiative and adiabatic gradients), QDTMI (a thermodynamic derivative), EPSILM (the energy generation rate).
SETUPV uses them to build the circulation coefficients; VCIRC reads RM, HGM, ALM; GETFC reads RM, DM; SECULR reads RM.

Design

A blanket USE INTVAR_MOD. At rp_1, CODIFF also declared the block without using it.

Definitions

moduleINTVAR_MOD

lines 19–40 · GitHub

Uses modules DIMENSIONS_MOD

Used by (4 units)

GETFC, SECULR, SETUPV, VCIRC

Change log

The changes since rp_1

Changed since rp_1. New file: COMMON/INTVAR/ as a module (CODIFF's unused declaration dropped).

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C INTVAR_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/INTVAR/ into this module. 8C 9C Imported with a blanket USE, the convention for the conversion. 10C All 4 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 INTVAR_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 :: ALM(JSON) 29 REAL*8 :: DELAMI(JSON) 30 REAL*8 :: DELMI(JSON) 31 REAL*8 :: DM(JSON) 32 REAL*8 :: EPSILM(JSON) 33 REAL*8 :: HGM(JSON) 34 REAL*8 :: HS3(JSON) 35 REAL*8 :: PM(JSON) 36 REAL*8 :: QDTMI(JSON) 37 REAL*8 :: RM(JSON) 38 REAL*8 :: TM(JSON) 39 40 END MODULE INTVAR_MOD

1"""Python version of intvar_mod.f (YREC 7be4133): the structure at the 2interfaces between shells, formerly COMMON/INTVAR/. 3 4Translates: intvar_mod.f#INTVAR_MOD (7be4133) 5 6Routines that USE INTVAR_MOD import this module and write 7``intvar_mod.<name>``. All set by setupv; index i is the interface between 8shells i-1 and i. 9""" 10 11import numpy as np 12 13from pyrec import dimensions_mod 14 15_JSON = dimensions_mod.json 16 17# The structure at the interfaces (Fortran lines 28–38). 18alm = np.zeros(_JSON + 1) # luminosity (cgs) 19delami = np.zeros(_JSON + 1) # adiabatic gradient 20delmi = np.zeros(_JSON + 1) # radiative gradient 21dm = np.zeros(_JSON + 1) # density 22epsilm = np.zeros(_JSON + 1) # energy generation rate 23hgm = np.zeros(_JSON + 1) # gravity 24hs3 = np.zeros(_JSON + 1) # mass 25pm = np.zeros(_JSON + 1) # pressure 26qdtmi = np.zeros(_JSON + 1) # thermodynamic derivative 27rm = np.zeros(_JSON + 1) # radius (mean of the two shells') 28tm = np.zeros(_JSON + 1) # temperature