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
|
function [lon,lat,var,var2,lon2,lat2,d,m] = extract (y,m,d)
isOk=1;
isok = 1;
root = '/home9/begmeil/cersat/ftp/products/gridded/sst-pathfinder/netcdf/v5-all';
rootIma = '/home/typhon/stagiaires/jaccart/SST';
%----TEST----
if~exist(root,'dir')
isOk=0;
error(sprintf('ERROR, invalid directory: %s',root))
elseif~exist(rootIma,'dir')
isOk=0;
error(sprintf('ERROR, invalid directory: %s',rootIma))
end
%--NOM-DU-FICHIER--
while m<=10
while d<=31
%--nom-du-fichier--
fileName = sprintf('%.4d%.2d%.2d.04d2pfv50-all.nc.bz2',y,m,d);
fileName = fullfile(root,sprintf('%.4d',y),sprintf('%.2d',m),fileName);
%--Existance-du-fichier--
if~exist(fileName,'file')
isok = 0
end
%----Lecture----
if isok
disp(fileName);
nc = cl_netcdf_writer;
nc = set_file_name (nc,fileName);
s = read (nc);
var = get_var_value (nc,s,'sst',1,1,1);
lat = get_var_value (nc,s,'lat',1,1,1);
lon = get_var_value (nc,s,'lon',1,1,1);
var2 = var(lat>=-15 & lat<=30, lon>=-80 & lon<=20);
end
%----Sauvegarde-des-variables----
if isok
Name = fullfile(rootIma,sprintf('%.4d',y),sprintf('%.4d.%.2d.%.2d.mat',y,m,d));
save(Name, 'var2');
d = d+1;
end
end
m = m+1
d=1
end
end |
