This dataset contains the input files (main_input and custom_reference_binary), basic simulation info (grid_info, jobinfo.txt, and equation_coefficients) and the last Rayleigh checkpoints for the two cases published in the ApJ Letter, Confinement of the Solar Tachocline by Dynamo Action in the Radiative Interior. The data for each case (the "HD case" and "MHD case" as referred to in the paper) are given in .zip files. NOTE: The MHD case has custom diagnostic outputs and the main_input file (as directly given here) can only be run using Loren Matilsky's personal branch of the code: https://github.com/illorenzo7/Rayleigh/tree/me_prod_subcurv (for full reproducibility, commit fcf1a7b81624a0af7866d52d5b7836549268044e). Otherwise, the MHD case should run using the current main branch of the Rayleigh code (https://github.com/geodynamics/Rayleigh), as long as the custom quantity codes (>2200) are removed.