function [th U] = IncFlow(t)

U = .2;

th = (20*pi/180)*sin(2*pi*t/20);

%if(t<20)
%    th = 0;
%else
%    th = 20*pi/180;
%end