***** New Yohkoh analysis software for Te-EM analysis ***** (/disk/yla/www/ystorage/documents/sxt_resp/sre_files, as of 28-Jan-2018) ---------------------------------------------------- 1) sxt_teem.pro ---------------------------------------------------- [Recent mods.] ; 9-nov-2011, Aki Takeda, added /corona, /hybrid, /photospheric, ; /mewe, and sre_file keywords (in response to the ; creation of the new response functions). [Examples] --- use default sre file : 'sre_ch601_hybrid_chianti.genx' --- IDL> sxt_teem,index1,data1,index2,data2,te,em --- use 'sre_ch601_corona_chianti.genx' --- IDL> sxt_teem,index1,data1,index2,data2,te,em,/corona --- use 'sre_ch601_photos_chianti.genx' --- IDL> sxt_teem,index1,data1,index2,data2,te,em,/photos --- use former standard sre file, 'sre950419_02.genx' --- IDL> sxt_teem,index1,data1,index2,data2,te,em,/mewe --- specify user specified sre file --- IDL> sxt_teem,index1,data1,index2,data2,te,em, $ sre_file='./sre_ch601_corona_mazzotta.genx' ---------------------------------------------------- 2) sxt_teem1.pro ---------------------------------------------------- [Recent mods.] ; 9-nov-2011, Aki Takeda, Added /mewe, /corona, /hybrid, /photospheric, ; and sre_file keywords. This procedure is called in sxt_teem.pro and thus rare to be used independantly. (In other words, sxt_teem.pro is a wrapper of sxt_teem1.pro) Usage (calling sequence) is the same as sxt_teem.pro ---------------------------------------------------- 3) sxt_teem2.pro (changed name from sxt_teem3.pro) ---------------------------------------------------- [Recent mods.] ; 20-Feb-2010 Aki Takeda, - change the default sre file ; (Mewe to Chianti with coronal abundance). ; - added SRE_FILE keyword to accept user-specified ; response functions. ; - setting hard-wired variable in a separate ; routine, teem_hardwired.pro ; 06-Aug-2010 Aki T, - added UNC_P1, UNC_P2 keyword to pass photon noise ; calculated outside this program. ; 28-Sep-2010 Aki T, - update the path of the default sre file ; (ydb_beta -> $SSW_SXT/response). ; 18-Mar-2011 Aki T, - added /corona, /photospheric, /hibrid ; and /mewe keywords. ; 27-Jan-2018 Aki T, - update response functions to chainti 800 version. ; also changed binding to abundance files. ; /hibrid : sun_coronal_2012_schmelz_ext.abund (from Fludra 1999) ; /photos : sun_photospheric_2011_caffau.abund (from Grevesse 2007) ; /corona : sun_coronal_1992_feldman_ext.abund (no change) [Examples] --- default sre file : 'sre_ch800_hybrid_chianti.genx' --- IDL> sxt_teem2, filt1, img1, filt2, img2, Te , EM --- use 'sre_ch800_corona_chianti.genx' --- IDL> sxt_teem2, filt1, img1, filt2, img2, Te , EM, /corona --- use 'sre_ch800_photos_chianti.genx' --- IDL> sxt_teem2, filt1, img1, filt2, img2, Te , EM, /photospheric --- use former standard sre file, 'sre950419_02.genx' --- IDL> sxt_teem2, filt1, img1, filt2, img2, Te , EM, /mewe --- specify user specified sre file --- IDL> sxt_teem2, filt1, img1, filt2, img2, Te , EM, $ sre_file='./sre_ch601_corona_mazzotta.genx' --------------------------------- 5) sxt_flux.pro --------------------------------- [Recent mods.] ; 2-jul-2010, Aki Takeda - Update default sre-file, Mewe --> Chianti. ; 14-jul-2010, AkT, Added the /mewe keyword to go back easily to Mewe sre. ; Name changed from SXT_FLUX2 to SXT_FLUX. ; 26-mar-2011, AkT, Change default CHIANTI sre-file (corona --> hybrid). ; 6-feb-2015, AkT, Change sre-files to those with CHIANTI ver. 7.14. ; 28-jan-2018, AkT, Change sre-files to those with CHIANTI ver. 8.00. [Examples] --- to get flux (DN/s) with the sre-file (with hybrid abundance: default) --- IDL> flux = sxt_flux(LogTe, Filt) --- to get flux (DN/s) with sre-file with coronal abund. --- IDL> flux = sxt_flux(LogTe, Filt, $ sre_file='$SSW_SXT/response/sre_ch800_corona_chianti.genx') --- to get flux (DN/s) with sre-file with photospheric abund. --- IDL> flux = sxt_flux(LogTe, Filt, $ sre_file='$SSW_SXT/response/sre_ch800_photos_chianti.genx') --- to get flux (DN/s) with the former sre-file (Mewe) --- IDL> flux = sxt_flux(LogTe, Filt, /mewe) --------------------------------- 6) chianti_spec.pro --------------------------------- A modified version of mewe_spec.pro and mewe_spec_lwa.pro. Compute a Chianti isothermal spectrum (line + continuum) in the specified wavelength bins. Calculation is made for EM=1.e44 cm^-3. [Recent mods.] ; 1-Jul-2010, Aki Takeda : Modified version of mewe_spec and ; mewe_spec_lwa.pro ; 3-Apr-2011, Aki T.: Default sre file, corona --> hybrid. ; 28-Jan-2018, Aki T.: Change abundance keyword binding. ; /hybrid : sun_coronal_2012_schmelz_ext.abund (from Fludra 1999) ; /photos : sun_photospheric_2011_caffau.abund (from Grevesse 2007) ; /corona : sun_coronal_1992_feldman_ext.abund (change only file name) Optional input keywords available: Ab_ion : 1 : sun_coronal_1992_feldman_ext.abund & chianti.ioneq : 2 (default): sun_coronal_2012_schmelz_ext.abund & chianti.ioneq : 3 : sun_photospheric_2011_caffau.abund & chianti.ioneq Abun : Full pass name for the abundance file. ex) Abun='$SSW/packages/chianti/dbase/abundance/'+ $ 'sun_coronal_1999_fludra_ext.abund' Ion : Ful pass name for the ionization equilbrium file. (If Ab_ion is set, Abun and Ioneq will be neglected.) [Examples] --- calculation with the default abund + ioneq set --- Flux = chianti_spec(temp,wave_edges,/photon) ; ph s-1 Flux = chianti_spec(temp,wave_edges,/erg) ; erg s-1 Flux = chianti_spec(5.0e6,[1,8,20,30],/erg) ; erg s-1 --- calculation with coronal abund. --- Flux = chianti_spec(temp,wave_edges,/erg, ab_ion=1) --- calculation with photospheric abund. --- Flux = chianti_spec(temp,wave_edges,/erg, ab_ion=3) --- calculation with hybrid abund.+ Mazzotta's ioneq. --- Flux = chianti_spec(temp,wave_edges,/erg, $ ion=!xuvtop+'/ioneq/mazzotta_etal.ioneq') --------------------------------- 7) teem_hardwired.pro --------------------------------- A batch file called from sxt_teem3, etc.,to set hardwired constants. ; 20-Feb-2010 (Aki Takeda) separated hard-wired part used to be defined ; in more than couple of different routines. ; 21-Jan-2015 (AkT) CCD gain updated according to Acton 2015 SolarPhys. ; ;----------------------------------------------------------------- em_assumed = 44. + 3 ; for /msec and EM=1.e44; HARD-WIRED *** label_id = [1, 2, 3, 6, 5, 4] ; HARD-WIRED *** thick_order = [1, 2, 3, 6, 5, 4] ; HARD-WIRED *** ; ;gain_default = 100. ; e/DN, gain at launch gain_default = 90. ; e/DN, gain at launch ;-----------------------------------------------------------------