figure(1)
[a1,c1,n1] = variogramfit(dist_uv_F,gamma_uv_F,a01_F,c01_F,[],...
'solver','fminsearchbnd',...
'nugget',0,...
'plotit',true)
hold on
[a1,c1,n1] = variogramfit(dist_uvV,gamma_uvV,a01_V,c01_V,[],...
'solver','fminsearchbnd',...
'nugget',0,...
'plotit',true)
hold on
ylabel('gamma (uv)')
xlabel('lag distance dataset1-2')
hold on
ax1 = gca;
set(ax1, 'XColor','k', 'YColor','k')
ax2 = axes('Position', get(ax1,'Position'),...
'XAxisLocation', 'top',...
'YAxisLocation', 'right',...
'Color','none',...
'XColor', 'k', 'YColor', 'k');
hold on
h = variogramfit(dist_uv,gamma_uv,a01,c01,[],...
'solver','fminsearchbnd',...
'nugget',0,...
'plotit',true)
hold on
xlabel('lag distance dataset3')