Computation/
curl.pro
Calculate the vertical component of the curl of a vectors field
located on Arakawa C-grid.
curl
Calculation
result = curl(uu, vv, DIREC=scalar string, MILLION=scalar: 0 or 1, _EXTRA=_EXTRA)
Return value
the vertical component of the curl of the input data (with the same
size) located at F point (see restrictions)
Parameters
uu
Matrix representing the zonal coordinates (at U point) of a field of vectors
A 2D (xy), 3D (xyz or yt) or a structure readable by litchamp and containing
a 2D (xy), 3D (xyz or yt) array (4D case is not coded yet).
Note that the dimension of the array must suit the domain dimension.
vv
Matrix representing the meridional coordinates (at V point) of a field of vectors
A 2D (xy), 3D (xyz or yt) or a structure readable by litchamp and containing
a 2D (xy), 3D (xyz or yt) array (4D case is not coded yet).
Note that the dimension of the array must suit the domain dimension.
Keywords
DIREC
type: scalar string
Use if you want to call moyenne or
grossemoyenne after the div computation
with a mean done in the DIREC direction
MILLION
type: scalar: 0 or 1
default: 0
Activate to multiply returned array by 1.e6
_EXTRA
Used to declare that this routine accepts inherited keywords
Examples
Version history
Version
$Id: curl.pro 371 2008-08-07 09:32:02Z pinsard $
History
Guillaume Roullet (grlod@ipsl.jussieu.fr)
Sebastien Masson (smasson@lodyc.jussieu.fr)
adaptation to work with a reduce domain
21/5/1999: missing values at !values.f_nan
Known issues
Todo items
code the 4D case
Restrictions
- Works only for Arakawa C-grid.
- UU must be on U grid, VV must be on V grid
- 4D case is not coded yet
- the common variable jpt is used to differ xyz (jpt=1) and xyt (jpt=1) cases.
- U and V arrays are cut in the same geographic domain. Because of the shift between
T, U, V and F grids, it is possible that these two arrays do not have the same
size and refer to different indexes. In this case, arrays are re-cut on
common indexes. To avoid these re-cuts, use the keyword /memeindice in
domdef
- When computing the divergence, we update, vargrid, varname, varunits and the
grid position parameters (firstxf, lastxf, nxf, firstyf, lastyf, nyf).
- points that cannot be computed (domain boundaries, coastline) are set to NaN
Other attributes
Uses routines
cm_4cal
cm_4data
cm_4mesh
Produced by IDLdoc 2.0.