mdphy_mod.f

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

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

Variables on this page (13)

amum, cpm, delam, delm, delrm, epsm, esumm, json, om, qdtm, thdifm, velm, viscm

Purpose

MDPHY_MOD holds physical quantities at each shell at an intermediate time of the model step, for the rotational-mixing and angular-momentum transport loop, shared by 10 routines. It replaces COMMON/MDPHY/ (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. MIDMOD sets them by linear interpolation, with weight FX, between the old model (…O) and the new one:

DELRM, DELAM, DELM: the radiative, adiabatic and actual temperature gradients; CPM ($c_P$), QDTM (a thermodynamic derivative), OM (opacity), ESUMM (energy generation, SESUM) and EPSM (the same plus SEG(6,J) and SEG(7,J)), THDIFM (thermal diffusivity), VISCM (viscosity). Read by SETUPV (which interpolates most of them to the interfaces, into INTVAR_MOD), VCIRC (THDIFM, VISCM), ROTMIX (DELM), ROTGRID and MIXGRID (EPSM), and LIBURN, LIBURN2 (DELAM, DELRM, to find the convective boundary); MIDMOD itself uses DELAM − DELRM to decide whether a shell that changes state over the step is convective.
VELM: the convective velocity, interpolated where the shell keeps its state, the larger of the two where it becomes convective, 0 where it becomes radiative (midmod.f, lines 178–188 ↗); read nowhere: unused?
AMUM: the mean molecular weight, set by MIDMOD, by CHECKC (from the new composition) and GETW (corrected for the change over the step); read by SETUPV and VCIRC (the μ-gradient terms).

Design

A blanket USE MDPHY_MOD. At rp_1, CODIFF and MIXCZ also declared the block; they no longer do.

Definitions

moduleMDPHY_MOD

lines 19–41 · GitHub

Uses modules DIMENSIONS_MOD

Used by (10 units)

CHECKC, GETW, LIBURN, LIBURN2, MIDMOD, MIXGRID, ROTGRID, ROTMIX, SETUPV, VCIRC

Change log

The changes since rp_1

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

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C MDPHY_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/MDPHY/ into this module. 8C 9C Imported with a blanket USE, the convention for the conversion. 10C All 10 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 MDPHY_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 :: AMUM(JSON) 29 REAL*8 :: CPM(JSON) 30 REAL*8 :: DELM(JSON) 31 REAL*8 :: DELAM(JSON) 32 REAL*8 :: DELRM(JSON) 33 REAL*8 :: ESUMM(JSON) 34 REAL*8 :: OM(JSON) 35 REAL*8 :: QDTM(JSON) 36 REAL*8 :: THDIFM(JSON) 37 REAL*8 :: VELM(JSON) 38 REAL*8 :: VISCM(JSON) 39 REAL*8 :: EPSM(JSON) 40 41 END MODULE MDPHY_MOD

1"""Python version of mdphy_mod.f (YREC 7be4133): the physics at each shell 2at an intermediate time of the model step, formerly COMMON/MDPHY/. 3 4Translates: mdphy_mod.f#MDPHY_MOD (7be4133) 5 6Routines that USE MDPHY_MOD import this module and write 7``mdphy_mod.<name>``. Set by midmod (amum also by checkc and getw). 8""" 9 10import numpy as np 11 12from pyrec import dimensions_mod 13 14_JSON = dimensions_mod.json 15 16# Fortran lines 28–39. 17amum = np.zeros(_JSON + 1) # mean molecular weight 18cpm = np.zeros(_JSON + 1) # c_P 19delm = np.zeros(_JSON + 1) # actual temperature gradient 20delam = np.zeros(_JSON + 1) # adiabatic gradient 21delrm = np.zeros(_JSON + 1) # radiative gradient 22esumm = np.zeros(_JSON + 1) # energy generation 23om = np.zeros(_JSON + 1) # opacity 24qdtm = np.zeros(_JSON + 1) # thermodynamic derivative 25thdifm = np.zeros(_JSON + 1) # thermal diffusivity 26velm = np.zeros(_JSON + 1) # convective velocity; read nowhere 27viscm = np.zeros(_JSON + 1) # viscosity 28epsm = np.zeros(_JSON + 1) # energy generation plus seg(6), seg(7)