⏲️ Update confinement time docs

documentation/proc-pages/fusion-devices/stellarator.md

@@ -74,7 +74,7 @@ ixc = 10 * "ISS04 Renormalization Factor" (can also be fixed)
 ixc = 50 * Coil current density, aka winding pack thickness (required!) 
 ixc = 59 * Winding Pack copper fraction
 ixc = 56 * Exponential Quench Dumping Time
-ixc = 109 * Thermal alpha particle pressure (iterated to consistency, use together with `taulimit`)
+ixc = 109 * Thermal alpha particle pressure (iterated to consistency, use together with `f_alpha_energy_confinement_min`)
 ixc = 169 * Achievable Temperature of the ECRH at the ignition point
 ```
 
@@ -142,13 +142,13 @@ The constraint equation `icc=91` together with `ixc = 169` enforces that the fou
 
 ### $\tau_E$ scaling laws
 
-Five confinement time scaling laws relevant to stellarators are present within `PROCESS`. The value of switch isc` determines which of these in the plasma energy balance calculation.
+Five confinement time scaling laws relevant to stellarators are present within `PROCESS`. The value of switch i_confinement_time` determines which of these in the plasma energy balance calculation.
 
-$\tau_E$ (Large Helical Device[^8]: `isc=21`) = $0.17 \, R^{0.75}_0 \, a^2_p \, \bar{n}^{0.69}_{20} \, B^{0.84}_0 \, P^{-0.58}$  
-$\tau_E$ (Gyro-reduced Bohm[^9]: `isc=22`) = $0.25 \, R^{0.6}_0 \, a^{2.4}_p \, \bar{n}^{0.6}_{20} \, B^{0.8}_0 \, P^{-0.6}$  
-$\tau_E$ (Lackner-Gottardi[^10]: `isc=23`) = $0.17 \, R_0 \, a^2_p \, \bar{n}^{0.6}_{20} \, B^{0.8}_0 \, P^{-0.6} \, \iota^{0.4}$  
-$\tau_E$ (ISS95[^11]: `isc=37`) = $0.079 \, R^{0.65}_0 \, a^{2.21}_p \, \bar{n}^{0.51}_{20} \, B^{0.83}_0 \, P^{-0.59} \, \bar{\iota}^{0.4}$  
-$\tau_E$ (ISS04[^12]: `isc=38`) = $0.134 \, R^{0.64}_0 \, a^{2.28}_p \, \bar{n}^{0.54}_{20} \, B^{0.84}_0 \, P^{-0.61} \, \bar{\iota}^{0.41}$
+$\tau_E$ (Large Helical Device[^8]: `i_confinement_time=21`) = $0.17 \, R^{0.75}_0 \, a^2_p \, \bar{n}^{0.69}_{20} \, B^{0.84}_0 \, P^{-0.58}$  
+$\tau_E$ (Gyro-reduced Bohm[^9]: `i_confinement_time=22`) = $0.25 \, R^{0.6}_0 \, a^{2.4}_p \, \bar{n}^{0.6}_{20} \, B^{0.8}_0 \, P^{-0.6}$  
+$\tau_E$ (Lackner-Gottardi[^10]: `i_confinement_time=23`) = $0.17 \, R_0 \, a^2_p \, \bar{n}^{0.6}_{20} \, B^{0.8}_0 \, P^{-0.6} \, \iota^{0.4}$  
+$\tau_E$ (ISS95[^11]: `i_confinement_time=37`) = $0.079 \, R^{0.65}_0 \, a^{2.21}_p \, \bar{n}^{0.51}_{20} \, B^{0.83}_0 \, P^{-0.59} \, \bar{\iota}^{0.4}$  
+$\tau_E$ (ISS04[^12]: `i_confinement_time=38`) = $0.134 \, R^{0.64}_0 \, a^{2.28}_p \, \bar{n}^{0.54}_{20} \, B^{0.84}_0 \, P^{-0.61} \, \bar{\iota}^{0.41}$
 
 Here $\bar{\iota}$ is the rotational transform, which is equivalent to the reciprocal of the tokamak safety factor $q$.
 

documentation/proc-pages/physics-models/error.txt

@@ -596,15 +596,15 @@ energy confinement times \(\tau_E\) are equal.
 
 Many energy confinement time scaling laws are present within PROCESS,
 for tokamaks, RFPs or stellarators. These are calculated in routine
-\texttt{pcond}. The value of \texttt{isc} determines which of the
+\texttt{pcond}. The value of \texttt{i_confinement_time} determines which of the
 scalings is used in the plasma energy balance calculation. The table
 below summarises the available scaling laws. The most commonly used is
 the so-called IPB98(y,2) scaling.
 
 \begin{longtable}[]{@{}cll@{}}
 \toprule
 \begin{minipage}[b]{0.05\columnwidth}\centering\strut
-\texttt{isc}\strut
+\texttt{i_confinement_time}\strut
 \end{minipage} & \begin{minipage}[b]{0.03\columnwidth}\raggedright\strut
 scaling law\strut
 \end{minipage} & \begin{minipage}[b]{0.03\columnwidth}\raggedright\strut
@@ -967,7 +967,7 @@ J. Menard 2019, Phil. Trans. R. Soc. A 377:201704401\strut
 \begin{minipage}[t]{0.05\columnwidth}\centering\strut
 48\strut
 \end{minipage} & \begin{minipage}[t]{0.03\columnwidth}\raggedright\strut
-Use input \texttt{tauee\_in}\strut
+Use input \texttt{t_electron_energy_confinement\_in}\strut
 \end{minipage} & \begin{minipage}[t]{0.03\columnwidth}\raggedright\strut
 \strut
 \end{minipage}\tabularnewline
@@ -987,21 +987,21 @@ additional radiation causes an almost equal drop in power transported by
 ions and electrons, leaving the confinement nearly unchanged.
 
 To allow for these uncertainties, three options are available, using the
-switch \texttt{iradloss}. In each case, the particle transport loss
+switch \texttt{i_rad_loss}. In each case, the particle transport loss
 power \(P_{\mbox{loss}}\), referred to in the code as \texttt{pscaling},
 is derived directly from the energy confinement scaling law.
 
-\texttt{iradloss\ =\ 0} -- Total power lost is scaling power plus
+\texttt{i_rad_loss\ =\ 0} -- Total power lost is scaling power plus
 radiation
 
 \texttt{pscaling\ +\ pden_plasma_rad_mw\ =\ f_alpha_plasma*alpha_power_density\ +\ charged_power_density\ +\ pden_plasma_ohmic_mw\ +\ pinjmw/vol_plasma}
 
-\texttt{iradloss\ =\ 1} -- Total power lost is scaling power plus core
+\texttt{i_rad_loss\ =\ 1} -- Total power lost is scaling power plus core
 radiation only
 
 \texttt{pscaling\ +\ pcoreradpv\ =\ f_alpha_plasma*alpha_power_density\ +\ charged_power_density\ +\ pden_plasma_ohmic_mw\ +\ pinjmw/vol_plasma}
 
-\texttt{iradloss\ =\ 2} -- Total power lost is scaling power only, with
+\texttt{i_rad_loss\ =\ 2} -- Total power lost is scaling power only, with
 no additional allowance for radiation. This is not recommended for power
 plant models.
 
@@ -1376,9 +1376,9 @@ too high.
 Laws}{Inverse Quadrature in \textbackslash{}tau\_E Scaling Laws}}\label{inverse-quadrature-in-tau_e-scaling-laws}
 
 Switch \texttt{iinvqd} determines whether the energy confinement time
-scaling laws due to Kaye-Goldston (\texttt{isc\ =\ 5}) and Goldston
-(\texttt{isc\ =\ 9}) should include an inverse quadrature scaling with
-the Neo-Alcator result (\texttt{isc\ =\ 1}). A value
+scaling laws due to Kaye-Goldston (\texttt{i_confinement_time\ =\ 5}) and Goldston
+(\texttt{i_confinement_time\ =\ 9}) should include an inverse quadrature scaling with
+the Neo-Alcator result (\texttt{i_confinement_time\ =\ 1}). A value
 \texttt{iinvqd\ =\ 1}includes this scaling.
 
 \subsubsection{Plasma / First Wall Gap}\label{plasma-first-wall-gap}