prevcz_mod.f

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

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

Variables on this page (6)

bbe9old, bli6old, bli7old, jenvos, lbliold, rbase0

Purpose

PREVCZ_MOD carries the state of the surface convection zone and of light-element burning from one step to the next, for LIBURN and LIBURN2. Shared by 3 routines. It replaces COMMON/PREVCZ/ (MHP via Claude, 10/14/26).

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

RBASE0: the radius of the base of the convection zone at the start of the step (with overshoot). LIBURN, LIBURN2 compute it from the starting model when it is 0, and store the step's ending base in it for the next step; the main program resets it to 0 when LEXCOM is set.
JENVOS: the shell of the envelope's base at the end of the step, stored by LIBURN, LIBURN2 and used as the next step's starting JENV0 when RBASE0 is already known.
BLI6OLD, BLI7OLD, BBE9OLD: the mass-averaged burning rates of ⁶Li, ⁷Li and ⁹Be in the convection zone at the end of the step (their logarithms, by the code's comments), stored by LIBURN for the start of the next.
LBLIOLD (.FALSE. at the start): whether those three hold LIBURN's end-of-step rates. LIBURN sets it when it stores them, clears it when it skips the convection zone, and reuses the stored rates as the next step's starting rates only when it is set (otherwise it recomputes them); LIBURN2 always computes its starting rates and clears it, and the main program clears it at each kind card when LEXCOM is set.

Design

A blanket USE PREVCZ_MOD in each routine.

Definitions

modulePREVCZ_MOD

lines 19–35 · GitHub

Used by (3 units)

(file level), LIBURN, LIBURN2

Change log

The changes since rp_1

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

Changed since a8a5950. The new flag LBLIOLD says when BLI6OLD, BLI7OLD, and BBE9OLD hold LIBURN's end-of-step rates. The old test, BLI6OLD ≤ 0, was always true for the logarithm of a small rate, so the stored rates were never reused; the main program now clears the flag instead of zeroing the three rates.

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C PREVCZ_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/PREVCZ/ 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 PREVCZ_MOD 20 21 IMPLICIT NONE 22 SAVE 23 24 REAL*8 :: RBASE0 25 REAL*8 :: BLI6OLD 26 REAL*8 :: BLI7OLD 27 REAL*8 :: BBE9OLD 28C 11/02/26 MHP via Claude 29C code changed to add LBLIOLD, set when BLI6OLD, BLI7OLD and BBE9OLD hold 30C LIBURN's end-of-step surface CZ rates. The old test BLI6OLD.LE.0 was 31C always true for these logs of small rates, so they were never reused. 32 LOGICAL*4 :: LBLIOLD = .FALSE. 33 INTEGER*4 :: JENVOS 34 35 END MODULE PREVCZ_MOD

1"""Python version of prevcz_mod.f (YREC 7be4133): the convection zone and 2light-element burning rates carried from one step to the next, formerly 3COMMON/PREVCZ/. 4 5Translates: prevcz_mod.f#PREVCZ_MOD (7be4133) 6 7Module-level names, as in the Fortran module, in lower case. Routines that 8USE PREVCZ_MOD import this module and write ``prevcz_mod.<name>``. 9Set and read by liburn and liburn2; main clears lbliold when lexcom is set. 10""" 11 12# Fortran lines 24–33. 13rbase0 = 0.0 # radius of the convection zone's base at the step's start 14bli6old = 0.0 # CZ-averaged 6Li, 7Li, 9Be burning rates at the step's end 15bli7old = 0.0 16bbe9old = 0.0 17lbliold = False # True when the three hold liburn's end-of-step rates 18jenvos = 0 # shell of the envelope's base at the step's end