Name:
helicity_br.pro
Purpose:
Calculates a proxy of magnetic helicity injection rate through the photosphere of active regions.
Calling sequence:
output_helicity=helicity_br(br_t0, idx_t0, br_t1, idx_t1, binxy_helicity=binxy_helicity, set_mpil=set_mpil, binxy_mpil=binxy_mpil, $
extend=extend, fwhm_align=fwhm_align, fwhm_dave=fwhm_dave, mpil_info=mpil_info, mpil_xy=mpil_xy)
Input parameters:
br_t0 - hmi_sharp_cea_720s Br image at t0 (note that t1 > t0)
idx_t0 - hmi_sharp_cea_720s meta data at t0
br_t1 - hmi_sharp_cea_720s Br image at t1
idx_t1 - hmi_sharp_cea_720s meta data at t1
Input Keywords:
binxy_helicity - number for xy spatial binning of the sharp magnetic field input data related to calculation of magnetic helicity injection rates [in pixels] ; default=1 (i.e., no binning)
set_mpil - if set, calculate magnetic helicity injection rates only near magnetic polarity inversion lines (MPILs) ; default=0
binxy_mpil - number for xy spatial binning of sharp magnetic field input data related to MPIL detection [pixels] ; default=1
extend - length of a square box surrounding each of MPIL pixels [Mm]. The box region will be considered for calculating magnetic helicity injection rates therein; default=5
fwhm_align - full width at half maximum for sharp image alignment [arcsecs] ; default=10
fwhm_dave - full width at half maximum for flow field calculation [arcsecs] ; default=10
mpil_info - information on types of MPILs
mpil_xy - x and y coordinates of MPILs [pixels]
Output:
output_helicity - a structure having 11 tags related to helicity injection properties
pos_dhdt_s: postive dhdt_s [10^35 Mx^2/sec]
abs_neg_dhdt_s: absolute of negative dhdt_s [10^35 Mx^2/sec]
abs_tot_dhdt_s: absolute of total dhdt_s [10^35 Mx^2/sec]
tot_uns_dhdt_s: total unsigned dhdt_s [10^35 Mx^2/sec]
pos_dhdt_m: postive dhdt_m [10^35 Mx^2/sec]
abs_neg_dhdt_m: absolute of negative dhdt_m [10^35 Mx^2/sec]
abs_tot_dhdt_m: absolute of total dhdt_m [10^35 Mx^2/sec]
tot_uns_dhdt_m: total unsigned dhdt_m [10^35 Mx^2/sec]
abs_tot_dhdt: absolute of total dhdt [10^35 Mx^2/sec]
abs_tot_dhdt_s_plus_m: abs_tot_dhdt_s + abs_tot_dhdt_m [10^35 Mx^2/sec]
tot_uns_dhdt: total unsigned dhdt [10^35 Mx^2/sec]
Property name in DB:
helicity_br.pos_dhdt_s, helicity_br.abs_neg_dhdt_s, helicity_br.abs_tot_dhdt_s, helicity_br.tot_uns_dhdt_s, helicity_br.pos_dhdt_m, helicity_br.abs_neg_dhdt_m, helicity_br.abs_tot_dhdt_m, helicity_br.tot_uns_dhdt_m, helicity_br.abs_tot_dhdt, helicity_br.abs_tot_dhdt_s_plus_m, helicity_br.tot_uns_dhdt
Repository
https://dev.flarecast.eu/stash/projects/FE/repos/helicity_br-idl/browse