{ "yrec_commit": "7be4133", "checked": "2026-11-09/10", "null_change": "nearly every opacity (3 <= log T <= 9, -6 <= log rho <= 9.75) times 1 + 1e-10", "bugs": [ { "file": "acdpth_mod.f", "number": 1, "result": "cannot be reached", "status": "unconfirmed", "title": "A switch that cannot be turned on", "runs": "`LACOUT` is only read in `src/` and is in no namelist group, so the acoustic-depth branch needs it set true in the source; none of 34 runs reached it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/acdpth_mod.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/acdpth_mod.f#L20-L25", "https://github.com/yreclab/yrec/blob/7be4133/src/acdpth_mod.f#L55" ] }, { "file": "alex06tab.f", "number": 1, "result": "harmless in the runs", "status": "unconfirmed", "title": "The 4th node in $X$ is tested on IZ, not on IX", "runs": "In `testsuite/Test_m1p4feh+0p0GS98_dif_TAMS`, case (ii) held in 363 calls; the fix moves log L by up to 0.0062 dex, against 0.057 for the null change; case (i) needs a surface Z above 0.06.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/alex06tab.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/alex06tab.f#L56-L68", "https://github.com/yreclab/yrec/blob/7be4133/src/alex06tab.f#L91-L104", "https://github.com/yreclab/yrec/blob/7be4133/src/alex06tab.f#L117-L121" ] }, { "file": "alfilein.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "Two tolerances for [Fe/H] and $\\alpha$", "runs": "In `testsuite/Test_solarGS98_allard`, the only Allard run, the test ran 316 times and never dropped a record; it needs a requested [Fe/H] between 1e-6 and 1e-5 off a table value.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/alfilein.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/alfilein.f#L209-L210", "https://github.com/yreclab/yrec/blob/7be4133/src/alfilein.f#L296-L297", "https://github.com/yreclab/yrec/blob/7be4133/src/alfilein.f#L326-L327" ] }, { "file": "alfilein.f", "number": 2, "result": "never occurred", "status": "unconfirmed", "title": "Loop bounds for the new [Fe/H] and $\\alpha$ values", "runs": "In `testsuite/Test_solarGS98_allard`, the loops ran for 315 records, all with [Fe/H] = 0 and alpha = 0, so no unset entry was compared; it needs records of several [Fe/H] or alpha.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/alfilein.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/alfilein.f#L248-L262" ] }, { "file": "altabinit.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "A format mismatch stops the run with a run-time error", "runs": "In `testsuite/Test_solarGS98_allard`, the test ran once, on a table of 62 rows, and the message was never written; a requested `ALATM_FeH` not in the table, such as -0.25, would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/altabinit.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/altabinit.f#L119-L124" ] }, { "file": "altabinit.f", "number": 2, "result": "never occurred", "status": "unconfirmed", "title": "A row with no valid entry is read out of bounds", "runs": "In `testsuite/Test_solarGS98_allard`, no row of the 62 had fewer than two valid entries, and every row of the three shipped Allard tables holds at least four; only a bad table would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/altabinit.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/altabinit.f#L68-L84", "https://github.com/yreclab/yrec/blob/7be4133/src/altabinit.f#L92-L98" ] }, { "file": "altabinit.f", "number": 3, "result": "not reached", "status": "unconfirmed", "title": "An empty table is read out of bounds", "runs": "ALTABINIT ran once, in `testsuite/Test_solarGS98_allard`, with 62 rows and 6 columns; no case requests an [Fe/H] or alpha that its table lacks.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/altabinit.f.html#b3", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/altabinit.f#L61-L64", "https://github.com/yreclab/yrec/blob/7be4133/src/altabinit.f#L119-L124", "https://github.com/yreclab/yrec/blob/7be4133/src/alfilein.f#L209-L213" ] }, { "file": "bsrotmix.f", "number": 1, "result": "cannot be reached", "status": "unconfirmed", "title": "$^3$He is left out of the convergence test", "runs": "GETW sets `LBURS` false on every call, and BURSMIX, the only caller of BSROTMIX, runs only when it is true; no namelist changes that, and no run reached it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/bsrotmix.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/bsrotmix.f#L88-L96", "https://github.com/yreclab/yrec/blob/7be4133/src/bsrotmix.f#L213-L220" ] }, { "file": "bsrotmix.f", "number": 2, "result": "cannot be reached", "status": "unconfirmed", "title": "STOP999 when no species is present", "runs": "GETW sets `LBURS` false on every call, and BURSMIX, the only caller of BSROTMIX, runs only when it is true; no namelist changes that, and no run reached it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/bsrotmix.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/bsrotmix.f#L76-L77" ] }, { "file": "bsstep.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "The step-size underflow test uses the hydrogen abundance", "runs": "In 33 runs, failed steps were shrunk up to 25,983 times (the 0.04 M☉ brown dwarf), but neither `X+H.EQ.X` nor `X0+H.EQ.X0` ever held; it would take an integration that fails completely.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/bsstep.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/bsstep.f#L64" ] }, { "file": "bsstep.f", "number": 2, "result": "never occurred", "status": "unconfirmed", "title": "No check of IMAX and NUSE against the sizes of the arrays", "runs": "BSSTEP ran in 33 runs, up to 4,437,788 times each, always with `IMAX` = 11 and `NUSE` = 7, the defaults that every case keeps; a run with other values would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/bsstep.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/bsstep.f#L27-L28", "https://github.com/yreclab/yrec/blob/7be4133/src/bsstep.f#L38-L42", "https://github.com/yreclab/yrec/blob/7be4133/src/bsstep.f#L57", "https://github.com/yreclab/yrec/blob/7be4133/src/ratext.f#L11-L12" ] }, { "file": "bursmix.f", "number": 1, "result": "cannot be reached", "status": "unconfirmed", "title": "NDIFCOM gets a region never set when LOK is false", "runs": "GETW sets `LBURS` false on every call and calls BURSMIX only when it is true; no namelist changes that, and no run reached it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/bursmix.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/bursmix.f#L28-L31", "https://github.com/yreclab/yrec/blob/7be4133/src/bursmix.f#L46-L49", "https://github.com/yreclab/yrec/blob/7be4133/src/bursmix.f#L103-L106" ] }, { "file": "calcad.f", "number": 1, "result": "cannot be reached", "status": "unconfirmed", "title": "Step sizes not restored", "runs": "CALCAD runs only when `LACOUT` is true, and `LACOUT` is false in the source and in no namelist group; not reached in 34 runs.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/calcad.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/calcad.f#L196-L207", "https://github.com/yreclab/yrec/blob/7be4133/src/calcad.f#L233-L235", "https://github.com/yreclab/yrec/blob/7be4133/src/wrtmod.f#L237-L242" ] }, { "file": "calcad.f", "number": 2, "result": "cannot be reached", "status": "unconfirmed", "title": "The interior's $Z$", "runs": "CALCAD runs only when `LACOUT` is true, and `LACOUT` is false in the source and in no namelist group; not reached in 34 runs.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/calcad.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/calcad.f#L246-L247", "https://github.com/yreclab/yrec/blob/7be4133/src/calcad.f#L283" ] }, { "file": "calcad.f", "number": 3, "result": "cannot be reached", "status": "unconfirmed", "title": "ICZ may be stale", "runs": "CALCAD runs only when `LACOUT` is true, and `LACOUT` is false in the source and in no namelist group; not reached in 34 runs.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/calcad.f.html#b3", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/calcad.f#L333-L351", "https://github.com/yreclab/yrec/blob/7be4133/src/calcad.f#L370-L373" ] }, { "file": "calcad.f", "number": 4, "result": "cannot be reached", "status": "unconfirmed", "title": "FXIONJ is a scalar", "runs": "CALCAD runs only when `LACOUT` is true, and `LACOUT` is false in the source and in no namelist group; not reached in 34 runs.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/calcad.f.html#b4", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/calcad.f#L82", "https://github.com/yreclab/yrec/blob/7be4133/src/calcad.f#L292", "https://github.com/yreclab/yrec/blob/7be4133/src/eqstat2.f#L54", "https://github.com/yreclab/yrec/blob/7be4133/src/eqstat2.f#L171-L173" ] }, { "file": "calcad.f", "number": 5, "result": "cannot be reached", "status": "unconfirmed", "title": "EOS errors ignored", "runs": "CALCAD runs only when `LACOUT` is true, and `LACOUT` is false in the source and in no namelist group; not reached in 34 runs.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/calcad.f.html#b5", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/calcad.f#L269-L273" ] }, { "file": "calcad.f", "number": 6, "result": "cannot be reached", "status": "unconfirmed", "title": "Output format", "runs": "CALCAD runs only when `LACOUT` is true, and `LACOUT` is false in the source and in no namelist group; not reached in 34 runs.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/calcad.f.html#b6", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/calcad.f#L445-L446" ] }, { "file": "calcad.f", "number": 7, "result": "cannot be reached", "status": "unconfirmed", "title": "A short profile", "runs": "CALCAD runs only when `LACOUT` is true, and `LACOUT` is false in the source and in no namelist group; not reached in 34 runs.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/calcad.f.html#b7", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/calcad.f#L317-L322", "https://github.com/yreclab/yrec/blob/7be4133/src/calcad.f#L441", "https://github.com/yreclab/yrec/blob/7be4133/src/boole.f#L56-L63" ] }, { "file": "calcad.f", "number": 8, "result": "cannot be reached", "status": "unconfirmed", "title": "The surface point is left out of the integrals", "runs": "Cannot be reached at 7be4133: CALCAD runs only with LACOUT, false and set nowhere; a counting line at L310 never fired in five of YREC's cases. The fix was not run.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/calcad.f.html#b8", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/calcad.f#L310-L322", "https://github.com/yreclab/yrec/blob/7be4133/src/calcad.f#L241-L256", "https://github.com/yreclab/yrec/blob/7be4133/src/boole.f#L41-L52" ] }, { "file": "chkcal.f", "number": 1, "result": "wrong output", "status": "checked by KC (Nov 10, 2026)", "title": "The first cycle prints an \"Old BL\" and \"Old RL\" that were never set", "runs": "In `examples/run_standard_solar_model/Test_solarcal_noGS_norot`, the first cycle printed the old luminosity and radius before `BLP` and `RLP` were set; nothing computed reads these two printed lines.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/chkcal.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/chkcal.f#L167-L170", "https://github.com/yreclab/yrec/blob/7be4133/src/chkcal.f#L179-L180" ] }, { "file": "chkcal.f", "number": 2, "result": "never occurred", "status": "unconfirmed", "title": "A target $Z/X$ of 0 sets $Z = 0$", "runs": "In `examples/run_standard_solar_model/Test_solarcal_noGS_norot`, `CALSOLZX` was 0.02292 both times; it needs `LCALS` and `LCALSOLZX` on with `CALSOLZX` at or below 0.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/chkcal.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/chkcal.f#L61", "https://github.com/yreclab/yrec/blob/7be4133/src/chkcal.f#L156-L157" ] }, { "file": "chkscal.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "A false crossing at the start of a run that begins above $R_*$", "runs": "`LCALST` is off by default and in every case; `LCALST` on, with a track whose second model lies above R*, such as a pre-main-sequence start, would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/chkscal.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/chkscal.f#L97-L103", "https://github.com/yreclab/yrec/blob/7be4133/src/setscal.f#L107", "https://github.com/yreclab/yrec/blob/7be4133/src/main.f#L1119-L1123" ] }, { "file": "coefft.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "The nuclear terms at exactly TCUT(1)", "runs": "The store ran in 33 runs, up to about 72 million times each, at log T from 5.5 to 8.40, but log T never equaled `TCUT(1)` exactly.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/coefft.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/coefft.f#L281", "https://github.com/yreclab/yrec/blob/7be4133/src/coefft.f#L416", "https://github.com/yreclab/yrec/blob/7be4133/src/coefft.f#L425-L427" ] }, { "file": "condopacp.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "A point on the table's last temperature or density stops the run", "runs": "The lookup ran in 33 runs, up to 223 million times each, at log T from 3 to 8.40 and log rho from -6 to 7.45; no point fell on log T = 9 or log rho = 9.75.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/condopacp.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/condopacp.f#L99-L109" ] }, { "file": "condopacpint.f", "number": 1, "result": "changes results, small", "status": "checked by KC (Nov 10, 2026)", "title": "The mixture's derivatives are multiplied by $\\log K$", "runs": "In `testsuite/Test_m0p3feh+0p0GS98_scv_ZAMS`, the fix moves log L by up to 1.4e-5 dex, against 1e-9 for the null change; on the rotating 3 M☉ giant the change is within the noise.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/condopacpint.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/condopacpint.f#L68", "https://github.com/yreclab/yrec/blob/7be4133/src/condopacpint.f#L72" ] }, { "file": "condopacpint.f", "number": 2, "result": "changes results, small", "status": "checked by KC (Nov 10, 2026)", "title": "The temperature derivative is taken from the density derivatives", "runs": "With this fix alone, `testsuite/Test_m0p3feh+0p0GS98_scv_ZAMS` moves log L by up to 5.0e-6 dex, against 1e-9 for the null change; on the rotating 3 M☉ giant the change is within the noise.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/condopacpint.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/condopacpint.f#L72" ] }, { "file": "condopacpint.f", "number": 3, "result": "changes results, small", "status": "checked by KC (Nov 10, 2026)", "title": "A sign in the extrapolation below $\\log\\rho = -6$", "runs": "In `examples/run_from_dbl_to_zams/Test_m0040_feh0_noGS_ZAMS`, the fix moves log L by up to 2.4e-4 dex, against 2.5e-5 for the null change, and the end age shifts from 8,681 to 8,692 Gyr.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/condopacpint.f.html#b3", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/condopacpint.f#L92", "https://github.com/yreclab/yrec/blob/7be4133/src/condopacpint.f#L64" ] }, { "file": "condopacpint.f", "number": 4, "result": "harmless in the runs", "status": "unconfirmed", "title": "The step of the extrapolation is reversed", "runs": "The extrapolation went the wrong way in nearly every call, but in `examples/run_from_dbl_to_zams/Test_m0040_feh0_noGS_ZAMS` log L moves by up to 9.5e-5 dex, against 2.5e-5 for the null change.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/condopacpint.f.html#b4", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/condopacpint.f#L99" ] }, { "file": "convec.f", "number": 1, "result": "harmless in the runs", "status": "unconfirmed", "title": "The upper edge of the last radiative region is written to the first", "runs": "In `testsuite/Test_m9p0feh+0p0GS98_base_TAHB`, the wrong element was written 3,344 times, but with the fix the track and `.store` file are identical, since later calls overwrite it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/convec.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/convec.f#L200-L205" ] }, { "file": "convec.f", "number": 2, "result": "changes results, small", "status": "checked by KC (Nov 10, 2026)", "title": "A merger keeps the upper zone's top", "runs": "In `testsuite/Test_m9p0feh+0p0GS98_base_TAHB`, 22 of the 130 mergers cut the overshot core back by one or two shells; the fix moves log L by up to 1.45e-3 dex, against 7.0e-6 for the null change.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/convec.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/convec.f#L161-L173" ] }, { "file": "crrect.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "JENV is unset at levels 1 and 2", "runs": "Two mass-loss runs did not reach it; it needs `LMDOT` with a mass gain (`DMDT0` > 0, `LREIMER` off), which no case or document sets.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/crrect.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/crrect.f#L57", "https://github.com/yreclab/yrec/blob/7be4133/src/crrect.f#L150-L161", "https://github.com/yreclab/yrec/blob/7be4133/src/crrect.f#L176-L179" ] }, { "file": "dadcoeft.f", "number": 1, "result": "cannot be reached", "status": "unconfirmed", "title": "The coefficients on the even grid are never set", "runs": "PARMIN stops any run that sets `LDIFAD` with rotation, and DADCOEFT is called only with `LDIFAD`, from SECULR, which runs only with rotation.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/dadcoeft.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/dadcoeft.f#L244", "https://github.com/yreclab/yrec/blob/7be4133/src/dadcoeft.f#L265-L269", "https://github.com/yreclab/yrec/blob/7be4133/src/dadcoeft.f#L592-L606" ] }, { "file": "dadcoeft.f", "number": 2, "result": "cannot be reached", "status": "unconfirmed", "title": "With ITDIF2 = 1, the default, DADCOEFT always stops the run", "runs": "PARMIN stops any run that sets `LDIFAD` with rotation, and DADCOEFT is called only with `LDIFAD`, from SECULR, which runs only with rotation.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/dadcoeft.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/dadcoeft.f#L179", "https://github.com/yreclab/yrec/blob/7be4133/src/dadcoeft.f#L558-L562", "https://github.com/yreclab/yrec/blob/7be4133/src/dadcoeft.f#L575-L585" ] }, { "file": "dadcoeft.f", "number": 3, "result": "cannot be reached", "status": "unconfirmed", "title": "Array bounds not checked", "runs": "PARMIN stops any run that sets `LDIFAD` with rotation, and DADCOEFT is called only with `LDIFAD`, from SECULR, which runs only with rotation.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/dadcoeft.f.html#b3", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/dadcoeft.f#L39", "https://github.com/yreclab/yrec/blob/7be4133/src/dadcoeft.f#L92", "https://github.com/yreclab/yrec/blob/7be4133/src/dadcoeft.f#L155", "https://github.com/yreclab/yrec/blob/7be4133/src/dadcoeft.f#L534-L535" ] }, { "file": "dadcoeft.f", "number": 4, "result": "cannot be reached", "status": "unconfirmed", "title": "A failed substep is not cut", "runs": "PARMIN stops any run that sets `LDIFAD` with rotation, and DADCOEFT is called only with `LDIFAD`, from SECULR, which runs only with rotation.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/dadcoeft.f.html#b4", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/dadcoeft.f#L183", "https://github.com/yreclab/yrec/blob/7be4133/src/dadcoeft.f#L558-L564", "https://github.com/yreclab/yrec/blob/7be4133/src/dadcoeft.f#L575" ] }, { "file": "dadcoeft.f", "number": 5, "result": "cannot be reached", "status": "unconfirmed", "title": "An unset entry in the printed history", "runs": "PARMIN stops any run that sets `LDIFAD` with rotation, and DADCOEFT is called only with `LDIFAD`, from SECULR, which runs only with rotation.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/dadcoeft.f.html#b5", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/dadcoeft.f#L212", "https://github.com/yreclab/yrec/blob/7be4133/src/dadcoeft.f#L534", "https://github.com/yreclab/yrec/blob/7be4133/src/dadcoeft.f#L544" ] }, { "file": "dburn.f", "number": 1, "result": "harmless in the runs", "status": "unconfirmed", "title": "FMASSACC is read under mass loss", "runs": "In a made-up case, `testsuite/Test_m3p0feh+0p0GS98_std25_TAMS` with `LMDOT` and `LREIMER` on, `FMASSACC` was read unset but 0; log L moves by up to 6.2e-7 dex, against 2e-7 for the null change.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/dburn.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/dburn.f#L75", "https://github.com/yreclab/yrec/blob/7be4133/src/dburn.f#L106", "https://github.com/yreclab/yrec/blob/7be4133/src/dburn.f#L128-L132" ] }, { "file": "dburnm.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "FMASSACC is read under mass loss", "runs": "DBURNM runs only in MIDMOD, called by GETW in rotating runs, and the mass-loss runs did not rotate; it needs `LMDOT` with `LROT`, which no case sets.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/dburnm.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/dburnm.f#L109", "https://github.com/yreclab/yrec/blob/7be4133/src/dburnm.f#L140-L145" ] }, { "file": "engeb.f", "number": 1, "result": "changes results, small", "status": "checked by KC (Nov 10, 2026)", "title": "The Itoh branch replaces the temperature derivative by the density derivative", "runs": "In `testsuite/Test_m2p0_feh+0p0_vzams150_IES2_TAMS`, the fix moves log L by up to 1.83e-5 dex, against 1.58e-8 for the null change, with the same models and end age.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/engeb.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/engeb.f#L1255-L1256" ] }, { "file": "engeb.f", "number": 2, "result": "changes results, small", "status": "checked by KC (Nov 10, 2026)", "title": "Relative derivatives added to absolute ones in the Itoh branch", "runs": "With this fix alone, `testsuite/Test_m2p0_feh+0p0_vzams150_IES2_TAMS` moves log T_c by 3.5e-7 dex, against 1.72e-8 for the null change; in `IES1` the change is within the noise.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/engeb.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/engeb.f#L1247-L1249", "https://github.com/yreclab/yrec/blob/7be4133/src/engeb.f#L1255-L1256" ] }, { "file": "engeb.f", "number": 3, "result": "changes results, small", "status": "checked by KC (Nov 10, 2026)", "title": "The alpha captures' temperature derivatives", "runs": "In `testsuite/Test_m0p95feh+0p0_GN93_TAHB`, the fix moves log L by up to 2.76e-5 dex, against 1.67e-7 for the null change; in the rotating 3 M☉ giant the change is within the noise.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/engeb.f.html#b3", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/engeb.f#L896-L898", "https://github.com/yreclab/yrec/blob/7be4133/src/engeb.f#L909-L912", "https://github.com/yreclab/yrec/blob/7be4133/src/engeb.f#L925-L927" ] }, { "file": "engeb.f", "number": 4, "result": "not reached", "status": "unconfirmed", "title": "The low-temperature path returns stale values", "runs": "Only YTIME could reach it, at the hottest point of a model with log rho_c of 5 or more and a temperature at or below `TCUT(1)`; no case gives such a model.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/engeb.f.html#b4", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/engeb.f#L537-L550", "https://github.com/yreclab/yrec/blob/7be4133/src/engeb.f#L1081-L1095", "https://github.com/yreclab/yrec/blob/7be4133/src/engeb.f#L1132-L1219" ] }, { "file": "envint.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "ENVG never set", "runs": "`LNEWTCZ` is true by default and in every case, so the call to TAUCAL is skipped; a rotating run with `LNEWTCZ` false would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/envint.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/envint.f#L82", "https://github.com/yreclab/yrec/blob/7be4133/src/envint.f#L969-L970", "https://github.com/yreclab/yrec/blob/7be4133/src/taucal.f#L130-L142" ] }, { "file": "envint.f", "number": 2, "result": "wrong output", "status": "checked by KC (Nov 10, 2026)", "title": "Wrong $c_P$ in a print", "runs": "In a made-up case, `testsuite/Test_m1p0feh+0p0GS98_base_TAMS` with `LSTCH` off, the line was written 626 times, each with `QCPP` = 0; with the fix the track is identical, and only that line of the `.store` file changes.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/envint.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/envint.f#L299-L313", "https://github.com/yreclab/yrec/blob/7be4133/src/envint.f#L403-L412" ] }, { "file": "envint.f", "number": 3, "result": "harmless in the runs", "status": "unconfirmed", "title": "The previous point's values on other paths", "runs": "The condition held in all 15 runs with `LSTCH` on and an integrated atmosphere; in `testsuite/Test_m0p3feh+0p0GS98_scv_ZAMS` the fix leaves the track identical and the `.store` file within 2%.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/envint.f.html#b3", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/envint.f#L315-L322", "https://github.com/yreclab/yrec/blob/7be4133/src/envint.f#L433-L474" ] }, { "file": "envint.f", "number": 4, "result": "wrong output", "status": "checked by KC (Nov 10, 2026)", "title": "A stale summary when the envelope is skipped", "runs": "The envelope was skipped once in four runs; in the rotating 3 M☉ case and `testsuite/Test_m2p0_feh+0p0_vzams150_IES2_TAMS`, the fix changes only the `.short` file's summary.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/envint.f.html#b4", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/envint.f#L530-L540", "https://github.com/yreclab/yrec/blob/7be4133/src/envint.f#L925-L930" ] }, { "file": "envint.f", "number": 5, "result": "never occurred", "status": "unconfirmed", "title": "No limit on the retries of the atmosphere's start", "runs": "In `testsuite/Test_m9p0feh+0p0GS98_base_TAHB`, the start was retried 68 times, but one division by 10 always sufficed; a model near the Eddington limit would test the missing limit.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/envint.f.html#b5", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/envint.f#L273-L281", "https://github.com/yreclab/yrec/blob/7be4133/src/envint.f#L367-L372" ] }, { "file": "eqburn.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "Division by a rate that can be zero", "runs": "The explicit burning ran in 30 runs, up to 4,009,747 times, and `GR2`, `GR4`, and `GR5` were never zero; it needs `TCUT(2)` or `TCUT(3)` lowered, or an `SSTANDARD` of 0.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/eqburn.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/eqburn.f#L174", "https://github.com/yreclab/yrec/blob/7be4133/src/eqburn.f#L190-L191" ] }, { "file": "eqscve.f", "number": 1, "result": "harmless in the runs", "status": "unconfirmed", "title": "H$_2$ counted as one hydrogen nucleus", "runs": "In `testsuite/Test_m0p95feh+0p0_GN93_TAHB`, the fix moves log L by at most 3e-8 dex, against 1.67e-7 for the null change, with the same 966 models and end age.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/eqscve.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/eqscve.f#L463", "https://github.com/yreclab/yrec/blob/7be4133/src/setscv.f#L58-L59" ] }, { "file": "eqscvg.f", "number": 1, "result": "changes results, small", "status": "checked by KC (Nov 10, 2026)", "title": "H$_2$ counted as one hydrogen nucleus", "runs": "In `examples/run_from_dbl_to_zams/Test_m0040_feh0_noGS_ZAMS`, the fix moves log L by up to 0.0035 dex, against 2.5e-5 for the null change, and the run ends at 7,928 Gyr instead of 8,681 Gyr.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/eqscvg.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/eqscvg.f#L326" ] }, { "file": "eqscvg.f", "number": 2, "result": "never occurred", "status": "unconfirmed", "title": "The logarithm of a zero mixing entropy", "runs": "The path ran in 28 runs, but no interpolated `SMIX` was zero or negative; an envelope with `XENV` = 0 or 1, which no case has, would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/eqscvg.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/eqscvg.f#L231-L233", "https://github.com/yreclab/yrec/blob/7be4133/src/setscv.f#L70-L71" ] }, { "file": "eqstat.f", "number": 1, "result": "harmless in the runs", "status": "unconfirmed", "title": "A misspelled counter", "runs": "In a made-up case, `testsuite/Test_m1p0feh+0p0GS98_base_TAMS` with `LNumDeriv` on, the counter was wrong in nearly every call, but with the fix the `.track` and `.store` files are identical.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/eqstat.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/eqstat.f#L31", "https://github.com/yreclab/yrec/blob/7be4133/src/eqstat.f#L85", "https://github.com/yreclab/yrec/blob/7be4133/src/eqstat.f#L111-L113" ] }, { "file": "eqstat2.f", "number": 1, "result": "cannot be reached", "status": "unconfirmed", "title": "A dead block with unset variables", "runs": "`LSCVD` is set false just before the test, so the block never runs.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/eqstat2.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/eqstat2.f#L189-L241" ] }, { "file": "equal_to_model.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "Division by $Z$", "runs": "These lines belong to the old settling routines; it takes `LDIFY` or `LDIFZ` with `LNEWDIF` false and `LDIFLI` off, and a shell with Z = 0, which no case has.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/equal_to_model.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/equal_to_model.f#L52", "https://github.com/yreclab/yrec/blob/7be4133/src/equal_to_model.f#L103", "https://github.com/yreclab/yrec/blob/7be4133/src/equal_to_model.f#L122" ] }, { "file": "esac06.f", "number": 1, "result": "harmless in the runs", "status": "unconfirmed", "title": "EOS(2) and EOS(5) rescaled again on every call with IORDER < 5", "runs": "The rescaling held in most calls of 31 runs, with `EOS(5)` up to 6.30e42, far below overflow; with the fix, `testsuite/Test_m1p0feh+0p0GS98_base_TAMS` gives identical `.track` and `.store` files. Kept (b): A comment in RHOOFP06 records crashes from this growth before its tolerance was loosened.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/esac06.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/esac06.f#L408", "https://github.com/yreclab/yrec/blob/7be4133/src/esac06.f#L477-L485" ] }, { "file": "esac06.f", "number": 2, "result": "never occurred", "status": "unconfirmed", "title": "The second $X$ quadratic's warning tests the wrong array", "runs": "Loop 47 ran in 26 runs, up to 84,320,569 times, but the second quadratic never exceeded 1e20; an X stencil reaching the table's empty region would test it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/esac06.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/esac06.f#L452-L458" ] }, { "file": "esac06.f", "number": 3, "result": "not reached", "status": "unconfirmed", "title": "The message at label 66 swaps IP and IQ", "runs": "Never reached in the 31 runs with the OPAL 2006 tables, whose shipped tables share one Z; only a file whose tables differ in Z reaches it, and the run then stops.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/esac06.f.html#b3", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/esac06.f#L498-L499" ] }, { "file": "esac06.f", "number": 4, "result": "cannot be reached", "status": "unconfirmed", "title": "The check on IORDER warns and goes on", "runs": "Every call passes `IORDER` = 1 or 9, so the warning cannot fire.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/esac06.f.html#b4", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/esac06.f#L230-L232" ] }, { "file": "getalex06.f", "number": 1, "result": "harmless in the runs", "status": "unconfirmed", "title": "The derivative is not extrapolated with the value", "runs": "In `examples/run_from_dbl_to_zams/Test_m0040_feh0_noGS_ZAMS`, the derivative disagreed at 45 million of 49 million points, but log L moves by up to 9.4e-5 dex, against 2.5e-5 for the null change.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/getalex06.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/getalex06.f#L90-L96", "https://github.com/yreclab/yrec/blob/7be4133/src/getalex06.f#L108-L111", "https://github.com/yreclab/yrec/blob/7be4133/src/getalex06.f#L113-L120", "https://github.com/yreclab/yrec/blob/7be4133/src/getalex06.f#L126-L129" ] }, { "file": "getfc.f", "number": 1, "result": "changes results, small", "status": "checked by KC (Nov 10, 2026)", "title": "The factor ½ in $\\alpha$", "runs": "In a made-up case, `testsuite/Test_m2p0_feh+0p0_vzams150_IES1_TAMS` with `LVFC` on, surface N-14 moves by up to 0.0039 dex, against 7e-6 for the null change, and the main sequence lasts 1.357 Gyr instead of 1.344.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/getfc.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/getfc.f#L63", "https://github.com/yreclab/yrec/blob/7be4133/src/getfc.f#L75-L79", "https://github.com/yreclab/yrec/blob/7be4133/src/getfc.f#L112" ] }, { "file": "getmodel2.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "The number of shells is not checked", "runs": "The header was read in all 33 runs, with `M` from 441 to 2,118 shells, always within 1 to `JSON`; a damaged model file, or one of more than 5,000 shells, would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/getmodel2.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/getmodel2.f#L57-L58", "https://github.com/yreclab/yrec/blob/7be4133/src/getmodel2.f#L103-L106" ] }, { "file": "getnewenv.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "OMEGA is a local array", "runs": "The routine runs only for a mass gain (`DMDT` > 0), which no case sets; a rotating run with `LMDOT` and `DMDT0` > 0 would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/getnewenv.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/getnewenv.f#L74-L75", "https://github.com/yreclab/yrec/blob/7be4133/src/getnewenv.f#L281-L299", "https://github.com/yreclab/yrec/blob/7be4133/src/getnewenv.f#L300" ] }, { "file": "getnewenv.f", "number": 2, "result": "not reached", "status": "unconfirmed", "title": "A sign error in the interpolation of $X$ and $Z$", "runs": "The routine runs only for a mass gain (`DMDT` > 0), which no case sets; a run with `LMDOT` and `DMDT0` > 0 would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/getnewenv.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/getnewenv.f#L217-L218" ] }, { "file": "getopac.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "A stale temperature in the pure-$Z$ trap's message", "runs": "Every case sets `LPUREZ` on; `LPUREZ` off, with a helium-burning core whose Z passes 0.1 above log T = 7.0, would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/getopac.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/getopac.f#L104-L114" ] }, { "file": "getopac.f", "number": 2, "result": "not reached", "status": "unconfirmed", "title": "The single-table test now only changes the stop", "runs": "Every case sets `LOPAL95` on; `LOPAL95` off would reach it, and stops the run at the first interior point either way.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/getopac.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/getopac.f#L140-L141", "https://github.com/yreclab/yrec/blob/7be4133/src/getopac.f#L162-L166" ] }, { "file": "getopal95.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "LSHIFT and LEXTH left over from the last call", "runs": "No case or document holds a hydrogen-poor surface, which it needs.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/getopal95.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/getopal95.f#L167-L176", "https://github.com/yreclab/yrec/blob/7be4133/src/getopal95.f#L214" ] }, { "file": "getopal95.f", "number": 2, "result": "changes results, large", "status": "checked by KC (Nov 10, 2026)", "title": "Empty temperatures in the OP tables above $\\log T = 8.0$", "runs": "In a made-up case, `examples/run_from_zahb_to_tahb/ZAHB_m1p0feh+0p0_TAHB` with `FLIV95` set to `op05/AAG2021zx0.0187.tab`, log L moves by up to 0.27 dex, against 1.4e-6 for the null change.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/getopal95.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/getopal95.f#L113-L131", "https://github.com/yreclab/yrec/blob/7be4133/src/getopal95.f#L177-L209", "https://github.com/yreclab/yrec/blob/7be4133/src/op952d.f#L39-L67" ] }, { "file": "getopal95.f", "number": 3, "result": "never occurred", "status": "unconfirmed", "title": "The wrong array in the empty-corner search", "runs": "The search ran in five runs, up to 113 million times, but `CAPPA2` never gave a different answer from `CAPPA`.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/getopal95.f.html#b3", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/getopal95.f#L193-L199" ] }, { "file": "getqua.f", "number": 1, "result": "changes results, large", "status": "checked by KC (Nov 10, 2026)", "title": "A wrong index in $d\\rho/dr$", "runs": "In `testsuite/Test_3Msolar_GS_rot_gray_OPALSCV_SF2_GS98_OP_CF10_RGp2diffrotI3cz2`, the fix moves log L by up to 0.42 dex, against 0.0076 for the null change, with 1,919 models instead of 1,853.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/getqua.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/getqua.f#L76-L77", "https://github.com/yreclab/yrec/blob/7be4133/src/getqua.f#L82-L84" ] }, { "file": "getqua.f", "number": 2, "result": "never occurred", "status": "unconfirmed", "title": "The central $d\\rho/dr$ is divided by the wrong interval", "runs": "Never occurred: in 2,252 calls in the 0.8 M☉ rotating case, r1/(r2-r1) was 16.0 to 17.9, so DR was r2-r1 and the central d rho/dr correct; the fix computes the same, and was not run.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/getqua.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/getqua.f#L54-L56", "https://github.com/yreclab/yrec/blob/7be4133/src/getqua.f#L61", "https://github.com/yreclab/yrec/blob/7be4133/src/getqua.f#L72-L77" ] }, { "file": "getrot.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "A one-shell convective surface is never treated", "runs": "A convection zone started in all five rotating test cases, but never at the last shell; it needs a rotating model whose surface convection zone is a single shell.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/getrot.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/getrot.f#L33-L37", "https://github.com/yreclab/yrec/blob/7be4133/src/getrot.f#L51-L58" ] }, { "file": "getw.f", "number": 1, "result": "changes results, small", "status": "checked by KC (Nov 10, 2026)", "title": "THP and WMST recomputed from a stale OMEGA after the first substep", "runs": "In `testsuite/Test_m0p8feh+0p0GS98_rotp8p0d0p1mag1e4std_TAMS`, the fix moves central log T by up to 0.0015 dex, against 0.0002 for the null change; the solar case stays within the noise.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/getw.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/getw.f#L280-L290", "https://github.com/yreclab/yrec/blob/7be4133/src/getw.f#L465-L467" ] }, { "file": "getyrec7.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "The flags record is read one column late", "runs": "No run started from a YREC7 file; a YREC7 file with a four-digit `JENV` would reach the misreading.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/getyrec7.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/getyrec7.f#L104-L105", "https://github.com/yreclab/yrec/blob/7be4133/src/putyrec7.f#L65-L66" ] }, { "file": "grsett.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "The coefficient correction is one midpoint off", "runs": "These lines belong to the old settling routines; `LDIFY` or `LDIFZ` with `LNEWDIF` false and `LDIFLI` off would reach them.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/grsett.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/grsett.f#L213-L215", "https://github.com/yreclab/yrec/blob/7be4133/src/get_imp_diffco.f#L19-L21" ] }, { "file": "grsett.f", "number": 2, "result": "not reached", "status": "unconfirmed", "title": "The convergence test ignores negative changes", "runs": "These lines belong to the old settling routines; `LDIFY` or `LDIFZ` with `LNEWDIF` false and `LDIFLI` off would reach them.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/grsett.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/grsett.f#L233-L246" ] }, { "file": "gtpurz.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "A density beyond the table's margin is reported as a temperature", "runs": "No run stopped here; a density beyond the table's margin at one of the four temperatures would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/gtpurz.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/gtpurz.f#L81-L105", "https://github.com/yreclab/yrec/blob/7be4133/src/gtpurz.f#L120-L126" ] }, { "file": "gtpurz.f", "number": 2, "result": "never occurred", "status": "unconfirmed", "title": "The extrapolation margin is a natural log", "runs": "GTPURZ's density test ran 57 million times in `testsuite/Test_m0p95feh+0p0_GN93_TAHB` and 54 million in `testsuite/Test_m9p0feh+0p0GS98_base_TAHB`; the density stayed at least 0.48 dex inside the table, so neither margin was used.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/gtpurz.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/gtpurz.f#L66", "https://github.com/yreclab/yrec/blob/7be4133/src/gtpurz.f#L87", "https://github.com/yreclab/yrec/blob/7be4133/src/gtpurz.f#L95-L96" ] }, { "file": "hpoint.f", "number": 1, "result": "harmless in the runs", "status": "unconfirmed", "title": "The start-of-step energy generation is interpolated against the new mass points", "runs": "In `testsuite/Test_m0p8feh+0p0GS98_rotp8p0d0p1mag1e4std_TAMS`, the fix moves log L by up to 6.4e-5 dex, against 1.4e-4 for the null change; the rotating 3 M☉ giant's track is identical.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/hpoint.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/hpoint.f#L602", "https://github.com/yreclab/yrec/blob/7be4133/src/hpoint.f#L733-L740", "https://github.com/yreclab/yrec/blob/7be4133/src/hpoint.f#L751-L758" ] }, { "file": "hpoint.f", "number": 2, "result": "never occurred", "status": "unconfirmed", "title": "The flag array can overflow", "runs": "The base was flagged in 26 runs, up to 7,196 times, with at most 21 flags set before it, far from the 98 needed; a run with `LOVSTE` and about a hundred jumps would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/hpoint.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/hpoint.f#L73", "https://github.com/yreclab/yrec/blob/7be4133/src/hpoint.f#L223-L227", "https://github.com/yreclab/yrec/blob/7be4133/src/hpoint.f#L264-L267", "https://github.com/yreclab/yrec/blob/7be4133/src/hpoint.f#L287-L290", "https://github.com/yreclab/yrec/blob/7be4133/src/hpoint.f#L294" ] }, { "file": "hpoint.f", "number": 3, "result": "changes results, large", "status": "checked by KC (Nov 10, 2026)", "title": "ALPHAE taken as a depth in $\\log_{10}P$", "runs": "In `examples/run_from_zahb_to_tahb/ZAHB_m0p8feh+0p0_TAHB`, an approximate version of the fix (`ALPHAE/CLN`) moves log L by up to 0.334 dex, against 2.0e-7 for the null change, and the run ends at 0.2356 Gyr instead of 0.2297.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/hpoint.f.html#b3", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/hpoint.f#L252-L255", "https://github.com/yreclab/yrec/blob/7be4133/src/main.f#L803-L804" ] }, { "file": "hpoint.f", "number": 4, "result": "not reached", "status": "unconfirmed", "title": "The debug table prints new values under \"old points.\"", "runs": "It needs `LDEBUG` with `NPOINT` below 9999, which no case sets.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/hpoint.f.html#b4", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/hpoint.f#L436-L441", "https://github.com/yreclab/yrec/blob/7be4133/src/hpoint.f#L493-L499", "https://github.com/yreclab/yrec/blob/7be4133/src/hpoint.f#L648-L655", "https://github.com/yreclab/yrec/blob/7be4133/src/hpoint.f#L713-L720" ] }, { "file": "hpoint.f", "number": 5, "result": "never occurred", "status": "unconfirmed", "title": "A third equal flag survives", "runs": "A duplicate flag was removed in 6 runs, up to 194 times, but a third equal flag never followed.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/hpoint.f.html#b5", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/hpoint.f#L309-L320" ] }, { "file": "htimer.f", "number": 1, "result": "harmless in the runs", "status": "unconfirmed", "title": "A zero previous step passes through the envelope-triangle criterion", "runs": "In a made-up case, `testsuite/Test_m1p0feh+0p0GS98_base_HeIgnite` with `LMDOT` and `LREIMER` on, the step fell below 1 s 919 times; the run stops the same way with or without the fix, with identical outputs.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/htimer.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/htimer.f#L158-L162", "https://github.com/yreclab/yrec/blob/7be4133/src/htimer.f#L171-L182", "https://github.com/yreclab/yrec/blob/7be4133/src/entime.f#L92-L95" ] }, { "file": "iounits_mod.f", "number": 1, "result": "cannot be reached", "status": "unconfirmed", "title": "IJLAST is never assigned", "runs": "`LJLAST` is set only in the main program's `LACOUT` branch, and `LACOUT` is in no namelist group and set nowhere in `src/`.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/iounits_mod.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/iounits_mod.f#L42", "https://github.com/yreclab/yrec/blob/7be4133/src/wrtout.f#L710-L715", "https://github.com/yreclab/yrec/blob/7be4133/src/parmin.f#L847-L849" ] }, { "file": "kcsurfp.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "Missing pressures at the highest gravities", "runs": "In a made-up case, `testsuite/Test_solarGS98_kurucz` with `KTTAU` = 5 and `atmk2004p00.tab`, the branch ran 10,156 times and never met a row missing its highest gravities.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/kcsurfp.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/kcsurfp.f#L97-L111", "https://github.com/yreclab/yrec/blob/7be4133/src/kcsurfp.f#L112-L132" ] }, { "file": "kemcom.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "The error message gives a fixed limit of 50", "runs": "The network always converged within `KEMMAX` = 50, the value in every case; the message is wrong only with another `KEMMAX`, and then only when the network fails.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/kemcom.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/kemcom.f#L377-L385" ] }, { "file": "kemcom.f", "number": 2, "result": "never occurred", "status": "unconfirmed", "title": "¹⁵N and ¹⁷O are not mixed in burning zones", "runs": "Never occurred: in 11,559 passes of burning convection zones over five cases, every shell's 15N and 17O already equaled the zone's mean within a relative 1e-10, so the fix was not run.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/kemcom.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/kemcom.f#L400-L410", "https://github.com/yreclab/yrec/blob/7be4133/src/kemcom.f#L121-L127" ] }, { "file": "liburn.f", "number": 1, "result": "harmless in the runs", "status": "unconfirmed", "title": "The sign of FRAD", "runs": "`FRAD` lay outside 0 to 1 in 23 runs; with the fix, `testsuite/Test_solarM22M_difkur` keeps its compared track columns unchanged, where the null change moves log L by up to 0.000167 dex.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/liburn.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/liburn.f#L509-L528" ] }, { "file": "liburn.f", "number": 2, "result": "harmless in the runs", "status": "unconfirmed", "title": "Points newly in a deepened zone are not burned", "runs": "The zone deepened in 23 runs; with the fix, `testsuite/Test_3Msolar_GS_rot_gray_OPALSCV_SF2_GS98_OP_CF10_RGp2diffrotI3cz2` keeps its compared track columns unchanged, where the null change moves log L by up to 0.00758 dex.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/liburn.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/liburn.f#L191-L193", "https://github.com/yreclab/yrec/blob/7be4133/src/liburn.f#L373-L379", "https://github.com/yreclab/yrec/blob/7be4133/src/liburn.f#L496-L500" ] }, { "file": "liburn.f", "number": 3, "result": "never occurred", "status": "unconfirmed", "title": "⁹Be is not guarded against overflow", "runs": "The division ran in 23 runs, up to 11,399 times, with `DBE9` at most 201, far below 709.78, so `EXP` never overflowed.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/liburn.f.html#b3", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/liburn.f#L529-L538" ] }, { "file": "liburn.f", "number": 4, "result": "harmless in the runs", "status": "unconfirmed", "title": "The bottom point's rates are left adjusted", "runs": "The rates were left adjusted in 24 runs; with the fix, `testsuite/Test_m0p8feh+0p0GS98_rotp8p0d0p1mag1e4std_TAMS` keeps an identical track and `.store` file.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/liburn.f.html#b4", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/liburn.f#L334-L343", "https://github.com/yreclab/yrec/blob/7be4133/src/liburn.f#L371-L372" ] }, { "file": "liburn2.f", "number": 1, "result": "cannot be reached", "status": "unconfirmed", "title": "The sign of FRAD", "runs": "LIBURN2 is called only from GETW and BURSMIX, both only when `LBURS` is true, and GETW sets `LBURS` false on every call, which no namelist changes.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/liburn2.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/liburn2.f#L319-L339" ] }, { "file": "liburn2.f", "number": 2, "result": "cannot be reached", "status": "unconfirmed", "title": "Points newly in a deepened zone are not burned", "runs": "LIBURN2 is called only from GETW and BURSMIX, both only when `LBURS` is true, and GETW sets `LBURS` false on every call, which no namelist changes.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/liburn2.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/liburn2.f#L155-L157", "https://github.com/yreclab/yrec/blob/7be4133/src/liburn2.f#L274-L278", "https://github.com/yreclab/yrec/blob/7be4133/src/liburn2.f#L308-L312" ] }, { "file": "liburn2.f", "number": 3, "result": "cannot be reached", "status": "unconfirmed", "title": "The overshoot zone's width mixes two models", "runs": "LIBURN2 is called only from GETW and BURSMIX, both only when `LBURS` is true, and GETW sets `LBURS` false on every call, which no namelist changes.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/liburn2.f.html#b3", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/liburn2.f#L135-L142" ] }, { "file": "liburn2.f", "number": 4, "result": "cannot be reached", "status": "unconfirmed", "title": "HRO(0) read for a fully convective start", "runs": "LIBURN2 is called only from GETW and BURSMIX, both only when `LBURS` is true, and GETW sets `LBURS` false on every call, which no namelist changes.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/liburn2.f.html#b4", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/liburn2.f#L113-L115" ] }, { "file": "liburn2.f", "number": 5, "result": "cannot be reached", "status": "unconfirmed", "title": "No overflow guard in the last loop", "runs": "LIBURN2 is called only from GETW and BURSMIX, both only when `LBURS` is true, and GETW sets `LBURS` false on every call, which no namelist changes.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/liburn2.f.html#b5", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/liburn2.f#L340-L346" ] }, { "file": "liburn2.f", "number": 6, "result": "cannot be reached", "status": "unconfirmed", "title": "The bottom point's rates are left adjusted", "runs": "LIBURN2 is called only from GETW and BURSMIX, both only when `LBURS` is true, and GETW sets `LBURS` false on every call, which no namelist changes.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/liburn2.f.html#b6", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/liburn2.f#L233-L242", "https://github.com/yreclab/yrec/blob/7be4133/src/liburn2.f#L270-L271" ] }, { "file": "ll95tbl.f", "number": 1, "result": "changes results, large", "status": "checked by KC (Nov 10, 2026)", "title": "The $X = 1 - Z$ tables at $Z$ = 0.06 and 0.08 overwrite the $X = 0$ tables at $Z$ = 0.08 and 0.10", "runs": "In `testsuite/Test_m0p95feh+0p0_GN93_TAHB`, the fix moves log L by up to 0.36 dex, against 1.7e-7 for the null change, and the track ends at 0.0943 Gyr instead of 0.1293.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/ll95tbl.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/ll95tbl.f#L117", "https://github.com/yreclab/yrec/blob/7be4133/src/ll95tbl.f#L176-L191", "https://github.com/yreclab/yrec/blob/7be4133/src/llot95a_mod.f#L46-L47" ] }, { "file": "main.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "The stop conditions are ignored on rescale-and-evolve runs", "runs": "In `sample_cases/UpperMS/Sample_m9p0feh+0p0GS98_none_TAHB`, the stop test was skipped at both models, but central Y stayed above its target, so no stop was missed.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/main.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/parmin.f#L1396-L1405", "https://github.com/yreclab/yrec/blob/7be4133/src/parmin.f#L1463", "https://github.com/yreclab/yrec/blob/7be4133/src/main.f#L1087-L1108" ] }, { "file": "massloss.f", "number": 1, "result": "changes results, large", "status": "checked by KC (Nov 10, 2026)", "title": "The Reimers rate is in g/s, read as $M_\\odot$/yr", "runs": "A made-up case: `testsuite/Test_m3p0feh+0p0GS98_std25_TAMS` with `LMDOT` and `LREIMER` on; unfixed, the 3 M☉ star falls to 0.11 M☉ and fails to converge, and with the fix it keeps its mass.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/massloss.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/massloss.f#L100-L103", "https://github.com/yreclab/yrec/blob/7be4133/src/mdot.f#L111", "https://github.com/yreclab/yrec/blob/7be4133/src/mdot.f#L159" ] }, { "file": "mdot.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "Accretion past the disk's lifetime gives a negative step", "runs": "No run had `LMDOT` with `LDISK` and accretion; it needs `LMDOT`, `LDISK`, and `DMDT0` > 0, from a starting model older than `TDISK`.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/mdot.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/mdot.f#L146-L157", "https://github.com/yreclab/yrec/blob/7be4133/src/main.f#L981" ] }, { "file": "mdot.f", "number": 2, "result": "stops the run", "status": "checked by KC (Nov 10, 2026)", "title": "A stop without a message when there is no surface convection zone", "runs": "A made-up case: `testsuite/Test_m3p0feh+0p0GS98_std25_TAMS` with Reimers mass loss on and `massloss.f` bug 1 fixed ends here at model 129, at `STOP 911` with no message; the 1 M☉ giant-branch case does too, at model 6,599.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/mdot.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/mdot.f#L118", "https://github.com/yreclab/yrec/blob/7be4133/src/massloss.f#L106-L122" ] }, { "file": "meval.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "No point of evaluation within the data", "runs": "In all 33 runs that got past their start, label 80 was reached up to 342,520 times, always with `START` = 1; a regridding with no point inside the old grid would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/meval.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/meval.f#L107-L121", "https://github.com/yreclab/yrec/blob/7be4133/src/meval.f#L138-L144", "https://github.com/yreclab/yrec/blob/7be4133/src/meval.f#L242-L243" ] }, { "file": "microdiff_cod.f", "number": 1, "result": "cannot be reached", "status": "unconfirmed", "title": "Without hydrogen settling, the hydrogen gradient is never computed", "runs": "`LDIFZ` and `LDIFLI` force `LDIFY` on, and `MIX` and `ROTMIX` call `MICRODIFF` only while `LDIFY` is true, so the hydrogen call always runs first, as in all seven runs with settling.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/microdiff_cod.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/microdiff_cod.f#L96", "https://github.com/yreclab/yrec/blob/7be4133/src/microdiff_cod.f#L191-L212" ] }, { "file": "microdiff_cod.f", "number": 2, "result": "harmless in the runs", "status": "unconfirmed", "title": "Skipped points keep old terms", "runs": "Points were skipped up to 21.9 million times per run, as in `testsuite/Test_solarM22M_difkur`, but the stored terms were finite and the abundance there zero, so they were harmless; the fix was not run.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/microdiff_cod.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/microdiff_cod.f#L100-L105", "https://github.com/yreclab/yrec/blob/7be4133/src/microdiff_cod.f#L163-L174" ] }, { "file": "microdiff_etm.f", "number": 1, "result": "harmless in the runs", "status": "unconfirmed", "title": "Light elements in the interior and the envelope are not floored at zero", "runs": "Light elements went negative by at most 3.8e-13; with the fix, `testsuite/Test_m1p4feh+0p0GS98_dif_TAMS` and `testsuite/Test_solarM22M_difkur` give identical tracks, since `MIX` floors them at zero.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/microdiff_etm.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/microdiff_etm.f#L66-L74", "https://github.com/yreclab/yrec/blob/7be4133/src/microdiff_etm.f#L122-L131", "https://github.com/yreclab/yrec/blob/7be4133/src/microdiff_etm.f#L156-L163" ] }, { "file": "microdiff_etm.f", "number": 2, "result": "never occurred", "status": "unconfirmed", "title": "Division by zero where $Z = 0$", "runs": "The line ran in all seven runs with settling, but Z was never below 0.0020; it needs `LDIFZ` and a model with Z = 0 somewhere.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/microdiff_etm.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/microdiff_etm.f#L53", "https://github.com/yreclab/yrec/blob/7be4133/src/microdiff_etm.f#L112", "https://github.com/yreclab/yrec/blob/7be4133/src/microdiff_etm.f#L143" ] }, { "file": "microdiff_mte.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "The first midpoint can be interpolated from the wrong model points", "runs": "The search ran in all seven runs with settling, such as `testsuite/Test_solarGS98_dif`, and never failed; it needs a first point beyond `NPT`, as for a large convective core under few points.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/microdiff_mte.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/microdiff_mte.f#L100-L113" ] }, { "file": "microdiff_mte.f", "number": 2, "result": "harmless in the runs", "status": "unconfirmed", "title": "The grid is coarser than the comment says", "runs": "In the rotating solar model with settling, DR exceeded DRMIN at all 3,702 grid setups (by up to 1.020); the fix's effect stays within the null change (log L 0.00021 against 0.00029 dex).", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/microdiff_mte.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/microdiff_mte.f#L61-L64", "https://github.com/yreclab/yrec/blob/7be4133/src/microdiff_mte.f#L83" ] }, { "file": "microdiff_run.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "With NITER_GS = 0, the implicit term is skipped and the warning prints old values", "runs": "The loop ran in all seven runs with settling, such as `testsuite/Test_solarGS98_dif`, always with `NITER_GS` = 10; every case sets 10, and it needs `NITER_GS` = 0.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/microdiff_run.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/microdiff_run.f#L152", "https://github.com/yreclab/yrec/blob/7be4133/src/microdiff_run.f#L178-L196" ] }, { "file": "microdiff_run.f", "number": 2, "result": "harmless in the runs", "status": "unconfirmed", "title": "The convergence test is one-sided", "runs": "The loop ended one iteration early 14 to 259 times per run; with the fix, `testsuite/Test_m1p4feh+0p0GS98_dif_TAMS` and `testsuite/Test_solarM22M_difkur` give identical tracks and stored models.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/microdiff_run.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/microdiff_run.f#L178-L190" ] }, { "file": "microdiff_run.f", "number": 3, "result": "cannot be reached", "status": "unconfirmed", "title": "The midpoint change is shifted by one point", "runs": "`EQCOD2X_H` is set to zero just before every use, and it was zero every time in all seven runs with settling.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/microdiff_run.f.html#b3", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/microdiff_run.f#L155-L157" ] }, { "file": "microdiff_setup.f", "number": 1, "result": "cannot be reached", "status": "unconfirmed", "title": "No check that IEND > IBEGIN", "runs": "`MIX` and `ROTMIX` switch settling off once central X falls below `XMIN`, so the searches give `IEND` at least `IBEGIN` + 2; only Y below `YMIN` outside the core could move `IEND` lower.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/microdiff_setup.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/microdiff_setup.f#L107-L116", "https://github.com/yreclab/yrec/blob/7be4133/src/microdiff_setup.f#L130-L139" ] }, { "file": "mixcz.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "The sentinel LCZ(M+1) is written past the model", "runs": "In the 33 runs that got past their start, models had 349 to 3,455 shells, never 5,000, so the sentinel never fell outside the array; it needs a model of 5,000 shells.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/mixcz.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/mixcz.f#L65-L67" ] }, { "file": "mixcz.f", "number": 2, "result": "never occurred", "status": "unconfirmed", "title": "Helium is not recomputed after the clean-up", "runs": "Never occurred: in 12,145 calls over five cases neither clamp acted; the fix, which also resets Y everywhere, changes results by about the size of the null change (log L 0.071 against 0.057 dex in the 1.4 M☉ case).", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/mixcz.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/mixcz.f#L117-L122", "https://github.com/yreclab/yrec/blob/7be4133/src/mix.f#L473-L474" ] }, { "file": "model_to_equal.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "The first midpoint's search runs over the wrong indices", "runs": "These are the old settling routines, which no case uses; reaching them takes `LDIFY` or `LDIFZ` with `LNEWDIF` false and `LDIFLI` off, and a model whose `IBEGIN` exceeds `NPT`.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/model_to_equal.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/model_to_equal.f#L91-L97" ] }, { "file": "mwind.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "WNEW overwrites the namelist's solid-body rate", "runs": "The overwrite happened at every step in `testsuite/Test_m0p8feh+0p0GS98_rotp8p0d0p1mag1e4std_TAMS`, but no later run read `WNEW`; a later run with `LFIRST(k)` true, k > 1, `LROT`, `LWNEW`, and `LMWIND` would.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/mwind.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/mwind.f#L186", "https://github.com/yreclab/yrec/blob/7be4133/src/mwind.f#L199", "https://github.com/yreclab/yrec/blob/7be4133/src/mwind.f#L208" ] }, { "file": "neutrino.f", "number": 1, "result": "changes results, small", "status": "checked by KC (Nov 10, 2026)", "title": "Metals other than C and O are left out", "runs": "In the 1 M☉ run to helium ignition the fractions summed to 0.984 to 0.991 in 33,460,185 calls; the fix changes the last neutrino luminosity from 0.845 to 0.876 L☉ (null change: 0.846), with log L within the noise and central log T a little beyond it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/neutrino.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/neutrino.f#L52-L58", "https://github.com/yreclab/yrec/blob/7be4133/src/nulosses.f#L87-L95" ] }, { "file": "op95xtab.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "Table 9 at $Z = 0.04$", "runs": "The high-X branch ran 350 times in `testsuite/Test_m1p4feh+0p0GS98_dif_TAMS`, at a table Z of 0.0188, never at Z = 0.04; `ZOPAL951` = 0.04 with a surface X above 0.8 would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/op95xtab.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/op95xtab.f#L53-L65", "https://github.com/yreclab/yrec/blob/7be4133/src/op95ztab.f#L68-L70" ] }, { "file": "op95xtab.f", "number": 2, "result": "harmless in the runs", "status": "unconfirmed", "title": "The last $X$ node taken as 1 instead of $1 - Z$", "runs": "In `testsuite/Test_m1p4feh+0p0GS98_dif_TAMS`, the misplaced node was in the stencil 203 times; with the fix, the track and the `.store` file are identical.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/op95xtab.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/op95xtab.f#L63-L65", "https://github.com/yreclab/yrec/blob/7be4133/src/op95xtab.f#L76-L79" ] }, { "file": "op95ztab.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "The run stops for $0.08 < Z \\le 0.1$", "runs": "No case sets `ZOPAL951` above 0.08; `ZOPAL951` between 0.08 and 0.1 would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/op95ztab.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/op95ztab.f#L34-L46" ] }, { "file": "op95ztab.f", "number": 2, "result": "cannot be reached", "status": "unconfirmed", "title": "Label 30 in the wrong place", "runs": "Any Z above 0.08, and so Z of 0.1 or more, stops the run before label 30, since the search for the four Z values runs only to `NUMZ-1` (bug 1).", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/op95ztab.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/op95ztab.f#L97-L124" ] }, { "file": "osplin.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "MEVAL's error code is ignored", "runs": "In all 33 runs that got past their start, MEVAL never returned `ERR` = 2; points out of order would reach it, and no caller gave them.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/osplin.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/osplin.f#L47", "https://github.com/yreclab/yrec/blob/7be4133/src/meval.f#L93-L98" ] }, { "file": "parmin.f", "number": 1, "result": "wrong output", "status": "checked by KC (Nov 10, 2026)", "title": "The wrong message for a rescale-and-evolve card", "runs": "In `examples/run_from_seed_to_start/m10000gs98z001000_prems`, the wrong message was printed once; with format 551 restored, the track and `.store` file are identical, and only the `.short` file changes.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/parmin.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/parmin.f#L1476-L1480" ] }, { "file": "parmin.f", "number": 2, "result": "changes results, small", "status": "checked by KC (Nov 10, 2026)", "title": "ISETMIX = 2 and ISETISO = 2 do nothing", "runs": "A made-up case: `testsuite/Test_solarAAG21_dif` with `ISETMIX` = `ISETISO` = 2 and its custom mixture on; the fix moves log L by up to 0.0025 dex on the pre-main sequence, against 0.00023 for the null change.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/parmin.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/parmin.f#L1141", "https://github.com/yreclab/yrec/blob/7be4133/src/parmin.f#L1149", "https://github.com/yreclab/yrec/blob/7be4133/src/parmin.f#L1193", "https://github.com/yreclab/yrec/blob/7be4133/src/parmin.f#L1196" ] }, { "file": "parmin.f", "number": 3, "result": "harmless in the runs", "status": "unconfirmed", "title": "ZENV set from every run card", "runs": "The line set `ZENV` = 0 in 32 runs; with it dropped, `testsuite/Test_m0p3feh+0p0GS98_scv_ZAMS` and `examples/run_standard_solar_model/Test_solarcal_noGS_norot` give identical tracks and `.store` files.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/parmin.f.html#b3", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/parmin.f#L1486" ] }, { "file": "parmin.f", "number": 4, "result": "never occurred", "status": "unconfirmed", "title": "No / in FTRACK leaves LAST_SLASH_IDX unset", "runs": "Every `FTRACK` in the 33 runs named a folder; four namelists in `sample_cases/Gyrochrones`, such as `rotevolin_m0p2per1p4`, name none and were not run.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/parmin.f.html#b4", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/parmin.f#L815-L828" ] }, { "file": "parmin.f", "number": 5, "result": "withdrawn", "status": "unconfirmed", "title": "EXPAND_VALUE leaves the end of TEMP undefined", "runs": "Every run set all three variables; the withdrawal rests on gfortran's behavior, and a compiler that leaves the argument unchanged for an unset variable would bring the item back.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/parmin.f.html#b5", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/parmin.f#L1527", "https://github.com/yreclab/yrec/blob/7be4133/src/parmin.f#L1533-L1538" ] }, { "file": "parmin.f", "number": 6, "result": "wrong output", "status": "checked by KC (Nov 10, 2026)", "title": "Long paths are cut to 256 characters without a word", "runs": "Checked by reading the code: a 201-character output folder cuts a case's names to 256 characters without a message; in the round of Nov 9, two fix runs lost their files this way.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/parmin.f.html#b6", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/parmin.f#L1527-L1538", "https://github.com/yreclab/yrec/blob/7be4133/src/parmin.f#L825-L828", "https://github.com/yreclab/yrec/blob/7be4133/src/pdist.f#L86-L111" ] }, { "file": "pdist.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "Two different lines under the same labels", "runs": "PDIST runs only with `LPOUT` on, which no case sets; any run with `LPOUT` on reaches this line at every model.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/pdist.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/pdist.f#L119-L122" ] }, { "file": "pindex.f", "number": 1, "result": "cannot be reached", "status": "unconfirmed", "title": "IEND undefined", "runs": "PINDEX's only caller, PUTSTORE, passes its own `LSHELL`, which nothing sets (PUTSTORE's bug 3); it was false at every call.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/pindex.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/pindex.f#L48-L56" ] }, { "file": "pindex.f", "number": 2, "result": "cannot be reached", "status": "unconfirmed", "title": "M listed twice", "runs": "PINDEX's only caller, PUTSTORE, passes its own `LSHELL`, which nothing sets (PUTSTORE's bug 3); it was false at every call.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/pindex.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/pindex.f#L65-L72", "https://github.com/yreclab/yrec/blob/7be4133/src/pindex.f#L83" ] }, { "file": "putmodel2.f", "number": 1, "result": "wrong output", "status": "checked by KC (Nov 10, 2026)", "title": "The caller's luminosities are left in erg s⁻¹", "runs": "Luminosities in erg/s reached HTIMER and the `ISNU` fractions in 27 runs; with the fix, `testsuite/Test_m0p95feh+0p0_GN93_TAHB` gives an identical track, and only its `.short` file differs.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/putmodel2.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/putmodel2.f#L90-L97" ] }, { "file": "putstore.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "The SCV + OPAL 2006 label is lost", "runs": "The `LSCV` branch ran in 30 runs, but YREC's cases with SCV and OPAL 2006 all set `LDH`; any of them with `LDH` off would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/putstore.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/putstore.f#L149-L157" ] }, { "file": "putstore.f", "number": 2, "result": "wrong output", "status": "checked by KC (Nov 10, 2026)", "title": "Atmosphere 5 has no label", "runs": "A made-up case: `testsuite/Test_solarGS98_kurucz` with `KTTAU` = 5 and the Castelli-Kurucz table writes `ATM` without its own value at both stores; with the fix, the track is identical.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/putstore.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/putstore.f#L135-L145" ] }, { "file": "putstore.f", "number": 3, "result": "never occurred", "status": "checked by KC (Nov 10, 2026)", "title": "PINDEX gets arguments that are never set", "runs": "In 10 runs the shell's flag was set at a store, but no shell went unlisted: every old-format case keeps NPRTPT = 1, which lists every shell; the fix's .store change comes from PINDEX's own bug 2 (M listed twice), which the fix reaches.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/putstore.f.html#b3", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/putstore.f#L51-L52", "https://github.com/yreclab/yrec/blob/7be4133/src/putstore.f#L262" ] }, { "file": "putstore.f", "number": 4, "result": "wrong output", "status": "checked by KC (Nov 10, 2026)", "title": "The old header calls the triple-α column \"3HE.\"", "runs": "The header was written in 14 runs; with the fix, `examples/run_from_dbl_to_zams/Test_m0040_feh0_noGS_ZAMS` gives an identical track and a `.store` file within 2%, as expected for a label.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/putstore.f.html#b4", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/putstore.f#L275", "https://github.com/yreclab/yrec/blob/7be4133/src/putstore.f#L309" ] }, { "file": "putyrec7.f", "number": 1, "result": "cannot be reached", "status": "unconfirmed", "title": "The flags record does not match its reader", "runs": "Nothing in `src/` calls PUTYREC7.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/putyrec7.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/putyrec7.f#L65-L66", "https://github.com/yreclab/yrec/blob/7be4133/src/getyrec7.f#L104-L105" ] }, { "file": "qenv.f", "number": 1, "result": "harmless in the runs", "status": "unconfirmed", "title": "The distortion factor in $d\\log r/d\\log P$", "runs": "In `testsuite/Test_m0p8feh+0p0GS98_rotp8p0d0p1mag1e4std_TAMS`, the fix moves log Teff by up to 0.000105, against 0.0000272 for the null change, while log L stays at the noise: not clearly above it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/qenv.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/qenv.f#L72", "https://github.com/yreclab/yrec/blob/7be4133/src/qenv.f#L74" ] }, { "file": "readalex06.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "A check that does not stop", "runs": "The check ran 128 times in each of the 33 runs and found no mismatch in `alexmol06gs98.tab`; a table file in another order or with other compositions would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/readalex06.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/readalex06.f#L63-L70" ] }, { "file": "readalex06.f", "number": 2, "result": "never occurred", "status": "unconfirmed", "title": "No STATUS='OLD'", "runs": "The `OPEN` ran once in each of the 33 runs, and the file always existed; a wrong `FALEX06` would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/readalex06.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/readalex06.f#L54" ] }, { "file": "readcoeos06.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "The 2006 tables called \"OEOS96\" in a message", "runs": "In the 31 runs with the OPAL 2006 tables, every density row was long enough; only a damaged or different file reaches the message, and the run then stops.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/readcoeos06.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/readcoeos06.f#L114-L118" ] }, { "file": "readcoeos06.f", "number": 2, "result": "harmless in the runs", "status": "unconfirmed", "title": "ALOGNE read into one element", "runs": "Each of the 31 runs with the OPAL 2006 tables put 116,595 of 117,440 rows in the wrong element, but `ALOGNE` is read nowhere, so no output can change; the fix was not run.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/readcoeos06.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/readcoeos06.f#L124-L126", "https://github.com/yreclab/yrec/blob/7be4133/src/readcoeos06.f#L81", "https://github.com/yreclab/yrec/blob/7be4133/src/readcoeos06.f#L89" ] }, { "file": "rotmix.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "The outer point for the settling step is not moved down", "runs": "In `testsuite/Test_m0p8feh+0p0GS98_rotp8p0d0p1mag1e4std_TAMS` and `testsuite/Test_solarGS98_difrotmixper18p4d1p0`, `J` always equaled `JMAX`; it matters only where Y at the top of the last radiative zone has fallen to `YMIN`.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/rotmix.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/rotmix.f#L156-L160", "https://github.com/yreclab/yrec/blob/7be4133/src/rotmix.f#L174-L180", "https://github.com/yreclab/yrec/blob/7be4133/src/rotmix.f#L191-L194" ] }, { "file": "rotmix.f", "number": 2, "result": "not reached", "status": "unconfirmed", "title": "The stop message names the wrong routine", "runs": "No run stopped here; a rotating run with `LDIFY` for a star without a surface convection zone would reach it, and stop.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/rotmix.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/rotmix.f#L163-L169" ] }, { "file": "rscale.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "The error message for one species names no species", "runs": "Every YREC case that sets `LNEWCP` has a target of 1.5D-5; a rescale card on a model with central hydrogen, with `LNEWCP` on and a target of 1 or more, would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/rscale.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/rscale.f#L175-L178", "https://github.com/yreclab/yrec/blob/7be4133/src/rscale.f#L38" ] }, { "file": "rscale.f", "number": 2, "result": "not reached", "status": "unconfirmed", "title": "The \"new method\" does not check the new total mass against the core", "runs": "The new method was not taken, the horizontal-branch cases included; a horizontal-branch rescale card with less than 0.01 M☉ outside the core and a target mass at or below the core's would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/rscale.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/rscale.f#L300-L305", "https://github.com/yreclab/yrec/blob/7be4133/src/rscale.f#L310-L311", "https://github.com/yreclab/yrec/blob/7be4133/src/rscale.f#L320-L322" ] }, { "file": "sconvec.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "A merger keeps the upper zone's top", "runs": "In `sample_cases/UpperMS/Sample_m9p0feh+0p0GS98_yaleeos_sc_TAMS`, the one run with semiconvection, no two zones merged; it needs `LSEMIC` and a second zone within the walk's reach.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/sconvec.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/sconvec.f#L278-L291" ] }, { "file": "sconvec.f", "number": 2, "result": "harmless in the runs", "status": "unconfirmed", "title": "The radius is not updated in the walk", "runs": "In `sample_cases/UpperMS/Sample_m9p0feh+0p0GS98_yaleeos_sc_TAMS`, the walk used a wrong radius once; with the fix, the track and `.store` file are identical.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/sconvec.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/sconvec.f#L144", "https://github.com/yreclab/yrec/blob/7be4133/src/sconvec.f#L227-L251" ] }, { "file": "setcal.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "The new stop criteria are not copied to the repeated runs", "runs": "Never occurred: SETCAL ran once, in YREC's one calibration case, with no END_DCEN, END_XCEN, or END_YCEN set for run 2 or 3, so copying them changes nothing; the fix was not run.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/setcal.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/setcal.f#L83-L100", "https://github.com/yreclab/yrec/blob/7be4133/src/nmlpar_mod.f#L261-L263", "https://github.com/yreclab/yrec/blob/7be4133/src/main.f#L1089-L1115" ] }, { "file": "setscv.f", "number": 1, "result": "harmless in the runs", "status": "unconfirmed", "title": "QSMIXT reads the entropy of mixing beyond the lower row's pressures", "runs": "In all 32 runs with `LSCV`, 46 of 3,720 table points read an unset entry, but no counted stencil read one; with the fix, `testsuite/Test_m0p3feh+0p0GS98_scv_ZAMS` gives identical `.track` and `.store` files.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/setscv.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/setscv.f#L150-L151" ] }, { "file": "setscv.f", "number": 2, "result": "harmless in the runs", "status": "unconfirmed", "title": "The same reads in QDTT and QAT, where a fix is undone", "runs": "In all 32 runs with `LSCV`, 46 of 3,720 points read unset entries; with the fix, `testsuite/Test_m0p3feh+0p0GS98_scv_ZAMS` gives identical files, since the two columns are read nowhere.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/setscv.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/setscv.f#L212-L213", "https://github.com/yreclab/yrec/blob/7be4133/src/setscv.f#L226-L231" ] }, { "file": "setup_grsett.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "The derivatives miss the mixing factors", "runs": "These are the old settling routines, which no case uses; reaching them takes `LDIFY` or `LDIFZ` with `LNEWDIF` false, `LDIFLI` off, `LTHOUL`, and `CSTMIXING` or `CSTDIFFMIX` other than 1.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/setup_grsett.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/setup_grsett.f#L288-L291", "https://github.com/yreclab/yrec/blob/7be4133/src/setup_grsett.f#L311", "https://github.com/yreclab/yrec/blob/7be4133/src/setup_grsett.f#L316" ] }, { "file": "setups.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "A wrong bound for the Kurucz/Castelli tables", "runs": "A made-up case: `testsuite/Test_solarGS98_kurucz` with `KTTAU` = 5 and the Castelli-Kurucz table; all 76 temperatures had a valid pressure, so the line never assigned `NG`.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/setups.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/setups.f#L355", "https://github.com/yreclab/yrec/blob/7be4133/src/setups.f#L59-L60" ] }, { "file": "shape.f", "number": 1, "result": "harmless in the runs", "status": "unconfirmed", "title": "The first shell's $\\eta_2$ is seven times too large", "runs": "In `testsuite/Test_3Msolar_GS_rot_gray_OPALSCV_SF2_GS98_OP_CF10_RGp2diffrotI3cz2`, the fix moves log L by up to 0.0012 dex, against 0.0076 for the null change; the faster 2 M☉ cases were not run with it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/shape.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/shape.f#L60", "https://github.com/yreclab/yrec/blob/7be4133/src/shape.f#L93" ] }, { "file": "shape.f", "number": 2, "result": "harmless in the runs", "status": "unconfirmed", "title": "The Newton test for $r_0$ compares with ACFPFT, not its cube root", "runs": "In `testsuite/Test_m2p0_feh+0p0_vzams150_IES1_TAMS`, 15,909 solves ran all 20 steps; with the fix, every track difference is at most 1.6e-8, against 1.1e-7 for the null change.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/shape.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/shape.f#L71", "https://github.com/yreclab/yrec/blob/7be4133/src/shape.f#L122", "https://github.com/yreclab/yrec/blob/7be4133/src/shape.f#L134" ] }, { "file": "simeqc.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "The error message goes to unit 5", "runs": "In the 30 runs with the pivot test, the largest pivot was never below 0.229, so the matrix was never singular; a singular system would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/simeqc.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/simeqc.f#L31-L33" ] }, { "file": "sneut.f", "number": 1, "result": "changes results, small", "status": "checked by KC (Nov 10, 2026)", "title": "A typo in tfac2", "runs": "In `testsuite/Test_m2p0_feh+0p0_vzams150_IES2_TAMS`, the fix moves log L by up to 7.95e-5 dex, against 1.58e-8 for the null change, with the same 546 models and end age.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/sneut.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/sneut.f#L124-L125" ] }, { "file": "solid.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "A zone without angular momentum divides by zero", "runs": "In all five rotating cases the zone's angular momentum was never zero; an initial rotation of zero (`LWNEW` with `WNEW` = 0), or a start model without rotation, would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/solid.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/solid.f#L49-L57", "https://github.com/yreclab/yrec/blob/7be4133/src/solid.f#L86-L88" ] }, { "file": "solid.f", "number": 2, "result": "harmless in the runs", "status": "unconfirmed", "title": "QIW in the Newton step", "runs": "In `testsuite/Test_m0p8feh+0p0GS98_rotp8p0d0p1mag1e4std_TAMS`, the fix moves log L by up to 0.00017 dex, against 0.00014 for the null change; in `testsuite/Test_m2p0_feh+0p0_vzams150_IES1_TAMS`, every difference equals the null change's.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/solid.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/solid.f#L91", "https://github.com/yreclab/yrec/blob/7be4133/src/momi.f#L84-L86", "https://github.com/yreclab/yrec/blob/7be4133/src/intmom.f#L28-L31" ] }, { "file": "starin.f", "number": 1, "result": "harmless in the runs", "status": "unconfirmed", "title": "The densities are not shifted in the core extension", "runs": "Wrong starting densities in 12 runs; with the fix, `testsuite/Test_solarGS98_allard` moves log L by up to 1.1e-4 dex, against 1.4e-4 for the null change, with the same models and end ages.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/starin.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/starin.f#L499-L510", "https://github.com/yreclab/yrec/blob/7be4133/src/starin.f#L518" ] }, { "file": "starin.f", "number": 2, "result": "not reached", "status": "unconfirmed", "title": "A reversed sign in the envelope interpolation", "runs": "No run, such as `examples/run_from_dbl_to_zams/Test_m0040_feh0_noGS_ZAMS`, placed the new point in the first envelope interval; a shallower fitting point there, with a last shell whose X or Z differs from the envelope's, would.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/starin.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/starin.f#L790-L791" ] }, { "file": "starin.f", "number": 3, "result": "harmless in the runs", "status": "unconfirmed", "title": "The defaults of R16_18 and R16_17 look swapped", "runs": "With 471.4 in place of 2693, `testsuite/Test_solarAAG21_dif` moves log L by up to 2.2e-5 dex, against 7.6e-5 for the null change; the starting oxygen-18 itself, the harm named, is not on the track.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/starin.f.html#b3", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/starin.f#L448-L449", "https://github.com/yreclab/yrec/blob/7be4133/src/nmlpar_mod.f#L232", "https://github.com/yreclab/yrec/blob/7be4133/src/nmlbypass_mod.f#L105" ] }, { "file": "starin.f", "number": 4, "result": "stops the run", "status": "checked by KC (Nov 10, 2026)", "title": "0/0 with one core point", "runs": "A made-up case: `testsuite/Test_solarGS98_kurucz` with `FCORE` = 1.0 stops in its first second in TPGRAD, at `STOP 912`; with the fix, it runs to its end.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/starin.f.html#b4", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/starin.f#L480", "https://github.com/yreclab/yrec/blob/7be4133/src/starin.f#L515", "https://github.com/yreclab/yrec/blob/7be4133/src/starin.f#L534" ] }, { "file": "starin.f", "number": 5, "result": "wrong output", "status": "checked by KC (Nov 10, 2026)", "title": "The warning prints an unset variable", "runs": "In `examples/class_exercises/exercise_4/Grid_m4p0feh+0p0GS98_std25_TAMS`, a mixing length of 1.853 against the run's 1.708 printed no warning; with the test fixed it prints, and only the `.short` file differs.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/starin.f.html#b5", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/starin.f#L328-L333", "https://github.com/yreclab/yrec/blob/7be4133/src/starin.f#L86" ] }, { "file": "stitch.f", "number": 1, "result": "changes results, small", "status": "checked by KC (Nov 10, 2026)", "title": "The integration steps are not restored", "runs": "In `testsuite/Test_m0p1feh+0p0GS98_spot25_ZAMS`, the fix moves log L by up to 5.8e-5 dex, against 9.4e-9 for the null change, and the last model is at 16.677 Gyr instead of 16.714.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/stitch.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/stitch.f#L140-L151" ] }, { "file": "stitch.f", "number": 2, "result": "wrong output", "status": "checked by KC (Nov 10, 2026)", "title": "β of the interior is read from the envelope's array", "runs": "In 17 runs the interior's written beta came from the envelope's array; with the fix, `testsuite/Test_m0p3feh+0p0GS98_scv_ZAMS` and `testsuite/Test_solarM22M_difkur` have identical tracks, and only their `.store` files change.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/stitch.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/stitch.f#L117-L120" ] }, { "file": "stitch.f", "number": 3, "result": "wrong output", "status": "checked by KC (Nov 10, 2026)", "title": "The atmosphere's gravity is the last point's", "runs": "Every atmosphere line writes the last envelope point's gravity (all 113 at 5.950E+02 in `testsuite/Test_m0p3feh+0p0GS98_scv_ZAMS`, model 7); with the fix, that case's track is identical and its `.store` file within 2%.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/stitch.f.html#b3", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/stitch.f#L212-L223" ] }, { "file": "stitch.f", "number": 4, "result": "never occurred", "status": "unconfirmed", "title": "The atmosphere's luminosity is set only in the envelope branch", "runs": "In the 13 runs that wrote the atmosphere's lines, `LSTENV` was on, so `B` was always set; `LSTCH` and `LSTATM` on with `LSTENV` off would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/stitch.f.html#b4", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/stitch.f#L153", "https://github.com/yreclab/yrec/blob/7be4133/src/stitch.f#L215" ] }, { "file": "stitch.f", "number": 5, "result": "never occurred", "status": "unconfirmed", "title": "LPULPT is never set", "runs": "In 17 runs `LPULPT` was false every time, as gfortran's zero-filled storage gives, so the results are right with that compiler.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/stitch.f.html#b5", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/stitch.f#L174-L176" ] }, { "file": "stitch.f", "number": 6, "result": "not reached", "status": "unconfirmed", "title": "Three routines apply the spots only without a surface convection zone", "runs": "Not reached in the spotted 1 M☉ case: STARIN's fitting point needed no move, and WRTMOD is not called with LSTCH, which every spotted YREC case sets; the fix left the outputs identical.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/stitch.f.html#b6", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/stitch.f#L167-L173", "https://github.com/yreclab/yrec/blob/7be4133/src/starin.f#L723", "https://github.com/yreclab/yrec/blob/7be4133/src/calcad.f#L228", "https://github.com/yreclab/yrec/blob/7be4133/src/wrtmod.f#L228" ] }, { "file": "surfp.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "Missing gravities used at the top of a row", "runs": "In `testsuite/Test_solarGS98_kurucz` and `testsuite/Test_solarM22M_difkur`, with `atmk1990p00.tab`, no row with missing highest gravities was met; a cool dwarf with a table of [Fe/H] at or below -2.3 would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/surfp.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/surfp.f#L88-L102", "https://github.com/yreclab/yrec/blob/7be4133/src/surfp.f#L104-L123" ] }, { "file": "taucal.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "ICZTOP is not set for a single zone", "runs": "TAUCAL runs only with `LNEWTCZ` off, which no case sets; `LNEWTCZ` = .FALSE., with the base of a single surface convection zone inside the envelope, would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/taucal.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/taucal.f#L87-L98", "https://github.com/yreclab/yrec/blob/7be4133/src/taucal.f#L152-L153", "https://github.com/yreclab/yrec/blob/7be4133/src/taucal.f#L159-L163" ] }, { "file": "tauintnew.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "The overturn time at IMAX is overwritten", "runs": "No model in the 33 runs took the branch; a single surface convection zone with few shells near its base would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/tauintnew.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/tauintnew.f#L176-L181", "https://github.com/yreclab/yrec/blob/7be4133/src/tauintnew.f#L229-L231", "https://github.com/yreclab/yrec/blob/7be4133/src/tauintnew.f#L235-L272" ] }, { "file": "tauintnew.f", "number": 2, "result": "never occurred", "status": "unconfirmed", "title": "The restarts can read past the profile, and can loop forever", "runs": "Never occurred: 120 restarts in the 1.4 M☉ diffusion case (IBEG 298 to 352, M 407 to 427), none in four other cases; IBEG never passed M, nor IMAX+2. The fix was not run.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/tauintnew.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/tauintnew.f#L75", "https://github.com/yreclab/yrec/blob/7be4133/src/tauintnew.f#L115-L131", "https://github.com/yreclab/yrec/blob/7be4133/src/tauintnew.f#L284-L287", "https://github.com/yreclab/yrec/blob/7be4133/src/tauintnew.f#L92-L98" ] }, { "file": "thdiff.f", "number": 1, "result": "cannot be reached", "status": "unconfirmed", "title": "No check of the array sizes", "runs": "Both callers pass `NSPEC` = 4.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/thdiff.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/thdiff.f#L48-L53", "https://github.com/yreclab/yrec/blob/7be4133/src/thdiff.f#L72" ] }, { "file": "tpgrad.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "Spots can stop the run on a point that is convective without them", "runs": "In `testsuite/Test_m0p1feh+0p0GS98_spot25_ZAMS`, the spot correction ran 1,418,129 times and `DELDEL` stayed at 0.0788 or more; a spotted point near a zone's edge where the correction makes `DELDEL` vanish would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/tpgrad.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/tpgrad.f#L82-L83", "https://github.com/yreclab/yrec/blob/7be4133/src/tpgrad.f#L127-L135", "https://github.com/yreclab/yrec/blob/7be4133/src/tpgrad.f#L145-L156" ] }, { "file": "tpgrad.f", "number": 2, "result": "changes results, small", "status": "checked by KC (Nov 10, 2026)", "title": "With spots, the derivatives of the gradient are those of the unspotted one", "runs": "In the spotted 1 M☉ case the derivatives were taken 974,865 times with S = 0.878 to 0.99982; the fix moves log L by up to 0.00055 dex, against 0.000051 for a null change.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/tpgrad.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/tpgrad.f#L127-L135", "https://github.com/yreclab/yrec/blob/7be4133/src/tpgrad.f#L201-L202" ] }, { "file": "trapzd.f", "number": 1, "result": "harmless in the runs", "status": "unconfirmed", "title": "The interpolation does not move with the midpoint", "runs": "In `testsuite/Test_m0p8feh+0p0GS98_rotp8p0d0p1mag1e4std_TAMS`, the fix moves central X by 1.7e-4 at one model, against at most 6.1e-5 for the null change, and log L by 0.00018 dex, against 0.00019.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/trapzd.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/trapzd.f#L34-L35", "https://github.com/yreclab/yrec/blob/7be4133/src/trapzd.f#L44-L47" ] }, { "file": "vcirc.f", "number": 1, "result": "cannot be reached", "status": "unconfirmed", "title": "The $\\mu$ factor compounds on FESTIME", "runs": "`FESTIME` is read only in the `LDIFAD` branch, and PARMIN stops any run that sets `LDIFAD` with rotation.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/vcirc.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/vcirc.f#L256", "https://github.com/yreclab/yrec/blob/7be4133/src/vcirc.f#L304-L317" ] }, { "file": "vcirc.f", "number": 2, "result": "changes results, small", "status": "checked by KC (Nov 10, 2026)", "title": "CC is never set", "runs": "In `testsuite/Test_m2p0_feh+0p0_vzams150_IES2_TAMS`, the fix moves log L by up to 0.0056 dex and central log T by 0.0055, against about 2e-8 for the null change, with the same 546 models.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/vcirc.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/vcirc.f#L121-L122", "https://github.com/yreclab/yrec/blob/7be4133/src/vcirc.f#L286-L288" ] }, { "file": "vcirc.f", "number": 3, "result": "never occurred", "status": "unconfirmed", "title": "FES for the GSF velocity", "runs": "In `testsuite/Test_m2p0_feh+0p0_vzams150_IES2_TAMS`, the fallback ran with `VGSF` = 0 every time, and `FES` = `FGSF` in every namelist; `IMU` = 3 with `FES` and `FGSF` different and a positive GSF velocity would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/vcirc.f.html#b3", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/vcirc.f#L564-L568" ] }, { "file": "vcirc.f", "number": 4, "result": "cannot be reached", "status": "unconfirmed", "title": "IES = 2 with LDIFAD reuses the previous call's parts", "runs": "PARMIN stops any run that sets `LDIFAD` with rotation.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/vcirc.f.html#b4", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/vcirc.f#L202-L226", "https://github.com/yreclab/yrec/blob/7be4133/src/vcirc.f#L602-L607", "https://github.com/yreclab/yrec/blob/7be4133/src/vcirc.f#L629-L632" ] }, { "file": "viscos.f", "number": 1, "result": "harmless in the runs", "status": "unconfirmed", "title": "ELECN is not the electron density", "runs": "In `testsuite/Test_3Msolar_GS_rot_gray_OPALSCV_SF2_GS98_OP_CF10_RGp2diffrotI3cz2`, the fix moves log L by up to 0.0022 dex, against 0.0076 for the null change; the other rotating cases were not run with it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/viscos.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/viscos.f#L48", "https://github.com/yreclab/yrec/blob/7be4133/src/viscos.f#L57", "https://github.com/yreclab/yrec/blob/7be4133/src/viscos.f#L61", "https://github.com/yreclab/yrec/blob/7be4133/src/viscos.f#L73" ] }, { "file": "viscos.f", "number": 2, "result": "harmless in the runs", "status": "unconfirmed", "title": "The total $Z$ counts as an ion", "runs": "In the rotating solar model species 3 entered COMSUM at 1,789,275 shells; the fix's effect is about the size of the null change, not clearly beyond it (log L 0.00021 against 0.00029 dex, surface rotation 0.0089 against 0.0068 km/s).", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/viscos.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/viscos.f#L38-L42", "https://github.com/yreclab/yrec/blob/7be4133/src/viscos.f#L52-L57", "https://github.com/yreclab/yrec/blob/7be4133/src/viscos.f#L73" ] }, { "file": "viscos.f", "number": 3, "result": "harmless in the runs", "status": "unconfirmed", "title": "Three isotope masses are mistyped", "runs": "In the rotating solar model the three mistyped masses were used at 1,789,275 shells; the fix's effect is about the size of the null change, not clearly beyond it (surface rotation 0.0095 against 0.0068 km/s, log L 0.00027 against 0.00029 dex).", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/viscos.f.html#b3", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/viscos.f#L38-L40" ] }, { "file": "wind.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "WNEW overwrites the namelist's solid-body rate", "runs": "Every rotating case uses the new law, so WIND never ran; `LROT` with `LMWIND` off, and a later run that reads its model from file with `LWNEW`, would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/wind.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/wind.f#L122", "https://github.com/yreclab/yrec/blob/7be4133/src/wind.f#L133", "https://github.com/yreclab/yrec/blob/7be4133/src/wind.f#L142" ] }, { "file": "wrthead.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "The headers of track versions 1, 2, and 3 do not match the data", "runs": "Every case sets `ITRVER` = 0; the default `ITRVER` = 1, which a namelist without the line gets, or 2 or 3, would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/wrthead.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/wrthead.f#L152-L176", "https://github.com/yreclab/yrec/blob/7be4133/src/wrtout.f#L570-L641" ] }, { "file": "wrtlst.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "A lost label", "runs": "The SCV branch ran at every model in 32 runs, but never with `LOPALE06` on and `LDH` off; YREC's cases use SCV only with `LDH`.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/wrtlst.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/wrtlst.f#L110-L118" ] }, { "file": "wrtlst.f", "number": 2, "result": "never occurred", "status": "unconfirmed", "title": "Unset atmosphere label", "runs": "In 33 runs `KTTAU` was never outside 0 to 5; only an invalid `KTTAU`, which nothing in YREC offers, would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/wrtlst.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/wrtlst.f#L93-L106" ] }, { "file": "wrtmil.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "The file is closed after each model", "runs": "`LMILNE` is off by default and set in no case; `LMILNE` on, with two or more models written, would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/wrtmil.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/wrtmil.f#L109", "https://github.com/yreclab/yrec/blob/7be4133/src/parmin.f#L956-L959" ] }, { "file": "wrtmod.f", "number": 1, "result": "wrong output", "status": "checked by KC (Nov 10, 2026)", "title": "A typo in $d\\rho/dr$", "runs": "In seven solar cases with `LPULSE`, such as `testsuite/Test_solarGS98_dif`, `QPR1` differed from `QPR2` by -5.6% to +10.3% on every shell; with the fix, `testsuite/Test_solarGS98_allard` and `testsuite/Test_solarM22M_difkur` have identical tracks and `.store` files.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/wrtmod.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/wrtmod.f#L165-L168" ] }, { "file": "wrtmod.f", "number": 2, "result": "wrong output", "status": "checked by KC (Nov 10, 2026)", "title": "The pulsation model's shell list is never set", "runs": "In all nine runs with `LPULSE`, such as `testsuite/Test_solarGS98_dif`, `IDM` was 0 every time, so the pulsation model was written with 0 shells; no fix was run.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/wrtmod.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/wrtmod.f#L51", "https://github.com/yreclab/yrec/blob/7be4133/src/wrtmod.f#L255-L256", "https://github.com/yreclab/yrec/blob/7be4133/src/wrtmod.f#L265-L266" ] }, { "file": "wrtmod.f", "number": 3, "result": "wrong output", "status": "checked by KC (Nov 10, 2026)", "title": "Debugging prints", "runs": "The prints ran at every call in 10 runs; with them removed, `testsuite/Test_solarAAG21_dif` and `testsuite/Test_solarGS98_allard` have identical tracks and `.store` files.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/wrtmod.f.html#b3", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/wrtmod.f#L127", "https://github.com/yreclab/yrec/blob/7be4133/src/wrtmod.f#L131" ] }, { "file": "wrtmod.f", "number": 4, "result": "never occurred", "status": "unconfirmed", "title": "LPRT may be unset", "runs": "In 10 runs the line always ran with `LSTORE` on, so `LPRT` was set; `LPULSE` with `LSTORE` off would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/wrtmod.f.html#b4", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/wrtmod.f#L186", "https://github.com/yreclab/yrec/blob/7be4133/src/envint.f#L206", "https://github.com/yreclab/yrec/blob/7be4133/src/envint.f#L222" ] }, { "file": "wrtout.f", "number": 1, "result": "not reached", "status": "unconfirmed", "title": "Track version 2 without rotation writes rotation values never set", "runs": "Every case sets `ITRVER` = 0; `ITRVER` = 2 without `LROT` would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/wrtout.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/wrtout.f#L594-L608", "https://github.com/yreclab/yrec/blob/7be4133/src/wrtout.f#L630-L634" ] }, { "file": "wrtout.f", "number": 2, "result": "wrong output", "status": "checked by KC (Nov 10, 2026)", "title": "The summary gives the neutrino fluxes and capture rates in erg", "runs": "With the labels corrected, `testsuite/Test_m0p3feh+0p0GS98_scv_ZAMS` and `testsuite/Test_3Msolar_GS_rot_gray_OPALSCV_SF2_GS98_OP_CF10_RGp2diffrotI3cz2` have identical tracks; only their `.short` files differ.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/wrtout.f.html#b2", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/wrtout.f#L401-L403", "https://github.com/yreclab/yrec/blob/7be4133/src/wrtout.f#L423-L426" ] }, { "file": "wtime.f", "number": 1, "result": "never occurred", "status": "unconfirmed", "title": "0/0 where $\\Omega$ is zero", "runs": "In all five rotating runs, the central rotation rate was never zero; a rotating run with no rotation at the center would reach it.", "write_up": "https://kailicao.github.io/pyrec-map/yrec-map/wtime.f.html#b1", "code": [ "https://github.com/yreclab/yrec/blob/7be4133/src/wtime.f#L53-L62" ] } ] }