% sens_1.m
% to draw a set of plots which display sensitivity of the hindcast
% skill to the EOF amplitude at depth (dh*).
%    dh1   dh2   dh3   dh4   dh5   dh7    
a=[  0.1   0.2   0.3   0.4   0.5   0.7 ];   % amplitude in a deep water
a1=[ 32.3  27.9  27.6  27.8  27.9  28.1     % NRMS (Predicted & OA)
     7.2   6.5   6.5   6.5   6.5   6.5      % RMS  (Predicted & ERS)
     0.86  0.89  0.89  0.89  0.89  0.89     % corr    -.-
     0.69  0.71  0.70  0.70  0.70  0.70     % b       -.-
     6.4   6.1   6.1   6.1   6.2   6.2      % RMS  (Predicted & CTD)
     0.71  0.73  0.73  0.72  0.72  0.71     % corr    -.-
     0.72  0.71  0.69  0.68  0.68  0.67  ]; % b       -.-
% 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',7,'Linewidth',0.3);
hold on
plot(a(3),a1(1,3),'sk','Markersize',7,'Markerfacecolor','k');
plot(a(3),un(1),'^k','Markersize',5);
hold off
axis([0.05 0.75 22  37]);
ylabel('NMRS Diff (%)')
legend('EOF sensitivity ','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',7,'Linewidth',0.3);
hold on
plot(a(3),a1(2,3),'sk','Markersize',7,'Markerfacecolor','k');
plot(a(3),un(2),'^k','Markersize',5);
hold off
axis([0.05 0.75 6.0  7.5]);
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',7,'Linewidth',0.3);
hold on
plot(a(3),a1(5,3),'sk','Markersize',7,'Markerfacecolor','k');
plot(a(3),un(5),'^k','Markersize',5);
hold off
axis([0.05 0.75 5.7  6.7]);
% 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',7,'Linewidth',0.3);
hold on
plot(a(3),a1(3,3),'sk','Markersize',7,'Markerfacecolor','k');
plot(a(3),un(3),'^k','Markersize',5);
hold off
axis([0.05 0.75 0.8  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',7,'Linewidth',0.3);
hold on
plot(a(3),a1(6,3),'sk','Markersize',7,'Markerfacecolor','k');
plot(a(3),un(6),'^k','Markersize',5);
hold off
axis([0.05 0.75 0.65  0.75]);
% 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',7,'Linewidth',0.3);
hold on
plot(a(3),a1(4,3),'sk','Markersize',7,'Markerfacecolor','k');
plot(a(3),un(4),'^k','Markersize',5);
hold off
axis([0.05 .75 0.6  0.8]);
ylabel('B (inclnation)')
%title('Predicted vs ERS')
%legend('beta=0.25','beta=0.5','beta=1.0')
xlabel(['Deep EOF Amplitude'],'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',7,'Linewidth',0.3);
hold on
plot(a(3),a1(7,3),'sk','Markersize',7,'Markerfacecolor','k');
plot(a(3),un(7),'^k','Markersize',5);
hold off
axis([0.05 0.75 0.6  0.8]);
% ylabel('RMS (cm) ')
%title('Predicted vs CTD')
%legend('beta=0.25','beta=0.5','beta=1.0')
xlabel(['Deep EOF Amplitude'],'fontsize',9)

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