%bellowsFillVeri.m

[tDat,~,cyl,bel,~,~,~] = efficiencyBEwrapper('Y:\Data\TestTankCPF2013Aug12T111737.xls');
set(0,'DefaultAxesFontSize',14)

%constants
cylDia = 2; %in
cylRad = .0254*cylDia/2; %meters
mtsConv = 34.798; %mm


%convert cyl to vol
cVol = (cyl/1000)*(pi*cylRad^2); %cubic meters

%convert bellows to vol
bVol = zeros(size(bel));
for bi = 1:length(bel)
    bVol(bi) = bellowsFill((bel(bi)-mtsConv)/1000);
end

figure('color','w');
plot(tDat,cVol-cVol(1),'k',tDat,bVol-bVol(1),'r:')
xlabel('Elapsed Time (sec)')
ylabel('Volume (m^3)')
legend('Measured Piston Volume','Predicted Bellows Volume')
