

function viscosity = visc(viscair, oil, visc2, visc1, T2, T1, Tactual)
%visc (****) calculates the actual viscosity of the thruster unit
%fill fluid, using theactual temperature

if oil == 1
    %calculate oil viscosity at a specified temperature
    viscosity = ((visc2-visc1)/(log(T2)-log(T1)))* (log(Tactual)-log(T1)) +visc1;
else
    viscosity = viscair
end







