setscal.f

Prepares the list of runs for the automatic calibration of a star to a given luminosity and effective temperature (or radius): 50 runs, 25 repetitions of a two-run pattern. The adjustment between repetitions is done in CHKSCAL. It is the analogue, for any star, of SETCAL for the Sun.

• Who calls it: the main program (main.f), when LCALST is true, before the loop over runs.
• Arguments: the state of the search, which the main program keeps and passes on to CHKSCAL: ALRI (named "log(R/Rsun) of the previous model," though CHKSCAL stores $R/R_\odot$; set to 0, line 107↗, which makes a track that starts above $R_*$ cross at once: chkscal.f, suspected bug 1) and the flags LSTAR, LPASSR (set to false); the comments in the routine define each.
• Through modules: NMLPAR_MOD, the target (XLS the luminosity in solar units, XLSTOL its tolerance, LTEFF, STEFF, SR) and the run list's arrays, and RUNLOOP_MOD, as in SETCAL; CONSTANTS_MOD; IOUNITS_MOD (ITRACK).
• Outputs: the run list; the target's radius or temperature; a line "Evolve to R, L =" on the terminal and in the track file (ITRACK).

From $L = 4\pi R^2\sigma T_\mathrm{eff}^4$, the missing one of $R$ and $T_\mathrm{eff}$ is computed:

$$\frac{R}{R_\odot} = \frac{\sqrt{L/(4\pi\sigma)}}{T_\mathrm{eff}^2\,R_\odot}\quad(\texttt{LTEFF}\ \text{true}),\qquad T_\mathrm{eff} = \left(\frac{L}{4\pi\sigma R^2}\right)^{1/4}\quad(\text{otherwise}),$$

with $L$ = XLS·CLSUN (erg s$^{-1}$), $\sigma$ = CSIG, $R_\odot$ = CRSUN (cm). Then NUMRUN = 50: runs 3, 5, …, 49 repeat run 1 (a new starting model), and runs 4, 6, …, 50 repeat run 2 (evolution from the previous model).

The comment's "24" counts pairs of runs. The comment speaks of "THE NEXT SERIES OF 24 CALIBRATING RUNS" (lines 108–112↗); the loops over runs 3 to 50 in steps of 2 (line 120↗; line 131↗) fill 48 runs after the first two, that is, 24 repetitions of the two-run pattern, so the comment could say 24 pairs.

The new stop criteria are not copied. END_DCEN, END_XCEN, and END_YCEN ("MHP 10/24 ADDED STOP CRITERIA," line 57↗; per-run arrays in NMLPAR_MOD, nmlpar_mod.f, lines 261–263↗) are not copied to the repeated runs, unlike ENDAGE and SETDT (lines 131–139↗), so a value set for run 2 stops run 2 itself but not its repetitions, runs 4 to 50: the same omission as SETCAL's suspected bug 1. That bug's fix patches SETCAL only; here it would add END_DCEN(I) = END_DCEN(2), and the same for END_XCEN and END_YCEN, to the loop at lines 131–139↗.