addpath /data/mgoes/AX8/Time_transform
addpath /data/mgoes/matlab/timeplt

ncload bc_transp_adt.nc

t0 = julian(1950,1,1,0);
time = gregorian(TIME/24+t0);

%load ../AX97/Data/trp_all_manual_abs.mat
%load ../AX97/Data/trp_all_manual_abs_ypos45_30.mat
load ../AX97/Data/trp_all_manual_abs_ypos45_36.mat
tempo = nanmean( ymd_dec_ax97');
bc = [tsp_bc.transp_400_minus]/1e6;
%bc = [tsp_bc.transp_800_minus]/1e6;
time2 = month2yfraction(gregorian2month(time));

bc(bc==0)=nan;
figure(2),clf
plot(time2,BC)
hold on
plot(tempo,bc,'-or','markerfacecolor','k')
xlim([2004 2015])

ylabel('Sv')
!save2eps('ADT_XBT.eps',2,300)

