oldphy_mod.f

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

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

Variables on this page (13)

amuo, cpo, delao, delo, delro, epso, esumo, json, oo, qdto, thdifo, velo, visco

Purpose

OLDPHY_MOD holds physical quantities of the model at the start of the step, which MIDMOD interpolates against the new model for the rotational-mixing substeps, shared by HPOINT, MIDMOD and GETW. It replaces COMMON/OLDPHY/ (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. HPOINT saves them after it has remapped the model ("for diffusion store the auxiliary quantities needed to calculate velocities at the start of the timestep with the new point distribution"):

DELRO, DELO, DELAO: the radiative, actual and adiabatic gradients (SDEL(1..3)); AMUO (mean molecular weight), OO (opacity), CPO ($c_P$), QDTO, VELO (convective velocity), VISCO (viscosity), THDIFO (thermal diffusivity).
ESUMO, EPSO: the energy generation, without and with the SEG(6), SEG(7) terms, interpolated by HPOINT onto the new grid ("MHP 6/00").
All are read by MIDMOD, which forms the mid-step values of MDPHY_MOD from them; GETW reads AMUO (the change of μ over the step).

Design

A blanket USE OLDPHY_MOD. HPOINT also uses DELAO − DELO for DDELM in OLDROT2_MOD.

Definitions

moduleOLDPHY_MOD

lines 19–41 · GitHub

Uses modules DIMENSIONS_MOD

Used by (3 units)

GETW, HPOINT, MIDMOD

Change log

The changes since rp_1

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

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C OLDPHY_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/OLDPHY/ into this module. 8C 9C Imported with a blanket USE, the convention for the conversion. 10C All 3 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 OLDPHY_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 :: DELO(JSON) 29 REAL*8 :: DELAO(JSON) 30 REAL*8 :: AMUO(JSON) 31 REAL*8 :: OO(JSON) 32 REAL*8 :: CPO(JSON) 33 REAL*8 :: QDTO(JSON) 34 REAL*8 :: VELO(JSON) 35 REAL*8 :: VISCO(JSON) 36 REAL*8 :: THDIFO(JSON) 37 REAL*8 :: ESUMO(JSON) 38 REAL*8 :: DELRO(JSON) 39 REAL*8 :: EPSO(JSON) 40 41 END MODULE OLDPHY_MOD

1"""Python version of oldphy_mod.f (YREC 7be4133): the physics of the model 2at the start of the step, formerly COMMON/OLDPHY/. 3 4Translates: oldphy_mod.f#OLDPHY_MOD (7be4133) 5 6Routines that USE OLDPHY_MOD import this module and write 7``oldphy_mod.<name>``. Saved by hpoint, read by midmod (amuo also 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. 17delo = np.zeros(_JSON + 1) # actual temperature gradient 18delao = np.zeros(_JSON + 1) # adiabatic gradient 19amuo = np.zeros(_JSON + 1) # mean molecular weight 20oo = np.zeros(_JSON + 1) # opacity 21cpo = np.zeros(_JSON + 1) # c_P 22qdto = np.zeros(_JSON + 1) # thermodynamic derivative 23velo = np.zeros(_JSON + 1) # convective velocity 24visco = np.zeros(_JSON + 1) # viscosity 25thdifo = np.zeros(_JSON + 1) # thermal diffusivity 26esumo = np.zeros(_JSON + 1) # energy generation 27delro = np.zeros(_JSON + 1) # radiative gradient 28epso = np.zeros(_JSON + 1) # energy generation plus seg(6), seg(7)