deuter_mod.f

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

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

Variables on this page (5)

drate, drate0, fmassacc, jcz, json

Purpose

DEUTER_MOD holds the deuterium burning rates, the mass accreted in the step and the base of the surface convection zone, shared by 11 routines. It replaces COMMON/DEUTER/ (MHP via Claude, 10/14/26); JSON comes from DIMENSIONS_MOD.

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

DRATE(JSON): the deuterium burning rate at each shell, set by DEUTRATE (zeroed in MIX); read by DBURN and MIDMOD.
DRATE0(JSON): the rate at the start of the step, copied from DRATE by DEUTRATE on the first iteration (ITLVL = 1) and remapped by HPOINT when points are added; read by DBURN and MIDMOD.
FMASSACC: the mass accreted in the step (DELM), set by MDOT; read by DBURN, DBURNM and LIBURN, which mix the accreted material's deuterium and lithium into the convection zone.
JCZ: the shell at the base of the surface convection zone (JENV), set by the main program and MASSLOSS; read by DEUTRATE, which limits burning in the zone to the convective turnover time scale, and ENGEB.

Design

A blanket USE DEUTER_MOD.

Definitions

moduleDEUTER_MOD

lines 19–33 · GitHub

Uses modules DIMENSIONS_MOD

Used by (11 units)

(file level), DBURN, DBURNM, DEUTRATE, ENGEB, HPOINT, LIBURN, MASSLOSS, MDOT, MIDMOD, MIX

Change log

The changes since rp_1

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

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C DEUTER_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/DEUTER/ into this module. 8C 9C Imported with a blanket USE, the convention for the conversion. 10C All 11 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 DEUTER_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 :: DRATE(JSON) 29 REAL*8 :: DRATE0(JSON) 30 REAL*8 :: FMASSACC 31 INTEGER*4 :: JCZ 32 33 END MODULE DEUTER_MOD

1"""Python version of deuter_mod.f (YREC 7be4133): the deuterium burning 2rates, the accreted mass and the base of the convection zone, formerly 3COMMON/DEUTER/. 4 5Translates: deuter_mod.f#DEUTER_MOD (7be4133) 6 7Module-level names, as in the Fortran module, in lower case. Routines that 8USE DEUTER_MOD import this module and write ``deuter_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. FMASSACC has no initial value in the Fortran (0 18# here); MDOT sets it only when mass is gained (see dburn's page). 19drate = np.zeros(_JSON + 1) # deuterium burning rate per shell (deutrate) 20drate0 = np.zeros(_JSON + 1) # the same at the start of the step 21fmassacc = 0.0 # mass accreted in the step, g (mdot) 22jcz = 0 # shell at the base of the surface CZ (main, massloss)