gridarea1d: grid cell area along a vector of latitudes

In MagicForrest/DGVMTools: DGVM Processing, Analysis and Plotting Tools

grid cell area along a vector of latitudes

Description

Returns the area in square meters along a vetor of latitudes by equal longitude distance.

Usage

gridarea1d(lat, dlon, scale = 1, ellipse = FALSE)

Arguments

lat	vetor of latitudes
dlon	longitudinal extent
scale	multiplicator. If 1 (default) unit m^2
ellipse	TRUE (polar and equatorial radius differ) or FALSE (default, polar and equatorial radius are the same)

Details

The function returns a vector of area in square meters along a vector of latitudes. These must not be of equal distance. However, for the longitude will be equal along the given latitude vector. The latitude is assumed to be the gridcell midpoint and the northern and southern edges are calculated as half the distance to the next element in the latitude vector.

Value

vector of gridcell area is m^2

Author(s)

Joerg Steinkamp [e-mail]

Examples

dlon <- 0.5
lat <- seq(-89.75,89.75,0.5)
sum(gridarea1d(lat,dlon))*720*1.e-12

MagicForrest/DGVMTools documentation built on Aug. 23, 2024, 8:05 a.m.