close all
clear all

fname='../CROCO_FILES_3/croco_avg_Y2022M08.nc';

nc=netcdf(fname);

u=squeeze(mean(nc{'ubar'}(:)));
v=squeeze(mean(nc{'vbar'}(:)));
lon=nc{'lon_rho'}(:);
lat=nc{'lat_rho'}(:);
angle=nc{'angle'}(:);
h=nc{'h'}(:);
zeta=squeeze(mean(nc{'zeta'}(:)));
close(nc)
[M L]=size(lon);

cosa = cos(angle);
sina = sin(angle);
ur=u2rho_2d(u);
vr=v2rho_2d(v);
u = ur.*cosa - vr.*sina;
v = vr.*cosa + ur.*sina;

minlon=108.0; maxlon=110.0;
minlat=  12; maxlat= 14;
res=0.03;
[LON LAT]=meshgrid(minlon:res:maxlon,minlat:res:maxlat);
U=griddata(lon,lat,u,LON,LAT);
V=griddata(lon,lat,v,LON,LAT);

minplot=0.5; %0.4;
speed=sqrt(u.^2+v.^2);
D=h+zeta;
speed(D<minplot)=NaN;
ur(D<minplot)=NaN;
ur(D<minplot)=NaN;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

hf = figure('pos',[500 500 500 500]);
axis tight; set(hf,'DoubleBuffer','on');
%m_proj('stereo','lat',13.5,'lon',109.5,'rad',0.7);
%m_proj('albers equal-area','lat',[13.5 13.8],'long',[108.5 109.5],'rect','on');
m_proj('equidistant','lat',[13  13.8],'lon',[109.15 109.5]);

cmin=0; cmax=0.1; nbcol=10;
cint=(cmax-cmin)/nbcol;
map=AdvancedColormap('wyor',nbcol);
colormap(map);

m_pcolor(lon,lat,speed); hold on
%m_contourf(lon,lat,speed,[cmin:cint:cmax]);
shading flat; colorbar; 
hold on
%m_contour(lon,lat,h,[-0.4:100:100.2],'color',[0.7 0.7 0.7],'linewidth',3); 
Vmax=max(max(sqrt(u.^2+v.^2)));
m_quiver(LON,LAT,U,V,5.*Vmax)
%m_streamslice(lon,lat,u,v,2)
caxis([cmin cmax])
title(['CROCO residual currents'],'fontsize',11)
m_gshhs_f('patch','k')
m_grid('box','fancy','tickdir','in','fontsize',12);
hold off
export_fig -transparent currents.pdf

return

figure
pcolor(lon,lat,speed); hold on
colormap('cool')
shading flat
colorbar
quiver(lon,lat,u,v,1)
axis([109.15 109.5 13 13.8])
caxis([0 1])