Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
function f = sw_f(lat)
% SW_F Coriolis factor "f"
%===========================================================================
% SW_F $Revision: 1.3 $ $Date: 1994/10/10 04:57:08 $
% Copyright (C) CSIRO, Phil Morgan 1993.
%
% USAGE: f = sw_f(lat)
%
% DESCRIPTION:
% Calculates the Coriolis factor "f" defined by
% f = 2*Omega*Sin(lat) where Omega = 7.292e-5 radians/sec
%
% INPUT:
% lat = Latitude in decimal degress north [-90..+90]
%
% OUTPUT:
% f = Coriolis Factor "f" [s-1]
%
% AUTHOR: Phil Morgan 93-04-20 (morgan@ml.csiro.au)
%
% DISCLAIMER:
% This software is provided "as is" without warranty of any kind.
% See the file sw_copy.m for conditions of use and licence.
%
% REFERENCE:
% S. Pond & G.Pickard 2nd Edition 1986
% Introductory Dynamical Oceanogrpahy
% Pergamon Press Sydney. ISBN 0-08-028728-X
%
% A.E. Gill 1982. p.597
% "Atmosphere-Ocean Dynamics"
% Academic Press: New York. ISBN: 0-12-283522-0
%=========================================================================
% CALLER: general purpose
% CALLEE: none
%-------------
% CHECK INPUTS
%-------------
if nargin ~= 1
error('sw_f.m: Requires one input argument')
end %if
%-------------
% BEGIN
%-------------
% Eqn p27. Unesco 1983.
DEG2RAD = pi/180;
OMEGA = 7.292e-5; %s-1 A.E.Gill p.597
f = 2*OMEGA*sin(lat*DEG2RAD);
return
%===========================================================================