
%%
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 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

%%

%%%%   S1B
%clear EALmonthly EALmonthly_DR
%EALmonthly=nansum(gridedlaywest(:,10:end),2);
%EALmonthly_DR=nanmedian(gridedDRwest(:,10:end),2);

% 
% 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')
    

 %%
figure(3)
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

%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(4)

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(profiwhma,snowalbedohgt.h/1000,sdeviwhma'./sqrt(countprofiwhma),'horizontal','o-','Linewidth',3,'Color',[255/255,20/255,147/255])
hold on



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])


 