zsulaol.f

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

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

Variables on this page (9)

numrhoz, numrsz, numtz, zolaol, zorho, zot, zsdorl, zslaoll, zsrhol

Purpose

ZSULAOL ("spline pure Z LAOL89") prepares the pure-metal LAOL89 opacity table (DBG 12/95) for interpolation in density, as SULAOL did for the main LAOL89 tables (sulaol.f, removed since rp_1 with those tables): for each temperature, a natural cubic spline of $\log\kappa$ against $\log\rho$, through the nonzero entries, if there are at least four.

How it fits

Inputs: ZLAOL_MOD (formerly COMMON/ZLAOL/: ZOLAOL, the grids ZOT, ZORHO, sizes NUMRHOZ, NUMTZ), filled by RDZLAOL when LPUREZ is set.
Who calls it: SETUPOPAC.
What it calls: CSPLINE (end slopes $10^{30}$, a natural spline).
Outputs: ZSLAOL_MOD (formerly COMMON/ZSLAOL/: $\log\kappa$, $\log\rho$, second derivatives, points per temperature), read by GTPURZ; a temperature with fewer than four nonzero entries gets NUMRSZ = 0, and GTPURZ stops the run if it needs that temperature (gtpurz.f, line 75 ↗, gtpurz.f, lines 106–111 ↗); the table of the test suite, PURECO.DBGLAOL, has at least 23 at every temperature. ZOT is converted to $\log_{10} T$ in place (lines 23–25 ↗); this is safe because the routine runs once per job (main program → SETUPS → SETUPOPAC), and its only call follows a fresh read of the table by RDZLAOL, so even a second call would not take the logarithm twice (settled from the code, Nov 9, 2026: setupopac.f, lines 79–81 ↗).

Definitions

subroutineZSULAOL

lines 6–51 · GitHub

Natural cubic splines of $\log\kappa$ against $\log\rho$ for each temperature of the pure-Z LAOL89 table.

Uses modules ZLAOL_MOD, ZSLAOL_MOD

Called by SETUPOPAC

Flowchart of its calls, in order (1)

for IT=1, NUMTZ, lines 26–49

if (NUMS .GE. 4), lines 38–45

CSPLINE 40

Compiler warnings (1)

gfortran -Wall -Wextra, compiling only; nothing was run.

line 32: Inequality comparison for REAL(8) [-Wcompare-reals]

Change log

The changes since rp_1

Changed since rp_1. Declarations only: COMMON/ZLAOL/ and /ZSLAOL/ to ZLAOL_MOD and ZSLAOL_MOD, IMPLICIT NONE; the main LAOL89 tables (SULAOL) are removed, the pure-Z one kept.

1C 2C 3CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 4C SPLINE PURE Z LAOL89 5CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 6 SUBROUTINE ZSULAOL 7C 10/14/26 MHP via Claude 8C code changed to convert COMMON/ZSLAOL/ into ZSLAOL_MOD. 9 USE ZSLAOL_MOD 10C 10/14/26 MHP via Claude 11C code changed to convert COMMON/ZLAOL/ into ZLAOL_MOD. 12 USE ZLAOL_MOD 13C 10/20/26 MHP via Claude 14C code changed to replace implicit typing with IMPLICIT NONE and 15C explicit declarations. 16 IMPLICIT NONE 17 REAL*8 TDOR,TOP,TR 18 INTEGER IR,IT,NUMS 19C DBG 4/94 Modified to do ZRAMP stuff. 20 DIMENSION TOP(104),TR(104),TDOR(104) 21C DBG 12/95 ARRAYS FOR PURE Z TABLE 22 SAVE 23 DO IT=1, NUMTZ 24 ZOT(IT) = LOG10(ZOT(IT)) 25 END DO 26 DO IT=1, NUMTZ 27 NUMS=0 28 DO IR=1, NUMRHOZ 29 ZSLAOLL(IR,IT) = 0.0D0 30 ZSRHOL(IR,IT) = 0.0D0 31 ZSDORL(IR,IT) = 0.0D0 32 IF(ZOLAOL(IR,IT) .NE. 0.0D0) THEN 33 NUMS = NUMS+1 34 TOP(NUMS) = LOG10(ZOLAOL(IR,IT)) 35 TR(NUMS) = LOG10(ZORHO(IR)) 36 END IF 37 END DO 38 IF (NUMS .GE. 4) THEN 39 NUMRSZ(IT)=NUMS 40 CALL CSPLINE(TR,TOP,NUMS,1.0D30,1.0D30,TDOR) 41 DO IR=1,NUMS 42 ZSLAOLL(IR,IT) = TOP(IR) 43 ZSRHOL(IR,IT) = TR(IR) 44 ZSDORL(IR,IT) = TDOR(IR) 45 END DO 46 ELSE 47 NUMRSZ(IT) = 0 48 END IF 49 END DO 50 RETURN 51 END

1"""Python version of zsulaol.f (YREC 7be4133): splines of the pure-Z LAOL89 2opacity table in density.""" 3 4import math 5 6import numpy as np 7 8from pyrec import zlaol_mod 9from pyrec import zslaol_mod 10from pyrec.cspline import cspline 11 12# The Fortran's local work arrays (SAVE). 13top = np.zeros(105) 14tr = np.zeros(105) 15tdor = np.zeros(105) 16 17 18def zsulaol() -> None: 19 """Convert the pure-Z table's temperatures to log10 T in place, and 20 for each temperature fit a natural cubic spline of log kappa against 21 log rho through the nonzero entries, if there are at least four 22 (otherwise numrsz is 0 and gtpurz skips that temperature). 23 24 Translates: zsulaol.f#ZSULAOL (7be4133) 25 26 No arguments; reads zlaol_mod (from rdzlaol), sets zslaol_mod. None is 27 returned. 28 29 """ 30 31 z = zlaol_mod 32 zs = zslaol_mod 33 34 # The temperatures to log10 T, in place (Fortran lines 23–25). 35 for it in range(1, z.numtz + 1): 36 z.zot[it] = math.log10(z.zot[it]) 37 38 # For each temperature, a spline in density (Fortran lines 26–49). 39 for it in range(1, z.numtz + 1): 40 # The nonzero entries, in logs (Fortran lines 27–37). 41 nums = 0 42 for ir in range(1, z.numrhoz + 1): 43 zs.zslaoll[ir, it] = 0.0 44 zs.zsrhol[ir, it] = 0.0 45 zs.zsdorl[ir, it] = 0.0 46 if z.zolaol[ir, it] != 0.0: 47 nums = nums + 1 48 top[nums] = math.log10(z.zolaol[ir, it]) 49 tr[nums] = math.log10(z.zorho[ir]) 50 51 # A natural spline through them, if there are at least four (Fortran lines 38–48). 52 if nums >= 4: 53 zs.numrsz[it] = nums 54 cspline(tr, top, nums, 1.0e30, 1.0e30, tdor) 55 for ir in range(1, nums + 1): 56 zs.zslaoll[ir, it] = top[ir] 57 zs.zsrhol[ir, it] = tr[ir] 58 zs.zsdorl[ir, it] = tdor[ir] 59 else: 60 zs.numrsz[it] = 0