clear;
load APG3/eps300
load APG3/uv300
load APG3/uy300

ui=0.6408;
ut=0.0215;
re=400;

ypl=uv300(:,1);
y=ypl./ut./400;

%plot production and dissipation in simulation coordinates
eps=eps300(:,2).*ut^4.*400;
load APG3/uv300
load APG3/uy300
pr=-uv300(:,2).*uy300(:,2).*ut^4.*400;
plot(y,eps,'--',y,pr,'-')
axis([0 65 0.0 0.0001 ])
zoom on
