%%% Script needs:
% oetsettings
% Vlugtenburg scripts folder added to path!


%%% Aanpassing Matthieu 7/9/2009
% Locaal path D:\Documents\Checkouts weggehaald, Script nu draaien in de vlugtenburg root directory



% Survey 1
walks = Read_XYZ('walking\raw\2009_06_10\Products\XYZ_RDNAP.txt');
jetski = Read_XYZ('jetski\raw\2009_06_09\Products\XYZ_jetski_RD_NAP.txt');

% Survey 2
% walks = Read_XYZ('vlugtenburg\walking\raw\2009_07_17\Products\XYZ_RDNAP.txt');
% jetski = Read_XYZ('vlugtenburg\jetski\raw\2009_07_16\Products\XYZ_jetski_RD_NAP.txt');

% Survey 3
% walks = Read_XYZ('walking\raw\2009_08_05\Products\XYZ_RDNAP.txt');
% jetski = Read_XYZ('jetski\raw\2009_08_06\Products\XYZ_jetski_RD_NAP.txt');


figure
scatter(jetski.data(:,1),jetski.data(:,2),3,jetski.data(:,3))
hold on
scatter(walks.data(:,1),walks.data(:,2),3,walks.data(:,3))
add_surveygrid
axis image
xlim([6.69e4 6.946e4])
ylim([4.459e5 4.484e5])
xlabel('RD coordinates [10^3 m]');
ylabel('RD coordinates [10^3 m]');
colorbar
title('Survey Grid');

%% add jarkus transects
if 0  %% toggle between jarkus lines off (value =0) and on (value =1)
jarkus_file = 'http://opendap.deltares.nl:8080/thredds/dodsC/opendap/rijkswaterstaat/jarkus/profiles/transect.nc';
profile = 1425;
profile_count = 8;
time = 43; % for 2008

%RD coordinates of profile
X_rd = nc_varget(jarkus_file,'x',[profile 0],[profile_count -1]);
Y_rd = nc_varget(jarkus_file,'y',[profile 0],[profile_count -1]);

% figure
plot(X_rd',Y_rd','m','linewidth',2);
end
