zlaol_mod.f

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

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

Variables on this page (5)

numrhoz, numtz, zolaol, zorho, zot

Purpose

ZLAOL_MOD holds the LAOL89 opacity table for pure metals ("pure Z," "DBG 12/95"), as read, used at high $Z$ and $T$ when LPUREZ is set. Shared by 3 routines. It replaces COMMON/ZLAOL/ (MHP via Claude, 10/14/26). The dimensions (104 densities, 52 temperatures) are literal, as in the block.

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

ZOLAOL(104,52): the opacity $\kappa$ at each density and temperature (0 where the table has no entry); ZORHO(104), ZOT(52): the density and temperature grids; NUMRHOZ, NUMTZ: their sizes. All read from the file by RDZLAOL.
ZSULAOL reads them to build the splines of ZSLAOL_MOD, and replaces ZOT by $\log_{10}$ ZOT in place; GTPURZ then searches ZOT (as $\log T$) with NUMTZ.

Design

A blanket USE ZLAOL_MOD in each routine. Since ZSULAOL takes the log of ZOT in place, it must run exactly once after each RDZLAOL.

Definitions

moduleZLAOL_MOD

lines 19–30 · GitHub

Used by (3 units)

GTPURZ, RDZLAOL, ZSULAOL

Change log

The changes since rp_1

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

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C ZLAOL_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/ZLAOL/ into this module. 8C 9C Imported with a blanket USE, the convention for the conversion. 10C All 3 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 ZLAOL_MOD 20 21 IMPLICIT NONE 22 SAVE 23 24 REAL*8 :: ZOLAOL(104,52) 25 REAL*8 :: ZOT(52) 26 REAL*8 :: ZORHO(104) 27 INTEGER*4 :: NUMRHOZ 28 INTEGER*4 :: NUMTZ 29 30 END MODULE ZLAOL_MOD

1"""Python version of zlaol_mod.f (YREC 7be4133): the LAOL89 opacity table for pure metals, 2as read, formerly COMMON/ZLAOL/. 3 4Translates: zlaol_mod.f#ZLAOL_MOD (7be4133) 5 6Module-level names, as in the Fortran module, in lower case. Routines that 7USE ZLAOL_MOD import this module and write ``zlaol_mod.<name>``. 8Read from the file by rdzlaol. Index 0 unused in every dimension. 9""" 10 11import numpy as np 12 13# The table and its grids (Fortran lines 24–28). 14zolaol = np.zeros((105, 53)) # opacity at (density, temperature); 0 = no entry 15zot = np.zeros(53) # temperature grid (zsulaol takes its log10 in place) 16zorho = np.zeros(105) # density grid 17numrhoz = 0 # number of densities 18numtz = 0 # number of temperatures