errmom_mod.f

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

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

Variables on this page (1)

toleri

Purpose

ERRMOM_MOD holds TOLERI, the relative tolerance of the iteration for the moment of inertia and the surface rotation rate, shared by CHECKJ and SOLID. It replaces COMMON/ERRMOM/ (MHP via Claude, 10/14/26).

What its variable holds, who sets it, who reads it.

TOLERI (REAL*8): $10^{-12}$, its initial value in the module (no namelist sets it). CHECKJ relaxes it to at least DJOK/100 around its call to GETROT in the early iterations and then restores it. SOLID iterates the surface $\Omega$ until the angular momentum of the convection zone matches within TOLERI (at most 20 tries).

Design

A blanket USE ERRMOM_MOD. The tolerance reaches GETROT's callees through the module rather than as an argument.

Definitions

moduleERRMOM_MOD

lines 19–30 · GitHub

Used by (2 units)

CHECKJ, SOLID

Change log

The changes since rp_1

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

Changed since a8a5950. TOLERI is given a value, $10^{-12}$: Nothing set it before, so SOLID's tolerance was 0 and every call ran all 20 Newton passes. This settles former question 1.

1C 2C 3C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 4C ERRMOM_MOD 5C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6C 10/14/26 MHP via Claude 7C code changed to convert COMMON/ERRMOM/ into this module. 8C 9C Imported with a blanket USE, the convention for the conversion. 10C All 2 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 ERRMOM_MOD 20 21 IMPLICIT NONE 22 SAVE 23 24C 10/30/26 MHP via Claude 25C code changed to give TOLERI a value. Nothing set it, so SOLID's 26C relative J tolerance was 0 and every call took all 20 Newton passes; 27C CHECKJ only loosens it for the transport's early iterations. 28 REAL*8 :: TOLERI = 1.0D-12 29 30 END MODULE ERRMOM_MOD

1"""Python version of errmom_mod.f (YREC 7be4133): the tolerance of the 2moment-of-inertia iteration, formerly COMMON/ERRMOM/. 3 4Translates: errmom_mod.f#ERRMOM_MOD (7be4133) 5 6Routines that USE ERRMOM_MOD import this module and write 7``errmom_mod.toleri``. 8""" 9 10# Given a value in the module; checkj only relaxes it temporarily (Fortran 11# line 28). 12toleri = 1.0e-12