% sens_3.m
%
% sensitivity of the hindcast skill to variations in nondimensional
% parameter GAMMA^2. The rest of the parameters corresponds to central 
% run: DEA=0.3, BETA=0.5, ALPHA=1.5
%   
% Gamma^2        0.15    0.3    0.6
%         ----------------------------
          a1=[  33.3    27.6   32.6        % NRMS (Predicted & OA)
                7.2     6.5    6.9         % RMS  (Predicted & ERS)
                0.87    0.89   0.89        % corr    -.-
                0.83    0.70   0.61        % b       -.-
                7.9     6.1    5.9         % RMS  (Predicted & CTD)
                0.72    0.73   0.71        % corr    -.-
                0.98    0.69   0.53 ];     % b       -.-
   a=[.15  .305  0.6];
% central run: unfiltered alongboundary vorticity
un=[  33.4 6.9 0.87 0.71 6.4 0.73 0.78];
set(0,'DefaultAxesFontsize',6) 
hh=subplot(421)
set(hh,'Position',[0.1700    0.7808    0.715    0.1642])
set(gca,'fontsize',7)
plot(a(:),a1(1,:),'sk-','Markersize',6,'Linewidth',0.3);
hold on
plot(a(2),a1(1,2),'sk-','Markersize',6,'Markerfacecolor','k')
plot(a(2),un(1),'^k','Markersize',5);
hold off
axis([0.1 0.65 23  45]);
ylabel('NMRS Diff (%)')
legend('Gamma^2 sens','central run','unfilt vorticity')

h3=subplot(423)
set(h3,'Position',[0.1700    0.5439    0.3270    0.1642])
set(gca,'fontsize',7)
plot(a(:),a1(2,:),'sk-','Markersize',6,'Linewidth',0.3);
hold on
plot(a(2),a1(2,2),'sk-','Markersize',6,'Markerfacecolor','k')
plot(a(2),un(2),'^k','Markersize',5);
hold off
axis([0.1 .65 6.0  8.0]);
ylabel('RMS (cm) ')
title(['Predicted vs ERS'],'fontsize',9)
%legend('beta=0.25','beta=0.5','beta=1.0')

h4=subplot(424)
set(h4,'Position',[0.5680    0.5439    0.3270    0.1642])
set(gca,'fontsize',7)
plot(a(:),a1(5,:),'sk-','Markersize',6,'Linewidth',0.3);
hold on
plot(a(2),a1(5,2),'sk-','Markersize',6,'Markerfacecolor','k')
plot(a(2),un(5),'^k','Markersize',5);
hold off
axis([0.1 .65 5.5  9.5]);
% ylabel('RMS (cm) ')
title(['Predicted vs CTD'],'fontsize',9)
% legend('beta=0.25','beta=0.5','beta=1.0')

h5=subplot(425)
set(h5,'Position',[0.1700    0.3269    0.3270    0.1642])
set(gca,'fontsize',7)
plot(a(:),a1(3,:),'sk-','Markersize',6,'Linewidth',0.3);
hold on
plot(a(2),a1(3,2),'sk-','Markersize',6,'Markerfacecolor','k')
plot(a(2),un(3),'^k','Markersize',5);
hold off
axis([0.1 .65 0.85  0.95]);
ylabel('CORR COEF')
%title('Predicted vs ERS')
%legend('beta=0.25','beta=0.5','beta=1.0')

h6=subplot(426)
set(h6,'Position',[0.5680    0.3269    0.3270    0.1642])
set(gca,'fontsize',7)
plot(a(:),a1(6,:),'sk-','Markersize',6,'Linewidth',0.3);
hold on
plot(a(2),a1(6,2),'sk-','Markersize',6,'Markerfacecolor','k')
plot(a(2),un(6),'^k','Markersize',5);
hold off
axis([0.1 .65 0.6  0.8]);
% ylabel('RMS (cm) ')
%title('Predicted vs CTD')
%legend('beta=0.25','beta=0.5','beta=1.0')

h7=subplot(427)
set(h7,'Position',[0.1700    0.11    0.3270    0.1642])
set(gca,'fontsize',7)
plot(a(:),a1(4,:),'sk-','Markersize',6,'Linewidth',0.3);
hold on
plot(a(2),a1(4,2),'sk-','Markersize',6,'Markerfacecolor','k')
plot(a(2),un(4),'^k','Markersize',5);
hold off
axis([0.1 .65 0.5  .9]);
ylabel('B (inclnation)')
%title('Predicted vs ERS')
%legend('beta=0.25','beta=0.5','beta=1.0')
xlabel(['Gamma^2'],'fontsize',9)

 h8=subplot(428)
set(h8,'Position',[0.5680    0.11   0.3270    0.1642])
set(gca,'fontsize',7)
plot(a(:),a1(7,:),'sk-','Markersize',6,'Linewidth',0.3);
hold on
plot(a(2),a1(7,2),'sk-','Markersize',6,'Markerfacecolor','k')
plot(a(2),un(7),'^k','Markersize',5);
hold off
axis([0.1 .65 0.5  1.05]);
% ylabel('RMS (cm) ')
%title('Predicted vs CTD')
%legend('beta=0.25','beta=0.5','beta=1.0')
xlabel(['Gamma^2'],'fontsize',9)

set(gcf,'PaperPosition',[1.0 2.0 6.5 6.5])
