wrthead.f

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

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

Variables on this page (9)

iowr, iresca, itrack, itrver, lfirst, ltrack, nmodls, xenv0, zenv0

Purpose

WRTHEAD writes the header lines of the output files at the start of each run of a job: a line describing the run in the main output, and the track file's column header.

How it fits

Inputs: the arguments, the stellar mass SMASS (solar masses), the mixing length CMIXL and the run index NK; and through modules: the run's kind IRESCA(NK) (RUNLOOP_MOD, formerly CKIND), its number of models NMODLS(NK) and LFIRST(NK) (NMLPAR_MOD, namelist settings), the initial $X$ and $Z$ (XENV0, ZENV0 in LABEL_MOD), the namelist switch LTRACK and the track format's version ITRVER (NMLPAR_MOD, ITRVER 1 by default; read by PARMIN), and the units IOWR, ITRACK (IOUNITS_MOD).
Who calls it: the main program, main.f, once per run.
Outputs: a line on IOWR ("RUN= … EVOLVE," "RESCALE" or "RESCALE&EVOLVE"); with LTRACK, for a run that reads its starting model (LFIRST(NK)), a version line and the column names on ITRACK. The data lines under these headers are written by WRTOUT, which must match them.

Design

The column names depend on ITRVER: 0, one long line of names (G. Somers 11/14); 1, a description in comment lines (except two lines of format 1505, " Neutrinos …" and ' " " cont: …', which do not start with #); 2, the same plus the rotation columns (whose line, format 1510, does not start with # either); 3, a short list. For ITRVER = 0 the 83 names of format 1504 match, in number and order, the 83 values that WRTOUT writes with format 1499; nothing in the code checks this, and the comment asks that ITRVER change whenever the track file does (lines 89–90 ↗) (settled from the code, Nov 9, 2026: lines 131–148 ↗; wrtout.f, lines 558–568 ↗). The headers of versions 1–3 do not match their data (bug 1).

For users (from YREC's User Guide). The track file is the run's history, one line per model with ITRVER = 0, the format recommended. The two lines of DESCRIP in the control namelist, the run's description, head the track file too, but PARMIN writes them, before this header.

From YREC's documentation

Control Options: ITRVER chooses the format of the .track file; 0 (one line per model) is strongly recommended, though the default is 1.

Suspected bugs, checked by running the code

Each bug is marked by what running the code showed: (changes results, large) or (small), the fix moving a run clearly beyond a null change; (stops the run); (wrong output), the fix changing only what is written; (harmless in the runs), within the noise; (never occurred), the code ran but the fault did not; (not reached); (cannot be reached), a guard prevents it; (withdrawn), the suspicion was wrong.

1. (not reached) The headers of track versions 1, 2, and 3 do not match the data.

(lines 152–176 ↗; wrtout.f, lines 570–641 ↗) For ITRVER = 1 (the default) and 2, WRTOUT writes for each model: a first line of 14 values whose last five are the base of the convection zone (CZsurf_R, CZsurf_T, CZsurf_D, CZsurf_P, CZsurf_O; wrtout.f, lines 570–572 ↗), which format 1505 does not name; a central line whose ninth value is the total moment of inertia HITOT (wrtout.f, lines 573–575 ↗), not named; a neutrino line of the Cl and Ga rates and 10 fluxes (wrtout.f, lines 587–588 ↗), where the header names 8 fluxes and no rates; and a second surface line of H2, Li6, Li7, Be9, X, Y, Z, and Z/X (wrtout.f, lines 591–593 ↗), where the header names the first four. For version 2 the rotation line has nine values (SJTOT, SKEROT, HITOT, HICZ, Ω at the surface and center, ROTP, ROTV, TAUCZ; wrtout.f, lines 632–634 ↗), where format 1510 names four. For version 3 the tenth value is CZsurf_R (wrtout.f, lines 639–641 ↗), where format 1515 says "% Grav Energy" (the line that computed it is commented out, wrtout.f, lines 637–638 ↗). And three header lines do not start with # (lines 158–159 ↗; line 164 ↗; line 171 ↗), so a reader that skips # lines takes them for data. What follows: a user who reads a track of version 1, 2, or 3 by its header assigns columns wrongly; the run itself is not affected. On by default (LTRACK true and ITRVER 1: nmlpar_mod.f, line 59 ↗; nmlpar_mod.f, line 142 ↗), though the documentation recommends version 0.

A possible fix (the comment at lines 89–90 ↗ may also call for a new ITRVER):

wrthead.f

@@ lines 149–179 @@
 C G Somers END.
          ELSE IF(ITRVER .EQ. 1) THEN
             WRITE(ITRACK, 1505)
+C 11/09/26 KC via Claude
+C code changed to make the headers of ITRVER = 1, 2 and 3 name the
+C columns WRTOUT writes, each line starting with '#'.
  1505       FORMAT(
      1'# Model #, shells, AGE, log(L/Lsun), log(R/Rsun), log(g),',
-     1' log(Teff), Mconv. core, Mconv. env.)' ,/,
-     2'# Central: log(T), log(RHO), log(P), BETA, ETA, X, Y, Z',/,
+     1' log(Teff), Mconv. core, Mconv. env.; CZ base: R, T, Rho,',
+     1' P, kappa',/,
+     2'# Central: log(T), log(RHO), log(P), BETA, ETA, X, Y, Z,',
+     2' total I',/,
      3'# Luminosity: ppI, ppII, ppIII, CNO, triple-alpha,',
      3' He-C, gravity, neutrinos (old)',/,
-     3' Neutrinos (1E10 erg/CM^2 at earth): pp, pep, hep, Be7,',
-     3' B8, N13, O15, F17',/,
+     3'# Cl SNU, Ga SNU, Neutrinos (1E10/CM^2/S at earth): pp,',
+     3' pep, hep, Be7, B8, N13, O15, F17, diag1, diag2',/,
      4'# Central Abundances: He3, C12, C13, N14, N15, O16,',
      4' O17, O18',/,
      5'# Surface Abundances: He3, C12, C13, N14, N15, O16,',
      5' O17,O18',/,
-     6'    "        " cont: H2, Li6, Li7, Be9',/,
+     6'#    "        " cont: H2, Li6, Li7, Be9, X, Y, Z, Z/X',/,
      7'# H shell loc: mass frac-base, midpoint, top; radius frac-',
      7'base, midpoint, top')
          ELSE IF(ITRVER .EQ. 2) THEN
             WRITE(ITRACK, 1505)
             WRITE(ITRACK, 1510)
  1510       FORMAT(
-     1' Jtot, K.E. Rotation, OMEGAsurf, OMEGAcenter')
+     1'# Jtot, K.E. Rotation, total I, CZ I, OMEGAsurf,',
+     1' OMEGAcenter, Prot (days), Vrot (km/s), TauCZ (s)')
          ELSE IF(ITRVER .EQ. 3) THEN
           WRITE(ITRACK, 1515)
  1515       FORMAT(
      1'#Model #, shells, AGE, log(L/Lsun), log(R/Rsun), log(g),',
-     1' log(Teff), Mconv. core, Mconv. env., % Grav Energy, X env')
+     1' log(Teff), Mconv. core, Mconv. env., R CZ base, X env')
          END IF
       END IF

Checked by running YREC (Nov 9, 2026; a line added at each suspected path counted how often it was taken, the outputs otherwise unchanged): not reached in any of 34 runs (26 of YREC's cases, from a 0.04 $M_\odot$ brown dwarf and the solar models to 9 $M_\odot$, rotation, and helium burning, and 8 made-up ones). The header was written in the 33 runs that got past their start (once each, three times in examples/run_standard_solar_model/Test_solarcal_noGS_norot, which has several runs), always with ITRVER = 0: Every one of YREC's cases sets it, and the made-up ones keep their case's value. Not run: the default ITRVER = 1, which a namelist without the line gets, or 2 or 3. Intent: The comment at lines 89–90 ↗ asks for a new version number whenever the track file changes, so the lines seem to have grown without their headers.

Definitions

subroutineWRTHEAD

lines 12–181 · GitHub

Writes the run line and the track file's header (by version ITRVER).

Uses modules IOUNITS_MOD, LABEL_MOD, NMLPAR_MOD, RUNLOOP_MOD

Called by (file level)

Change log

The changes since rp_1

Changed since rp_1. CMIXL and NK are now arguments (from the main program) instead of COMMON members; the isochrone header (the LISO option, on IISO) is removed with the option; declarations: COMMON blocks to modules, IMPLICIT NONE.

1C 2C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 3C WRTHEAD 4C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 5C write the headers for all the appropriate output files 6 7C 10/9/26 MHP via Claude 8C code changed to pass NK as an explicit argument instead of via 9C COMMON/ZRAMP/, which this file otherwise never used (NK was its only 10C member ever referenced here). See the CALL WRTHEAD site in main.f, 11C which already has NK as its own local run-index variable. 12 SUBROUTINE WRTHEAD (SMASS,CMIXL,NK) 13C 10/14/26 MHP via Claude 14C code changed to convert COMMON/LABEL/ into LABEL_MOD. 15 USE LABEL_MOD 16C 10/14/26 MHP via Claude 17C code changed to convert COMMON/NMLPAR/ into NMLPAR_MOD. 18 USE NMLPAR_MOD 19C 10/14/26 MHP via Claude 20C code changed to convert COMMON/IOUNITS/ into IOUNITS_MOD. 21 USE IOUNITS_MOD 22C 10/14/26 MHP via Claude 23C code changed to convert COMMON/RUNLOOP/ into RUNLOOP_MOD. A blanket 24C USE imports exactly the names the COMMON declaration made global here, 25C so it cannot change meaning. See runloop_mod.f for the block history 26C and for why the module declares the types the implicit rules supplied. 27 USE RUNLOOP_MOD 28C 10/20/26 MHP via Claude 29C code changed to replace implicit typing with IMPLICIT NONE and 30C explicit declarations. 31 IMPLICIT NONE 32 REAL*8 CMIXL,SMASS 33 INTEGER NK 34 35C MHP 8/25 Removed character file names from common block 36C MHP 10/02 added proper dimensions for flaol2, fopal2 37C CHARACTER*256 FISO,FLAOL2, FOPAL2 38C 10/12/26 MHP via Claude 39C code changed to extend COMMON/NMLPAR/ to 234 members, adding the 40C namelist variables that the physics routines read from 27 COMMON 41C blocks which also held derived or working members. Those blocks keep 42C only their non-namelist members; five of them -- NEWMX, VMULT, 43C PMMWIND, CCOUT1 and ENVGEN -- had nothing left and are gone. 44C Derivation-only namelist variables from the same blocks went to 45C COMMON/NMLBYPASS/ instead, not here. 46C 47C Declared identically in all 121 files that hold any of these members. 48C Generated from one member list, so every copy is byte-identical; the 49C linker's 'size of symbol changed' warning is the check that they are. 50C 10/13/26 MHP via Claude: removed the unused COMMON/CONSTANTS/ declaration 51C (no member referenced here; residue of the obsolete-option removals). 52C 9/24/26 MHP via Claude 53C code changed to pass CMIXL as an explicit argument instead of via 54C common block, to reduce implicit global state ahead of the F90 55C module conversion 56C 9/25/26 MHP via Claude 57C code changed to consolidate all I/O logical unit numbers, previously 58C scattered across COMMON/LUOUT/, COMMON/LUNUM/, COMMON/ALEX06/, 59C COMMON/ACDPTH/, COMMON/NEWOPAC/, COMMON/LOPAL95/, COMMON/ATMOS2/, 60C COMMON/ALATM03/, COMMON/OPALEOS/, and COMMON/SCV2/, into a single 61C COMMON/IOUNITS/ 62C 9/24/26 MHP via Claude 63C code changed to remove dead ILLDAT, IKUR from COMMON/LUNUM/ 64C 9/24/26 MHP via Claude 65C code changed to remove dead IDYN from COMMON/LUNUM/ 66C 9/25/26 MHP via Claude 67C code changed to remove dead COMMON/LUNUM/ISCOMP (declared, 68C assigned a unit number, but never opened or read anywhere) 69C MHP 8/25 Removed character file names from common block 70C 9/25/26 MHP via Claude 71C code changed to remove the LISO isochrone-output option 72 SAVE 73 74 IF (IRESCA(NK) .EQ. 1) THEN 75 WRITE(IOWR, 47) NK, XENV0, ZENV0, CMIXL, NMODLS(NK) 76 ELSE IF (IRESCA(NK) .EQ. 2) THEN 77 WRITE(IOWR, 48) NK, XENV0, ZENV0, CMIXL, NMODLS(NK) 78 ELSE IF (IRESCA(NK) .EQ. 3) THEN 79 WRITE(IOWR, 49) NK, XENV0, ZENV0, CMIXL, NMODLS(NK) 80 END IF 81 47 47 FORMAT(/, ' RUN=',I2,' EVOLVE ', ' X=',F8.6, 82 * * ' Z=',F8.6,' CMIXL=', F8.6, ' NO.MODS=', I5) 83 48 48 FORMAT(/, ' RUN=',I2,' RESCALE ', ' X=',F8.6, 84 * * ' Z=',F8.6,' CMIXL=', F8.6, ' NO.MODS=', I5) 85 49 49 FORMAT(/, ' RUN=',I2,' RESCALE&EVOLVE ', ' X=',F8.6, 86 * * ' Z=',F8.6,' CMIXL=', F8.6, ' NO.MODS=', I5) 87 88 IF (LTRACK .AND. LFIRST(NK)) THEN 89C ITRVER identifies version of track out file. If you change 90C the track out file then change this version number. 91 WRITE(ITRACK, 1500)ITRVER,SMASS,XENV0,ZENV0,CMIXL 92 1500 1500 FORMAT('#Version=',i3,' Mtot/Msun =',1PE16.8, 93 * * ' Initial: X =',1PE16.8,' Z =',1PE16.8, 94 * * ' Mix. length =', 1PE16.8) 95 IF(ITRVER .EQ. 0) THEN 96c WRITE(ITRACK, 1503) 97c 1503 FORMAT( 98c 1'# Model #, shells, AGE(Gyr), log(L/Lsun), log(R/Rsun), log(g),', 99c 1' log(Teff), Mconv. core (Msun), Mconv. env, R,T,Rho,P,cappa env',/, 100c 2'# Central: log(T), log(RHO), log(P), BETA, ETA, X, Y, Z',/, 101c 3'# Luminosity: ppI, ppII, ppIII, CNO, triple-alpha,', 102c 3' He-C, gravity, neutrinos (old)',/, 103c 3'# Cl SNU, Ga SNU, Neutrinos (1E10 erg/CM^2 at earth): pp, pep, hep, Be7,', 104c 3' B8, N13, O15, F17 2xdiagnostic',/, 105c 4'# Central Abundances: He3, C12, C13, N14, N15, O16,', 106c 4' O17, O18',/, 107c 5'# Surface Abundances: He3, C12, C13, N14, N15, O16,', 108c 5' O17,O18 H2, Li6, Li7, Be9 X Y Z Z/X',/, 109c 6'# Jtot, KE rot tot, total I, CZ I, Omega center, surface, Prot (day), Vrot (km/s), TauCZ (sec) ',/, 110c 7'# H shell loc: mass frac-base, midpoint, top; radius frac-', 111c 7'base, midpoint, top, Pphot, mass (msun)') 112C G Somers 11/14; Added option to create .track file header. Uncomment the following 113C block if desired. 114 WRITE(ITRACK, 1504) 115C 1504 FORMAT( 116C 1'# ',/, 117C 2'# Step Shls Age (Gyr) log(L/Lsun) log(R/Rsun) log(g) log(Teff)', 118C 2' Mconv.core Mconv.env Rcore Tcore Rho_core P_core kappa_env ', 119C 2' log(T)_cen log(Rho)_cen log(P)_cen BETA ETA X_cen ', 120C 2' Y_cen Z_cen ppI_lum ppII_lum ppIII_lum CNO_lum ', 121C 2' 3-alpha_lum He-C_lum gravity OLD NUTRINOS Cl SNU Ga SNU **pp** ', 122C 2' **pep** **hep** **Be7** **B8** **N13** **O15** **F17** **diag1** **diag2** ', 123C 2' He3_cen C12_cen C13_cen N14_cen B10_cen O16_cen ', 124C 2' B11_cen O18_cen He3_surf C12_surf C13_surf N14_surf ', 125C 2' B10_surf O16_surf B11_surf O18_surf H2_surf Li6_surf', 126C 2' Li7_surf Be9_surf X_surf Y_surf Z_surf Z/X_su', 127C 2'rf Jtot KE_rot_tot total I CZ I Omega_surf Omega', 128C 2'_cen Prot (days) Vrot (km/s) TauCZ (s) Mfrac_base Mfrac_midp Mfr', 129C 2'ac_top Rfrac_base Rfrac_midp Rfrac_top P_phot Mass ',/, 130C 2'# ') 131 1504 1504 FORMAT( 132 1 1'# ',/, 133 2 2' Step Shls Age_gyr LogL_lsun LogR_rsun Log_g', 134 2 2' log_Teff Mco_core Mco_env Rco_env Tco_env Dco_env', 135 2 2' Pco_env Oco_env LogT_cen LogD_cen logP_cen Beta_cen', 136 2 2' Eta_cen X_cen Y_cen Z_cen ppI_lsun', 137 2 2' ppII_lsun ppIII_lsun CNO_lsun 3a_lsun HeC_lsun', 138 2 2' Egrav_lsun Neut_lsun Cl_snu Ga_snu pp_neut pep_neut hep_neut', 139 2 2' Be7_neut B8_neut N13_neut O15_neut F17_neut diag1 diag2 He3_cen', 140 2 2' C12_cen C13_cen N14_cen N15_cen O16_cen', 141 2 2' O17_cen O18_cen He3_sur C12_sur C13_sur', 142 2 2' N14_sur N15_sur O16_sur O17_sur O18_sur', 143 2 2' H2_sur Li6_sur Li7_sur Be9_sur X_sur', 144 2 2' Y_sur Z_sur Z_X_sur Jtot KE_rot_tot', 145 2 2' I_tot I_cz Omega_sur Omega_cen Prot_sur_d', 146 2 2' Vrot_kms TauCZ_s MHshell_base MHshell_mid MHshell_top', 147 2 2' RHshell_base RHShell_mid RHshell_top logP_phot Mass_msun',/, 148 2 2'# ') 149C G Somers END. 150 ELSE IF(ITRVER .EQ. 1) THEN 151 WRITE(ITRACK, 1505) 152 1505 1505 FORMAT( 153 1 1'# Model #, shells, AGE, log(L/Lsun), log(R/Rsun), log(g),', 154 1 1' log(Teff), Mconv. core, Mconv. env.)' ,/, 155 2 2'# Central: log(T), log(RHO), log(P), BETA, ETA, X, Y, Z',/, 156 3 3'# Luminosity: ppI, ppII, ppIII, CNO, triple-alpha,', 157 3 3' He-C, gravity, neutrinos (old)',/, 158 3 3' Neutrinos (1E10 erg/CM^2 at earth): pp, pep, hep, Be7,', 159 3 3' B8, N13, O15, F17',/, 160 4 4'# Central Abundances: He3, C12, C13, N14, N15, O16,', 161 4 4' O17, O18',/, 162 5 5'# Surface Abundances: He3, C12, C13, N14, N15, O16,', 163 5 5' O17,O18',/, 164 6 6' " " cont: H2, Li6, Li7, Be9',/, 165 7 7'# H shell loc: mass frac-base, midpoint, top; radius frac-', 166 7 7'base, midpoint, top') 167 ELSE IF(ITRVER .EQ. 2) THEN 168 WRITE(ITRACK, 1505) 169 WRITE(ITRACK, 1510) 170 1510 1510 FORMAT( 171 1 1' Jtot, K.E. Rotation, OMEGAsurf, OMEGAcenter') 172 ELSE IF(ITRVER .EQ. 3) THEN 173 WRITE(ITRACK, 1515) 174 1515 1515 FORMAT( 175 1 1'#Model #, shells, AGE, log(L/Lsun), log(R/Rsun), log(g),', 176 1 1' log(Teff), Mconv. core, Mconv. env., % Grav Energy, X env') 177 END IF 178 END IF 179 180 RETURN 181 END

1"""Python version of wrthead.f (YREC 7be4133): the header lines of the 2output files at the start of each run.""" 3 4from pyrec import iounits_mod 5from pyrec import label_mod 6from pyrec import nmlpar_mod 7from pyrec import runloop_mod 8 9_RUN_KIND = {1: 'EVOLVE ', 2: 'RESCALE ', 3: 'RESCALE&EVOLVE '} 10 11# FORMAT 1504 (ITRVER = 0): one line of column names, 16 records of the Fortran 12# joined into one; it does not start with '#' (G. Somers 11/14) (Fortran 13# lines 131–148). 14_COLUMNS_1504 = ( 15 ' Step Shls Age_gyr LogL_lsun LogR_rsun Log_g' 16 ' log_Teff Mco_core Mco_env Rco_env Tco_env Dco_env' 17 ' Pco_env Oco_env LogT_cen LogD_cen logP_cen Beta_cen' 18 ' Eta_cen X_cen Y_cen Z_cen ppI_lsun' 19 ' ppII_lsun ppIII_lsun CNO_lsun 3a_lsun HeC_lsun' 20 ' Egrav_lsun Neut_lsun Cl_snu Ga_snu pp_neut pep_neut hep_neut' 21 ' Be7_neut B8_neut N13_neut O15_neut F17_neut diag1 diag2 He3_cen' 22 ' C12_cen C13_cen N14_cen N15_cen O16_cen' 23 ' O17_cen O18_cen He3_sur C12_sur C13_sur' 24 ' N14_sur N15_sur O16_sur O17_sur O18_sur' 25 ' H2_sur Li6_sur Li7_sur Be9_sur X_sur' 26 ' Y_sur Z_sur Z_X_sur Jtot KE_rot_tot' 27 ' I_tot I_cz Omega_sur Omega_cen Prot_sur_d' 28 ' Vrot_kms TauCZ_s MHshell_base MHshell_mid MHshell_top' 29 ' RHshell_base RHShell_mid RHshell_top logP_phot Mass_msun' 30) 31_FMT_1504 = '# \n' + _COLUMNS_1504 + '\n# ' 32 33# FORMAT 1505 (ITRVER = 1 and 2). Two of its lines do not start with '#' 34# (Fortran lines 152–166). Suspected bug (see the page): the headers of 35# versions 1-3 do not name all the columns WRTOUT writes. 36_FMT_1505 = '\n'.join([ 37 '# Model #, shells, AGE, log(L/Lsun), log(R/Rsun), log(g),' 38 ' log(Teff), Mconv. core, Mconv. env.)', 39 '# Central: log(T), log(RHO), log(P), BETA, ETA, X, Y, Z', 40 '# Luminosity: ppI, ppII, ppIII, CNO, triple-alpha,' 41 ' He-C, gravity, neutrinos (old)', 42 ' Neutrinos (1E10 erg/CM^2 at earth): pp, pep, hep, Be7,' 43 ' B8, N13, O15, F17', 44 '# Central Abundances: He3, C12, C13, N14, N15, O16,' 45 ' O17, O18', 46 '# Surface Abundances: He3, C12, C13, N14, N15, O16,' 47 ' O17,O18', 48 ' " " cont: H2, Li6, Li7, Be9', 49 '# H shell loc: mass frac-base, midpoint, top; radius frac-' 50 'base, midpoint, top', 51]) 52 53# FORMAT 1510 (ITRVER = 2): the rotation columns, without '#' (Fortran lines 54# 170–171). 55_FMT_1510 = ' Jtot, K.E. Rotation, OMEGAsurf, OMEGAcenter' 56 57# FORMAT 1515 (ITRVER = 3) (Fortran lines 174–176). 58_FMT_1515 = ('#Model #, shells, AGE, log(L/Lsun), log(R/Rsun), log(g),' 59 ' log(Teff), Mconv. core, Mconv. env., % Grav Energy, X env') 60 61 62def wrthead(smass: float, cmixl: float, nk: int) -> None: 63 """Write the run's line on IOWR and, for a run that reads its starting 64 model (LFIRST(NK)) with LTRACK on, the track file's version line and 65 column header on ITRACK, in the layout chosen by ITRVER. 66 67 Translates: wrthead.f#WRTHEAD (7be4133) 68 69 Arguments: the stellar mass smass (solar masses), the mixing length 70 cmixl and the run index nk. Nothing is returned. 71 72 """ 73 74 n = nmlpar_mod 75 r = runloop_mod 76 xenv0 = label_mod.xenv0 77 zenv0 = label_mod.zenv0 78 79 # FORMATs 47, 48, 49: the same line but for the run's kind (Fortran lines 80 # 74–86). 81 kind = r.iresca[nk] 82 if kind in _RUN_KIND: 83 iounits_mod.write( 84 iounits_mod.iowr, 85 f"\n RUN={nk:2d} {_RUN_KIND[kind]} X={xenv0:8.6f}" 86 f" Z={zenv0:8.6f} CMIXL={cmixl:8.6f}" 87 f" NO.MODS={r.nmodls[nk]:5d}") 88 89 if n.ltrack and r.lfirst[nk]: 90 # ITRVER identifies the version of the track file: change it when 91 # the track file changes (and WRTOUT with it) (Fortran lines 88–94). 92 iounits_mod.write( 93 iounits_mod.itrack, 94 f"#Version={n.itrver:3d} Mtot/Msun ={smass:16.8E}" 95 f" Initial: X ={xenv0:16.8E} Z ={zenv0:16.8E}" 96 f" Mix. length ={cmixl:16.8E}") 97 98 # The column header by version (Fortran lines 95–177). 99 if n.itrver == 0: 100 # Question (see the page): the comment says to uncomment a 101 # block, but FORMAT 1504's second version is active; its line 102 # of column names does not start with '#'. 103 iounits_mod.write(iounits_mod.itrack, _FMT_1504) 104 elif n.itrver == 1: 105 iounits_mod.write(iounits_mod.itrack, _FMT_1505) 106 elif n.itrver == 2: 107 iounits_mod.write(iounits_mod.itrack, _FMT_1505) 108 iounits_mod.write(iounits_mod.itrack, _FMT_1510) 109 elif n.itrver == 3: 110 iounits_mod.write(iounits_mod.itrack, _FMT_1515)