cross_mod.f

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

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

Variables on this page (4)

fc12, fo16, qqs0ee, qs0e

Purpose

CROSS_MOD holds the factors by which the nuclear reaction rates are adjusted for user-changed cross sections: the S-factor derivative factors and the N15+p branching factors, shared by 3 routines. It replaces COMMON/CROSS/ (MHP via Claude, 10/14/26), less SSTANDARD (now in NMLPAR_MOD) and LNEWNUC (now in NMLBYPASS_MOD).

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

QS0E(8), QQS0EE(8): factors on the first-derivative ($Q_2$, $Q_3$) and second-derivative ($Q_4$, $Q_5$) terms of the eight pp and CNO rates, proportional to $S'/S$ and $S''/S$ relative to the reference values (REMAP's comments). REMAP sets them to 1, or, with LNEWNUC, from the user's S-factors and their derivatives.
FC12, FO16: the factors on the N15+p → C12+α and → O16+γ branches, 1 or, with LNEWNUC, the ratio of the user's S-factor to the reference; set in REMAP.
All four are read by ENGEB and RATES, which multiply the rate brackets and the branching ratio by them.

Design

A blanket USE CROSS_MOD in each of the three routines.

Definitions

moduleCROSS_MOD

lines 19–29 · GitHub

Used by (3 units)

ENGEB, RATES, REMAP

Change log

The changes since rp_1

Changed since rp_1. New file: COMMON/CROSS/ as a module, with SSTANDARD and LNEWNUC moved to the namelist modules.

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C CROSS_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/CROSS/ 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 CROSS_MOD 20 21 IMPLICIT NONE 22 SAVE 23 24 REAL*8 :: QS0E(8) 25 REAL*8 :: QQS0EE(8) 26 REAL*8 :: FO16 27 REAL*8 :: FC12 28 29 END MODULE CROSS_MOD

1"""Python version of cross_mod.f (YREC 7be4133): the cross-section 2adjustment factors of the nuclear rates, formerly COMMON/CROSS/. 3 4Translates: cross_mod.f#CROSS_MOD (7be4133) 5 6Module-level names, as in the Fortran module, in lower case (REAL*8). 7Routines that USE CROSS_MOD import this module and write 8``cross_mod.<name>``. 9""" 10 11import numpy as np 12 13# Fortran lines 24–27. 14qs0e = np.zeros(9) # S'/S factors of the 8 pp/CNO rates (remap), index 0 unused 15qqs0ee = np.zeros(9) # S''/S factors (remap), index 0 unused 16fo16 = 0.0 # factor on N15+p -> O16+gamma (remap) 17fc12 = 0.0 # factor on N15+p -> C12+alpha (remap)