egrid_mod.f

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

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

Variables on this page (6)

chi, dchi, echi, es1, json, ntot

Purpose

EGRID_MOD holds the evenly spaced grid on which angular momentum transport and rotational mixing are solved: the grid variable $\chi$ at the model's shells and at the even points, the even points' masses, the spacing and their number. Shared by 7 routines. It replaces COMMON/EGRID/ (MHP via Claude, 10/14/26); JSON comes from DIMENSIONS_MOD.

What its variables hold, who sets them, who reads them. All are set by GETGRID for the region IBEG…IEND:

CHI(JSON): $\chi = \log m/\Delta_m + L/\Delta_L - \log P/\Delta_P$ at each shell of the region, where $\Delta_m$ = HPTTOL(2), $\Delta_L$ = HPTTOL(9)$\cdot L_\mathrm{surf}$ and $\Delta_P$ = HPTTOL(11), the mesh tolerances (GETGRID).
NTOT (INTEGER*4): the number of even points, the number of shells in the region but at least 3.
DCHI: the spacing, $(\chi_\mathrm{last} - \chi_1)/(\mathtt{NTOT}-1)$.
ECHI(JSON): $\chi$ at the even points, from CHI(1) in steps of DCHI.
ES1(JSON): the mass (g) at the even points, interpolated in $\log m$ by OSPLIN.
Read by ROTGRID and MIXGRID, which interpolate quantities between the two grids, MIXCOM, EQ2MOD, SECULR (which passes NTOT to the solvers) and CHECKJ (printing).

Design

A blanket USE EGRID_MOD. At rp_1, SECULR named the third member EQS1, harmless with COMMON's storage association.

Definitions

moduleEGRID_MOD

lines 19–34 · GitHub

Uses modules DIMENSIONS_MOD

Used by (7 units)

CHECKJ, EQ2MOD, GETGRID, MIXCOM, MIXGRID, ROTGRID, SECULR

Change log

The changes since rp_1

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

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C EGRID_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/EGRID/ into this module. 8C 9C Imported with a blanket USE, the convention for the conversion. 10C All 7 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 EGRID_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 :: CHI(JSON) 29 REAL*8 :: ECHI(JSON) 30 REAL*8 :: ES1(JSON) 31 REAL*8 :: DCHI 32 INTEGER*4 :: NTOT 33 34 END MODULE EGRID_MOD

1"""Python version of egrid_mod.f (YREC 7be4133): the evenly spaced grid of 2the rotation solvers, formerly COMMON/EGRID/. 3 4Translates: egrid_mod.f#EGRID_MOD (7be4133) 5 6Routines that USE EGRID_MOD import this module and write 7``egrid_mod.<name>``. All set by getgrid. 8""" 9 10import numpy as np 11 12from pyrec import dimensions_mod 13 14_JSON = dimensions_mod.json 15 16# Fortran lines 28–32. 17chi = np.zeros(_JSON + 1) # grid variable at the model's shells 18echi = np.zeros(_JSON + 1) # grid variable at the even points 19es1 = np.zeros(_JSON + 1) # mass (g) at the even points 20dchi = 0.0 # spacing of the even points in chi 21ntot = 0 # number of even points (at least 3)