fac_mod.f

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

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

Variables on this page (8)

imin, imin2, iminmax, iminmax2, jp, kp, nt, ntc

Purpose

FAC_MOD holds the range of valid gravities at each temperature of the two surface-pressure tables, and the interpolation indices of the last lookup, shared by SETUPS, SURFP and KCSURFP. It replaces COMMON/FAC/ (MHP via Claude, 10/14/26). The table dimensions NT = 57 and NTC = 76 are PRIVATE, since the importing files declare their own.

What its variables hold, who sets them, who reads them. All are INTEGER*4.

IMIN(NT), IMINMAX(NT): for each temperature of the surface-pressure table, the lowest and highest gravity index with a valid pressure (entries ≤ 0 are missing, −999), set by SETUPS when it reads the table ("G Somers 5/15"); read by SURFP. When no pressure is valid at a temperature, IMIN is the number of gravities, and SETUPS's comment expects the code to fail there.
IMIN2(NTC), IMINMAX2(NTC): the same for the 76-temperature table, read by KCSURFP.
JP, KP(4): the temperature index and the four gravity indices of the last interpolation, set by SURFP and KCSURFP; read nowhere.

Design

A blanket USE FAC_MOD.

Definitions

moduleFAC_MOD

lines 23–38 · GitHub

Used by (3 units)

KCSURFP, SETUPS, SURFP

Change log

The changes since rp_1

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

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C FAC_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/FAC/ 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. 19C 20C The dimension parameters are PRIVATE: the importing files 21C declare their own PARAMETER statements for these, and an 22C exported copy would collide with them. 23 MODULE FAC_MOD 24 25 IMPLICIT NONE 26 SAVE 27 28 INTEGER*4, PARAMETER, PRIVATE :: NT = 57 29 INTEGER*4, PARAMETER, PRIVATE :: NTC = 76 30 31 INTEGER*4 :: IMIN(NT) 32 INTEGER*4 :: IMINMAX(NT) 33 INTEGER*4 :: JP 34 INTEGER*4 :: KP(4) 35 INTEGER*4 :: IMIN2(NTC) 36 INTEGER*4 :: IMINMAX2(NTC) 37 38 END MODULE FAC_MOD

1"""Python version of fac_mod.f (YREC 7be4133): the valid gravity ranges of 2the surface-pressure tables, formerly COMMON/FAC/. 3 4Translates: fac_mod.f#FAC_MOD (7be4133) 5 6Routines that USE FAC_MOD import this module and write ``fac_mod.<name>``. 7All INTEGER*4; arrays with index 0 unused. 8""" 9 10import numpy as np 11 12# Table dimensions, PRIVATE in the Fortran (Fortran lines 28–29). 13_NT = 57 14_NTC = 76 15 16# Fortran lines 31–36. 17imin = np.zeros(_NT + 1, dtype=int) # lowest valid gravity per T (setups) 18iminmax = np.zeros(_NT + 1, dtype=int) # highest valid gravity per T (setups) 19jp = 0 # last T index (surfp, kcsurfp); read nowhere 20kp = np.zeros(5, dtype=int) # last gravity indices; read nowhere 21imin2 = np.zeros(_NTC + 1, dtype=int) # the same for the 76-T table (kcsurfp) 22iminmax2 = np.zeros(_NTC + 1, dtype=int)