splin_mod.f

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

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

Variables on this page (5)

json, xtab, xval, ytab, yval

Purpose

SPLIN_MOD holds the work arrays of the interpolations by OSPLIN (the shape-preserving spline) in the routines that move quantities between the model's grid and the even grid in $\chi$ of rotational mixing: the table (XTAB, YTAB) and the points of evaluation (XVAL) with the values found (YVAL). Shared by 4 routines. It replaces COMMON/SPLIN/ (MHP via Claude, 10/14/26); JSON comes from DIMENSIONS_MOD.

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

XTAB(JSON), YTAB(JSON): the table to interpolate; XVAL(JSON): the abscissas at which to evaluate; YVAL(JSON): the values returned by OSPLIN. Each user fills them, calls OSPLIN with them as arguments and reads the result before returning: EQ2MOD, MIXCOM, MIXGRID and ROTGRID. OSPLIN itself does not use the module (its arguments have the same names).

Design

A blanket USE SPLIN_MOD in each routine. The arrays are scratch space: no value is passed from one routine to another through them, so they could be local arrays of each user.

Definitions

moduleSPLIN_MOD

lines 19–33 · GitHub

Uses modules DIMENSIONS_MOD

Used by (4 units)

EQ2MOD, MIXCOM, MIXGRID, ROTGRID

Change log

The changes since rp_1

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

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C SPLIN_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/SPLIN/ into this module. 8C 9C Imported with a blanket USE, the convention for the conversion. 10C All 4 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 SPLIN_MOD 20C 10/20/26 MHP via Claude 21C code changed to take JSON from DIMENSIONS_MOD instead of declaring 22C a PRIVATE copy here. 23 USE DIMENSIONS_MOD 24 25 IMPLICIT NONE 26 SAVE 27 28 REAL*8 :: XVAL(JSON) 29 REAL*8 :: YVAL(JSON) 30 REAL*8 :: XTAB(JSON) 31 REAL*8 :: YTAB(JSON) 32 33 END MODULE SPLIN_MOD

1"""Python version of splin_mod.f (YREC 7be4133): work arrays for interpolation with osplin, formerly COMMON/SPLIN/. 2 3Translates: splin_mod.f#SPLIN_MOD (7be4133) 4 5Module-level names, as in the Fortran module, in lower case. Routines that 6USE SPLIN_MOD import this module and write ``splin_mod.<name>``. 7Scratch space of eq2mod, mixcom, mixgrid and rotgrid. Index 0 unused. 8""" 9 10import numpy as np 11 12from pyrec import dimensions_mod 13 14_JSON = dimensions_mod.json 15 16xval = np.zeros(_JSON + 1) # abscissas of evaluation 17yval = np.zeros(_JSON + 1) # values returned by osplin 18xtab = np.zeros(_JSON + 1) # the table's abscissas 19ytab = np.zeros(_JSON + 1) # the table's values