egridder_mod.f

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

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

Variables on this page (5)

eqqcod, eqqqcod, json, qqcod, qqqcod

Purpose

EGRIDDER_MOD holds the coefficients of the higher-derivative terms of the diffusion-plus-advection treatment of angular momentum (LDIFAD), shared by SETUPV, VCIRC and DADCOEFT. It replaces COMMON/EGRIDDER/ (MHP via Claude, 10/14/26); JSON comes from DIMENSIONS_MOD.

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

QQCOD(JSON): $\tfrac23\,r^4/(G\,m)$-type factors at each shell (CC23*RM**4/CG/HS3), set by SETUPV; VCIRC multiplies them by $\Omega$ (WM) into EQQCOD.
EQQCOD(JSON): read by DADCOEFT as the coefficients of its $d^2\Omega/dr^2$ rows.
QQQCOD(JSON): factors set by SETUPV (HTM, HTP, setupv.f, line 410 ↗, setupv.f, line 435 ↗); read nowhere: unused?
EQQQCOD(JSON): read by DADCOEFT as the coefficients of its $d^3\Omega/dr^3$ rows (dadcoeft.f, lines 371–373 ↗), but set by no routine (so at rp_1 too), so those coefficients are zero. This gap is part of DADCOEFT's suspected bug 1, whose possible fix names QQQCOD as the model-grid version to map onto the even grid (only with LDIFAD, off by default).

Design

A blanket USE EGRIDDER_MOD. The pairs mirror each other: VCIRC turns QQCOD into EQQCOD, but nothing turns QQQCOD into EQQQCOD, the same gap as in DIFAD3_MOD.

Definitions

moduleEGRIDDER_MOD

lines 19–33 · GitHub

Uses modules DIMENSIONS_MOD

Used by (3 units)

DADCOEFT, SETUPV, VCIRC

Change log

The changes since rp_1

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

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C EGRIDDER_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/EGRIDDER/ 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 EGRIDDER_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 :: EQQCOD(JSON) 29 REAL*8 :: EQQQCOD(JSON) 30 REAL*8 :: QQCOD(JSON) 31 REAL*8 :: QQQCOD(JSON) 32 33 END MODULE EGRIDDER_MOD

1"""Python version of egridder_mod.f (YREC 7be4133): the higher-derivative 2coefficients of the diffusion-plus-advection treatment, formerly 3COMMON/EGRIDDER/. 4 5Translates: egridder_mod.f#EGRIDDER_MOD (7be4133) 6 7Routines that USE EGRIDDER_MOD import this module and write 8``egridder_mod.<name>``. 9""" 10 11import numpy as np 12 13from pyrec import dimensions_mod 14 15_JSON = dimensions_mod.json 16 17# Fortran lines 28–31. 18eqqcod = np.zeros(_JSON + 1) # qqcod * Omega (vcirc); read by dadcoeft 19# Question (see the page): eqqqcod is read by dadcoeft but set nowhere. 20eqqqcod = np.zeros(_JSON + 1) 21qqcod = np.zeros(_JSON + 1) # setupv; read by vcirc 22qqqcod = np.zeros(_JSON + 1) # setupv; read nowhere