dwmax_mod.f

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

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

Variables on this page (3)

json, qwrmax, qwrmx0

Purpose

DWMAX_MOD holds the critical rotation gradient for the shear instability at each shell, and its value at the start of the step, shared by PHYSIC and HPOINT. It replaces COMMON/DWMAX/ (MHP via Claude, 10/14/26); JSON comes from DIMENSIONS_MOD.

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

QWRMAX(JSON): the critical $d\Omega/d\ln r$ (Endal and Sofia 1978, as PHYSIC's comments cite), set by PHYSIC for the interface between shells IM-1 and IM: 0 where both shells are convective or the root is not real (physic.f, lines 155–186 ↗). Read only by HPOINT, to copy it (hpoint.f, line 855 ↗), so PHYSIC's computation reaches nothing (so also at rp_1): unused?
QWRMX0(JSON): the start-of-step copy, made by HPOINT when it saves the old model; read nowhere: unused?

Design

A blanket USE DWMAX_MOD. SECULR has a local array of the same name, QWRMAX, which SETUPV fills for the angular momentum transport; it is not this module's. At rp_1, MIDMOD also declared the block without using it.

Definitions

moduleDWMAX_MOD

lines 19–31 · GitHub

Uses modules DIMENSIONS_MOD

Used by (2 units)

HPOINT, PHYSIC

Change log

The changes since rp_1

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

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C DWMAX_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/DWMAX/ 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 DWMAX_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 :: QWRMAX(JSON) 29 REAL*8 :: QWRMX0(JSON) 30 31 END MODULE DWMAX_MOD

1"""Python version of dwmax_mod.f (YREC 7be4133): the critical rotation 2gradient of the shear instability, formerly COMMON/DWMAX/. 3 4Translates: dwmax_mod.f#DWMAX_MOD (7be4133) 5 6Routines that USE DWMAX_MOD import this module and write 7``dwmax_mod.<name>``. 8""" 9 10import numpy as np 11 12from pyrec import dimensions_mod 13 14_JSON = dimensions_mod.json 15 16# Fortran lines 28–29. 17qwrmax = np.zeros(_JSON + 1) # critical dOmega/dln r per interface (physic) 18qwrmx0 = np.zeros(_JSON + 1) # its start-of-step copy (hpoint; read nowhere)