clc;clear all;close all

addpath('/home/yiwen/matlab_test/tight_subplot')
addpath('/data3/yiwen/cesm_final_exp/matlab_toolbox/colorbrewer')
addpath('/data3/yiwen/cesm_final_exp/matlab_toolbox/m_map')

%=============================================
%EXP(1)region (SH) POS_[51S,53S,55S]  NEG[47S,43S]
%=============================================
m_proj('miller','lat',[25,26],'lon',[121,122]);
m_gshhs_f('color','k');              % Coastline...
m_gshhs_f('speckle','color','k');    % with speckle added

grid on;box on;
m_proj('miller','lon',[0 360],'lat',[-90 0]); 
m_coast('patch',[.9 .9 .9]);
m_grid('box','fancy','tickdir','out');
h1=m_line([0 360],[-51 -51],'linewi',1.5,'color','k');     
h2=m_line([0 360],[-53 -53],'linewi',1.5,'color','b','linestyle', '--');     
h3=m_line([0 360],[-55 -55],'linewi',1.5,'color',[0.4940 0.1840 0.5560],'linestyle', '--'); 
h4=m_line([0 360],[-47 -47],'linewi',1.5,'color',[0.9290 0.6940 0.1250],'linestyle', '--');     
h5=m_line([0 360],[-43 -43],'linewi',1.5,'color',[0.8500 0.3250 0.0980],'linestyle', '--');     
legend([h1(1),h2(1),h3(1),h4(1),h5(1)],...
      'control','POS 1X','POS 3X','NEG 1X','NEG 3X',...
      'location','northeast')
%=============================================
%EXP(2)region (NH) Wind[50-70N] HOS[50-70N]
%=============================================
grid on;box on;
m_proj('miller','lon',[121 122],'lat',[25 26]); 
m_coast('patch',[.9 .9 .9]);
m_grid('box','fancy','tickdir','out');
h1=m_line([0 360],[30 30],'linewi',1.5,'color',[0.6350 0.0780 0.1840]);   
h2=m_line([0 360],[85 85],'linewi',1.5,'color',[0.6350 0.0780 0.1840]); 
h3=m_line([0 0],[30 85],'linewi',1.5,'color',[0.6350 0.0780 0.1840]); 
h4=m_line([360 360],[30 85],'linewi',1.5,'color',[0.6350 0.0780 0.1840]); 








