-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathread0
63 lines (52 loc) · 1.81 KB
/
read0
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
%% read all data in files....
filenum=29 %constant... lazy
for i=1:filenum
g=i-1;
file = ['Yukon_Jason2_' num2str(g) '.nc'];
finfo= ncinfo(file);
Global(i).lat=ncread(file,'lat');
Global(i).lon=ncread(file,'lon');
Global(i).ID=ncread(file,'ID');
Global(i).sat=ncread(file,'sat');
Global(i).pass=ncread(file,'pass');
Timeseries(i).time=ncread(file,'/Timeseries/time');
Timeseries(i).cycle=ncread(file,'/Timeseries/cycle');
Timeseries(i).hbar=ncread(file,'/Timeseries/hbar');
Timeseries(i).hwbar=ncread(file,'/Timeseries/hwbar');
Timeseries(i).sig0bar=ncread(file,'/Timeseries/sig0bar');
Timeseries(i).pkbar=ncread(file,'/Timeseries/pkbar');
% Sampling.scene
% Sampling.lonbox
% Sampling.latbox
% Sampling.island
Level3(i).lon=ncread(file,'/Level3/lon');
Level3(i).lat=ncread(file,'/Level3/lat');
Level3(i).h=ncread(file,'/Level3/h');
Level3(i).sig0=ncread(file,'/Level3/sig0');
Level3(i).pk=ncread(file,'/Level3/pk');
Level3(i).cycle=ncread(file,'/Level3/cycle');
Level3(i).time=ncread(file,'/Level3/time');
Level3(i).heightfilter=ncread(file,'/Level3/heightfilter');
Level3(i).icefilter=ncread(file,'/Level3/icefilter');
Level3(i).allfilter=ncread(file,'/Level3/allfilter');
Filter(i).nNODATA=ncread(file,'Filter/nNODATA');
Filter(i).riverh=ncread(file,'Filter/riverh');
Filter(i).maxh=ncread(file,'Filter/maxh');
Filter(i).minh=ncread(file,'Filter/minh');
Filter(i).icethaw=ncread(file,'Filter/icethaw');
Filter(i).icefreeze=ncread(file,'Filter/icefreeze');
end
Global
Timeseries
Level3
Filter
%% usgs data
file2 = ['USGS.nc'];
Ufinfo= ncinfo(file2);
title=Ufinfo.Attributes.Value
USGS.datum=ncread(file2,'datum');
USGS.rmse=ncread(file2,'rmse');
USGS.bias=ncread(file2,'bias');
USGS.SD=ncread(file2,'SD');
USGS.height=ncread(file2,'height');
USGS.time=ncread(file2,'time');