clc ;clear all;close all

pathi='/data3/yiwen/Timcom_exp/matlab_code/timeseries/'
pathi2='/data3/yiwen/Timcom_exp/matlab_code/timeseries_renew/'
load([pathi,'test01_6th']);SV_control=AMOC_JRAc_yr(41:end);
load([pathi2,'change_wind_45']);SV_45=SV_yr;
load([pathi2,'change_FWF_45']);SV_FWF45=SV_yr;
load([pathi2,'change_wind_medius']);SV_medius=SV_yr;
% load([pathi,'change_windNS']);SV_NS=SV_yr;
% load([pathi,'change_FWFNS']);SV_FWFNS=SV_yr;


pathi3='/data3/yiwen/cesm_final_exp/plot_test/old_data_test/timeseries_result/'
load([pathi3,'MOC_OBS_time'])

%%
set(0,'defaultfigurecolor','w')
plot([1998:1:2018],SV_control,'linewidth',1.5);hold on
plot([1998:1:2018],SV_medius,'linewidth',1.5);hold on
plot([1998:1:2018],SV_45,'linewidth',1.5);hold on
% plot([1998:1:2018],SV_NS,'linewidth',1.5);hold on
plot([1998:1:2018],SV_FWF45,'linewidth',1.5);hold on
% plot([1998:1:2018],SV_FWFNS,'linewidth',1.5);hold on
plot([2004:1:2018],OBS_MOC_year,'linewidth',1.5,'color','k')
legend({'control','WIND0.1:medius','WIND0.1:45','HOST0.1:45','RAPID'},'location','SouthWest','Fontsize',12)
xlabel('year','fontsize',16)
ylabel('SV','fontsize',16)
title('TIMCOM:GJRA55v2.10 case 26.5 maximum transport','fontsize',18)
set(gca,'Linewidth',2)
axis([1998 2018 8 20])

