% sens_2.m
% 
% sensetivity of the hindcast skill to variations in nondimensional
% parameters ALPHA and BETA.
% EOF amplitude at depth = 0.3
%
% (i) BETA=0.25 
  a=[    0.75    1.5    2.0     3.0];

   a1=[  27.6    29.2   30.0    30.8      % NRMS (Predicted & OA)
         6.5     6.5    6.6     6.6       % RMS  (Predicted & ERS)
         0.89    0.89   0.88    0.88      % corr    -.-
         0.71    0.71   0.70    0.71      % b       -.-
         6.1     6.4    6.6     6.9       % RMS  (Predicted & CTD)
         0.73    0.70   0.69    0.67      % corr    -.-
         0.70    0.68   0.68    0.67  ];  % b       -.-

% (ii) BETA=0.5
%
% alpha  0.75    1.5    2.0     3.0  
a2=[     33.3    27.6   28.2    29.3      % NRMS (Predicted & OA)
         7.1     6.5    6.5     6.5       % RMS  (Predicted & ERS)
         0.87    0.89   0.89    0.89      % corr    -.-
         0.69    0.70   0.70    0.70      % b       -.-
         7.0     6.1    6.1     6.4       % RMS  (Predicted & CTD)
         0.66    0.73   0.72    0.7       % corr    -.-
         0.65    0.69   0.68    0.68  ];  % b       -.-
% (iii) BETA=1.0
%
% alpha  0.75    1.5    2.0     3.0  
a3=[     41.6    33.4   29.2    27.7      % NRMS (Predicted & OA)
         7.7     7.1    6.6     6.5       % RMS  (Predicted & ERS)
         0.83    0.87   0.89    0.89      % corr    -.-
         0.67    0.69   0.70    0.70      % b       -.-
         9.2     6.9    6.4     6.1       % RMS  (Predicted & CTD)
         0.5     0.66   0.70    0.73      % corr    -.-
         0.58    0.65   0.68    0.69  ];  % b       -.-

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,:),'^k-','Markersize',5,'Linewidth',0.3);
hold on
plot(a(:),a2(1,:),'sk-','Markersize',6,'Linewidth',0.3);
plot(a(:),a3(1,:),'ok-','Markersize',5,'Linewidth',0.3);
plot(a(2),a2(1,2),'sk','Markersize',6,'Markerfacecolor','k')
hold off
axis([0.25 3.5 23  45]);
ylabel('NMRS Diff (%)')
legend('beta=0.25','beta=0.5','beta=1.0','central run')

h3=subplot(423)
set(h3,'Position',[0.1700    0.5439    0.3270    0.1642])
set(gca,'fontsize',7)
plot(a(:),a1(2,:),'^k-','Markersize',5,'Linewidth',0.3);
hold on
plot(a(:),a2(2,:),'sk-','Markersize',6,'Linewidth',0.3);
plot(a(:),a3(2,:),'ok-','Markersize',5,'Linewidth',0.3);
plot(a(2),a2(2,2),'sk','Markersize',6,'Markerfacecolor','k')
hold off
axis([0.25 3.5 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,:),'^k-','Markersize',5,'Linewidth',0.3);
hold on
plot(a(:),a2(5,:),'sk-','Markersize',6,'Linewidth',0.3);
plot(a(:),a3(5,:),'ok-','Markersize',5,'Linewidth',0.3);
plot(a(2),a2(5,2),'sk','Markersize',6,'Markerfacecolor','k')
hold off
axis([0.25 3.5 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,:),'^k-','Markersize',5,'Linewidth',0.3);
hold on
plot(a(:),a2(3,:),'sk-','Markersize',6,'Linewidth',0.3);
plot(a(:),a3(3,:),'ok-','Markersize',5,'Linewidth',0.3);
plot(a(2),a2(3,2),'sk','Markersize',6,'Markerfacecolor','k')
hold off
axis([0.25 3.5 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,:),'^k-','Markersize',5,'Linewidth',0.3);
hold on
plot(a(:),a2(6,:),'sk-','Markersize',6,'Linewidth',0.3);
plot(a(:),a3(6,:),'ok-','Markersize',5,'Linewidth',0.3);
plot(a(2),a2(6,2),'sk','Markersize',6,'Markerfacecolor','k')
hold off
axis([0.25 3.5 0.4  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,:),'^k-','Markersize',5,'Linewidth',0.3);
hold on
plot(a(:),a2(4,:),'sk-','Markersize',6,'Linewidth',0.3);
plot(a(:),a3(4,:),'ok-','Markersize',5,'Linewidth',0.3);
plot(a(2),a2(4,2),'sk','Markersize',6,'Markerfacecolor','k')
hold off
axis([0.25 3.5 0.6  0.8]);
ylabel('B (inclnation)')
%title('Predicted vs ERS')
%legend('beta=0.25','beta=0.5','beta=1.0')
xlabel(['Alpha'],'fontsize',9)

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

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