%%
clear dir
files=dir('/Volumes/Untitled/CALLIPSO_DULI/ALay_L2_05km_Himalya_*')

l=length(files)




A07=load('ALay_L2_05km_Himalya_2007_03_assembled_lat_lon_alt.dat.txt');
B07=load('ALay_L2_05km_Himalya_2007_04_assembled_lat_lon_alt.dat.txt');
C07=load('ALay_L2_05km_Himalya_2007_05_assembled_lat_lon_alt.dat.txt');
D07=load('ALay_L2_05km_Himalya_2007_06_assembled_lat_lon_alt.dat.txt');
E07=load('ALay_L2_05km_Himalya_2007_07_assembled_lat_lon_alt.dat.txt');
F07=load('ALay_L2_05km_Himalya_2007_08_assembled_lat_lon_alt.dat.txt');


A08=load('ALay_L2_05km_Himalya_2008_03_assembled_lat_lon_alt.dat.txt');
B08=load('ALay_L2_05km_Himalya_2008_04_assembled_lat_lon_alt.dat.txt');
C08=load('ALay_L2_05km_Himalya_2008_05_assembled_lat_lon_alt.dat.txt');
D08=load('ALay_L2_05km_Himalya_2008_06_assembled_lat_lon_alt.dat.txt');
E08=load('ALay_L2_05km_Himalya_2008_07_assembled_lat_lon_alt.dat.txt');
F08=load('ALay_L2_05km_Himalya_2008_08_assembled_lat_lon_alt.dat.txt');


A09=load('ALay_L2_05km_Himalya_2009_03_assembled_lat_lon_alt.dat.txt');
B09=load('ALay_L2_05km_Himalya_2009_04_assembled_lat_lon_alt.dat.txt');
C09=load('ALay_L2_05km_Himalya_2009_05_assembled_lat_lon_alt.dat.txt');
D09=load('ALay_L2_05km_Himalya_2009_06_assembled_lat_lon_alt.dat.txt');
E09=load('ALay_L2_05km_Himalya_2009_07_assembled_lat_lon_alt.dat.txt');
F09=load('ALay_L2_05km_Himalya_2009_08_assembled_lat_lon_alt.dat.txt');


A10=load('ALay_L2_05km_Himalya_2010_03_assembled_lat_lon_alt.dat.txt');
B10=load('ALay_L2_05km_Himalya_2010_04_assembled_lat_lon_alt.dat.txt');
C10=load('ALay_L2_05km_Himalya_2010_05_assembled_lat_lon_alt.dat.txt');
D10=load('ALay_L2_05km_Himalya_2010_06_assembled_lat_lon_alt.dat.txt');
E10=load('ALay_L2_05km_Himalya_2010_07_assembled_lat_lon_alt.dat.txt');
F10=load('ALay_L2_05km_Himalya_2010_08_assembled_lat_lon_alt.dat.txt');


A11=load('ALay_L2_05km_Himalya_2011_03_assembled_lat_lon_alt.dat.txt');
B11=load('ALay_L2_05km_Himalya_2011_04_assembled_lat_lon_alt.dat.txt');
C11=load('ALay_L2_05km_Himalya_2011_05_assembled_lat_lon_alt.dat.txt');
D11=load('ALay_L2_05km_Himalya_2011_06_assembled_lat_lon_alt.dat.txt');
E11=load('ALay_L2_05km_Himalya_2011_07_assembled_lat_lon_alt.dat.txt');
F11=load('ALay_L2_05km_Himalya_2011_08_assembled_lat_lon_alt.dat.txt');

A12=load('ALay_L2_05km_Himalya_2012_03_assembled_lat_lon_alt.dat.txt');
B12=load('ALay_L2_05km_Himalya_2012_04_assembled_lat_lon_alt.dat.txt');
C12=load('ALay_L2_05km_Himalya_2012_05_assembled_lat_lon_alt.dat.txt');
D12=load('ALay_L2_05km_Himalya_2012_06_assembled_lat_lon_alt.dat.txt');
E12=load('ALay_L2_05km_Himalya_2012_07_assembled_lat_lon_alt.dat.txt');
F12=load('ALay_L2_05km_Himalya_2012_08_assembled_lat_lon_alt.dat.txt');


A13=load('ALay_L2_05km_Himalya_2013_03_assembled_lat_lon_alt.dat.txt');
B13=load('ALay_L2_05km_Himalya_2013_04_assembled_lat_lon_alt.dat.txt');
C13=load('ALay_L2_05km_Himalya_2013_05_assembled_lat_lon_alt.dat.txt');
D13=load('ALay_L2_05km_Himalya_2013_06_assembled_lat_lon_alt.dat.txt');
E13=load('ALay_L2_05km_Himalya_2013_07_assembled_lat_lon_alt.dat.txt');
F13=load('ALay_L2_05km_Himalya_2013_08_assembled_lat_lon_alt.dat.txt');


A14=load('ALay_L2_05km_Himalya_2014_03_assembled_lat_lon_alt.dat.txt');
B14=load('ALay_L2_05km_Himalya_2014_04_assembled_lat_lon_alt.dat.txt');
C14=load('ALay_L2_05km_Himalya_2014_05_assembled_lat_lon_alt.dat.txt');
D14=load('ALay_L2_05km_Himalya_2014_06_assembled_lat_lon_alt.dat.txt');
E14=load('ALay_L2_05km_Himalya_2014_07_assembled_lat_lon_alt.dat.txt');
F14=load('ALay_L2_05km_Himalya_2014_08_assembled_lat_lon_alt.dat.txt');


A15=load('ALay_L2_05km_Himalya_2015_03_assembled_lat_lon_alt.dat.txt');
B15=load('ALay_L2_05km_Himalya_2015_04_assembled_lat_lon_alt.dat.txt');
C15=load('ALay_L2_05km_Himalya_2015_05_assembled_lat_lon_alt.dat.txt');
D15=load('ALay_L2_05km_Himalya_2015_06_assembled_lat_lon_alt.dat.txt');
E15=load('ALay_L2_05km_Himalya_2015_07_assembled_lat_lon_alt.dat.txt');
F15=load('ALay_L2_05km_Himalya_2015_08_assembled_lat_lon_alt.dat.txt');


A16=load('ALay_L2_05km_Himalya_2016_03_assembled_lat_lon_alt.dat.txt');
B16=load('ALay_L2_05km_Himalya_2016_04_assembled_lat_lon_alt.dat.txt');
C16=load('ALay_L2_05km_Himalya_2016_05_assembled_lat_lon_alt.dat.txt');
D16=load('ALay_L2_05km_Himalya_2016_06_assembled_lat_lon_alt.dat.txt');
E16=load('ALay_L2_05km_Himalya_2016_07_assembled_lat_lon_alt.dat.txt');
F16=load('ALay_L2_05km_Himalya_2016_08_assembled_lat_lon_alt.dat.txt');







dataspat=[A11; B11; C11; D11; E11; F11; A12; B12; C12; D12; E12; F12; A13; B13; C13; D13; E13; F13; A14; B14; C14; D14; E14; F14; A15; B15; C15; D15; E15; F15; A16; B16; C16; D16; E16; F16]


clear LT LN gridedlay gridedDR lati longi
LT=[25:2:40];
LN=[60:2:90];


for i= 1:7
    for j=1:14
        
 id=find(dataspat(:,2)<LN(j+1) & dataspat(:,2)>LN(j) & dataspat(:,1)>LT(i) & dataspat(:,1)<LT(i+1)& dataspat(:,3)>2 );
  
    gridedlay(i,j)=length(id);
    gridedDR(i,j)=nanmedian(dataspat(id,4));
    lati(i,j)=mean(LT(i:i+1));
    longi(i,j)=mean(LN(j:j+1));
    
    end
end

load('/Volumes/Untitled/HMA/paper1_analysis/mycolormap.mat')

figure(1)
colormap(parula)
imagesc(longi(1,:),lati(:,1),gridedlay)
hold on
set(gca,'YDir','normal')
caxis([500 7000])
colormap('jet')
S=shaperead('/Volumes/Untitled/HMA/paper1_analysis/00_rgi60_regions/00_rgi60_O2Regions.shp');

plot([S.X],[S.Y],'LineWidth',2,'color','k','Linestyle','-')
xlim([67 88])
ylim([25 39])
set(gca,'fontsize',24)
hold on
xlabel('Longitude (^oE)')
ylabel('Latitude (^oN)')


figure(2)
id=gridedlay<2500
gridedDR(id)=NaN
colormap(parula)
imagesc(longi(1,:),lati(:,1),gridedDR)
hold on
set(gca,'YDir','normal')
caxis([.05 .3])
colormap('jet')

S1=shaperead('/Volumes/Untitled/EPA/HMS/world_countries_2017.shp');
 plot([S1.X],[S1.Y],'LineWidth',1,'color','k')
 
 hold on
 
 S1=shaperead('/Volumes/Untitled/HMA/India Shape/india_st.shp');
 plot([S1.X],[S1.Y],'LineWidth',2,'color','k')
 
 %S=shaperead('/Volumes/Untitled/HMA/paper1_analysis/00_rgi60_regions/00_rgi60_O2Regions.shp');

%plot([S.X],[S.Y],'LineWidth',2,'color','r','Linestyle','-')

 
%xlim([64 88])
%ylim([25 39])
set(gca,'fontsize',22)
hold on
xlabel('Longitude (^oE)')
ylabel('Latitude (^oN)')

%%

%%

y1w=load('2011_Dalb_profile_WH.mat','mati');
y2w=load('2012_Dalb_profile_WH.mat','mati');
y3w=load('2013_Dalb_profile_WH.mat','mati');
y4w=load('2015_Dalb_profile_WH.mat','mati');
y5w=load('2016_Dalb_profile_WH.mat','mati');
y6w=load('2014_Dalb_profile_WH.mat','mati');

y1h=load('2011_Dalb_profile_HK.mat','mati');
y2h=load('2012_Dalb_profile_HK.mat','mati');
y3h=load('2013_Dalb_profile_HK.mat','mati');
y4h=load('2015_Dalb_profile_HK.mat','mati');
y5h=load('2016_Dalb_profile_HK.mat','mati');
y6h=load('2014_Dalb_profile_HK.mat','mati');

y1k=load('2011_Dalb_profile_KK.mat','mati');
y2k=load('2012_Dalb_profile_KK.mat','mati');
y3k=load('2013_Dalb_profile_KK.mat','mati');
y4k=load('2015_Dalb_profile_KK.mat','mati');
y5k=load('2016_Dalb_profile_KK.mat','mati');
y6k=load('2014_Dalb_profile_KK.mat','mati');


snowalbedohgt=load('2016_Dalb_profile_KK.mat','h');


clear Ywhma profiwhma sumprofiwhma

Ywhma=[y1k.mati y2k.mati y3k.mati y4k.mati y5k.mati y6k.mati y1h.mati y2h.mati y3h.mati y4h.mati y5h.mati y6h.mati y1w.mati y2w.mati y3w.mati y4w.mati y5w.mati y6w.mati]

Ywhma(Ywhma==0 )= NaN
profiwhma=nanmedian(Ywhma,2);
sdeviwhma=nanstd(Ywhma');
countprofiwhma= sum(Ywhma>0 ,2)
sumprofiwhma=nansum(Ywhma,2);
id=(countprofiwhma<54)
    profiwhma(id)=NaN
    sumprofiwhma(id)=NaN
    
clear Yi Yii Yiii profi profii profiii
Yi=[ y1w.mati y2w.mati y3w.mati y4w.mati y5w.mati y6w.mati]

Yii=[y1k.mati y2k.mati y3k.mati y4k.mati y5k.mati y6k.mati ]

Yiii=[y1h.mati y2h.mati y3h.mati y4h.mati y5h.mati y6h.mati ]


Yi(Yi==0 )= NaN
profi=nanmedian(Yi,2);
sdevi=nanstd(Yi');
sumprofi=nansum(Yi,2);

Yii(Yii==0 )= NaN
profii=nanmedian(Yii,2);
sdevii=nanstd(Yii');
sumprofii=nansum(Yii,2);

Yiii(Yiii==0 )= NaN
profiii=nanmedian(Yiii,2);
sdeviii=nanstd(Yiii');
sumprofiii=nansum(Yiii,2);

countprofi= sum(Yi>0 ,2)
countprofii= sum(Yii>0 ,2)
countprofiii= sum(Yiii>0 ,2)

id=(countprofi<105)
    profi(id)=NaN

    id=(countprofii<36)
    profii(id)=NaN
    
    id=(countprofii<36)
    profiii(id)=NaN
    

DA_monthly=nanmean(Yi(20:end,:),1);

for i=1:48
DAmonthly(i)=(DA_monthly(i));%+DA_monthly(i+48)+DA_monthly(i+96))./3

end

id=find(~isnan(DAmonthly))
delta_albmonth=DAmonthly(id);




%%


clear ALb* bin1
bin1=[1:.05:7]
for j=1:length(bin1)-1
    
    id=find(snowalbedohgt.h/1000>bin1(j) & snowalbedohgt.h/1000<bin1(j+1)  );
   
    ALbwest(j)=(nanmean(sumprofiwhma(id)./countprofiwhma(id)));
    ALbweststd(j)=(nanstd(profiwhma(id)));
     %ALbeast(j)=nanmean(sumprof(id)./countprof(id));
  
        ALbwh(j)=(nanmedian(profi(id)));
        ALbwhstd(j)=nanmean(sdevi(id));
    
      
        ALbkk(j)=(nanmedian(profii(id)));
        ALbkkstd(j)=nanmean(sdevii(id));
    
        ALbhk(j)=nanmedian(profiii(id));
     ALbhkstd(j)=nanmean(sdeviii(id));
     
       
    
       
        
     
end


%%


clear h grided* bin
bin=[1:.15:7];

for i=1:6:60
    
    if i>=1 & i<=6
    y=2007
    
     for m=3:8
    
        
    A=load(sprintf('ALay_L2_05km_Himalya_%d_0%d_assembled_lat_lon_alt.dat.txt',y,m));
   
    data11=A;
    
for j=1:length(bin)-1
    
    id=find(data11(:,2)<95 & data11(:,2)>79 & data11(:,1)>25 & data11(:,1)<29 & data11(:,3)>bin(j) & data11(:,3)<bin(j+1));
    
    id1=find(data11(:,2)<79 & data11(:,2)>60 & data11(:,1)>25 & data11(:,1)<35 & data11(:,3)>bin(j) & data11(:,3)<bin(j+1));
    
    gridedlayeast(i,j)=length(id);
    gridedlaywest(i,j)=length(id1);
    gridedDReast(i,j)=nanmean(data11(id,4));
    gridedDRwest(i,j)=nanmean(data11(id1,4));
   h(j)=nanmean(bin(j:j+1));
   
end
i=i+1;
     end
    
     
    elseif i>=7 & i<=12
    y=2008
    
    
     for m=3:8
    
        
    A=load(sprintf('ALay_L2_05km_Himalya_%d_0%d_assembled_lat_lon_alt.dat.txt',y,m));
   
    data11=A;
    
for j=1:length(bin)-1
    
    id=find(data11(:,2)<95 & data11(:,2)>79 & data11(:,1)>25 & data11(:,1)<29 & data11(:,3)>bin(j) & data11(:,3)<bin(j+1));
    
    id1=find(data11(:,2)<79 & data11(:,2)>60 & data11(:,1)>25 & data11(:,1)<35 & data11(:,3)>bin(j) & data11(:,3)<bin(j+1));
    
    gridedlayeast(i,j)=length(id);
    gridedlaywest(i,j)=length(id1);
    gridedDReast(i,j)=nanmean(data11(id,4));
    gridedDRwest(i,j)=nanmean(data11(id1,4));
   h(j)=nanmean(bin(j:j+1));
   
end
i=i+1;
     end
    
     
    elseif i>=13 & i<=18
    y=2009
    
    
     for m=3:8
    
        
    A=load(sprintf('ALay_L2_05km_Himalya_%d_0%d_assembled_lat_lon_alt.dat.txt',y,m));
   
    data11=A;
    
for j=1:length(bin)-1
    
    id=find(data11(:,2)<95 & data11(:,2)>79 & data11(:,1)>25 & data11(:,1)<29 & data11(:,3)>bin(j) & data11(:,3)<bin(j+1));
    
    id1=find(data11(:,2)<79 & data11(:,2)>60 & data11(:,1)>25 & data11(:,1)<35 & data11(:,3)>bin(j) & data11(:,3)<bin(j+1));
    
    gridedlayeast(i,j)=length(id);
    gridedlaywest(i,j)=length(id1);
    gridedDReast(i,j)=nanmean(data11(id,4));
    gridedDRwest(i,j)=nanmean(data11(id1,4));
   h(j)=nanmean(bin(j:j+1));
   
end
i=i+1;
     end
    
     
    elseif i>=19 & i<=24
    y=2010
    
    
     for m=3:8
    
        
    A=load(sprintf('ALay_L2_05km_Himalya_%d_0%d_assembled_lat_lon_alt.dat.txt',y,m));
   
    data11=A;
    
for j=1:length(bin)-1
    
    id=find(data11(:,2)<95 & data11(:,2)>79 & data11(:,1)>25 & data11(:,1)<29 & data11(:,3)>bin(j) & data11(:,3)<bin(j+1));
    
    id1=find(data11(:,2)<79 & data11(:,2)>60 & data11(:,1)>25 & data11(:,1)<35 & data11(:,3)>bin(j) & data11(:,3)<bin(j+1));
    
    gridedlayeast(i,j)=length(id);
    gridedlaywest(i,j)=length(id1);
    gridedDReast(i,j)=nanmean(data11(id,4));
    gridedDRwest(i,j)=nanmean(data11(id1,4));
   h(j)=nanmean(bin(j:j+1));
   
end
i=i+1;
     end
    
     
    elseif i>=25 & i<=30
    y=2011
    
    
     for m=3:8
    
        
    A=load(sprintf('ALay_L2_05km_Himalya_%d_0%d_assembled_lat_lon_alt.dat.txt',y,m));
   
    data11=A;
    
for j=1:length(bin)-1
    
    id=find(data11(:,2)<95 & data11(:,2)>79 & data11(:,1)>25 & data11(:,1)<29 & data11(:,3)>bin(j) & data11(:,3)<bin(j+1));
    
    id1=find(data11(:,2)<79 & data11(:,2)>60 & data11(:,1)>25 & data11(:,1)<35 & data11(:,3)>bin(j) & data11(:,3)<bin(j+1));
    
    gridedlayeast(i,j)=length(id);
    gridedlaywest(i,j)=length(id1);
    gridedDReast(i,j)=nanmean(data11(id,4));
    gridedDRwest(i,j)=nanmean(data11(id1,4));
   h(j)=nanmean(bin(j:j+1));
   
end
i=i+1;
     end
    
     
    elseif i>=31 & i<=36
    y=2012
    
    
     for m=3:8
    
        
    A=load(sprintf('ALay_L2_05km_Himalya_%d_0%d_assembled_lat_lon_alt.dat.txt',y,m));
   
    data11=A;
    
for j=1:length(bin)-1
    
    id=find(data11(:,2)<95 & data11(:,2)>79 & data11(:,1)>25 & data11(:,1)<29 & data11(:,3)>bin(j) & data11(:,3)<bin(j+1));
    
    id1=find(data11(:,2)<79 & data11(:,2)>60 & data11(:,1)>25 & data11(:,1)<35 & data11(:,3)>bin(j) & data11(:,3)<bin(j+1));
    
    gridedlayeast(i,j)=length(id);
    gridedlaywest(i,j)=length(id1);
    gridedDReast(i,j)=nanmean(data11(id,4));
    gridedDRwest(i,j)=nanmean(data11(id1,4));
   h(j)=nanmean(bin(j:j+1));
   
end
i=i+1;
     end
    
     
    elseif i>=37 & i<=42
    y=2013
    
    
     for m=3:8
    
        
    A=load(sprintf('ALay_L2_05km_Himalya_%d_0%d_assembled_lat_lon_alt.dat.txt',y,m));
   
    data11=A;
    
for j=1:length(bin)-1
    
    id=find(data11(:,2)<95 & data11(:,2)>79 & data11(:,1)>25 & data11(:,1)<29 & data11(:,3)>bin(j) & data11(:,3)<bin(j+1));
    
    id1=find(data11(:,2)<79 & data11(:,2)>60 & data11(:,1)>25 & data11(:,1)<35 & data11(:,3)>bin(j) & data11(:,3)<bin(j+1));
    
    gridedlayeast(i,j)=length(id);
    gridedlaywest(i,j)=length(id1);
    gridedDReast(i,j)=nanmean(data11(id,4));
    gridedDRwest(i,j)=nanmean(data11(id1,4));
   h(j)=nanmean(bin(j:j+1));
   
end
i=i+1;
     end
    
     
    elseif i>=43 & i<=48
    y=2014
    
     for m=3:8
    
        
    A=load(sprintf('ALay_L2_05km_Himalya_%d_0%d_assembled_lat_lon_alt.dat.txt',y,m));
   
    data11=A;
    
for j=1:length(bin)-1
    
    id=find(data11(:,2)<95 & data11(:,2)>79 & data11(:,1)>25 & data11(:,1)<29 & data11(:,3)>bin(j) & data11(:,3)<bin(j+1));
    
    id1=find(data11(:,2)<79 & data11(:,2)>60 & data11(:,1)>25 & data11(:,1)<35 & data11(:,3)>bin(j) & data11(:,3)<bin(j+1));
    
    gridedlayeast(i,j)=length(id);
    gridedlaywest(i,j)=length(id1);
    gridedDReast(i,j)=nanmean(data11(id,4));
    gridedDRwest(i,j)=nanmean(data11(id1,4));
   h(j)=nanmean(bin(j:j+1));
   
end
i=i+1;
     end
    
        elseif i>=49 & i<=54
            
    y=2015
     for m=3:8
    
        
    A=load(sprintf('ALay_L2_05km_Himalya_%d_0%d_assembled_lat_lon_alt.dat.txt',y,m));
   
    data11=A;
    
for j=1:length(bin)-1
    
    id=find(data11(:,2)<95 & data11(:,2)>79 & data11(:,1)>25 & data11(:,1)<29 & data11(:,3)>bin(j) & data11(:,3)<bin(j+1));
    
    id1=find(data11(:,2)<79 & data11(:,2)>60 & data11(:,1)>25 & data11(:,1)<35 & data11(:,3)>bin(j) & data11(:,3)<bin(j+1));
    
    gridedlayeast(i,j)=length(id);
    gridedlaywest(i,j)=length(id1);
    gridedDReast(i,j)=nanmean(data11(id,4));
    gridedDRwest(i,j)=nanmean(data11(id1,4));
   h(j)=nanmean(bin(j:j+1));
   
end
i=i+1;
    end
    
        else
            y=2016
            
         for m=3:8
    
        
    A=load(sprintf('ALay_L2_05km_Himalya_%d_0%d_assembled_lat_lon_alt.dat.txt',y,m));
   
    data11=A;
    
for j=1:length(bin)-1
    
    id=find(data11(:,2)<95 & data11(:,2)>=79 & data11(:,1)>25 & data11(:,1)<=29 & data11(:,3)>bin(j) & data11(:,3)<=bin(j+1));
    
    id1=find(data11(:,2)<79 & data11(:,2)>=60 & data11(:,1)>25 & data11(:,1)<=35 & data11(:,3)>bin(j) & data11(:,3)<=bin(j+1));
    
    gridedlayeast(i,j)=length(id);
    gridedlaywest(i,j)=length(id1);
    gridedDReast(i,j)=nanmean(data11(id,4));
    gridedDRwest(i,j)=nanmean(data11(id1,4));
   h(j)=nanmean(bin(j:j+1));
   
end
i=i+1;
    end
        
        
    end
  
   
end
%%
errorbar(nanmean(gridedlaywest,1),bin(1:end-1),nanstd(gridedlaywest,1)./length(),'k--','horizontal','Linewidth',3,'color',[255/255,20/255,147/255])
hold on
%errorbar(nanmean(gridedlayeast,1),bin(1:end-1),nanstd(gridedlayeast,1),'b--','horizontal','Linewidth',3,'color',[215/255,183/255,64/255])
hold on
xlim([-50 400])
ylim([1.5 6.5])
set(gca,'fontsize',25)
hold on
xlabel(' EALs (#)')
    ylabel('Altitude (km)')
  %   text(.65,.12,'WH+KK+HK','Fontsize',25,'Units','normalized','Color',[255/255,20/255,147/255])

 % text(.75,.22,'EH+CH','Fontsize',25,'Units','normalized','color',[215/255,183/255,64/255])

%%

errorbar(nanmean(gridedDRwest,1),bin(2:end),nanstd(gridedDRwest,1),'k--','horizontal','Linewidth',3)
hold on
errorbar(nanmean(gridedDReast,1),bin(2:end),nanstd(gridedDReast,1),'b--','horizontal','Linewidth',3)
hold on
errorbar(gridedlay2,bin(2:end),1,'horizontal','S-')

%%
clear EALmonthly EALmonthly_DR
EALmonthly=nansum(gridedlaywest(:,10:end),2);
EALmonthly_DR=nanmedian(gridedDRwest(:,10:end),2);
%%
v=(EALmonthly(25:54))
Xd=prctile(v,[0:10:100])

for i=1:10
    
   id=find(v>Xd(i) & v<=Xd(i+1)) 
    
    EALs(i)=nanmean(v(id));
    SA(i)=nanmean((delta_albmonth(id)));
     EALssd(i)=nanstd(v(id));
    SAsd(i)=nanstd((delta_albmonth(id)));
end

errorbar((EALs),(SA),SAsd,SAsd,EALssd,EALssd)

%%
figure(3)
pos3 = [0.1 0.1 0.37 0.27];
pos1 = [0.1 0.5 0.49 0.49];

pos4 = [0.6 0.1 0.37 0.27];
pos2 = [0.55 0.5 0.4 0.52];


subplot('position',pos4)

%%
scatter((EALmonthly(25:60)),(EALmonthly_DR(25:60)),200,'MarkerFaceColor',[255/255,20/255,147/255],'MarkerEdgeColor',[255/255,20/255,147/255])

xlim([2000 8000])
ylim([0.12 .2])
set(gca,'fontsize',22)
hold on
[p2 S]= polyfit((EALmonthly(25:60)),(EALmonthly_DR(25:60)),1)
 [y delta] = polyval(p2,(EALmonthly(25:60)),S);
hold on
plot((EALmonthly(25:60)),y,'--k','Linewidth',5)
box on
[px2 rmse] = corrcoef((EALmonthly(25:60)),(EALmonthly_DR(25:60))); 
%caption = sprintf('Y = %.2f * X + %.2f', p2(1), p2(2));
%text(.5,.05, caption, 'Fontsize',15,'Units','normalized','color','b');

 text(.75,.12,[ 'R = ',num2str(px2(1,2),'%.2f') ],'Fontsize',20,'Units','normalized','color','k')

  xlabel(' EALs (#)')
    ylabel('Depolarization Ratio')
    
 %%

subplot('position',pos3)
%%
scatter((EALmonthly(25:60)),(delta_albmonth(1:36)),200,(EALmonthly_DR(25:60)),'MarkerFaceColor',[255/255,20/255,147/255],'MarkerEdgeColor',[255/255,20/255,147/255])
%xlim([1000 11000])
ylim([0.4 4])
set(gca,'fontsize',22)
hold on
[p2 S]= polyfit((EALmonthly(25:54)),(delta_albmonth(1:30)),1)
 [y delta] = polyval(p2,(EALmonthly(25:54)),S);
hold on
plot((EALmonthly(25:54)),y,'--k','Linewidth',5)
box on
[px2 rmse] = corrcoef((EALmonthly(25:54)),(delta_albmonth(1:30))); 
%caption = sprintf('Y = %.2f * X + %.2f', p2(1), p2(2));
%text(.5,.05, caption, 'Fontsize',15,'Units','normalized','color','b');

 text(.75,.12,[ 'R = ',num2str(px2(1,2),'%.2f') ],'Fontsize',20,'Units','normalized','color','k')

  xlabel(' EALs (#)')
    ylabel('\Delta \alpha (%)')
    
 %%   
 subplot('position',pos1)   
 %%
figure(1)
errorbar(nanmedian(gridedlaywest,1),bin(1:end-1),nanstd(gridedlaywest,1)./sqrt(60),'k-','horizontal','Linewidth',4,'color',[255/255,20/255,147/255])
hold on
%errorbar(nanmedian(gridedlayeast,1),bin(1:end-1),nanstd(gridedlayeast,1)./sqrt(60),'b-','horizontal','Linewidth',3,'color',[215/255,183/255,64/255])

%plot(nanmedian(gridedlaywest,1),bin(1:end-1)+.125,'o-','Linewidth',3,'color',[240/255,20/255,147/255])
%plot(prctile((gridedlaywest),50)./max(prctile((gridedlaywest),50)),bin(1:end-1)+.125,'o--','Linewidth',3,'color',[240/255,20/255,147/255])
hold on
hold on
%plot(prctile((gridedlaywest),25)./max(prctile((gridedlaywest),25)),bin(1:end-1)+.125,'o--','Linewidth',1,'color',[240/255,20/255,147/255])
hold on
%plot(prctile((gridedlaywest),75)./max(prctile((gridedlaywest),75)),bin(1:end-1)+.125,'x--','Linewidth',1,'color',[240/255,20/255,147/255])

%bar(nanmedian(gridedlaywest,1)./max(nanmedian(gridedlaywest,1)),'color',[240/255,20/255,147/255])

hold on
%xlim([0.2 1])
xlim([0 400])
ylim([1 6.5])
set(gca,'fontsize',22)
hold on
xlabel(' EALs (#)')
    ylabel('Altitude (km)')
grid on
text(.65,.2,'Indus basin box','Fontsize',20,'Units','normalized','Color',[255/255,20/255,147/255])
text(.65,.15,'(WH+HK+KK)','Fontsize',20,'Units','normalized','Color',[255/255,20/255,147/255])
%text(.1,.85,'Ganga basin','Fontsize',20,'Units','normalized','Color',[215/255,183/255,64/255])
%text(.1,.8,'(CH+EH)','Fontsize',20,'Units','normalized','Color',[215/255,183/255,64/255])


figure(2)

box on
  % errorbar((X2),bin1(1:end-1)+.125,(ALbchstd+ALbehstd)./sqrt(36),'horizontal','s-','Linewidth',3,'color',[215/255,183/255,64/255])
hold on
%errorbar((X),bin1(1:end-1)+.125,(ALbwhstd+ALbhkstd+ALbkkstd)./sqrt(36),'horizontal','s-','Linewidth',3,'color',[255/255,20/255,147/255])

%errorbar(ALbch*2,bin1(1:end-1)+.125,ALbchstd./sqrt(36),'horizontal','s-','Linewidth',3,'color','g')
%errorbar(ALbwh*3,bin1(1:end-1)+.125,ALbwhstd./sqrt(36),'horizontal','s-','Linewidth',3,'color','r')
%errorbar(ALbhk*4,bin1(1:end-1)+.125,ALbhkstd./sqrt(36),'horizontal','s-','Linewidth',3,'color','b')
%errorbar(ALbkk*5,bin1(1:end-1)+.125,ALbkkstd./sqrt(36),'horizontal','s-','Linewidth',3,'color','k')

errorbar(profiwhma,snowalbedohgt.h/1000,sdeviwhma'./sqrt(countprofiwhma),'horizontal','o-','Linewidth',3,'Color',[255/255,20/255,147/255])
hold on

%plot(profiwhma,snowalbedohgt.h/1000,'o-','Linewidth',3,'Color',[255/255,20/255,147/255])
hold on
%errorbar(profehma,snowalbedohgt.h/1000,sdevehma'./sqrt(36),'horizontal','o-','Linewidth',3,'Color',[215/255,183/255,64/255])
%errorbar((ALbwest),bin1(1:end-1)+.125,ALbweststd,'horizontal','o-','Linewidth',3,'Color',[255/255,20/255,147/255])
hold on

%plot((profiwhma),snowalbedohgt.h/1000,'o-','Linewidth',3,'Color',[255/255,20/255,147/255])
hold on
%plot(profehma./max(profehma),snowalbedohgt.h/1000,'Linewidth',3,'Color',[215/255,183/255,64/255])
%plot(nanmedian(gridedlaywest,1)./max(nanmedian(gridedlaywest,1)),bin(1:end-1)+.125,'x--','Linewidth',3,'color',[240/255,20/255,147/255])



%plot((profehma./max(profehma)),snowalbedohgt.h/1000,'o-','Linewidth',3,'Color',[255/255,20/255,147/255])
hold on
%plot(profehma./max(profehma),snowalbedohgt.h/1000,'Linewidth',3,'Color',[215/255,183/255,64/255])
%plot(nanmedian(gridedlayeast,1)./max(nanmedian(gridedlayeast,1)),bin(1:end-1)+.125,'x--','Linewidth',3,'color',[240/255,20/255,147/255])




    hold on
    %errorbar(ALbwest.*3,bin1(1:end-1)+.125,ALbweststd./sqrt(36*3),'horizontal','s-','Linewidth',3,'Color',[255/255,20/255,147/255])
  

xlim([0.1 4.5])
%xlim([0.2 1])
ylim([1 6.5])
set(gca,'fontsize',22)
hold on
hold on
ylabel('Altitude (km)')
xlabel('\Delta \alpha (%)')
%xlabel('\Delta \alpha (normalized)')
  grid on
text(.70,.15,'HK+WH+KK','Fontsize',20,'Units','normalized','Color',[255/255,20/255,147/255])
%text(.70,.25,'CH+EH','Fontsize',20,'Units','normalized','Color',[215/255,183/255,64/255])


  %%
 for i=1:length(profi)
    
    X(i)=sum(profi(i)+profii(i)+profiii(i));
    X2(i)=sum(prof(i)+profch(i));
    
 end

 
 figure(2)
    errorbar(X2,snowalbedohgt.h/1000,sdevch./sqrt(36),'horizontal','s-','Linewidth',2,'color',[215/255,183/255,64/255])
%errorbar(ALbeast,bin1(1:end-1)+.125,ALbeaststd./sqrt(36*2),'horizontal','s-','Linewidth',3,'color',[215/255,183/255,64/255])
hold on
errorbar(X,snowalbedohgt.h/1000,sdev./sqrt(36),'horizontal','s-','Linewidth',2,'color',[255/255,20/255,147/255])

    hold on
    %errorbar(ALbwest.*3,bin1(1:end-1)+.125,ALbweststd./sqrt(36*3),'horizontal','s-','Linewidth',3,'Color',[255/255,20/255,147/255])
  
%errorbar(sumprofi,snowalbedohgt.h/1000,sdevi./sqrt(36),'horizontal','o-','Linewidth',2,'Color','r')

%errorbar(sumprofii,snowalbedohgt.h/1000,sdevii./sqrt(36),'horizontal','o-','Linewidth',2,'Color','k')

%errorbar(sumprofiii,snowalbedohgt.h/1000,sdeviii./sqrt(36),'horizontal','o-','Linewidth',2,'Color','b')

xlim([0 15])
ylim([2.5 6.5])
set(gca,'fontsize',22)
hold on
hold on
ylabel('Altitude (km)')
xlabel('\Delta \alpha (%)')
  grid on

%%
figure(3)
   errorbar((sumprofiii+sumprofi+sumprofii),snowalbedohgt.h/1000,sdevi','horizontal','o-','Linewidth',3,'Color',[255/255,20/255,147/255])
hold on
errorbar(sumprofehma,snowalbedohgt.h/1000,sdevi','horizontal','o-','Linewidth',3,'Color',[215/255,183/255,64/255])
xlim([0 500])
ylim([1 6])
set(gca,'fontsize',22)
hold on
hold on
ylabel('Altitude (km)')
xlabel('\Delta \alpha (%)')
  
%text(.65,.95,'WH+KK+HK','Fontsize',20,'Units','normalized','Color',[255/255,20/255,147/255])
%text(.75,.9,'EH+CH','Fontsize',20,'Units','normalized','color',[215/255,183/255,64/255])


grid on
%%

    figure(4)
 %errorbar(nanmedian(gridedlaywest,1),bin(1:end-1),nanstd(gridedlaywest,1),'k--','horizontal','Linewidth',3,'color',[255/255,20/255,147/255])


    errorbar(countprofehma,snowalbedohgt.h/1000,sdevch./sqrt(36),'horizontal','s-','Linewidth',2,'color',[215/255,183/255,64/255])
%errorbar(ALbeast,bin1(1:end-1)+.125,ALbeaststd./sqrt(36*2),'horizontal','s-','Linewidth',3,'color',[215/255,183/255,64/255])
hold on
errorbar(countprofiwhma,snowalbedohgt.h/1000,sdev./sqrt(36),'horizontal','s-','Linewidth',2,'color',[255/255,20/255,147/255])
hold on
xlim([0 650])
ylim([1 5])
set(gca,'fontsize',22)
hold on
xlabel(' EALs (#)')
    ylabel('Altitude (km)')
grid on
text(.70,.15,'EH+CH','Fontsize',20,'Units','normalized','color',[215/255,183/255,64/255])
 % subplot('position',pos2) 
 

%%
annotation('textbox',[0.01 .95 0.25 .05],'String',{'A)'},'LineStyle','none','FontSize',20);
annotation('textbox',[0.52 .95 0.25 .05],'String',{'B)'},'LineStyle','none','FontSize',20);
annotation('textbox',[0.01 .4 0.25 .05],'String',{'C)'},'LineStyle','none','FontSize',20);
annotation('textbox',[0.52 .4 0.25 .05],'String',{'D)'},'LineStyle','none','FontSize',20);