From c9e150388a9cedf9cfacf7c74f8860f7b9adfdef Mon Sep 17 00:00:00 2001 From: nperez Date: Fri, 22 Nov 2019 16:53:05 +0100 Subject: [PATCH] accept positive values for west longitudes --- R/WeightedMean.R | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/R/WeightedMean.R b/R/WeightedMean.R index 70c33b4..2ab1a1a 100644 --- a/R/WeightedMean.R +++ b/R/WeightedMean.R @@ -120,7 +120,8 @@ WeightedMean <- function(data, lon, lat, region = NULL, mask = NULL, londim = NU cosphi <- t(array(cos(lat * pi / 180), dim = c(length(lat), length(lon)))) nblat <- length(lat) nblon <- length(lon) - dlon <- abs(c(lon[2 : nblon] - lon[1 : nblon - 1])) * pi / 180 + lon[lon > 180] = lon[lon > 180] - 360 + dlon <- abs(c(abs(lon[2 : nblon]) - abs(lon[1 : nblon - 1]))) * pi / 180 dlon <- c(dlon, dlon[1]) dlon <- array(dlon, dim = c(nblon, nblat)) dlat <- abs(c(lat[2 : nblat] - lat[1 : nblat - 1])) * pi / 180 -- GitLab