PRO xrt_teem_ch713, index1, data1, index2, data2, te, em, et, ee, $ solar_spectrum_file=solar_spectrum_file, $ ; for spectrum file input hybrid=hybrid, photospheric=photospheric, corona=corona, $ ; abundance switch ; bin = bin, trange = trange, no_threshold = no_threshold, $ te_err_threshold = te_err_threshold, $ photon_noise_threshold = photon_noise_threshold, $ mask1 = mask1, mask2 = mask2, eff_area = eff_area ; ========================================================================= ;+ ; PROJECT: ; ; Solar-B / XRT ; ; NAME: ; ; XRT_TEEM_CH713 ; ; CATEGORY: ; ; XRT CORONAL TEMPERATURE DIAGNOSTICS ; ; PURPOSE: ; ; Caliclate coronal temperature and emission measure using filter ; ratio method. Values are derived in log scale. ; ; This is a wrapper routine to . This routine allows ; to use a variety of solar spectrum files calculated with CHIANTI ; database version 7.1.3 with different combination of elementalry ; abundances and ionization equilibrium. The solar spectrum file ; to be used is passed to xrt_teem.pro (and propagated further to ; xrt_flux.pro) through the /spectrum keyword. ; ; As default, this program uses the solar spectrum, ; $SSW_XRT/idl/response/spectrum_files/solspec_ch713_corona_chianti.genx, ; (currently in MSU, /disk/yla/software/xrt_ssw/spectrum_files/) ; that is calculated with CHIANTI ver. 7.1.3 (density: 10^9 [cm^-3], ; ionization equilibrium: chianti.ioneq, ; abundance: sun_coronal_1992_feldman_ext). ; ; CALLING SEQUENCE: ; ; XRT_TEEM_CH713, index1, data1, index2, data2, te, em, et, ee, $ ; [solar_spectrum_file = solar_spectrum_file], $ ; [corona=corona], [hybrid=hybrid], [photospheric=photoshperic], $ ; [bin=bin], [te_err_threshold=te_err_threshold], $ ; [photon_noise_threshold=photon_noise_threshold], $ ; [/no_threshold], [/info], [eff_area = eff_area] ; ; INPUTS: ; ; INDEX1 - (structure array) XRT index. ; DATA1 - (2 or 3-dim float array, [x, y, t]) ; XRT *non-normalized* level1 data. ; INDEX2 - (structure array) XRT index. ; DATA2 - (2 or 3-dim float array, [x, y, t]) ; XRT *non-normalized* level1 data. ; ; KEYWORDS: ; ; SOLAR_SPECTRUM_FILE - [Optional input] (scalar string) ; specify the name of the solar sectrum files (with path). ; CORONA - [Optional input] (boolean) ; switch to use the spectrum file, ; $SSW_XRT/idl/response/spectrum_files/ ; solspec_ch713_corona_chianti.genx. ; HYBRID - [Optional input] (boolean) ; switch to use the spectrum file, ; $SSW_XRT/idl/response/spectrum_files/ ; solspec_ch713_hybrid_chianti.genx. ; PHOTOSPHERIC - [Optional input] (boolean) ; switch to use the spectrum file, ; $SSW_XRT/idl/response/spectrum_files/ ; solspec_ch713_photos_chianti.genx. ; BIN - [Optional input] same as xrt_teem.pro ; TE_ERR_THRESHOLD - [Optional input] same as xrt_teem.pro ; PHOTON_NOISE_THRESHOLD - [Optional input] same as xrt_teem.pro ; /NO_THRESHOLD - [Optional] same as xrt_teem.pro ; /INFO - [Optional] same as xrt_teem.pro ; EFF_AREA - [Optional] same as xrt_teem.pro ; ; OUTPUTS: ; ; TE - (2-dim float array, [x, y]) ; derived temperature in log scale [log K]. ; EM - (2-dim float array, [x, y]) ; derived volume emission measure in log scale [log cm^-3]. ; ET - (2-dim float array, [x, y]) ; error of temperature in log scale [log K]. ; EE - (2-dim float array, [x, y]) ; error of volume emission measure in log scale [log cm^-3]. ; ; EXAMPLES: ; ; Use default spectrum file, i.e., 'solspec_ch713_corona_chianti.genx' ; in $SSW_XRT/idl/response/spectrum_files. ; IDL> xrt_teem_ch713, index1, data1, index2, data2, te, em, et, ee ; ; Use spectrum file with hybrid abundance, i.e., ; $SSW_XRT/idl/response/spectrum_files/solspec_ch713_hybrid_chianti.genx. ; IDL> xrt_teem_ch713, index1, data1, index2, data2, te, em, et, ee, /hybrid ; ; Use other spectrum file (file path required). ; IDL> sp_file='$SSW_XRT/idl/response/spectrum_files/'+ $ ; 'solspec_ch713_cr2012schmelz_chianti.genx' ; IDL> xrt_teem_ch713, index1, data1, index2, data2, te, em, et, ee, $ ; solar_spectrum_file=sp_file ; ; ; COMMON BLOCKS: ; ; none ; ; NOTES: ; ; call xrt_teem.pro. ; ; CONTACT: ; ; MODIFICATION HISTORY: ; ; progver = 'v2014-Feb-22' Aki Takeda (MSU). ; ;- ; ========================================================================= ps = path_sep() dir_spec = get_logenv('$SSW_XRT') + ps + 'idl' + ps + 'response' + ps dir_spec = '/disk/yla/software/xrt_ssw/solar_spectrum/' ; abd_switch = 0 if keyword_set(corona) or keyword_set(hybrid) or keyword_set(photospheric) then $ abd_switch = 1 ; if not keyword_set(solar_spectrum_file) and abd_switch eq 0 then $ spfile = dir_spec+'solspec_ch713_corona_chianti' $ ; default file else begin if keyword_set(photospheric) then $ ; photospheric abundance spfile = dir_spec+'solspec_ch713_photos_chianti' if keyword_set(hybrid) then $ ; hybrid abundance spfile = dir_spec+'solspec_ch713_hybrid_chianti' if keyword_set(corona) then $ ; coronal abundance spfile = dir_spec+'solspec_ch713_corona_chianti' if abd_switch eq 0 then $ spfile = solar_spectrum_file ; input by users endelse ; restgen, te_spectrum, lambda_spectrum, spectrum, file = spfile ; help, spfile, te_spectrum, lambda_spectrum, spectrum ; xrt_teem, index1, data1, index2, data2, te, em, et, ee, $ lambda_spectrum = lambda_spectrum, te_spectrum = te_spectrum, $ spectrum = spectrum, $ ; bin = bin, trange = trange, no_threshold = no_threshold, $ te_err_threshold = te_err_threshold, $ photon_noise_threshold = photon_noise_threshold, $ mask1 = mask1, mask2 = mask2, eff_area = eff_area ; end