;+
; Batch File: THM_CRIB_SUPERPOSED_INTERPOL_HISTO
;
; Purpose: Demonstrates the application of the routines
; superpo_interpol and superpo_histo
;
; Note: The examples in this crib sheet call 'superpo_interpol'. Simply
; replace the calls with 'superpo_histo' to try the other routine.
; Both routines accept the same input parameters.
;
; Written by Andreas Keiling
;
; $LastChangedBy: $
; $LastChangedDate: $
; $LastChangedRevision: $
; $URL $
;-
;;;;;;;;;;;;;;;;;EXAMPLE #1;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
del_data,'*'
thm_init
timespan,'2007-03-23/00:00:00' ; one whole day by default
thm_load_gmag,site='ccnv'
thm_load_gmag,site='drby'
thm_load_gmag,site='fsim'
thm_load_gmag,site='fsmi'
thm_load_gmag,site='fykn'
split_vec,'thg_mag_ccnv'
split_vec,'thg_mag_drby'
split_vec,'thg_mag_fsim'
split_vec,'thg_mag_fsmi'
split_vec,'thg_mag_fykn'
superpo_interpol,'thg_mag_ccnv_x thg_mag_drby_x thg_mag_fsim_x thg_mag_fsmi_x thg_mag_fykn_x', $
min='thg_pseudoAL', $
max='thg_pseudoAU', $
dif='thg_pseudoAE', $
avg='thg_avg', $
med='thg_median', $
res=30.0
;superpo_histo,'thg_mag_ccnv_x thg_mag_drby_x thg_mag_fsim_x thg_mag_fsmi_x thg_mag_fykn_x', $
; min='thg_pseudoAL', $
; max='thg_pseudoAU', $
; dif='thg_pseudoAE', $
; avg='thg_avg', $
; med='thg_median', $
; res=30.0
options,'*_x thg_pseudoAU thg_pseudoAL thg_pseudoAE thg_avg thg_median', /ynozero
tplot,'*_x thg_pseudoAU thg_pseudoAL thg_pseudoAE thg_avg thg_median'
stop
;;;;;;;;;;;;;; EXAMPLE #2 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
del_data,'*'
thm_init
timespan,'2006-12-23/00:00:00' ; one whole day by default
thm_load_gmag,/subtract_median ; load all stations and subtract median from each station
split_vec,'thg_mag_????'
superpo_interpol,'thg_mag_????_x',min='thg_pseudoAL', res=600.0 ; use all available stations
;superpo_histo,'thg_mag_????_x',min='thg_pseudoAL', res=600.0 ; use all available stations
tplot,'*pseudo*'
stop
;;;;;;;;;;;;;;;;;EXAMPLE #3;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
del_data,'*'
thm_init
timespan,'2006-12-23/00:00:00',18,/hour
thm_load_gmag,/subtract_median ; load all stations and subtract median from each station
split_vec,'thg_mag_????'
superpo_interpol,'thg_mag_????_x', res=1.0 ; does default values for all keywords except res
;superpo_histo,'thg_mag_????_x', res=1.0 ; does default values for all keywords except res
tplot,'*arr*'
stop
;;;;;;;;;;;;;;;;;EXAMPLE #4;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
del_data,'*'
thm_init
timespan,'2006-12-23/00:00:00',18,/hour
thm_load_gmag ; load all stations
split_vec,'thg_mag_????'
superpo_interpol,'thg_mag_????_x' ; does default values for all keywords
;superpo_histo,'thg_mag_????_x' ; does default values for all keywords
options,'thg_mag_????_x *arr*', /ynozero
tplot,'thg_mag_????_x *arr*'
end