Function: thm_fgm_find_shift_1p1c Description: Returns the vector to be added to l1 data to correct the '390 nT' jumps for one probe and one component. Input: b = the uncorrected measurements threshold = threshold matrix (probe x component). A jump occurs when the field crosses this value. Note that the actual threshold depends on the field derivative. The value given here is for constant field (th0). th_slope = slope of the threshold: threshold=th0+th_slope*(db/dt) jump = jump matrix (probe x component). The expected values of the jumps. datatype= 'fgl' or -> low resolution (4Hz) 'fgh' or -> high resolution (128 Hz) 'fge' -> engeneering data (8 Hz) Keywords: /fixed: take the correction from the jump matrix. Default for fgh is to compute the correction for each jump such way that the first derivative is constant across the jump from low to high values. It has no effect for fge or fgl data for which the correction is taken by default from the jump matrix. /constant: average the correction over all jumps for fgh. It has no effect for fge or fgl. Return value: the correction vector Example: corr = thm_fgm_find_shift_1p1c(b, th_up, th_dn, jump, /constant, datatype='fgh') Notes - this function is used by fgm_find_shift Written by Dragos Constantinescu Minor modifications by Patrick Cruce for incorporation into TDAS. Specifically: 1. Changed Function Name 2. Disabled output of statistics to file. $LastChangedBy: aaflores $ $LastChangedDate: 2012-02-22 11:30:33 -0800 (Wed, 22 Feb 2012) $ $LastChangedRevision: 9811 $ $URL: svn+ssh://thmsvn@ambrosia.ssl.berkeley.edu/repos/spdsoft/trunk/projects/themis/spacecraft/fields/thm_fgm_find_shift_1p1c.pro $