function floatstr_edit_setup(h, options)
    % Takes the handle of an edit_text, and sets up callbacks such that it
    % becomes a float_edit box, with arrowup/arrowdown functionality.
    % Needs an options struct, which is attached to h as user_data. Must
    % contain fields:
    %  - .format       Format of the floating point number
    %  - .min          Minimum number
    %  - .max          Maximum Number
    %  - .default_num  Default Number (for bad input)
    %  - .default_str  Default String (for bad input)

    set(h, 'UserData', options) 
    set(h, 'FontName', 'Courier', 'FontSize', 10)
    % Sets callback, calls it to make sure string is in suitable format
    set(h,'Callback',@editfloat_Callback);
    editfloat_Callback(h, 0)

    % Note: 'persistent' fails when the GUI is closed and started again. You need to
    % restart Matlab in this case. I stopped using 'persistent'.
    % Sets up Callback of underlying Java Swing object.
    jh = findjobj(h, 'nomenu'); % external function by Yair Altman.
    set(jh,'KeyPressedCallback', ...
        {@editfloat_KeyPressedCallback_Java, h});

    function editfloat_Callback(hObject,~)
    % Called on "commit" event ("Enter" and "LostFocus" or so). 
    % Parses the element's input and updates the string with parsed input.

    edit_str = get(hObject, 'String');
    user_data = get(hObject, 'UserData');
    [~, clean_str] = editfloat_str_parser(edit_str, user_data);
    set(hObject, 'String', clean_str)

    function editfloat_KeyPressedCallback_Java(jhObject, event, hObject)
    % Called on each button press inside edit text. 
    % Increments/decrements the value at the current caret position.
    % Needs the following fields in the "user_data" struct:
    %  - .format       Format of the floating point number
    %  - .min          Minimum number
    %  - .max          Maximum Number
    %  - .default_num  Default Number (for bad input)
    %  - .default_str  Default String (for bad input)

    keynum = get(event, 'keyCode'); 

    if keynum==38 || keynum==40  % arrow up or arrow down
        user_data = get(hObject, 'UserData');
        % format = user_data.format;
        caret_position = get(jhObject, 'CaretPosition') + 1; %matlab indexing
        current_str = get(jhObject, 'Text'); %using jhObject, because during edition, hObject contains old stuff
        [old_num, parsed_str] = editfloat_str_parser(current_str, user_data);
        
        if caret_position > length(parsed_str)
            % we're at the end of the string. do nothing
            return
        end
        
        if length(parsed_str) ~= length(current_str)
            % User was during input. Update field, return
            % Update: In the current version, this is impossible, I think
            set(hObject, 'String', parsed_str) 
            drawnow;
            return
        end
        
        % Prepare a delta_str -> delta_num from the given string and the caret
        % position. 
        numeric_period = int16('.');
        delta_str = parsed_str;
        delta_str(parsed_str ~= numeric_period) = '0';  % -> 00000.00000 or so
        
        new_num = old_num;
        if ~strcmp(delta_str(caret_position), '.')
            delta_str(caret_position) = '1';          % -> 00100.00000 or so
            delta_num = str2double(delta_str);
            if keynum == 38 % up
                new_num = new_num + delta_num;
            elseif keynum == 40 %down
                new_num = new_num - delta_num;
            end
        end
        
        new_str = sprintf(user_data.format, new_num);
        [~, new_str] = editfloat_str_parser(new_str, user_data);
        set(hObject, 'String', new_str) %using hObject here, because jhObject led to weird errors
        drawnow;   % otherwise set() takes place after the setCaretPosition!
        
        % Set correct caret position again. Adjust for possibly new string
        % length.
        length_difference = length(new_str) - length(current_str);
        jhObject.setCaretPosition(caret_position - 1+length_difference);
        drawnow;  % I was getting more errors. Maybe this helps?
    end

function [num, clean_str] = editfloat_str_parser(float_str, options)
% Converts float_str to a number. Returns this number as
% numeric and as a string.
%
% The input "options" must be a struct that contains the fields:
%  - .format       Output format of the number string.
%  - .min          Minimum number
%  - .max          Maximum Number
%  - .default_num  Default Number (for bad input)
%  - .default_str  Default String (for bad input)

% Example output_format: '%016.6f'
%   - 0 padded to field length
%   - 16 field length (including decimal point! -> 15 digits)
%   - 6 digits after the decimal points
%   - floating point


no_commas = strrep(float_str, ',', '.'); % having "," as decimal is fine
num = str2double(no_commas);
if isnan(num)  % if contained text
    num = options.default_num;
    clean_str = options.default_str;
    return
end
num = real(num); % only keep the real part
if num < options.min     num = options.min; elseif num > options.max
    num = options.max;
end
clean_str = sprintf(options.format, num);

% --- Outputs from this function are returned to the command line.
function varargout = Advanced_THEs_Plot_GUI_OutputFcn(hObject, eventdata, handles) 
% varargout  cell array for returning output args (see VARARGOUT);

% This function is executed just after the GUI is made visible. All the
% java swing object stuff needs the object to be visible, which is why the
% following code bits are here.
% Note: 
% Makes edittext into editfloats by setting up the right callbacks (for
% arrowup/arrowdown functionality)

% FreqOffset field
editfloat_options.format = '%12.3f';
editfloat_options.min = 0;
editfloat_options.max = 100e6-1e-3;
editfloat_options.default_num = 0;
editfloat_options.default_str = '0.000';
editfloat_setup(handles.FreqOffset_edittext, editfloat_options)

%% Link two plots together so that the x-axis always shows the same range.
figure(1); ah(1) = gca; figure(2); ah(2) = gca; linkaxes(ah, 'x')

%% Copy lines from one figure to another
L = findobj(2,'type','line'); % first argument is the figure to copy from
for i = 1:length(L)
    copyobj(L(i),findobj(3,'type','axes')); %first arg is figure to copy to
end
% By hand: Copy/Paste works too

%% Make line invisible to "legend" command. Mark line first!
hasbehavior(gco, 'legend', false);

%% Do Allan Deviation
xdata = get(get(gca,'Children'),'xdata');
% xdata = (xdata - xdata(1))*24*3600; % in s from env data
ydata = get(get(gca,'Children'),'ydata');
ydata = ydata*200e-6/20/87.6; % relative deviation
delta_t = mean(diff(xdata));
[~,ADEV,~,~,~,TAU] = avar(ydata,delta_t); 
figure
loglog(TAU,ADEV,'x-')
set(gca,'YLim',[1e-8,1e-7])

%% Linear Fit with Fit uncertainty
% Fit *column* vectors x,y with linear polynomial y = a*x + b
% Remember to remove mean(x) first, if you want he errors to be
% uncorrelated!
cf = fit(x,y,'poly1'); 

cf_coeff = coeffvalues(cf);
cf_confint = confint(cf);
a = cf_coeff(1);
b = cf_coeff(2);
a_uncert = (cf_confint(2,1) - cf_confint(1,1))/2;
b_uncert = (cf_confint(2,2) - cf_confint(1,2))/2;

%% Smoothing with moving average
smoothing_length = 10;
filtered_data = filter(ones(1,smoothing_length)/smoothing_length,1,data);

%% Smoothing with butterworth filter
cutoff_freq = 0.02; % in Hz - this is a 50 second filter
sample_freq = 1/(mean(diff(time_data))*24*3600);
[A,B] = butter(2,cutoff_freq/(sample_freq/2));
filtered_data = filtfilt(A,B,data);