setupopac.f

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

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

Variables on this page (6)

l2z, lalex06, ldifz, lopal95, lpurez, lzramp

Purpose

Reads the opacity tables chosen for the run and sets up their interpolation. Its comment: "THIS SUBROUTINE READS IN SPECIFIED OPACITY TABLES AND SET UP SPLINES FOR THE TABLES. WHEN LZRAMP=T OR LDIFZ=T THEN READ IN SECOND SET OF OPACITY TABLES AT DIFFERENT Z."

How it fits

Who calls it: SETUPS, once per run; the file names are passed through (MHP 8/25).
What it calls, by the namelist switches in NMLPAR_MOD:
interior tables: OPAL95, LL95TBL then OP95XTAB (LOPAL95); the pure-Z table, RDZLAOL then ZSULAOL (LPUREZ);
low-temperature tables: Ferguson et al. 2005, READALEX06 (LALEX06).
Inputs: XENV (the envelope's $X$, for the tables built at fixed composition) and the file names FALEX06, FLIV95, FPUREZ. Sets L2Z = LZRAMP or LDIFZ (both namelist parameters) in NEWOPAC_MOD (line 64 ↗); despite the header's comment, no reader reads it, and only GETOPAC does, in its single-table test (getopac.f, lines 140–141 ↗; settled from the code, Nov 9, 2026).
Output: the tables, in the modules of the routines called. They are used at run time by GETOPAC.

From YREC's documentation

Control Options (Input files): FALEX06 is the molecular opacity table (Ferguson et al. 2005, with LALEX06, recommended), FLIV95 the atomic one (OP/OPAL format, with LOPAL95, recommended), and FPUREZ the pure Y/Z table (with LPUREZ).

Definitions

subroutineSETUPOPAC

lines 14–95 · GitHub

Uses modules NEWOPAC_MOD, NMLPAR_MOD

Called by SETUPS

Flowchart of its calls, in order (5)

if (LOPAL95), lines 69–72

LL95TBL 70→OP95XTAB 71

if (LPUREZ), lines 79–82

RDZLAOL 80→ZSULAOL 81

if (LALEX06), lines 87–93

READALEX06 88

Change log

The changes since rp_1

Changed since rp_1. The OPAL92 (SETLLO, LL4TH), LAOL89 (RDLAOL, SULAOL), Alexander 1995 (ALXTBL, ALX8TH) and Kurucz 1990 (SETKRZ) tables are removed (all removed since rp_1), with their file-name arguments and the argument V; this settles former question 1 (several interior tables at once), since OPAL95 is now the only interior table; COMMON blocks to modules, IMPLICIT NONE.

1C 2C 3CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 4C SETUPOPAC 5CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 6C 9/23/26 MHP via Claude 7C code changed to remove the OPAL92 opacity tables 8C 9/23/26 MHP via Claude 9C code changed to remove the Kurucz 1990 opacity tables 10C 9/23/26 MHP via Claude 11C code changed to remove the LAOL89 opacity tables 12C 9/24/26 MHP via Claude 13C code changed to remove the Alexander 1995 opacity tables 14 SUBROUTINE SETUPOPAC(XENV, FALEX06, 15 * *FLIV95,FPUREZ) 16C 10/14/26 MHP via Claude 17C code changed to convert COMMON/NEWOPAC/ into NEWOPAC_MOD. 18 USE NEWOPAC_MOD 19C 10/14/26 MHP via Claude 20C code changed to convert COMMON/NMLPAR/ into NMLPAR_MOD. 21 USE NMLPAR_MOD 22C 10/20/26 MHP via Claude 23C code changed to replace implicit typing with IMPLICIT NONE and 24C explicit declarations. 25 IMPLICIT NONE 26 REAL*8 XENV 27 28C MHP 8/25 Added pass-through file names 29 CHARACTER*256 FALEX06,FLIV95, 30 * * FPUREZ 31C 10/13/26 MHP via Claude: removed the unused COMMON/IOUNITS/ declaration 32C (no member referenced here; residue of the obsolete-option removals). 33C 10/12/26 MHP via Claude 34C code changed to extend COMMON/NMLPAR/ to 234 members, adding the 35C namelist variables that the physics routines read from 27 COMMON 36C blocks which also held derived or working members. Those blocks keep 37C only their non-namelist members; five of them -- NEWMX, VMULT, 38C PMMWIND, CCOUT1 and ENVGEN -- had nothing left and are gone. 39C Derivation-only namelist variables from the same blocks went to 40C COMMON/NMLBYPASS/ instead, not here. 41C 42C Declared identically in all 121 files that hold any of these members. 43C Generated from one member list, so every copy is byte-identical; the 44C linker's 'size of symbol changed' warning is the check that they are. 45C MHP 8/25 Removed character file names from common block 46C 9/24/26 MHP via Claude 47C code changed to remove dead IOOPAL2 from COMMON/ZRAMP/ 48C 10/9/26 MHP via Claude 49C code changed to remove NK from COMMON/ZRAMP/: it was only ever a 50C redundant copy of the run-index loop variable main.f already passes 51C explicitly to the subroutines that need it (STARIN/RSCALE), and was 52C never read by name in this file. 53C 10/13/26 MHP via Claude 54C code changed to remove COMMON/ZRAMP/: its sole remaining member 55C IOLAOL2 is neither written nor read anywhere. The namelist members 56C that shared this block moved to the consolidated block in stage 2b. 57 SAVE 58 59C 60C THIS SUBROUTINE READS IN SPECIFIED OPACITY TABLES AND 61C SET UP SPLINES FOR THE TABLES. 62C WHEN LZRAMP=T OR LDIFZ=T THEN READ IN SECOND SET OF 63C OPACITY TABLES AT DIFFERENT Z (E.G. ZOPAL952). 64 L2Z = LZRAMP.OR.LDIFZ 65C 66C INTERIOR TABLES 67C 68C READ IN OPAL95 TABLES 69 IF (LOPAL95) THEN 70 CALL LL95TBL(FLIV95) 71 CALL OP95XTAB(XENV) 72 END IF 73 74C 9/23/26 MHP via Claude 75C code changed to remove the LAOL89 opacity tables 76C 77C READ IN PURE Z TABLE 78C 79 IF(LPUREZ)THEN 80 CALL RDZLAOL(FPUREZ) 81 CALL ZSULAOL 82 END IF 83C 84C LOW TEMP TABLES 85C 86C READ IN ALEX 2006 TABLES 87 IF (LALEX06) THEN 88 CALL READALEX06(FALEX06) 89C 9/23/26 MHP via Claude 90C code changed to remove the Kurucz 1990 opacity tables 91C 9/24/26 MHP via Claude 92C code changed to remove the Alexander 1995 opacity tables 93 END IF 94 RETURN 95 END

1"""Python version of setupopac.f (YREC 7be4133): read the opacity tables 2chosen for the run.""" 3 4from pyrec import newopac_mod 5from pyrec import nmlpar_mod 6from pyrec.ll95tbl import ll95tbl 7from pyrec.op95xtab import op95xtab 8from pyrec.rdzlaol import rdzlaol 9from pyrec.readalex06 import readalex06 10from pyrec.zsulaol import zsulaol 11 12 13def setupopac(xenv: float, falex06: str, fliv95: str, fpurez: str) -> None: 14 """Read the interior (OPAL95, and the pure-Z table) and low-temperature 15 (Ferguson et al. 2005) opacity tables the namelist switches on, and set 16 L2Z, which tells the readers to read a second table at another Z. 17 18 Translates: setupopac.f#SETUPOPAC (7be4133) 19 20 Arguments: xenv, the envelope's X; the file names falex06, fliv95, 21 fpurez. None is set. 22 23 """ 24 25 n = nmlpar_mod 26 27 # With LZRAMP or LDIFZ, read a second set of tables at another Z (Fortran line 64). 28 newopac_mod.l2z = n.lzramp or n.ldifz 29 30 # Interior tables: OPAL95 (Fortran lines 69–72). 31 if n.lopal95: 32 ll95tbl(fliv95) 33 op95xtab(xenv) 34 35 # The pure-Z table (Fortran lines 79–82). 36 if n.lpurez: 37 rdzlaol(fpurez) 38 zsulaol() 39 40 # Low-temperature tables: Alexander/Ferguson 2006 (Fortran lines 87–93). 41 if n.lalex06: 42 readalex06(falex06)