This folder contains two types of files:
1. CSV (*.csv) files that list the critical values used in each figure.
2. DAT (*.dat) files that give radial profiles of u_r, w_r, and Θ at the onset of convection for a specific parameter set.
 
CSV format (columns in this order):
HT, MBC, E, Pr, Chi, Ra_c^0, Ra_c^1, Ra_c^2, Ra_c^3, Ra_c^4, m_c, omega_c

    • HT: heating type ("dif" for differential, "int" for internal)
    • MBC: mechanical boundary condition ("NSNS", "SFSF", or "NSSF")
    • E: Ekman number
    • Pr: Prandtl number
    • Chi: radius ratio
    • Ra_c^0 through Ra_c^4: critical Rayleigh numbers for formulations 0–3
    • m_c: critical azimuthal wavenumber
    • omega_c: critical frequency
    
Formulation 0 is the direct output from the code and is related to Formulation 1 by:

Differential heating: Ra_c^0 = Ra_c^1 * 1/(1 - Chi)^2
Internal heating: Ra_c^0 = Ra_c^1 * (1 - Chi)

DAT format: three space-separated columns:

1. u_r radial velocity
2. w_r radial vorticity
3. Theta temperature perturbation
    
DAT filenames follow: HT_MBC_Pr_Chi_E.dat
for example:

dif_NSNS_1.0_0.35_1e-4.dat
means differential heating, no-slip boundaries, Pr=1.0, Chi=0.35, E=1e-4.


See OnsetPlot.py and CriticalValues.py for examples of how to plot the equatorial radial velocity at onset and find the critical values for a specific point in parameter space.