zslaol_mod.f

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

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

Variables on this page (4)

numrsz, zsdorl, zslaoll, zsrhol

Purpose

ZSLAOL_MOD holds the density splines of the pure-metal LAOL89 opacity table (ZLAOL_MOD), one natural cubic spline of $\log\kappa$ against $\log\rho$ per temperature, through the table's nonzero entries. Shared by 2 routines. It replaces COMMON/ZSLAOL/ (MHP via Claude, 10/14/26).

What its variables hold, who sets them, who reads them. All are set by ZSULAOL and read by GTPURZ.

NUMRSZ(52): the number of nonzero entries at each temperature, or 0 when there are fewer than four (no spline).
ZSLAOLL(104,52), ZSRHOL(104,52): $\log\kappa$ and $\log\rho$ at those entries; ZSDORL(104,52): the spline's second derivatives, from CSPLINE.

Design

A blanket USE ZSLAOL_MOD in each routine.

Definitions

moduleZSLAOL_MOD

lines 19–29 · GitHub

Used by (2 units)

GTPURZ, ZSULAOL

Change log

The changes since rp_1

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

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C ZSLAOL_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/ZSLAOL/ 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 ZSLAOL_MOD 20 21 IMPLICIT NONE 22 SAVE 23 24 REAL*8 :: ZSLAOLL(104,52) 25 REAL*8 :: ZSRHOL(104,52) 26 REAL*8 :: ZSDORL(104,52) 27 INTEGER*4 :: NUMRSZ(52) 28 29 END MODULE ZSLAOL_MOD

1"""Python version of zslaol_mod.f (YREC 7be4133): the density splines of the pure-metal LAOL89 2opacity table, formerly COMMON/ZSLAOL/. 3 4Translates: zslaol_mod.f#ZSLAOL_MOD (7be4133) 5 6Module-level names, as in the Fortran module, in lower case. Routines that 7USE ZSLAOL_MOD import this module and write ``zslaol_mod.<name>``. 8Set by zsulaol, read by gtpurz. Index 0 unused in every dimension. 9""" 10 11import numpy as np 12 13# The splines, one per temperature (Fortran lines 24–27). 14zslaoll = np.zeros((105, 53)) # log kappa at the nonzero entries 15zsrhol = np.zeros((105, 53)) # log rho there 16zsdorl = np.zeros((105, 53)) # the spline's second derivatives (cspline) 17numrsz = np.zeros(53, dtype=np.int32) # entries per temperature (0: no spline)