function [f, dfdy] = cor(lat);
% Returns the coriolis and beta parameters at the input latitude (in degrees).

f = 2 * 7.292e-5 * sin((pi/180.0)*lat);
dfdy = 2 * 7.292e-5 * cos((pi/180.0)*lat) / 6370000.0;

