% Buoyancy check for SOLOPC
%Assume P,T,S profile already read in

rho= ( 1 +  density(T,S,P)/1000);   %density profile

%SOLO description
m0 = 29857.2 ;%mass(grams) of SOLOPC with dummy LOPC & standard SOLO SBE top cap 
alpha = (2.45e-6); %Pressco of SOLOPC hull
beta = (6.9e-5);  %Tempco of SOLO hull
v0=m0/rho(1);%guessed SOLO volume = SOLO mass estimate/(rho at surface)

vol=v0*(1 - alpha*P + beta*T);  %SOLO volume grams
buoy=vol.*rho-m0;%net buoy of SOLO

