Magnetic helicity and energy injection rates

Name:

helicity_energy_bvec.pro

Purpose:

Calculates magnetic helicity and energy injection rates through the photosphere of active regions.

Calling sequence:

output_helicity_energy=helicity_energy_bvec(br_t0, bt_t0, bp_t0, idx_t0, br_t1, bt_t1, bp_t1, idx_t1, binxy_helicity_energy=binxy_helicity_energy, set_mpil=set_mpil, binxy_mpil=binxy_mpil,

extend=extend, fwhm_align=fwhm_align, fwhm_dave4vm=fwhm_dave4vm, mpil_info=mpil_info, mpil_xy=mpil_xy)

Input parameters:

br_t0 - hmi_sharp_cea_720s Br image at t0 (note that t1 > t0)

bt_t0 - hmi_sharp_cea_720s Bt image at t0

bp_t0 - hmi_sharp_cea_720s Bp image at t0

idx_t0 - hmi_sharp_cea_720s meta data at t0

br_t1 - hmi_sharp_cea_720s Br image at t1

bt_t1 - hmi_sharp_cea_720s Bt image at t1

bp_t1 - hmi_sharp_cea_720s Bp image at t1

idx_t1 - hmi_sharp_cea_720s meta data at t1

Input Keywords:

binxy_helicity_energy - number for xy spatial binning of the sharp magnetic field input data related to calculation of magnetic helicity and energy injection rates [in pixels] ; default=1 (i.e., no binning)

set_mpil - if set, calculate magnetic helicity and energy 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 and energy injection rates therein; default=5

fwhm_align - full width at half maximum for sharp image alignment [arcsecs] ; default=10

fwhm_dave4vm - full width at half maximum for flow field calculation [arcsecs] ; default=4

mpil_info - information on types of MPILs

mpil_xy - x and y coordinates of MPILs [pixels]

Output:

output_helicity_energy - a structure having 22 tags related to helicity and energy injection properties

pos_dhdt_in: postive dhdt_in [10^35 Mx^2/sec]

abs_neg_dhdt_in: absolute of negative dhdt_in [10^35 Mx^2/sec]

abs_tot_dhdt_in: absolute of total dhdt_in [10^35 Mx^2/sec]

tot_uns_dhdt_in: total unsigned dhdt_in [10^35 Mx^2/sec]

pos_dhdt_sh: postive dhdt_sh [10^35 Mx^2/sec]

abs_neg_dhdt_sh: absolute of negative dhdt_sh [10^35 Mx^2/sec]

abs_tot_dhdt_sh: absolute of total dhdt_sh [10^35 Mx^2/sec]

tot_uns_dhdt_sh: total unsigned dhdt_sh [10^35 Mx^2/sec]

abs_tot_dhdt: absolute of total dhdt [10^35 Mx^2/sec]

abs_tot_dhdt_in_plus_sh: abs_tot_dhdt_in + abs_tot_dhdt_sh [10^35 Mx^2/sec]

tot_uns_dhdt: total unsigned dhdt [10^35 Mx^2/sec]

pos_dedt_in: postive dedt_in [10^24 erg/sec]

abs_neg_dedt_in: absolute of negative dedt_in [10^24 erg/sec]

abs_tot_dedt_in: absolute of total dedt_in [10^24 erg/sec]

tot_uns_dedt_in: total unsigned dedt_in [10^24 erg/sec]

pos_dedt_sh: postive dedt_sh [10^24 erg/sec]

abs_neg_dedt_sh: absolute of negative dedt_sh [10^24 erg/sec]

abs_tot_dedt_sh: absolute of total dedt_sh [10^24 erg/sec]

tot_uns_dedt_sh: total unsigned dedt_sh [10^24 erg/sec]

abs_tot_dedt: absolute of total dedt [10^24 erg/sec]

abs_tot_dhdt_in_plus_sh: abs_tot_dedt_in + abs_tot_dedt_sh [10^24 erg/sec]

tot_uns_dedt: total unsigned dedt [10^24 erg/sec]

Property name in DB:

helicity_energy_bvec.pos_dhdt_in, helicity_energy_bvec.abs_neg_dhdt_in, helicity_energy_bvec.abs_tot_dhdt_in, helicity_energy_bvec.tot_uns_dhdt_in, helicity_energy_bvec.pos_dhdt_sh, helicity_energy_bvec.abs_neg_dhdt_sh, helicity_energy_bvec.abs_tot_dhdt_sh, helicity_energy_bvec.tot_uns_dhdt_sh, helicity_energy_bvec.abs_tot_dhdt, helicity_energy_bvec.abs_tot_dhdt_in_plus_sh, helicity_energy_bvec.tot_uns_dhdt, helicity_energy_bvec.pos_dedt_in, helicity_energy_bvec.abs_neg_dedt_in, helicity_energy_bvec.abs_tot_dedt_in, helicity_energy_bvec.tot_uns_dedt_in, helicity_energy_bvec.pos_dedt_sh, helicity_energy_bvec.abs_neg_dedt_sh, helicity_energy_bvec.abs_tot_dedt_sh, helicity_energy_bvec.tot_uns_dedt_sh, helicity_energy_bvec.abs_tot_dedt, helicity_energy_bvec.abs_tot_dedt_in_plus_sh, helicity_energy_bvec.tot_uns_dedt