%clear;close all;
LATEXtools_param

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fname_SaralAltika='tracerefaltiKa_t429.txt'



fname_jason2_t9='J2_T009_C270_prev_20151026.asc'
fname_jason2_t222='J2_T222_C270_prev_20151026.asc'

fname_waypoints='tethys_waypoints.txt'; %files containing the positions

fname_RADAR='Mask_radar_DYFAMED.mat'

fname_TSG='20150926000042.COR'

fname_drifters='DRIFTERS/Oscahr_Drifter_Data_20151107T081357'%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%- Map boundaries -%
lonmin=6.7;
lonmax=9.1;
latmin=42.9;
latmax=44.1;
%
%- Geographic projection for m_map -%plot_
projection = 'mercator';
%
%- User coast in dirin -%
UserCoast='GOL.mat'; % in "sub_functions" directory
%
%- Plot the figures
m_proj(projection,'lon',[lonmin lonmax],'lat',[latmin latmax])

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%        SST L3       %%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

h=figure;%('Visible','off'); 
hold on; 

%plot_BATHY_ETOPO1
plot_BATHY(dirin_BATHY,[200 500 1000 2000]);

%plot_FSLE(fname_FSLE,0.2);

plot_SaralAltika(fname_SaralAltika,lonmin,lonmax,latmin,latmax,'r')
%plot_SaralAltika_modified('AL_T674_C28_prev_20151026_modifed.asc')
%plot_JASON2(fname_jason2_t9)
%plot_JASON2(fname_jason2_t222)

plot_MVPcasts('lonlatMVP.txt')

%plot_waypoints([dirin_Log_Position,fname_waypoints]);%,[7:11]);

%%%% Look for cast along satellite traces
%load('./MVP/L2/MVP_all_casts.mat')
%%jason2 #9	30/10 15:30 	30/10 23:30
%cast_jason=[24:89];
%m_plot(lon_mat(cast_jason),lat_mat(cast_jason),'bo');






%- User's High Resolution Coastline -%
m_usercoast(UserCoast,'patch',[.7 .7 .7],'edgecolor','k');
%m_coast('patch',[.7 .7 .7],'edgecolor','k');
%
%- Geographical Axis -%
m_grid('box','fancy','tickdir','in');
xlabel('Longtitude');
ylabel('Latitude');
%
%
print(h,'tmp','-dpng','-r200')
%print(h,[figdir,'route_SST_',datestr(now,'yyyymmddTHHMM')],'-dpng','-r200')





