% clear all
% 



% fileDuck= fullfile('/data/olvac/travail_en_cours/graffin/Coastline','Duck_shorelinePerDayFinal.mat');
% fileNarrabeen= fullfile('/data/olvac/travail_en_cours/graffin/Coastline','NarrabeenDataForRob.mat');
% fileTairua= fullfile('/data/olvac/travail_en_cours/graffin/Coastline','TairuaShorelineData.mat');

% 
% dataDuck=load(fileDuck);
% dataNarrabeen=load(fileNarrabeen);
% dataTairua=load(fileTairua);
% 
% timeDuck=dataDuck.toHat;
% timeTairua=dataTairua.Shorecast.time;
% timeNarrabeen=dataNarrabeen.shoreData.dates;
% 
% 
% shorelineTairua=dataTairua.Shorecast.average; %already averaged in space
% shorelineDuck=dataDuck.xoHat;
% shorelineNarrabeen=dataNarrabeen.shoreData.x;
% 
% for i=1:8235
%     shorelineDuckAverage(i)=mean(shorelineDuck(i,:),'omitnan');
% end
% 
% %shorelineDuckAverageDetrend=detrend(shorelineDuckAverage);
% 
% shorelines=[];
% time=725008:738156;
% 
% 
% % i1=725007;
% for i=725008:738156
%     shorelines(i-725007,1)=NaN;
% end
% 
% for i=1:length(timeTairua)
%     temp=datestr(timeTairua(i));
%     timeTairua(i)=datenum([temp(1:11) ' 00:00:00']);
% end
% 
% shorelines(timeDuck-725007,1)=shorelineDuckAverage;
% shorelines(timeNarrabeen-725007,2)=shorelineNarrabeen;
% shorelines(timeTairua-725007,3)=shorelineTairua;
% 
% % 
% shorelines(:,1)=interpolNaN(shorelines(:,1),timeDuck); %
% shorelines(:,2)=interpolNaN(shorelines(:,2),timeNarrabeen); %lon 151.30°E ; lat -33.72°N
% shorelines(:,3)=interpolNaN(shorelines(:,3),timeTairua);

% for i=1:7750
%     if sum(X(i,:),'omitnan')~=0
%     Tnan=find(isnan(X(i,:)));
%     Tval=find(~isnan(X(i,:)));
%     Tnancool=Tnan(find(Tnan>Tval(1) & Tnan<Tval(length(Tval))));
%     X(i,Tnancool)=interp1(Tval,X(i,Tval),Tnancool);
%     end
% en

% for i=1:14140
%     D=sqrt((lon(i)-LON_loc).^2+(lat(i)-LAT_loc).^2);
%     pos=find(D==min(D));
%     shorelineGlobal_2016_2018(i,:)=X(pos,:);
% end
    

% a='./in_situ_shoreline.nc';
% ncdisp(a)
% 
TrucVert=ncread(f3,'shorelineTrucVert');
time=ncread(f3,'time');

[TrucVert_mounthly mounths]=fromDaily2Mounthly(TrucVert,time);
TrucVert_mounthly(225:420)=detrend(TrucVert_mounthly(225:420));
% dt_TrucVert_mounthly=TrucVert_mounthly;
% dt_TrucVert_test=dt_TrucVert_mounthly;
% 
% for i=1:12
%     dt_TrucVert_test(i+12:12:432)=dt_TrucVert_mounthly(i+12:12:432)-mean(dt_TrucVert_mounthly(i+12:12:432),'omitnan');
% end
% 
% shoreline_norm_TrucVert=dt_TrucVert_test./std(dt_TrucVert_test,'omitnan');