Name: THM_MAKE_AE Purpose: This routine calculates the "pseudo" AE, AL, and AU geomagnetic indices by using ground magnetometer data from THEMIS GBOs. The names of all stations used for calculation are printed on the screen. In future, it will be possible to include ground data from other magnetometer networks. Note that currently the calculation of the "pseudo" indices does not subtract quiet day variation but simply the median. Syntax: THM_MAKE_AE [, RES = float] [, SITES = string ] [, /NO_LOAD ] Parameters: None. Keywords: res = sampling interval (by default 60 sec) sites = observatory name; default is to use high-latitude THEMIS sites plus a few more: ['atha', 'chbg', 'ekat', 'fsim', 'fsmi', 'fykn', $ 'gako', 'gbay', 'gill', 'inuv', 'kapu', 'kian', $ 'kuuj', 'mcgr', 'pgeo', 'pina', 'rank', 'snap', $ 'snkq', 'tpas', 'whit', 'yknf', 'fcc', 'cmo', $ 'naq', 'lrv'] ;made an array to facilitate the use of split_vec later If set to 'all', all available sites will be loaded and used. no_load = if set, use existing gmag (THEMIS) tplot variables which have already been loaded into the active TDAS environment if not set, load gmag data (either remotely or from computer disk) max_deviation = The maximum deviation that the magnetic field data can go from the median; points that exceed this value are omitted. The default value is plus or minus 1500 nT Example: see crib sheet "thm_crib_make_AE.pro" Notes: Written by Andreas Keiling, 15 May 2008 Modifications: Edited header, put subroutine before main body so that it would compile before being called by the main body, W.M.Feuerstein, 6/2/2008. Changed default sites to be only high-lat THEMIS stations: jmm, 25-nov-2009 Added FCC, CMO, NAQ, LRV as sites, 17-sep-2012, jmm Added max deviation, extra despike of magnetic field prior to index calculation, 4-nov-2013, jmm $LastChangedBy: egrimes $ $LastChangedDate: 2014-01-07 08:31:19 -0800 (Tue, 07 Jan 2014) $ $LastChangedRevision: 13810 $ $URL $