gravez_mod.f

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

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

Variables on this page (8)

ecod1z, ecod1z_h, ecod2z_h, eqcod1z_h, eqcod2z_h, ez, ez_h, json

Purpose

GRAVEZ_MOD holds the metals ($Z$) on the equally spaced radial grid of the older settling routine GRSETT, with their diffusion coefficients at the grid's points and midpoints, used with metal diffusion (LDIFZ). Shared by 5 routines. It replaces COMMON/GRAVEZ/ (MHP via Claude, 10/14/26); JSON comes from DIMENSIONS_MOD.

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

EZ(JSON): $Z$ at the grid points, interpolated from the model by MODEL_TO_EQUAL, advanced by LAX_WENDROF2, and turned into the change over the step by GRSETT (EZ - EZ_ORIG), which EQUAL_TO_MODEL applies to the model's points.
EZ_H(JSON): $Z$ at the midpoints, set by MODEL_TO_EQUAL and by LAX_WENDROF1 (the half step); read by GRSETT.
ECOD1Z(JSON): the settling coefficient at the grid points, set by MODEL_TO_EQUAL, read by LAX_WENDROF1.
ECOD1Z_H(JSON): the settling coefficient at the midpoints, set by MODEL_TO_EQUAL and updated by GRSETT with the half-step abundances (ECOD1Z_H + EZ_H*EQCOD1Z_H); read by LAX_WENDROF2.
EQCOD1Z_H(JSON): the abundance-dependent part of that coefficient, set by MODEL_TO_EQUAL, read by GRSETT.
ECOD2Z_H, EQCOD2Z_H (JSON each): the diffusion coefficient and its abundance-dependent part at the midpoints, set by MODEL_TO_EQUAL (model_to_equal.f, lines 105–107 ↗; model_to_equal.f, lines 180–187 ↗); read nowhere (unused?), since the metals get only the explicit settling term (see GRSETT's notes).

Design

A blanket USE GRAVEZ_MOD. The hydrogen counterparts (EX, ECOD1, …) are arguments of the same routines; only the metals go through the module.

Definitions

moduleGRAVEZ_MOD

lines 19–36 · GitHub

Uses modules DIMENSIONS_MOD

Used by (5 units)

EQUAL_TO_MODEL, GRSETT, LAX_WENDROF1, LAX_WENDROF2, MODEL_TO_EQUAL

Change log

The changes since rp_1

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

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C GRAVEZ_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/GRAVEZ/ into this module. 8C 9C Imported with a blanket USE, the convention for the conversion. 10C All 5 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 GRAVEZ_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 :: ECOD1Z(JSON) 29 REAL*8 :: ECOD1Z_H(JSON) 30 REAL*8 :: ECOD2Z_H(JSON) 31 REAL*8 :: EQCOD1Z_H(JSON) 32 REAL*8 :: EQCOD2Z_H(JSON) 33 REAL*8 :: EZ(JSON) 34 REAL*8 :: EZ_H(JSON) 35 36 END MODULE GRAVEZ_MOD

1"""Python version of gravez_mod.f (YREC 7be4133): the metals and their 2coefficients on grsett's equally spaced grid, formerly COMMON/GRAVEZ/. 3 4Translates: gravez_mod.f#GRAVEZ_MOD (7be4133) 5 6Routines that USE GRAVEZ_MOD import this module and write 7``gravez_mod.<name>``. "_h" marks the grid's midpoints. 8""" 9 10import numpy as np 11 12from pyrec import dimensions_mod 13 14_JSON = dimensions_mod.json 15 16ecod1z = np.zeros(_JSON + 1) # settling coefficient at the points (model_to_equal) 17ecod1z_h = np.zeros(_JSON + 1) # the same at the midpoints (model_to_equal, grsett) 18ecod2z_h = np.zeros(_JSON + 1) # diffusion coefficient (model_to_equal); read nowhere 19eqcod1z_h = np.zeros(_JSON + 1) # abundance-dependent part of ecod1z_h 20eqcod2z_h = np.zeros(_JSON + 1) # the same for ecod2z_h; read nowhere 21ez = np.zeros(_JSON + 1) # Z at the points, then its change (grsett) 22ez_h = np.zeros(_JSON + 1) # Z at the midpoints (model_to_equal, lax_wendrof1)