surfopac.f

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

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

Variables on this page (1)

lopal95

Purpose

SURFOPAC sets up the opacity tables for the surface hydrogen abundance $X$ (XENV): if the OPAL 1995 tables are switched on (LOPAL95), it calls the routine that interpolates them in $X$. Its comment says the tables must already have been read. The low-temperature tables of Ferguson et al. (2005) (LALEX06), now the only molecular tables, need no set-up here: GETALEX06 rebuilds their table itself, through ALEX06TAB, whenever the $X$ or $Z$ of a point differs from that of the current table by more than $10^{-8}$ (settled from the code, Nov 9, 2026: getalex06.f, lines 26–33 ↗).

How it fits

Who calls it: STARIN, when the starting model is ready, and HPOINT, when the surface $X$ (HCOMP(1,M)) has changed by more than $10^{-8}$.
What it calls: OP95XTAB (OPAL95, if LOPAL95).
Modules: NMLPAR_MOD, for the namelist switch LOPAL95, read by PARMIN.

Definitions

subroutineSURFOPAC

lines 6–44 · GitHub

Interpolates the OPAL 1995 tables to the surface $X$, if they are switched on.

Uses modules NMLPAR_MOD

Called by HPOINT, STARIN

Flowchart of its calls, in order (1)

if (LOPAL95), lines 39–41

OP95XTAB 40

Change log

The changes since rp_1

Changed since rp_1. The OPAL 1992 and Alexander 1995 tables were removed from the code, and with them the calls to LL4TH and ALX8TH; COMMON NEWOPAC gave way to NMLPAR_MOD; IMPLICIT NONE. Question 1 of rp_1 (why the Alexander 2006, LAOL 1989 and Kurucz 1990 tables need no surface-$X$ set-up) is settled in part: LAOL 1989 and Kurucz 1990 are gone, so only Alexander 2006 remains, and the question is narrowed to it.

1C 2C 3CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 4C SURFOPAC 5CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 6 SUBROUTINE SURFOPAC(XENV) 7C 10/14/26 MHP via Claude 8C code changed to convert COMMON/NMLPAR/ into NMLPAR_MOD. 9 USE NMLPAR_MOD 10C 10/20/26 MHP via Claude 11C code changed to replace implicit typing with IMPLICIT NONE and 12C explicit declarations. 13 IMPLICIT NONE 14 REAL*8 XENV 15C 10/13/26 MHP via Claude: removed the unused COMMON/IOUNITS/ declaration 16C (no member referenced here; residue of the obsolete-option removals). 17C 10/12/26 MHP via Claude 18C code changed to extend COMMON/NMLPAR/ to 234 members, adding the 19C namelist variables that the physics routines read from 27 COMMON 20C blocks which also held derived or working members. Those blocks keep 21C only their non-namelist members; five of them -- NEWMX, VMULT, 22C PMMWIND, CCOUT1 and ENVGEN -- had nothing left and are gone. 23C Derivation-only namelist variables from the same blocks went to 24C COMMON/NMLBYPASS/ instead, not here. 25C 26C Declared identically in all 121 files that hold any of these members. 27C Generated from one member list, so every copy is byte-identical; the 28C linker's 'size of symbol changed' warning is the check that they are. 29C 10/13/26 MHP via Claude: removed the unused COMMON/NEWOPAC/ declaration 30C (no member referenced here; residue of the obsolete-option removals). 31 SAVE 32 33C THIS SUBROUTINE SETS UP SURF OPACITY TABLES 34C ASSUMES TABLES HAVE ALREADY BEEN READ IN 35 36C INTERIOR TABLES 37 38C SETUP OPAL95 TABLES 39 IF (LOPAL95) THEN 40 CALL OP95XTAB(XENV) 41 END IF 42 43 RETURN 44 END

1"""Python version of surfopac.f (YREC development, 7be4133).""" 2 3from pyrec import nmlpar_mod 4from pyrec.op95xtab import op95xtab 5 6 7def surfopac(xenv: float) -> None: 8 """Set up the opacity tables for the surface hydrogen abundance. 9 10 Assumes the tables have already been read. Only the OPAL 1995 11 tables are interpolated in X here. 12 13 Translates: surfopac.f#SURFOPAC (7be4133) 14 15 Parameters 16 ---------- 17 xenv : float 18 Surface (envelope) hydrogen mass fraction X. 19 20 """ 21 # Interior tables: set up the OPAL95 tables (Fortran lines 36–41). 22 if nmlpar_mod.lopal95: 23 op95xtab(xenv)