pulse1_mod.f

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

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

Variables on this page (11)

json, lpumod, pemu, pqcp, pqdp, pqdt, pqed, pqet, pqod, pqot, prmu

Purpose

PULSE1_MOD holds the thermodynamic and opacity derivatives at every shell of the interior, stored for the pulsation output, and the switch that makes the envelope routines keep their derivatives. Shared by 6 routines. It replaces COMMON/PULSE1/ (MHP via Claude, 10/14/26); JSON comes from DIMENSIONS_MOD.

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

PQDP, PQDT, PQED, PQET, PQOD, PQOT, PQCP, PRMU, PEMU (JSON each): at each shell, the equation of state's $\partial\ln\rho/\partial\ln P$ and $\partial\ln\rho/\partial\ln T$ (QDP, QDT), the energy generation's and the opacity's logarithmic derivatives in $\rho$ and $T$ (QED, QET, QOD, QOT), $c_P$ (QCP), and RMU, EMU. Stored by COEFFT when LPULSE is set; PQDP and PQDT always, since GAM1 (SOUND_MOD) and the storage file need them ("JVS 01/11"). WRTMOD writes them in the pulsation model; STITCH writes PQDT to the stored model.
LPUMOD: set by ENVINT to LPULPT .AND. LPRT (pulsation output of the envelope); QENV and QATM then save their point's values in PULSE2_MOD.

Design

A blanket USE PULSE1_MOD in each routine.

Definitions

modulePULSE1_MOD

lines 19–39 · GitHub

Uses modules DIMENSIONS_MOD

Used by (6 units)

COEFFT, ENVINT, QATM, QENV, STITCH, WRTMOD

Change log

The changes since rp_1

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

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C PULSE1_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/PULSE1/ into this module. 8C 9C Imported with a blanket USE, the convention for the conversion. 10C All 6 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 PULSE1_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 :: PQDP(JSON) 29 REAL*8 :: PQED(JSON) 30 REAL*8 :: PQET(JSON) 31 REAL*8 :: PQOD(JSON) 32 REAL*8 :: PQOT(JSON) 33 REAL*8 :: PQCP(JSON) 34 REAL*8 :: PRMU(JSON) 35 REAL*8 :: PQDT(JSON) 36 REAL*8 :: PEMU(JSON) 37 LOGICAL*4 :: LPUMOD 38 39 END MODULE PULSE1_MOD

1"""Python version of pulse1_mod.f (YREC 7be4133): derivatives at every shell 2for the pulsation output, formerly COMMON/PULSE1/. 3 4Translates: pulse1_mod.f#PULSE1_MOD (7be4133) 5 6Module-level names, as in the Fortran module, in lower case. Routines that 7USE PULSE1_MOD import this module and write ``pulse1_mod.<name>``. 8""" 9 10import numpy as np 11 12from pyrec import dimensions_mod 13 14_JSON = dimensions_mod.json 15 16# At each shell, index 0 unused (coefft, when lpulse; pqdp and pqdt always) 17# (Fortran lines 28–36). 18pqdp = np.zeros(_JSON + 1) # dln rho / dln P 19pqed = np.zeros(_JSON + 1) # energy generation's derivative in rho 20pqet = np.zeros(_JSON + 1) # ... in T 21pqod = np.zeros(_JSON + 1) # opacity's derivative in rho 22pqot = np.zeros(_JSON + 1) # ... in T 23pqcp = np.zeros(_JSON + 1) # c_P 24prmu = np.zeros(_JSON + 1) 25pqdt = np.zeros(_JSON + 1) # dln rho / dln T 26pemu = np.zeros(_JSON + 1) 27 28# The envelope switch (Fortran line 37). 29lpumod = False # envelope pulsation output on (envint)