DOI: 10.5281/zenodo.16896515
title: Modeling of AR 12760 with GX Simulator and Evidence for the Extended Transition Region in Peripheral  Active Region Loops - data files

description: 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.