function results_bcs

% based on CRUISE_Q4_3_BCS.m

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%
% Load and plot climatological bcs
%
BCF=['M2';'Z0']
clpt=['b*- ';'bx- ']
%
tallfigure;

file=['../CLIM_BCS/unpert/BCS_M2.s2c']
[BCS,freq,gridname]=read_s2c(file);
plot(BCS(:,3),'b*-');

file=['../fcast_archives/EN302_FC1/CRUISE_T6.ALL/BCS_PRT3.s2c']
[BCS,freq,gridname]=read_s2c(file);
freq
plot(BCS(:,3),'b-');

file=['../CLIM_BCS/pert_FC1/BCS_M2_FC1.s2c']
[BCS,freq,gridname]=read_s2c(file);
freq
plot(BCS(:,3),'r-');

legend('M2 CLIM','Z0 CLIM','M2 Pert','Z0 Pert','M2 tot','Z0 tot')
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 EN302_FC1'])
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
