Commit 8686e5ff authored by Paul Nicot's avatar Paul Nicot

ajout de # de spot + flag dans les fichiers output

parent 862f313a
......@@ -9,19 +9,18 @@ function tr_sort(num_spot,mission,debut,h1,fin,h2)
% -» mission = diminutif de la mission (3 caractères en capitale)
% -» debut = année, mois, jour de début de dérive (ex: 20150315)
% -» h1 = heure,min de début de dérive (hhmm)
% -» fin = année, mois, jour de fin de dérive
% -» h2 = heure,min de fin de dérive (hhmm)
% (facultatif, si non précisé -» données compilées jusqu'à maintenant)
% -» fin* = année, mois, jour de fin de dérive
% -» h2* = heure,min de fin de dérive (hhmm)
% (*facultatif, si non précisé -» données compilées jusqu'à maintenant)
%
% exemple : sort_spot([23 24],'BDC',20150612,1255,20150623,1835)
%
% Le fichier brute de données est situé sur Brandypot:
% /sas/usagers/share_lasso/data/MEOPAR/spots/spots_pos_clean.csv
% Le fichier brute de données est situé sur Mingan:
% /share/archives/partage_lasso/spot/spots_pos_flag_clean.csv
%
% Enregistrements sous forme de structure sur Brandypot:
% sas/usagers/share_lasso/data/MEOPAR/drifter
% Enregistrements sous forme de structure sur Mingan:
% /share/archives/partage_lasso/spot/drifter
%rootdir = '/sas/usagers/share_lasso/data/MEOPAR';
rootdir = '/share/archives/partage_lasso/spot';
fid1 = fopen([ rootdir '/spots_pos_flag_clean.csv']);
......@@ -29,11 +28,8 @@ header = textscan(fid1,'%s;%s;%s;%s;%s;%s',1);
spot_clean = textscan(fid1,'%f %f %f %f %s %s','Delimiter',';');
fclose(fid1);
spot = spot_clean{1,1};
lat = spot_clean{1,2}; lat(find(lat == -99999)) = NaN;
lon = spot_clean{1,3}; lon(find(lon == -99999)) = NaN;
nU = spot_clean{1,4}; % temps UNIX
nM = nU/86400 + datenum(1970,1,1); % temps matlab
% temps matlab
nM = spot_clean{1,4}/86400 + datenum(1970,1,1);
% date de début de dérive (en heure UTC)
dat1 = num2str(debut); H1 = num2str(h1,'%04i');
......@@ -60,32 +56,18 @@ for i = 1:length(num_spot);
filename = ['s' num2str(num_spot(i),'%03i') '_' mission yr ];
I = find(spot(:) == num_spot(i));
tmp = [nM(I) lon(I) lat(I) nU(I)];
K = find(tmp(:,1) > dat1 & tmp(:,1) < dat2);
data = struct('timeM',{0},'lon',{''},'lat',{''},'timeU',{0},...
'start',{''},'end',{''},'mission',{''},...
'name',{''},'spot',{''});
data.timeM = tmp(K,1);
data.lon = tmp(K,2);
data.lat = tmp(K,3);
data.timeU = tmp(K,4);
data.start = [datestr(data.timeM(1)) ' UTC'];
if(~exist('fin','var'))
data.end = [''];
else
data.end = [datestr(data.timeM(end)) ' UTC'];
end
data.spot = [ num_spot(i) ];
data.mission = [ mission ];
data.name = [ filename ];
if(~exist([ rootdir '/drifter/' mission '/csv' ]));
K = find(spot_clean{1,1} == num_spot(i) & nM> dat1 & nM < dat2);
data = struct('spot',{''},'lon',{''},'lat',{''},'timeU',{0},...
'timeM',{0},'flag',{''});
data.timeM = nM(K);
data.lon = spot_clean{1,2}(K);
data.lat = spot_clean{1,3}(K);
data.timeU = spot_clean{1,4}(K);
data.flag = spot_clean{1,6}(K);
data.spot = spot_clean{1,1}(K);
if(~exist([ rootdir '/drifter/' mission ]));
mkdir([ rootdir '/drifter/' mission '/csv' ]);
mkdir([ rootdir '/drifter/' mission '/mat' ]);
end
......@@ -94,9 +76,7 @@ for i = 1:length(num_spot);
save(pathname,'data');
csvname = [ rootdir '/drifter/' mission '/csv/' filename '.csv'];
CSV = struct2cell(data);
CSV = [ CSV(1,1) CSV(2,1) CSV(3,1) CSV(4,1)];
CSV = cell2mat(CSV);
dlmwrite(csvname,CSV,'precision', 11);
CSV = struct2table(data);
writetable(CSV,csvname);
end
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment