;-----------------------------
fname='teem_amap_20230219'
;-----------------------------
;--- read XRT qsbp info structure ---
in_dir='/disk/yla/takeda/REU_2026/qsbp_teem/'
restgen,file=in_dir+fname,teem_xrt
yyyymmdd=strmid(fname,10,8)
file_am=teem_xrt.file_am
file_ap=teem_xrt.file_ap
ofrac=teem_xrt.ofrac
fctm_am=teem_xrt.fctm_am
fctm_ap=teem_xrt.fctm_ap
teem_fov=teem_xrt.teem_fov
teem_qs=teem_xrt.teem_qs
teem_bps=teem_xrt.teem_bps
nn_bp=n_elements(teem_bps(0,*))
;
te_arr=[teem_fov(0),teem_qs(0)]
em_arr=[teem_fov(1),teem_qs(1)]
for i=0,nn_bp-1 do begin
   te_arr=[te_arr,teem_bps(0,i)]
   em_arr=[em_arr,teem_bps(1,i)]
endfor
ndata=n_elements(te_arr)
irrad_arr=dblarr(ndata)
;
;*** calculate irradiance  ***
AU = 1.495979d13                  ; (cm)
fact1 = 1./(4.*!pi*AU^2)   ; turned out OK per consideration on 12-Mar-2019.
; AU_m = 1.495979d11              ; (m)
; fact2 = 1./(AU_m^2)
;
;--- irrad set ver. 7 ---
;wbins1=[1.,8.]     ; GOES/XRS
;wbins2=[0.3,10.]
wbins3=[2.8,36.6] ; Yohkoh/SXT with Pevtsov_2003
;wbins4=[3.,30.]    ; Yohkoh/SXT
;wbins5=[0.4,20.]   ; MinXSS
;wbins6=[5.,60.]    ; XRT narrower
;wbins7=[4.,200.]   ; XRT wider
ab_num=1   ; Chanti ver.11.02, chianti_2021 + adv ioneq(@3-oct-2025)
;
for i=0,ndata-1 do begin
  te=te_arr(i) & em=em_arr(i)
  rad3=ch1102_spec_fast(10d^te,wbins3,/erg,ab_ion=ab_num) ; erg/s per em=1e44
  work=alog10(rad3)-44.+em
  irrad_arr(i)=alog10(fact1)+work-3.        ;watt m^-2
endfor
;
help,irrad_arr
print,irrad_arr
;
;--- make final info structure ---
irrad_fov=irrad_arr(0)
irrad_qs=irrad_arr(1)
irrad_bps=irrad_arr(2:*)
;
irrad_xrt=create_struct('file_am',file_am,'file_ap',file_ap, $
              'ofrac',ofrac, 'fctm_am',fctm_am, 'fctm_ap',fctm_ap, $
                        'irrad_fov',irrad_fov, $
                        'irrad_qs',irrad_qs, $
                        'irrad_bps',irrad_bps)
help,irrad_xrt,/str
;
yn='n'
read,prompt='*** save info_structure? [y/n] ',yn
if yn eq 'y' then begin $
  outdir='/disk/yla/takeda/REU_2026/qsbp_irrad/'
  fname='irrad_amap_'+yyyymmdd
  text='Created by TEEM2IRRAD_QSBP.PRO. '+ $
       'A structure contains signals of '+ $
       'file_am(Al_mesh filename), file_ap(Al_poly filename), '+ $
       'ofrac(prefilter open fraction assumed), '+ $
       'fctm_am(Alms filter contam. ratio assumed), '+ $
       'fctm_ap(Apol filter contam. ratio assumed), '+ $
       'irrad_fov(Irradiance [2.8, 36.6]A for FOV data '+ $
       'irrad_qs(Irradiance [2.8, 36.6]A for QS data '+ $
       'irrad_bps(Irradiance [2.8, 36.6]A for BPs data '+ $
       'irradiance was calculated with Chianti ver. 11.02, '+ $
       ' Chianti_coronal elemental abundance in the unit of W/m^2.'
  name=['irrad_xrt']
  savegen,irrad_xrt,file=outdir+fname,text=text,name=name
endif
;
FIN:
END