-
Notifications
You must be signed in to change notification settings - Fork 2
/
HM_Step_01_02_Run_Pairs_dup.m
48 lines (40 loc) · 2.07 KB
/
HM_Step_01_02_Run_Pairs_dup.m
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
% A scipt that pick out ICOADS pairs and
% subset pairs by distance in each month
HM_load_package;
% *******************
% Fixed Parameters **
% *******************
varname = 'SST';
method = 'Bucket';
pick_limit = 1850;
env = 1;
[y_st,m_st] = ind2sub([50,12],num);
for yr = (1850+y_st):50:2014
for mon = m_st
% ######################################################################
% The following parameters are for results shown in Chan et al., 2019
% They can be updated in various ways to perform sensitivity tests
% ######################################################################
if 1, % Pick out all possible pairs
do_NpD = 1; % Group by nation and deck
do_rmdup = 0; % Do not remove DUPS != 0
do_other = 0; % Do not combine JP Kobe collections
% at this step, but they will be combined
% in later steps
do_rmsml = 0; % Do not remove small contributing groups
HM_pair_01_Raw_Pairs_dup...
(yr,mon,do_NpD,varname,method,do_rmdup,do_rmsml,env,do_other);
end
% ######################################################################
% The following parameters are for results shown in Chan et al., 2019
% They can be updated in various ways to perform sensitivity tests
% ######################################################################
if 1,
do_NpD = 1; % Group by nation and deck
P.do_rmdup = 0; % Do not remove DUPS != 0
P.do_rmsml = 0; % Do not remove small contributing groups
P.do_fewer_first = 0; % Screen pairs with ascending distance
HM_pair_02_Screen_Pairs_dup(yr,mon,do_NpD,varname,method,env,P);
end
end
end