Modeling of AR 12760 with GX Simulator and Evidence for the Extended Transition Region in Peripheral Active Region Loops - data files

Files to accompany "Modeling of AR~12760 with GX~Simulator and Evidence for the Extended Transition Region in Peripheral Active Region Loops" by Kucera, Nita, Klimchuck, & Fleishman, submitted to ApJ 2025 This work used the gx_simulator Solar Soft IDL package which can be found at https://github.com/Gelu-Nita/GX_SIMULATOR Used with Solar Soft https://www.lmsal.com/solarsoft/ Directories and files: SDO - directory containing time averaged SDO files in IDL Save format each map is a MAP structure array containing the map and the uncertainties (set to a value of 1 in this case). GXM - Files containing the magnetic field model as well as information about the transition region mask and EBTEL files to be used. hmi.M_720s.20200428_111042.E55S7CR.CEA.NAS.CHR_no_background_200TR_e_scl_1_alpha_-1.gxm - Magnetic Model 1 hmi.M_720s.20200428_111042.E55S7CR.CEA.NAS.CHR_180x180x200_dx2800_no_bckgrnd_TR200_eb1-1.gxm - Magnetic Model 2 can create and inspect with gx.pro routine, and used by gx_search4bestq.pro (see below) FitMasks - Directory containing fit masks specifying which pixel should be used in intensity model/data comparisons IDL save files containing arrays (named 'MASK')specifying what portions of the images should be used for image model comparison. Used by gx_search4bestq.pro mask_0percent211A_noedge_Nx30Ny30xr-195t45yr-170t70.sav - used with Magnetic Model 1, 8" resolution comparison mask_0percent211A_noedge_Nx120Ny120xr-195t45yr-170t70.sav - used with Magnetic Model 1, 2" resolution comparison mask_0percentI_noedge_Nx120Ny120xr-195t45yr-170t70_dx2800.sav - used with Magnetic Model 2, 2" resolution comparison mask_0percent211A_noedge_noTR_Nx120Ny120xr-195t45yr-170t70_dx1400.sav - used with Magnetic Model 1, 2" resolution comparison, removing pixels thought to contain TR emission. Results These files contain the results for the best q_0 fits for various a,b grids. Results for Magnetic Model 1, 2" model/data comparison, including both model corona and transition region emission results_94A_yesTR.sav results_131A_yesTR.sav results_171A_yesTR.sav results_193A_yesTR.sav results_211A_yesTR.sav results_335A_yesTR.sav Results for Magnetic Model 2, 2" model/data comparison, including both model corona and transition region emission results_171A_yesTR_MMod2.sav results_211A_yesTR_MMod2.sav Results for Magnetic Model 1, 8" model/data comparison, including both model corona and transition region emission results_211A_yesTR_resize30.sav results_335A_yesTR_resize30.sav Results for Magnetic Model 1, 2" model/data comparison, including only the model corona and not transition region emission results211A_noTR.sav The input SDO data, GMX model and FitMasks can be used to calculate the best q0 values for a given a, b pair with calls like this: IDL> gxmfile='GXM/hmi.M_720s.20200428_111042.E55S7CR.CEA.NAS.CHR_no_background_200TR_e_scl_1_alpha_-1.gxm' IDL> aiafile='AIA/AIA211_r20200428111500_s20200428091511_e20200428131011_const_uncert1.sav' IDL> Fitmaskfile='FitMasks/mask_0percent211A_noedge_noTR_Nx120Ny120xr-195t45yr-170t70_dx1400.sav' IDL> a=1 & b=2 IDL> result=gx_search4bestq(gxmpath=gxmfile,$ moddir='moddir',$ ;output directory for calculated model maps psDir='psdir',$ ;output directory for plot files with information about the fit to q0 tmpDir='tmpdir',$. ;output directory for temp files refdatapath=aiafile,$ ebtel_path='ebtel_scale=1_alpha=-1.sav',$. ;one of existing EBTEL files in GX simulator distribution renderer='aia_lib.pro',$. ;existing aia renderer file in GX Simulator distribution Nx=120, Ny=120, resize=[120,120],$ ;Size of input and output maps in pixels xc=-75, yc=-50, xfov=240, yfov=240,$ ; center and extent of field of view in arcsec a_arr= a, b_arr=b, q_start=[ 0.00100000, 0.0100000],$ ; a, b, and initial range for q0 values shift=[0,0],$ ;do not attempt to coalign the magnetic field and intensity data (not needed for SDO) addTR=1,$ ;Include the transition region (1 is default) all=1,$ ;make calculations parallel mask=FitMaskFile) In order to do fits to a full a,b grid it is recommended that the user use the GUI based gx_chmp.pro The results files are produced by the gx_chmp.pro routine and include and array of the results from a grid of gx_search4bestq .pro results. for example: IDL> ResultFileName='Results/results_211A_yesTR.sav' IDL> restore,ResultFilename,/ver IDL> tmp=min(result.chi2_best,w) ;find the index,w, of the best fit IDL> a=result[w].a ;best fit a IDL> b=result[w].b ;best fit b IDL> q0= result[w].q_chi2_best ;best fit q0 IDL> modmap=result[w].chi2_best_metrics->get(0,/map) ;model map. SolarSoft MAP structure IDL> AIAmap=result[w].chi2_best_metrics->get(1,/map) ;The AIA map compared to model map. SolarSoft MAP structure.