newrat_mod.f

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

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

Variables on this page (4)

json, rbe9, rli6, rli7

Purpose

NEWRAT_MOD holds the proton-capture rates of ⁶Li, ⁷Li and ⁹Be at each shell for the current model, shared by 5 routines. It replaces COMMON/NEWRAT/ (MHP via Claude, 10/14/26); JSON comes from DIMENSIONS_MOD. The rates at the start of the step are in OLDRAT_MOD.

What its variables hold, who sets them, who reads them. All are arrays of JSON.

RLI6, RLI7, RBE9: the rates (the Caughlan–Fowler fit times $\rho X$, scaled by the user factors SLI6, SLI7, SBE91–SBE93), set by LIRATE88 when called with J = 1, and zero above the burning region. LIBURN and LIBURN2 read them to burn the light elements, and adjust the value at the base of the convection zone by interpolation (0.5*FX toward the shell below).
GETW and BURSMIX copy them into the start-of-step rates of OLDRAT_MOD for the next step.

Design

A blanket USE NEWRAT_MOD.

Definitions

moduleNEWRAT_MOD

lines 19–32 · GitHub

Uses modules DIMENSIONS_MOD

Used by (5 units)

BURSMIX, GETW, LIBURN, LIBURN2, LIRATE88

Change log

The changes since rp_1

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

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C NEWRAT_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/NEWRAT/ 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 NEWRAT_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 :: RLI6(JSON) 29 REAL*8 :: RLI7(JSON) 30 REAL*8 :: RBE9(JSON) 31 32 END MODULE NEWRAT_MOD

1"""Python version of newrat_mod.f (YREC 7be4133): the light-element burning 2rates of the current model, formerly COMMON/NEWRAT/. 3 4Translates: newrat_mod.f#NEWRAT_MOD (7be4133) 5 6Routines that USE NEWRAT_MOD import this module and write 7``newrat_mod.<name>``. Set by lirate88 (j = 1). 8""" 9 10import numpy as np 11 12from pyrec import dimensions_mod 13 14_JSON = dimensions_mod.json 15 16# Fortran lines 28–30. 17rli6 = np.zeros(_JSON + 1) # Li6 + p rate 18rli7 = np.zeros(_JSON + 1) # Li7 + p rate 19rbe9 = np.zeros(_JSON + 1) # Be9 + p rate