In energyRt/merra2ools: R-tools to process MERRA-2 data (NASA) with focus on renewable energy

options(rmarkdown.html_vignette.check_title = FALSE)
library(knitr)
knitr::opts_chunk$set(
  collapse = TRUE, 
  out.height = "75%", out.width = "75%", fig.align = "center",
  comment = "#>"
)

library(merra2ools)
library(tidyverse)
# data_month <- 1
data_month <- 1:12
# get example of the data for all locations
x <- merra2_sample(data_month, add.coord = T)
x

Included variables

gif_file <- "images/merra_576x540_10fps.gif"
# if (file.exists(gif_file)) 
include_graphics(gif_file)

# gif_file <- "images/POA_map_720x544_5fps.gif"
# if (file.exists(gif_file)) include_graphics(gif_file)

Temperature at 10 meters hight (T10M), Kelvin degrees

Figure code

summary(x$T10M)
summary(x$T10M - 272) # Kelvin -> Celsius

plot_merra(x, name = "T10M", 
           scale = 1, intercept = -272, # Kelvin -> Celsius
           limits = c(-50, 50), 
           legend.name = "\u00B0C",
           datetime.position = c(144, 87))

gif_merra(x, name = "T10M", dirname = "images", fps = 10,
          limits = c(-50, 50), intercept = -272, 
          legend.position = c(0.9, 0.05), 
          legend.name = expression("\u00B0C"),
          datetime.position = c(144, 87))

gif_file <- "images/merra_T10M_576x360_10fps.gif"
if (file.exists(gif_file)) include_graphics(gif_file)

Wind speed at 10 and 50 meters hights (W10M and W50M), m/s

Figure code

summary(x$W50M)
# fig.W50M <- 
plot_merra(x, name = "W50M", limits = c(0, 30), 
           legend.name = expression("m/s"),
           datetime.position = c(144, 87))

# fig.W10M <- 
plot_merra(x, name = "W10M", limits = c(0, 30),
           legend.name = expression("m/s"),
           datetime.position = c(144, 87))

gif_merra(x, name = "W50M", dirname = "images", fps = 12,
          limits = c(0, 30), 
          legend.position = c(0.9, 0.05), 
          legend.name = expression("m/s"),
          datetime.position = c(144, 87))


gif_merra(x, name = "W10M", dirname = "images", fps = 12,
          limits = c(0, 30), 
          legend.position = c(0.9, 0.05), 
          legend.name = expression("m/s"),
          datetime.position = c(144, 87))

gif_file <- "images/merra_W10M_576x360_12fps.gif"
if (file.exists(gif_file)) include_graphics(gif_file)

gif_file <- "images/merra_W50M_576x360_12fps.gif"
if (file.exists(gif_file)) include_graphics(gif_file)

Global Horizontal Irradiance (SWGDN), Watt/m2

Figure code

summary(x$SWGDN)

plot_merra(x, name = "SWGDN", limits = c(0, 1200), 
           legend.position = c(0.88, 0.05), 
           legend.name = ("W/\u33A1"),
           datetime.position = c(144, 87))

plot_merra(x, name = "SWGDN", dirname = "images", fps = 5,
          limits = c(0, 1), scale = 1/1000,
          legend.position = c(0.88, 0.05), 
          legend.name = ("kW/\u33A1"),
          datetime.position = c(144, 87))

gif_merra(x, name = "SWGDN", dirname = "images", fps = 5,
          limits = c(0, 1200), scale = 1,
          legend.position = c(0.88, 0.05), 
          legend.name = ("W/\u33A1"),
          datetime.position = c(144, 87))

gif_file <- "images/merra_SWGDN_576x360_5fps.gif"
if (file.exists(gif_file)) include_graphics(gif_file)

Albedo (ALBEDO), index [0..1]

Figure code

summary(x$ALBEDO)

plot_merra(x, name = "ALBEDO", dirname = "images", fps = 5,
          limits = c(0, 1), scale = 1,
          legend.position = c(0.9, 0.05), 
          legend.name = (""),
          datetime.position = c(144, 87))

gif_merra(x, name = "ALBEDO", dirname = "images", fps = 5,
          limits = c(0, 1), scale = 1,
          legend.position = c(0.9, 0.05), 
          legend.name = (""),
          datetime.position = c(144, 87))

gif_file <- "images/merra_ALBEDO_576x360_5fps.gif"
if (file.exists(gif_file)) include_graphics(gif_file)

Surface atmospheric preassure (PS), hPa

Figure code

summary(x$PS_hPa)

plot_merra(x, name = "PS_hPa", dirname = "images", fps = 5,
          limits = c(500, 1000), scale = 1, direction = 1,
          legend.position = c(0.89, 0.05), palette = "YlGnBu",
          legend.name = ("hPa"),
          datetime.position = c(144, 87))

gif_merra(x, name = "PS_hPa", dirname = "images", fps = 24,
          limits = c(500, 1000), scale = 1, direction = 1,
          legend.position = c(0.89, 0.05), 
          legend.name = ("hPa"),
          datetime.position = c(144, 87))

gif_file <- "images/merra_PS_hPa_576x360_24fps.gif"
if (file.exists(gif_file)) include_graphics(gif_file)

Figure code

PS_mean <- x %>% 
  group_by(locid) %>%
  summarise(PS_mean = mean(PS_hPa, na.rm = TRUE), .groups = "drop")

x <- full_join(x, PS_mean) %>%
  mutate(dPS = PS_hPa - PS_mean)
summary(x$dPS)

plot_merra(x, name = "dPS", dirname = "images", fps = 5,
          limits = c(-30, 30), scale = 1, direction = -1,
          legend.position = c(0.9, 0.05), 
          legend.name = ("hPa"),
          datetime.position = c(144, 87))

gif_merra(x, name = "dPS", dirname = "images", fps = 5,
          limits = c(-30, 30), scale = 1, direction = -1,
          legend.position = c(0.89, 0.05), 
          legend.name = ("hPa"),
          datetime.position = c(144, 87))

gif_file <- "images/merra_dPS_576x360_10fps.gif"
if (file.exists(gif_file)) include_graphics(gif_file)

Combined figures code

fig.T10M <- function(x) {
  plot_merra(x, name = "T10M", intercept = -272, 
             limits = c(-50, 50), legend.name = "\u00B0C",
           # datetime.position = c(144, 87), 
           # legend.position = c(0.89, 0.05), 
           legend.position = NULL,
           datetime.position = NA,
           expand.x = c(0.002, 0.005), expand.y = c(0.002, 0.005))
}
# fig.T10M(x)

fig.W50M <- function(x) {
    plot_merra(x, name = "W50M", limits = c(0, 30), 
           legend.name = expression("m/s"),
           datetime.position = c(110, 87),
           # legend.position = c(0.89, 0.05), 
           legend.position = NULL,
           # datetime.position = NA,
           expand.x = c(0.002, 0.005), expand.y = c(0.002, 0.005)
           )
}
# fig.W50M(x)

fig.W10M <- function(x) {
    plot_merra(x, name = "W10M", limits = c(0, 30), 
           legend.name = expression("m/s"),
           # datetime.position = c(144, 87)
           # legend.position = c(0.89, 0.05), 
           legend.position = NULL,
           datetime.position = NA,
           expand.x = c(0.002, 0.005), expand.y = c(0.002, 0.005)
           )
}
# fig.W10M(x)

fig.ALBEDO <- function(x) {
  plot_merra(x, name = "ALBEDO", dirname = "images", fps = 5,
          limits = c(0, 1), scale = 1,
          # legend.position = c(0.9, 0.05), 
          legend.position = NULL,
          legend.name = (""), 
           datetime.position = NA,
           expand.x = c(0.002, 0.005), expand.y = c(0.002, 0.005)
           )
}
# fig.ALBEDO(x)

fig.SWGDN <- function(x) {
  plot_merra(x, name = "SWGDN", limits = c(0, 1200), 
           # legend.position = c(0.88, 0.05), 
           legend.position = NULL,
           legend.name = ("W/\u33A1"), 
           datetime.position = NULL,
           expand.x = c(0.002, 0.005), expand.y = c(0.002, 0.005)
           )
}
# fig.SWGDN(x)

fig.PS_hPa <- function(x) {
  plot_merra(x, name = "PS_hPa", dirname = "images", fps = 5,
          limits = c(500, 1000), scale = 1, direction = 1,
          # palette = "YlGnBu",
          # legend.position = c(0.89, 0.05), 
          legend.position = NULL,
          legend.name = ("hPa"),
           datetime.position = NA,
           expand.x = c(0.002, 0.005), expand.y = c(0.002, 0.005)
           )
}
# fig.PS_hPa(x)

fig.dPS <- function(x) {
  plot_merra(x, name = "dPS", dirname = "images", fps = 5,
          limits = c(-50, 50), scale = 1, direction = -1,
          # legend.position = c(0.9, 0.05), 
          legend.position = NULL,
          legend.name = ("hPa"),
           datetime.position = NA,
           expand.x = c(0.002, 0.005), expand.y = c(0.002, 0.005)
           )
}
# fig.dPS(x)

library("cowplot")
fig.grid <- function(x) {
  plot_grid(fig.W10M(x), fig.W50M(x),
            fig.SWGDN(x), fig.ALBEDO(x),
            fig.T10M(x), fig.PS_hPa(x),
            # fig.dPS(x),
            labels = c("W10M", "W50M", "SWGDN", "ALBEDO", "T10M", "PS"),
            ncol = 2, nrow = 3, hjust = -0.1, vjust = 1.1, label_size = 10)
}
# ii <- lubridate::month(x$UTC) == 2; summary(ii)
fig.grid(x)

gif_merra(x, FUN = "fig.grid", dirname = "images", fps = 10,
          gif.width = 576, gif.height = 360 * 1.5)


# gif_merra(x, FUN = "fig.grid", dirname = "images", fps = 10,
#           gif.width = 576 * 2, gif.height = 360 * 3)


# gif_merra(x, FUN = "fig.grid", dirname = "images", fps = 10,
#           gif.width = 576/1.5, gif.height = 360)


fig.grid2 <- function(x) {
  plot_grid(fig.W10M(x), fig.W50M(x),
            fig.SWGDN(x), fig.ALBEDO(x),
            fig.T10M(x), fig.dPS(x),
            # fig.dPS(x),
            labels = c("W10M", "W50M", "SWGDN", "ALBEDO", "T10M", "dPS"),
            ncol = 2, nrow = 3, hjust = -0.1, vjust = 1.1, label_size = 10)
}

gif_merra(x, FUN = "fig.grid2", dirname = "images", fps = 5,
          gif.width = 576, gif.height = 360 * 1.5)

energyRt/merra2ools documentation built on May 2, 2024, 4:53 a.m.