%EFD2 BLATT 6 AUFGABE 1	
function aufgabe1
%zu a)
display('zu a)')
%gegebene Werte und einfache Berechnungen
u_10=5;
Az=0.1;
C_10=1.0e-3;
deg=-40;
windrichtung=[sin(5/8*pi),cos(5/8*pi)]';
rho=999.975;
rho_luft=1.204;

omega=7.29e-5;
f=2*omega*sind(deg);

%windvektor
tau=u_10*(windrichtung(1)+i*windrichtung(2)); % - =??

%Lösung an der Oberfläche
z=0;

display('Geschwindigkeit der Strömung (komplex (u+iv)).')
zeta=(1-i)*(-tau/(sqrt(2*Az*f)*rho))*exp((1+i)*z*(sqrt(f/(2*Az))))

%Geschwindigkeit in u und v
display('Geschwindigkeit in u und v.')
u_o=real(zeta)
v_o=imag(zeta)

%Gesamtbetrag der Geschwindigkeit
V_oberflache=abs(zeta);

%Richtungsvektoren der Strömung
display('Richtung der Strömung')
u_d=u_o/V_oberflache
v_d=v_o/V_oberflache
display('dies entspricht (sin(3pi/8),cos(3pi/8)) also ENE');

%b)
display('zu b)')
display('Massentransport (komplex M_x+iM_y).')
Msc=-i*tau/f

display('Gesamtmassentransport:')
Msc_abs=abs(Msc)

display('Richtung des Massentransports')
Msc_dir=Msc/Msc_abs
display('dies entspricht (cos(3pi/8),sin(3pi/8)) also NEN');

%c
display('zu c)')
display('Tiefe in der die Strömung um 180Grad zur Oberfläche gedreht ist.')
t_180=sqrt(2*Az/abs(f))*pi

%moeglicher Test!
%z=sqrt(2*Az/abs(f))*pi
%zeta=(1-i)*(tau/(sqrt(2*Az*f)*rho))*exp((1+i)*z*(sqrt(f/(2*Az))))
%zeta/norm(zeta)

%d)
display('zu d)')

display('M_s=rho*int(zeta(z))')
% falls A_z(z) gegeben wäre, numerische Integration -> fertig.


end

