ccr_mod.f

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

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

Variables on this page (4)

dtable, dval, etat8, iden

Purpose

CCR_MOD holds the table of the electron (Fermi–Dirac) equation of state used for partial degeneracy. It replaces COMMON/CCR/ (MHP via Claude, 10/14/26).

What its variables hold, who sets them, who reads them. SETUPS reads from the Fermi-table file (unit IFERMI) the 43 abscissae in density DVAL, the degeneracies ETAT8 and the table DTABLE(5,43,20), and builds IDEN(261), which maps a density bin of width 0.05 to the table row; EQRELV reads them all.

Definitions

moduleCCR_MOD

lines 19–29 · GitHub

Used by (2 units)

EQRELV, SETUPS

Change log

The changes since rp_1

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

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C CCR_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/CCR/ into this module. 8C 9C Imported with a blanket USE, the convention for the conversion. 10C All 2 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 CCR_MOD 20 21 IMPLICIT NONE 22 SAVE 23 24 REAL*8 :: DVAL(43) 25 REAL*8 :: ETAT8(43) 26 REAL*8 :: DTABLE(5,43,20) 27 INTEGER*4 :: IDEN(261) 28 29 END MODULE CCR_MOD

1"""Python version of ccr_mod.f (YREC 7be4133): the Fermi-Dirac electron table, formerly COMMON/CCR/. 2 3Translates: ccr_mod.f#CCR_MOD (7be4133) 4 5Module-level names, as in the Fortran module, in lower case. Routines that 6USE CCR_MOD import this module and write ``ccr_mod.<name>``. 7""" 8 9import numpy as np 10 11# All set by setups (read from unit IFERMI; iden built there); read by eqrelv 12# (Fortran lines 24–27). 13dval = np.zeros(44) # density abscissae 14etat8 = np.zeros(44) # degeneracies 15dtable = np.zeros((6, 44, 21)) # dtable[k, i, j] for DTABLE(K, I, J) 16iden = np.zeros(262, dtype=int) # density bin (width 0.05) -> table row