GravityWaveAnalysis/miningRegionalForLearning.m at main · arthurMedforth/GravityWaveAnalysis · 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
function [ind_data,GW] = miningRegionalForLearning(region_data,data_included)

lats = [];
lons = [];
A = [];
wave_mag = [];
wave_angle = [];
u_137 = [];
u_100 = [];
u_50 = [];
u_23 = [];
v_137 = [];
v_100 = [];
v_50 = [];
v_23 = [];
t_137 = [];
precip = [];
switch data_included
    case 'zonal mean'
        for day = 1:length(region_data)
            lats = [lats;region_data(day).data(:,1)];

            A = [A;region_data(day).data(:,2)];
            wave_mag = [wave_mag;region_data(day).data(:,3)];
            wave_angle = [wave_angle;region_data(day).data(:,4)];

            u_137 = [u_137;region_data(day).data(:,6)];
            u_100 = [u_100;region_data(day).data(:,7)];
            u_50 = [u_50;region_data(day).data(:,8)];
            u_23 = [u_23;region_data(day).data(:,9)];

            v_137 = [v_137;region_data(day).data(:,10)];
            t_137 = [t_137;region_data(day).data(:,11)];
            v_100 = [v_100;region_data(day).data(:,12)];
            v_50 = [v_50;region_data(day).data(:,13)];
            v_23 = [v_23;region_data(day).data(:,14)];
            precip = [precip;region_data(day).data(:,15)];


        end
        GW = [A,wave_mag,wave_angle];
        ind_data = [lats,u_137,u_100,u_50,u_23,v_137,v_100,v_50,v_23,t_137,precip];

    case 'time average'
        for day =1:length(region_data)
            lats = [lats;region_data(day).data(:,1)];
            lons = [lons;region_data(day).data(:,2)];

            A = [A;region_data(day).data(:,3)];
            wave_mag = [wave_mag;region_data(day).data(:,4)];
            wave_angle = [wave_angle;region_data(day).data(:,5)];

            u_137 = [u_137;region_data(day).data(:,7)];
            u_100 = [u_100;region_data(day).data(:,8)];
            u_50 = [u_50;region_data(day).data(:,9)];
            u_23 = [u_23;region_data(day).data(:,10)];

            v_137 = [v_137;region_data(day).data(:,11)];
            t_137 = [t_137;region_data(day).data(:,12)];
            v_100 = [v_100;region_data(day).data(:,13)];
            v_50 = [v_50;region_data(day).data(:,14)];
            v_23 = [v_23;region_data(day).data(:,15)];
            precip = [precip;region_data(day).data(:,16)];

        end
        GW = [A,wave_mag,wave_angle];
        ind_data = [lats,lons,u_137,u_100,u_50,u_23,v_137,v_100,v_50,v_23,t_137,precip];

end