function rake = pitch(plunge, delta)
%Simple determination of pitch (rake) given plunge and dip
plungerads=plunge.*pi./180;
deltarads=delta.*pi./180;
rakerads=asin(sin(plungerads)./sin(deltarads));
rake=rakerads.*180./pi;


