%% Micromanometer velocity computations
level = [3.39,3.5,4.15,5.1,5.48]*1e-1;
deltaHmano = level -level(1)
deltPmano = 2*998*981*2.54/100*deltaHmano
velMano = sqrt(2*deltPmano/1.19)*100/2.54/12/10
%% dpt velocity comp
volts = [1.5,1.72,2.18,2.88,3.21];
deltaHdpt = (volts-volts(1))/4
delPdpt = 998*9.81*deltaHdpt*2.54/100
velDPT = sqrt(2*delPdpt/1.19)*100/2.54/12
%% comparisons
velWindWright = [0 15 25 35 45]
fitEqn = 'm*x';
fitOps = fitoptions('Method','NonlinearLeastSquares',...
     'Lower',[0],...
    'Upper',[inf],...
    'StartPoint',[1/2]);
fitType = fittype(fitEqn,'options',fitOps)
[fitInfoMano,GoFMano]=fit(velWindWright',velMano',fitType)
figure(9)
plot(fitInfoMano,velWindWright,velMano)
% manoAccuracy = (fitInfoMano*m-1/1*100)
%% comparisons
velWindWright = [0 15 25 35 45]
fitEqn = 'm*x';
fitOps = fitoptions('Method','NonlinearLeastSquares',...
     'Lower',[0],...
    'Upper',[inf],...
    'StartPoint',[1/2]);
fitType = fittype(fitEqn,'options',fitOps)
[fitInfoDPT,GoFDPT]=fit(velWindWright',velDPT',fitType)
figure(9)
plot(fitInfoDPT,velWindWright,velDPT)
% manoAccuracy = (fitInfoMano*m-1/1*100)