Untitled

clc; clear all;clf;close all;

EQualMomWithRTI=figure(1); set(EQualMomWithRTI, 'PaperType', 'a4'); set(EQualMomWithRTI, 'PaperOrientation', 'portrait');set(EQualMomWithRTI, 'Resize', 'on');set(EQualMomWithRTI, 'Position',[0 0 2480*(80/100) 3508*(80/100)]);

[A,B,C]=xlsread('Correlation_FullDayData.xlsx',1);

ht1= A(1,(1:44))';

nrgs = numel(ht1);

[A,B,C]=xlsread('Correlation_FullDayData.xlsx',1);

ht= A(1,(1:44))';

WindSpeedfreqNewComp_5 =A((2:10),(1:44));WindSpeedfreqNewComp_6 =A((11:63),(1:44));WindSpeedfreqNewComp_7 =A((64:115),(1:44));WindSpeedfreqNewComp_8 =A((116:167),(1:44));WindSpeedfreqNewComp_9 =A((168:220),(1:44));WindSpeedfreqNewComp_10 =A((221:272),(1:44));WindSpeedfreqNewComp_11 =A((273:324),(1:44));WindSpeedfreqNewComp_12 =A((325:376),(1:44));WindSpeedfreqNewComp_13 =A((377:429),(1:44));WindSpeedfreqNewComp_14 =A((430:481),(1:44));WindSpeedfreqNewComp_15 =A((482:533),(1:44));WindSpeedfreqNewComp_16 =A((534:586),(1:44));WindSpeedfreqNewComp_17 =A((587:611),(1:44));

[B,C,D]=xlsread('Correlation_FullDayData.xlsx',2);

ZonalfreqNewComp_5 =B((2:10),(1:44));ZonalfreqNewComp_6 =B((11:63),(1:44));ZonalfreqNewComp_7 =B((64:115),(1:44));ZonalfreqNewComp_8 =B((116:167),(1:44));ZonalfreqNewComp_9 =B((168:220),(1:44));ZonalfreqNewComp_10 =B((221:272),(1:44));ZonalfreqNewComp_11 =B((273:324),(1:44));ZonalfreqNewComp_12 =B((325:376),(1:44));ZonalfreqNewComp_13 =B((377:429),(1:44));ZonalfreqNewComp_14 =B((430:481),(1:44));ZonalfreqNewComp_15 =B((482:533),(1:44));ZonalfreqNewComp_16 =B((534:586),(1:44));ZonalfreqNewComp_17 =B((587:611),(1:44));

[C,D,E]=xlsread('Correlation_FullDayData.xlsx',3);

MeridionalfreqNewComp_5 =C((2:10),(1:44));MeridionalfreqNewComp_6 =C((11:63),(1:44));MeridionalfreqNewComp_7 =C((64:115),(1:44));MeridionalfreqNewComp_8 =C((116:167),(1:44));MeridionalfreqNewComp_9 =C((168:220),(1:44));MeridionalfreqNewComp_10 =C((221:272),(1:44));MeridionalfreqNewComp_11 =C((273:324),(1:44));MeridionalfreqNewComp_12 =C((325:376),(1:44));MeridionalfreqNewComp_13 =C((377:429),(1:44));MeridionalfreqNewComp_14 =C((430:481),(1:44));MeridionalfreqNewComp_15 =C((482:533),(1:44));MeridionalfreqNewComp_16 =C((534:586),(1:44));MeridionalfreqNewComp_17 =C((587:611),(1:44));

[A,B,C]=xlsread('Correlation_FullDayData.xlsx',4);

gpsWindSpeed_5 = A((2),(1:44))';gpsWindSpeed_6 = A((3),(1:44))';gpsWindSpeed_7 = A((4),(1:44))';gpsWindSpeed_8 = A((5),(1:44))';gpsWindSpeed_9 = A((6),(1:44))';gpsWindSpeed_10 = A((7),(1:44))';gpsWindSpeed_11 = A((8),(1:44))';gpsWindSpeed_12 = A((9),(1:44))';gpsWindSpeed_13 = A((10),(1:44))';gpsWindSpeed_14 = A((11),(1:44))';gpsWindSpeed_15 = A((12),(1:44))';gpsWindSpeed_16 = A((13),(1:44))';gpsWindSpeed_17 = A((14),(1:44))';

[A,B,C]=xlsread('Correlation_FullDayData.xlsx',5);

gpsZonal_5 = A((2),(1:44))';gpsZonal_6 = A((3),(1:44))';gpsZonal_7 = A((4),(1:44))';gpsZonal_8 = A((5),(1:44))';gpsZonal_9 = A((6),(1:44))';gpsZonal_10 = A((7),(1:44))';gpsZonal_11 = A((8),(1:44))';gpsZonal_12 = A((9),(1:44))';gpsZonal_13 = A((10),(1:44))';gpsZonal_14 = A((11),(1:44))';gpsZonal_15 = A((12),(1:44))';gpsZonal_16 = A((13),(1:44))';gpsZonal_17 = A((14),(1:44))';

[A,B,C]=xlsread('Correlation_FullDayData.xlsx',6);

gpsMeridional_5 = A((2),(1:44))';gpsMeridional_6 = A((3),(1:44))';gpsMeridional_7 = A((4),(1:44))';gpsMeridional_8 = A((5),(1:44))';gpsMeridional_9 = A((6),(1:44))';gpsMeridional_10 = A((7),(1:44))';gpsMeridional_11 = A((8),(1:44))';gpsMeridional_12 = A((9),(1:44))';gpsMeridional_13 = A((10),(1:44))';gpsMeridional_14 = A((11),(1:44))';gpsMeridional_15 = A((12),(1:44))';gpsMeridional_16 = A((13),(1:44))';gpsMeridional_17 = A((14),(1:44))';

mYWSarr1 = [WindSpeedfreqNewComp_5; WindSpeedfreqNewComp_6; WindSpeedfreqNewComp_7; WindSpeedfreqNewComp_8; WindSpeedfreqNewComp_9;WindSpeedfreqNewComp_10;WindSpeedfreqNewComp_11;WindSpeedfreqNewComp_12;WindSpeedfreqNewComp_13;WindSpeedfreqNewComp_14;WindSpeedfreqNewComp_15;WindSpeedfreqNewComp_16;WindSpeedfreqNewComp_17;];

mYUarr1 = [ZonalfreqNewComp_5; ZonalfreqNewComp_6; ZonalfreqNewComp_7; ZonalfreqNewComp_8; ZonalfreqNewComp_9;ZonalfreqNewComp_10;ZonalfreqNewComp_11;ZonalfreqNewComp_12;ZonalfreqNewComp_13;ZonalfreqNewComp_14;ZonalfreqNewComp_15;ZonalfreqNewComp_16;ZonalfreqNewComp_17;];

mYVarr1 = [MeridionalfreqNewComp_5; MeridionalfreqNewComp_6; MeridionalfreqNewComp_7; MeridionalfreqNewComp_8; MeridionalfreqNewComp_9;MeridionalfreqNewComp_10;MeridionalfreqNewComp_11;MeridionalfreqNewComp_12;MeridionalfreqNewComp_13;MeridionalfreqNewComp_14;MeridionalfreqNewComp_15;MeridionalfreqNewComp_16;MeridionalfreqNewComp_17;];

mYWSarr1 = mYWSarr1(:,(1:44));mYUarr1 = mYUarr1(:,(1:44));mYVarr1 = mYVarr1(:,(1:44));

WindSpeedfreqNewComp_Avg_5=mean(WindSpeedfreqNewComp_5,1)';

WindSpeedfreqNewComp_Avg_6=mean(WindSpeedfreqNewComp_6,1)';

WindSpeedfreqNewComp_Avg_7=mean(WindSpeedfreqNewComp_7,1)';

WindSpeedfreqNewComp_Avg_8=mean(WindSpeedfreqNewComp_8,1)';

WindSpeedfreqNewComp_Avg_9=mean(WindSpeedfreqNewComp_9,1)';

WindSpeedfreqNewComp_Avg_10=mean(WindSpeedfreqNewComp_10,1)';

WindSpeedfreqNewComp_Avg_11=mean(WindSpeedfreqNewComp_11,1)';

WindSpeedfreqNewComp_Avg_12=mean(WindSpeedfreqNewComp_12,1)';

WindSpeedfreqNewComp_Avg_13=mean(WindSpeedfreqNewComp_13,1)';

WindSpeedfreqNewComp_Avg_14=mean(WindSpeedfreqNewComp_14,1)';

WindSpeedfreqNewComp_Avg_15=mean(WindSpeedfreqNewComp_15,1)';

WindSpeedfreqNewComp_Avg_16=mean(WindSpeedfreqNewComp_16,1)';

WindSpeedfreqNewComp_Avg_17=mean(WindSpeedfreqNewComp_17,1)';

ZonalfreqNewComp_Avg_5=mean(ZonalfreqNewComp_5,1)';

ZonalfreqNewComp_Avg_6=mean(ZonalfreqNewComp_6,1)';

ZonalfreqNewComp_Avg_7=mean(ZonalfreqNewComp_7,1)';

ZonalfreqNewComp_Avg_8=mean(ZonalfreqNewComp_8,1)';

ZonalfreqNewComp_Avg_9=mean(ZonalfreqNewComp_9,1)';

ZonalfreqNewComp_Avg_10=mean(ZonalfreqNewComp_10,1)';

ZonalfreqNewComp_Avg_11=mean(ZonalfreqNewComp_11,1)';

ZonalfreqNewComp_Avg_12=mean(ZonalfreqNewComp_12,1)';

ZonalfreqNewComp_Avg_13=mean(ZonalfreqNewComp_13,1)';

ZonalfreqNewComp_Avg_14=mean(ZonalfreqNewComp_14,1)';

ZonalfreqNewComp_Avg_15=mean(ZonalfreqNewComp_15,1)';

ZonalfreqNewComp_Avg_16=mean(ZonalfreqNewComp_16,1)';

ZonalfreqNewComp_Avg_17=mean(ZonalfreqNewComp_17,1)';

MeridionalfreqNewComp_Avg_5=mean(MeridionalfreqNewComp_5,1)';

MeridionalfreqNewComp_Avg_6=mean(MeridionalfreqNewComp_6,1)';

MeridionalfreqNewComp_Avg_7=mean(MeridionalfreqNewComp_7,1)';

MeridionalfreqNewComp_Avg_8=mean(MeridionalfreqNewComp_8,1)';

MeridionalfreqNewComp_Avg_9=mean(MeridionalfreqNewComp_9,1)';

MeridionalfreqNewComp_Avg_10=mean(MeridionalfreqNewComp_10,1)';

MeridionalfreqNewComp_Avg_11=mean(MeridionalfreqNewComp_11,1)';

MeridionalfreqNewComp_Avg_12=mean(MeridionalfreqNewComp_12,1)';

MeridionalfreqNewComp_Avg_13=mean(MeridionalfreqNewComp_13,1)';

MeridionalfreqNewComp_Avg_14=mean(MeridionalfreqNewComp_14,1)';

MeridionalfreqNewComp_Avg_15=mean(MeridionalfreqNewComp_15,1)';

MeridionalfreqNewComp_Avg_16=mean(MeridionalfreqNewComp_16,1)';

MeridionalfreqNewComp_Avg_17=mean(MeridionalfreqNewComp_17,1)';

mYWSarr1_t = [WindSpeedfreqNewComp_Avg_5 WindSpeedfreqNewComp_Avg_6 WindSpeedfreqNewComp_Avg_7 WindSpeedfreqNewComp_Avg_8 WindSpeedfreqNewComp_Avg_9 WindSpeedfreqNewComp_Avg_10 WindSpeedfreqNewComp_Avg_11 WindSpeedfreqNewComp_Avg_12 WindSpeedfreqNewComp_Avg_13 WindSpeedfreqNewComp_Avg_14 WindSpeedfreqNewComp_Avg_15 WindSpeedfreqNewComp_Avg_16 WindSpeedfreqNewComp_Avg_17 ]';

mYUarr1_t = [ZonalfreqNewComp_Avg_5 ZonalfreqNewComp_Avg_6 ZonalfreqNewComp_Avg_7 ZonalfreqNewComp_Avg_8 ZonalfreqNewComp_Avg_9 ZonalfreqNewComp_Avg_10 ZonalfreqNewComp_Avg_11 ZonalfreqNewComp_Avg_12 ZonalfreqNewComp_Avg_13 ZonalfreqNewComp_Avg_14 ZonalfreqNewComp_Avg_15 ZonalfreqNewComp_Avg_16 ZonalfreqNewComp_Avg_17]';

mYVarr1_t = [MeridionalfreqNewComp_Avg_5 MeridionalfreqNewComp_Avg_6 MeridionalfreqNewComp_Avg_7 MeridionalfreqNewComp_Avg_8 MeridionalfreqNewComp_Avg_9 MeridionalfreqNewComp_Avg_10 MeridionalfreqNewComp_Avg_11 MeridionalfreqNewComp_Avg_12 MeridionalfreqNewComp_Avg_13 MeridionalfreqNewComp_Avg_14 MeridionalfreqNewComp_Avg_15 MeridionalfreqNewComp_Avg_16 MeridionalfreqNewComp_Avg_17]';

mYWSarr1_t = mYWSarr1_t(:,(1:44));mYUarr1_t = mYUarr1_t(:,(1:44));mYVarr1_t = mYVarr1_t(:,(1:44));

gpsWS1 = [gpsWindSpeed_5 gpsWindSpeed_6 gpsWindSpeed_7 gpsWindSpeed_8 gpsWindSpeed_9 gpsWindSpeed_10 gpsWindSpeed_11 gpsWindSpeed_12 gpsWindSpeed_13 gpsWindSpeed_14 gpsWindSpeed_15 gpsWindSpeed_16 gpsWindSpeed_17]';

gpsU1 = [gpsZonal_5 gpsZonal_6 gpsZonal_7 gpsZonal_8 gpsZonal_9 gpsZonal_10 gpsZonal_11 gpsZonal_12 gpsZonal_13 gpsZonal_14 gpsZonal_15 gpsZonal_16 gpsZonal_17]';

gpsV1 = [gpsMeridional_5 gpsMeridional_6 gpsMeridional_7 gpsMeridional_8 gpsMeridional_9 gpsMeridional_10 gpsMeridional_11 gpsMeridional_12 gpsMeridional_13 gpsMeridional_14 gpsMeridional_15 gpsMeridional_16 gpsMeridional_17]';

gpsWS1 = gpsWS1(:,(1:44));gpsU1 = gpsU1(:,(1:44));gpsV1 = gpsV1(:,(1:44));

gpsWS_Avg = nanmean(gpsWS1,1)';

gpsU_Avg = nanmean(gpsU1,1)';

gpsV_Avg = nanmean(gpsV1,1)';

mYWSarr1_t_Avg= nanmean(mYWSarr1_t,1)';

mYUarr1_t_Avg= nanmean(mYUarr1_t,1)';

mYVarr1_t_Avg= nanmean(mYVarr1_t,1)';

mYWSarr1_t_Avg_std = 1*nanstd(mYWSarr1_t,1)';

mYUarr1_t_Avg_std = 1*nanstd(mYUarr1_t,1)';

mYVarr1_t_Avg_std = 1*nanstd(mYVarr1_t,1)';

gpsWS_Avg_std = 1*nanstd(gpsWS1,1)';

gpsU_Avg_std = 1*nanstd(gpsU1,1)';

gpsV_Avg_std = 1*nanstd(gpsV1,1)';

subplot(3,3,4);

% xlim([0 25]);

arrr= mYWSarr1_t_Avg; barr= ht1; x = [arrr barr];

arrr= gpsWS_Avg; barr= ht1; y = [arrr barr];

out1=abs(mYWSarr1_t_Avg-gpsWS_Avg);

datas=[mYWSarr1_t_Avg ;gpsWS_Avg;out1;mYWSarr1_t_Avg_std;gpsWS_Avg_std];

AvgGain = mean(out1);

[MaxVal MaxInd]=max(mYWSarr1_t_Avg);x = [MaxVal ht1(MaxInd)];

[MaxVal MaxInd]=max(gpsWS_Avg);y = [MaxVal ht1(MaxInd)];

MaxGain = sqrt(((y(1,1)-x(1,1))^2)+((y(1,2)-x(1,2))^2));

[MinVal MinInd]=min(mYWSarr1_t_Avg);x = [MinVal ht1(MinInd)];

[MinVal MinInd]=min(gpsWS_Avg);y = [MinVal ht1(MinInd)];

MinGain = sqrt(((y(1,1)-x(1,1))^2)+((y(1,2)-x(1,2))^2));

hold on;

errorbar(gpsWS_Avg, ht1, gpsWS_Avg_std, 'horizontal', 'color', '#64F8FF', 'LineWidth', .01);

plot(gpsWS_Avg,ht1,'color','#01969B', 'LineWidth', 1,"LineStyle","-");

errorbar(mYWSarr1_t_Avg, ht1, mYWSarr1_t_Avg_std, 'horizontal', 'color', '#9696FA', 'LineWidth', .01);

plot(mYWSarr1_t_Avg,ht1,'color','#1E28FA', 'LineWidth', 1,"LineStyle","-");

[maxv maxi] = max(out1);

xxx=out1';yyy=ht1';

mYVal=AvgGain;

xind = find(xxx<=mYVal,1);

xxx = [xxx(1,1:xind-1 ) mYVal (xxx(1,xind:end))];

yyy = [yyy(1,1:xind-1 ) nan (yyy(1, xind: end))];

adash=[xxx(1,xind-1) yyy(1,xind-1)]'; bdash=[mYVal nan]';cdash =[xxx(1,xind+1) yyy(1,xind+1)]';

p=((cdash(1)-bdash(1))/(bdash(1)-adash(1)));

bdash(2)= (cdash(2) + (adash(2)*p))/(1+p);

yyy(1,xind) = bdash(2);

plot(out1,ht1,'Color','#D95319','LineWidth', 1,"LineStyle","-");

plot(zeros(size(gpsWS_Avg)),ht,'-.','Color','#2B7B2B', 'LineWidth', 1);

vals = [mYWSarr1_t_Avg; gpsWS_Avg;out1];mystds=[mYWSarr1_t_Avg_std;gpsWS_Avg_std;zeros(44,1)];

myacts1 = vals-mystds;

myacts2 = vals+mystds;

datas=[adash ;bdash;out1;xxx';vals((1:44),1);vals((45:88),1);];

plot(-10,(min(ht1)-1),'o','MarkerSize',0.1,'MarkerFaceColor', 'white');

plot((30),(max(ht1)+1),'o','MarkerSize',0.1,'MarkerFaceColor', 'white');

myTitle = strcat('(a) Horizontal Wind');

t= title(myTitle,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'bold','HorizontalAlignment',"center");

x = xlabel('Wind Speed (m/s)','FontSize',11,'FontName','Times New Roman', 'FontWeight', 'Normal');

y = ylabel('Range (km)','FontSize',11,'FontName','Times New Roman', 'FontWeight', 'Normal');

y.FontSize = 12;

set(gca, 'YTick',linspace( min(ht1), max(ht1), 7));

set(gca,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'Normal');

str = strcat('Avg Diff: ',num2str(AvgGain,4),'m/s');

ytickformat('%.1f');

set(gca, 'XTick',linspace(-10,30,5));

set(gca,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'Normal');

xtickformat('%.0f');

grid minor; grid on;

axis tight;

box on;

WindspeedMaxGain = MaxGain

WindspeedMaxGain = 4.0445

WindspeedAvgGain = AvgGain

WindspeedAvgGain = 1.7674

subplot(3,3,5);

arrr= mYUarr1_t_Avg; barr= ht1; x = [arrr barr];

arrr= gpsU_Avg; barr= ht1; y = [arrr barr];

out1=abs(mYUarr1_t_Avg-gpsU_Avg);

datas=[mYUarr1_t_Avg ;gpsU_Avg;out1;mYUarr1_t_Avg_std;gpsU_Avg_std];

AvgGain = mean(out1);

[MaxVal MaxInd]=max(mYUarr1_t_Avg);x = [MaxVal ht1(MaxInd)];

[MaxVal MaxInd]=max(gpsU_Avg);y = [MaxVal ht1(MaxInd)];

MaxGain = sqrt(((y(1,1)-x(1,1))^2)+((y(1,2)-x(1,2))^2));

[MinVal MinInd]=min(mYUarr1_t_Avg);x = [MinVal ht1(MinInd)];

[MinVal MinInd]=min(gpsU_Avg);y = [MinVal ht1(MinInd)];

MinGain = sqrt(((y(1,1)-x(1,1))^2)+((y(1,2)-x(1,2))^2));

set(gca, 'YTick',linspace( min(ht1), max(ht1), 7));

hold on;

errorbar(gpsU_Avg, ht1, gpsU_Avg_std, 'horizontal', 'color', '#64F8FF', 'LineWidth', .01);

plot(gpsU_Avg,ht1,'color','#01969B', 'LineWidth', 1,"LineStyle","-");

errorbar(mYUarr1_t_Avg, ht1, mYUarr1_t_Avg_std, 'horizontal', 'color', '#9696FA', 'LineWidth', .01);

plot(mYUarr1_t_Avg,ht1,'color','#1E28FA', 'LineWidth', 1,"LineStyle","-");

[maxv maxi] = max(out1);

xxx=out1';yyy=ht1';

mYVal=AvgGain;

xind = find(xxx<=mYVal,1);

xxx = [xxx(1,1:xind-1 ) mYVal (xxx(1,xind:end))];

yyy = [yyy(1,1:xind-1 ) nan (yyy(1, xind: end))];

plot(out1,ht1,'Color','#D95319','LineWidth', 1,"LineStyle","-");

plot(zeros(size(gpsU_Avg)),ht,'-.','Color','#2B7B2B', 'LineWidth', 1);

vals = [mYUarr1_t_Avg; gpsU_Avg;out1];mystds=[mYUarr1_t_Avg_std;gpsU_Avg_std;zeros(44,1)];

myacts1 = vals-mystds;

myacts2 = vals+mystds;

datas=[adash ;bdash;out1;xxx';vals((1:44),1);vals((45:88),1);];

plot(-30,(min(ht1)-1),'o','MarkerSize',0.1,'MarkerFaceColor', 'white');

plot((10),(max(ht1)+1),'o','MarkerSize',0.1,'MarkerFaceColor', 'white');

myTitle = strcat('(b) Zonal Wind');

t= title(myTitle,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'bold','HorizontalAlignment',"center");

x = xlabel('Zonal (m/s)','FontSize',11,'FontName','Times New Roman', 'FontWeight', 'Normal');

ytickformat('%.1f');

set(gca, 'XTick',linspace(-30,10,5));

set(gca, 'YTick',linspace( min(ht1), max(ht1), 7));

set(gca,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'Normal');

xtickformat('%.0f');

grid minor; grid on;

axis tight;

box on;

ZonalMaxGain = MaxGain

ZonalMaxGain = 4.4577

ZonalAvgGain = AvgGain

ZonalAvgGain = 2.4424

subplot(3,3,6);

arrr= mYVarr1_t_Avg; barr= ht1; x = [arrr barr];

arrr= gpsV_Avg; barr= ht1; y = [arrr barr];

out1 = abs(mYVarr1_t_Avg-gpsV_Avg);

datas=[mYVarr1_t_Avg ;gpsV_Avg;out1;mYVarr1_t_Avg_std;gpsV_Avg_std];

AvgGain = mean(out1);

[MaxVal MaxInd]=max(mYVarr1_t_Avg);x = [MaxVal ht1(MaxInd)];

[MaxVal MaxInd]=max(gpsV_Avg);y = [MaxVal ht1(MaxInd)];

MaxGain = sqrt(((y(1,1)-x(1,1))^2)+((y(1,2)-x(1,2))^2));

[MinVal MinInd]=min(mYVarr1_t_Avg);x = [MinVal ht1(MinInd)];

[MinVal MinInd]=min(gpsV_Avg);y = [MinVal ht1(MinInd)];

MinGain = sqrt(((y(1,1)-x(1,1))^2)+((y(1,2)-x(1,2))^2));

%plot(tSNR_t_avg,ht1,'-r', 'LineWidth', 1);

hold on;

errorbar(gpsV_Avg, ht1, gpsV_Avg_std, 'horizontal', 'color', '#64F8FF', 'LineWidth', .01);

plot(gpsV_Avg,ht1,'color','#01969B', 'LineWidth', 1,"LineStyle","-");

errorbar(mYVarr1_t_Avg, ht1, mYVarr1_t_Avg_std, 'horizontal', 'color', '#9696FA', 'LineWidth', .01);

plot(mYVarr1_t_Avg,ht1,'color','#1E28FA', 'LineWidth', 1,"LineStyle","-");

plot(zeros(size(gpsV_Avg)),ht,'-.','Color','green', 'LineWidth', 1);

[maxv maxi] = max(out1);

xxx=out1';yyy=ht1';

mYVal=AvgGain;

plot(out1,ht1,'Color','#D95319','LineWidth', 1,"LineStyle","-");

plot(zeros(size(gpsV_Avg)),ht,'-.','Color','#2B7B2B', 'LineWidth', 1);

xind = find(xxx<=mYVal,1);

xxx = [xxx(1,1:xind-1 ) mYVal (xxx(1,xind:end))];

yyy = [yyy(1,1:xind-1 ) nan (yyy(1, xind: end))];

str = strcat('Avg Diff: ',num2str(AvgGain,4),'m/s');

vals = [mYVarr1_t_Avg; gpsV_Avg;out1];mystds=[mYVarr1_t_Avg_std;gpsV_Avg_std;zeros(44,1)];

myacts1 = vals-mystds;

myacts2 = vals+mystds;

datas=[adash ;bdash;out1;xxx';vals((1:44),1);vals((45:88),1);];

plot(-10,(min(ht1)-1),'o','MarkerSize',0.1,'MarkerFaceColor', 'white');

plot((10),(max(ht1)+1),'o','MarkerSize',0.1,'MarkerFaceColor', 'white');

myTitle = strcat('(c) Meridional Wind');

t= title(myTitle,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'bold','HorizontalAlignment',"center");

x = xlabel('Meridional (m/s)','FontSize',11,'FontName','Times New Roman', 'FontWeight', 'Normal');

set(gca, 'YTick',linspace( min(ht1), max(ht1), 7));ylim([2.0, 30]);

ytickformat('%.1f');

set(gca, 'XTick',linspace(-10,10,5));

set(gca,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'Normal');

xtickformat('%.0f');

grid minor; grid on;

axis tight;

box on;

MeridionalMaxGain = MaxGain

MeridionalMaxGain = 5.0829

MeridionalAvgGain = AvgGain

MeridionalAvgGain = 1.0138

legend ({'GPS','','eBayes', '','','Difference'},"Location","none","Position",[0.487624118647838 0.646250897015525 0.417348135464659 0.0192592596301325],"Orientation","Horizontal","Box","off",'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'normal');