ind=find((~isnan(profile.Vh));
[locs pks]=peakseek(diff(ind),1,2);
hfig=figure('Position',[100 100 1000 500]);
haxes1=axes('Position',[0.1 0.1 .25 .9],'Parent',hfig);
hplt1=plot(profile.Vh,profile.heights,'Linewidth',.1);
hold on
hplt2=plot(profile.Vh-profile.dVh',profile.heights,':r');
hplt3=plot(profile.Vh+profile.dVh',profile.heights,':r');
haxes2=axes('Position',[0.40 0.1 .25 .9],'Parent',hfig);
hplt4=scatter(profile.theta,profile.heights,'MarkerSize',8);
