killingobservers/LoadOceanData.m at master · vccvisualization/killingobservers · 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
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Copyright (C) 2018 King Abdullah University Of Science and Technology
%
% Contact:
% Peter Rautek peter.rautek@kaust.edu.sa
% Matej Mlejnek matej.mlejnek@kaust.edu.sa
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function domain = LoadOceanData(filename)

dataV = load(filename);

domain.tmin = dataV.time(1);
domain.tmax = domain.tmin+91;

% domain sizes
domain.w = 390;
domain.h = 210;
domain.d = 1;
domain.timeSteps = 14;

% support other common variable names
domain.xsize = domain.w;
domain.ysize = domain.h;
domain.zsize = domain.d;
domain.tsize = domain.timeSteps;

% actual domain of ocean data set:
% lon: -12..29 lat: -45..-19

% actual grid size:
% w: 165, h: 105, t: 14

% paper parameters of Haller 2016:
domain.xmin = -4;
domain.xmax = 9;
domain.ymin = -35;
domain.ymax = -28;
domain.zmin = 1;
domain.zmax = 1;

u_interp = griddedInterpolant({dataV.lon,dataV.lat,dataV.time},permute(dataV.UT,[2,1,3]),'linear','none');
v_interp = griddedInterpolant({dataV.lon,dataV.lat,dataV.time},permute(dataV.VT,[2,1,3]),'linear','none');

domain.u = @(x,y,t) u_interp(x,y,t);
domain.v = @(x,y,t) v_interp(x,y,t);