Interpolation/
angle.pro
north stereographic polar projection
Compute angles between grid lines and direction of the North
(fom angle.F,v 2.2 in OPA8.2)
Routine summary
- result = fsnspp(plam, pphi, DOUBLE=DOUBLE)
-
- angle, fileocemesh, gcosu, gsinu, gcosv, gsinv, gcost, gsint, IODIRECTORY=IODIRECTORY, DOUBLE=DOUBLE
-
fsnspp
Interpolation
result = fsnspp(plam, pphi, DOUBLE=DOUBLE)
Return value
gsinu,gcosu : sinus and cosinus of the angle
gsinv,gcosv between north-south direction
gsint,gcost and the j-direction of the mesh
Parameters
plam
in
required
pphi
in
required
Keywords
DOUBLE
default: 0
use double precision (default is float)
Version history
Version
$Id: angle.pro 157 2006-08-21 09:01:50Z navarro $
History
Original : 96-07 (O. Marti)
98-06 (G. Madec)
Feb 2005: IDL adaptation S. Masson
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 periodicity rule, first column is
equal to the next to last column
=> gcosv[0, *] = gcosv[jpj-2, *]
=> gsinv[0, *] = gsinv[jpj-2, *]
angle
angle, fileocemesh, gcosu, gsinu, gcosv, gsinv, gcost, gsint, IODIRECTORY=IODIRECTORY, DOUBLE=DOUBLE
Parameters
fileocemesh
in
required
a netcdf file that contains (at least):
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
in
required
gsinu
in
required
gcosv
in
required
gsinv
in
required
gcost
in
required
gsint
in
required
Keywords
IODIRECTORY
the directory path where is located fileocemesh
DOUBLE
default: 0
use double precision (default is float)
Produced by IDLdoc 2.0 on Tue Aug 29 14:49:14 2006.