clear all;close all;
var_20C_L32=ncread('/Users/heyc/Micom/output/NOI20C_sr10m60d28/dz/NOI20C_sr10m60d28_micom_hm_dz_TS_LY32_YR_0557-0695_MM_03_wgtFLAG_Thickness.nc','var_ts');
var_20C_L33=ncread('/Users/heyc/Micom/output/NOI20C_sr10m60d28/dz/NOI20C_sr10m60d28_micom_hm_dz_TS_LY33_YR_0557-0695_MM_03_wgtFLAG_Thickness.nc','var_ts');
var_20C=var_20C_L32+var_20C_L33;

var_CORE_L32=ncread('/Users/heyc/output/NOIIA_T62_tn11_sr10m60d_01/dz/NOIIA_T62_tn11_sr10m60d_01_micom_hm_dz_TS_LY32_YR_0241-0300_MM_03_wgtFLAG_Thickness.nc','var_ts');
var_CORE_L33=ncread('/Users/heyc/output/NOIIA_T62_tn11_sr10m60d_01/dz/NOIIA_T62_tn11_sr10m60d_01_micom_hm_dz_TS_LY33_YR_0241-0300_MM_03_wgtFLAG_Thickness.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');


