function results_bcs(runid)

% based on CRUISE_Q4_3_BCS.m

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%
% Load and plot climatological bcs
%
BCF=['M2';'N2';'S2';'K1';'O1';'Z0']
clpt=['b*- ';'bx- ';'b+- ';'ro--';'rs--';'md-.']
%
tallfigure;
for i=1:6
   file=['../CLIM_BCS/BCS_',BCF(i,:),'.s2c']
   [BCS,freq,gridname]=read_s2c(file);
   plot(BCS(:,2),clpt(i,:));
end
clpt=['k- ';'k--';'k-.';'k: ';]
for i=1:4
   file=['../fcast_archives/',runid,'/CRUISE_T6.ALL/BCS_PRT',num2str(i),'.s2c']
   [BCS,freq,gridname]=read_s2c(file);
   freq
   plot(BCS(:,2),clpt(i,:));
end
legend('M2 CLIM','N2 CLIM','S2 CLIM','K1 CLIM','O1 CLIM',...
'Z0 CLIM',verbatim('M2 T6_1'),verbatim('Z0 T6_1'),...
verbatim('M2 T6_2'),verbatim('Z0 T6_2'))
ylabel('Amplitude of BCS (m) for Quoddy Simulation 3')
xlabel('Boundary node number (CCW around boundary, beginning at SW corner)')
title(['Boundary amplitudes for experiment ',runid])
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
