debhu_mod.f

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

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

Variables on this page (6)

cdh, dhnue, xxdh, yydh, zdh, zzdh

Purpose

DEBHU_MOD holds the inputs of the Debye–Hückel correction to the equation of state: the constant CDH, the charges of the 18 heavy elements and the local composition, shared by 11 routines. It replaces COMMON/DEBHU/ (MHP via Claude, 10/14/26), less the switch LDH and ETADH0, ETADH1 (now in NMLPAR_MOD).

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

CDH: $-\sqrt{\pi/(k\,m_H^3)}\,e^3/3$, set in SETUPS.
DHNUE(18): the charges of C, N, O, Ne, Na, Mg, Al, Si, P, S, Cl, Ar, Ca, Ti, Cr, Mn, Fe, Ni (SETUPS's comment), set in SETUPS.
XXDH, YYDH, ZZDH, ZDH(18): the local $X$, $Y$ (He3 + He4), $Z$ and the heavy-element abundances (C12 + C13, N14 + N15, O16 + O18 from HCOMP), copied from the shell's composition just before the equation of state is called, by nine routines: COEFFT, CRRECT, GETNEWENV, HSUBP, PHYSIC, SCONVEC, STARIN, WRTMOD and WRTOUT.
All are read only by EQRELV, which, with LDH, forms $\sum_i$ DHNUE(i) ZDH(i) and the Debye–Hückel pressure $P_\mathrm{DH} = $ CDH × (charge factor) × $\sqrt{\rho/T}\,\rho$.

Design

A blanket USE DEBHU_MOD. The composition reaches EQRELV through the module rather than as arguments, so each caller of the equation of state must copy it first; the copies are the same few lines in nine files. At rp_1, two files (EQSTAT, EQSTAT2) named the second member XXDY, harmless with COMMON's storage association.

Definitions

moduleDEBHU_MOD

lines 19–31 · GitHub

Used by (11 units)

COEFFT, CRRECT, EQRELV, GETNEWENV, HSUBP, PHYSIC, SCONVEC, SETUPS, STARIN, WRTMOD, WRTOUT

Change log

The changes since rp_1

Changed since rp_1. New file: COMMON/DEBHU/ as a module, with LDH, ETADH0, ETADH1 moved to NMLPAR_MOD.

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C DEBHU_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/DEBHU/ 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 DEBHU_MOD 20 21 IMPLICIT NONE 22 SAVE 23 24 REAL*8 :: CDH 25 REAL*8 :: ZDH(18) 26 REAL*8 :: XXDH 27 REAL*8 :: YYDH 28 REAL*8 :: ZZDH 29 REAL*8 :: DHNUE(18) 30 31 END MODULE DEBHU_MOD

1"""Python version of debhu_mod.f (YREC 7be4133): the inputs of the 2Debye-Hueckel correction, formerly COMMON/DEBHU/. 3 4Translates: debhu_mod.f#DEBHU_MOD (7be4133) 5 6Module-level names, as in the Fortran module, in lower case (REAL*8). 7Routines that USE DEBHU_MOD import this module and write 8``debhu_mod.<name>``. 9""" 10 11import numpy as np 12 13# Fortran lines 24–29. 14cdh = 0.0 # Debye-Hueckel constant (setups) 15zdh = np.zeros(19) # abundances of the 18 heavy elements, index 0 unused 16xxdh = 0.0 # local X (copied by the callers of the EOS) 17yydh = 0.0 # local Y 18zzdh = 0.0 # local Z 19dhnue = np.zeros(19) # charges of the 18 heavy elements (setups), index 0 unused