Easy_NetCDF/FUN_nc_gen_time_val_limit.m at main · HappySpring/Easy_NetCDF · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
function [ time_limit_in_ori, time_mtl ]= FUN_nc_gen_time_val_limit( filename0, time_var_str, timelimit, varargin )
% [ time_limit_in_ori ]= FUN_nc_gen_time_val_limit( filename0, time_var_str, timelimit )
% calculate time limit in the same unit as used in the netcdf file
% -------------------------------------------------------------------------
% INPUT
%      filename0: netcdf file path
%      time_var_str: name of time variable
%      timelimit: time limit in matlab format
% -------------------------------------------------------------------------
% OUTPUT
%       time_limit_in_ori
%       time_mtl
% -------------------------------------------------------------------------
%
% V1.10 by L. Chi: support loading time series from pre-saved .mat file
%                   [option: presaved_total_timeseries]
% V1.01 By L. Chi: Improve performance
% V1.00 By L. Chi (L.Chi.Ocean@outlook.com)
% -------------------------------------------------------------------------
%%

[presaved_total_timeseries, varargin] =  FUN_codetools_read_from_varargin( varargin, 'presaved_total_timeseries', [], true );

% This is only used to verify the early version of the script. please keep it off
[is_double_test,            varargin] =  FUN_codetools_read_from_varargin( varargin, 'is_double_test', false, true );

[calendar_in,               varargin] =  FUN_codetools_read_from_varargin( varargin, 'calendar_in', [], true );

if ~isempty(varargin)
    error('unknown input parameter!')
end


% load time series
is_load_from_local = false;

if ~isempty(presaved_total_timeseries)
    url_ind = FUN_struct_find_field_ind(presaved_total_timeseries,'filename0',filename0);

    if ~isnan(url_ind)
        is_load_from_local = true;
    end

end

if is_load_from_local
    time_in_ori = presaved_total_timeseries(url_ind).time_ori;
    time_mtl    = presaved_total_timeseries(url_ind).time_mtl;

else

    [time_mtl, time_in_ori] = FUN_nc_get_time_in_matlab_format( filename0,  time_var_str, 'calendar_in', calendar_in );


end

% calculate time limit in the original format of the given netcdf file
    time_ind = find( time_mtl >= timelimit(1) & time_mtl <= timelimit(2) );

    if isempty( time_ind )
        disp('Data are not available within the required range!')
        time_limit_in_ori = [nan nan];
        time_mtl = [];
        return
    end

    time_limit_in_ori = [ min( time_in_ori( time_ind) ),  max( time_in_ori( time_ind) ) ] ;
    time_mtl = time_mtl( time_ind );

% test results
% useless and to be deleted later
if is_double_test
    time_test = FUN_nc_varget_enhanced_region_2( filename0, time_var_str, {time_var_str}, {time_limit_in_ori} );

    if isequal(time_test.time, time_in_ori(time_ind))
    else
        error('Failed in test, Unexpected values!');
    end

end