% example showing the solution to the plane shear problem using 
% the Jinkins & Richman (1985) theory

clear all; close all;

param.e = 0.7;
param.phibar = 0.3;
param.u_bot = -0.3;
param.u_top = 0.3;
param.w_bot = 0.435;
param.w_top = 0.435;
param.Fr = 0;
param.H = 12;

xdisp = [0:0.01:1]';

sol_JR = planeShear_JR(param, xdisp);

figure;
subplot(311);
plot(xdisp, sol_JR.phi, 'b-');
title('solution to shear flow using J-R theory: <\nu>=0.3, H/\sigma=12');
xlabel('y/H');
ylabel('\nu(y)');

subplot(312);
plot(xdisp, sol_JR.u, 'g-');
title('solution to shear flow using J-R theory: <\nu>=0.3, H/\sigma=12');
xlabel('y/H');
ylabel('u(y)');

subplot(313);
plot(xdisp, sol_JR.w, 'r-');
title('solution to shear flow using J-R theory: <\nu>=0.3, H/\sigma=12');
xlabel('y/H');
ylabel('T^{1/2}(y)');

