Interpolation/
angle.pro
north stereographic polar projection
Compute angles between grid lines and direction of the North pole
(fom angle.F,v 2.2 in OPA8.2)
angle
Interpolation
angle, fileocemesh, gcosu, gsinu, gcosv, gsinv, gcost, gsint, gcosf, gsinf, IODIRECTORY=scalar string, DOUBLE=1 or 0
Parameters
fileocemesh
in
required
type: scalar string
a netcdf file that contains (at least) the following variables:
glamt, gphit: longitudes and latitudes at T-points
glamu, gphiu: longitudes and latitudes at U-points
glamv, gphiv: longitudes and latitudes at V-points
glamf, gphif: longitudes and latitudes at F-points
gcosu
out
type: 2d array
cosinus of the angle between grid lines at U points and direction of the North pole
gsinu
out
type: 2d array
sinus of the angle between grid lines at U points and direction of the North pole
gcosv
out
type: 2d array
cosinus of the angle between grid lines at V points and direction of the North pole
gsinv
out
type: 2d array
sinus of the angle between grid lines at V points and direction of the North pole
gcost
out
type: 2d array
cosinus of the angle between grid lines at T points and direction of the North pole
gsint
out
type: 2d array
sinus of the angle between grid lines at T points and direction of the North pole
gcosf
out
type: 2d array
cosinus of the angle between grid lines at F points and direction of the North pole
gsinf
out
type: 2d array
sinus of the angle between grid lines at F points and direction of the North pole
Keywords
IODIRECTORY
type: scalar string
default: ''
the directory path where is located fileocemesh
DOUBLE
type: 1 or 0
default: 0
put 1 to use double precision (default is float)
Version history
Version
$Id: angle.pro 327 2007-12-13 16:22:35Z pinsard $
History
Original : 96-07 (O. Marti)
98-06 (G. Madec)
Feb 2005: IDL adaptation S. Masson
May 2007: F points + call to lbcorca
Known issues
Restrictions
to compute the lateral boundary conditions, we assume that:
(1) the first line is similar to the second line
=> gcosu[*, 0] = gcosu[*, 1]
=> gsinu[*, 0] = gsinu[*, 1]
(2) the grid follows OPA x and north pole periodicity rule (see lbcorca)
Produced by IDLdoc 2.0.