%plot_GOL
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%% APetrenko 15 mars 2013
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Besoin : mettre aussi dans directoire de travail GoL.mat
%% et les fichiers de bathymétrie line0.txt et line200.txt
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%  This  file  is  free software ;  it  is distributed  in  %
%%  the hope that  it  will  be  useful,  but  without  any  %
%%  warranty.  You  can  redistribute  it and/or modify  it  %
%%  under  the  terms  of  the  GNU  General Public License  %
%%  as  published  by  the   Free  Software  Foundation  at  %
%%       http://www.gnu.org/copyleft/gpl.html                %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Si besoin mais a adapter par rapport a votre chemin de repertoire de travail
% Sinon a mettre dans le setpath
addpath(genpath('C:\Users\annep\Documents\Programs\m_map'))

%- area limits -% (in degres)
lon_min =  3; 
lon_max =  6;
lat_min = 42;
lat_max = 43.75;
%NB: for demo the Gulf of Lion is chosen
%
%- User high resolution coastline
UserCoast='GOL.mat';

%- Plot the figure
figure(1); hold on; 
%set(gcf,'position',[ScreenWidth/2 1 ScreenWidth/2  ScreenHeigth]);
m_proj('mercator','lon',[lon_min lon_max],'lat',[lat_min lat_max])

%-current vector - if data available
%m_quiver(lon,lat,u_rho,v_rho,10)

%- User's High Resolution Coastline -%
m_usercoast(UserCoast,'patch',[.7 .7 .7],'edgecolor','k');
%- Geographical Axis -%
m_grid('box','fancy','tickdir','in');
xlabel('Longtitude');
ylabel('Latitude');

% exemple de bathymétrie
load line100.txt
load line200.txt
hold on
m_plot(line100(:,2),line100(:,1),'k')
m_plot(line200(:,2),line200(:,1),'k')  %isoligne -200 m
m_plot(5.255,43.135,'kx')
m_text(5.26,43.135,'J','FontSize',15)
m_text(5.38,43.28,'M','FontSize',15)
m_text(5.48,43.1,'C','FontSize',15)
%m_text(5.4,43.4,'J = JULIO','FontSize',15)
%m_text(5.4,43.35,'M = MARSEILLE','FontSize',15)
%m_text(5.4,43.45,'C = Canyon de Cassidaigne','FontSize',15)


%%% limites du Parc National
%%% last point cote LA Ciotat to be checked done approximately
long_parc=[5+21.08/60 5+13.58/60 5+11.5/60  5+30.2/60 5+38/60];
lat_parc=[43+14.23/60 43+12.88/60 43+10.4/60 43+1.16/60 43+10/60];
m_plot(long_parc,lat_parc,'g--')