clear all;close all;
var_20C_L32=ncread('LSW_Thickness_L32_20CR.nc','var_ts');
var_20C_L33=ncread('LSW_Thickness_L33_20CR.nc','var_ts');
var_20C=var_20C_L32+var_20C_L33;

var_CORE_L32=ncread('LSW_Thickness_L32_CORE.nc','var_ts');
var_CORE_L33=ncread('LSW_Thickness_L33_CORE.nc','var_ts');
var_CORE=var_CORE_L32+var_CORE_L33;

figure;hold on;grid on;box on;
plot(1:139,var_20C,'r-');
plot(1:139,var_20C_L32,'g-');
plot(1:139,var_20C_L33,'b-');
set(gca,'xtick',0:10:139,'xticklabel',1870:10:2009,'ylim',[200 1600]);
xlabel('year');ylabel('tickness (m)');
legend('Layer32+33','Layer32','Layer33','location','northwest');
print('-painter','-depsc','-r300','LSW_Thickness_20CR_1871');

figure;hold on;grid on;box on;
plot(78:139,var_20C(78:139),'r-');
plot(78:139,var_20C_L32(78:139),'g-');
plot(78:139,var_20C_L33(78:139),'b-');
set(gca,'xtick',80:10:139,'xticklabel',1950:10:2009,'ylim',[200 1600]);
xlabel('year');ylabel('tickness (m)');
legend('Layer32+33','Layer32','Layer33','location','northwest');
print('-painter','-depsc','-r300','LSW_Thickness_20CR_1948');

figure;hold on;grid on;box on;
plot(78:137,var_CORE,'r-');
plot(78:137,var_CORE_L32,'g-');
plot(78:137,var_CORE_L33,'b-');
set(gca,'xtick',80:10:137,'xticklabel',1950:10:2009,'ylim',[200 1600]);
xlabel('year');ylabel('tickness (m)');
legend('Layer32+33','Layer32','Layer33','location','northwest');
print('-painter','-depsc','-r300','LSW_Thickness_CORE_1948');


