#!/usr/bin/env Rscript
args = commandArgs(trailingOnly=TRUE)
# how to use
# Rscript rwrf1.R wrfout_d01_RADM2_BASE
# Edit and change variables
library(raster)
library(wrftools)
data(cetesb)
cetesb <- cetesb[!is.na(cetesb$Station), "Station"]
t2 = c("T2", "o3", "co", "no", "no2", "U10", "V10", "PSFC", "PM10", "PM2_5_DRY")
# time series ####
wrfi <- args[1]
newdir <- substr(wrfi, start = 12, stop = 100)
newdir <- gsub(":", "_", newdir)
newdir <- gsub(".nc", "", newdir)
dir.create(newdir)
df <- xtractor(atmos = wrfi,
vars = t2,
points = cetesb,
# stations = cetesb$Station,
return_list = TRUE)
df$data$dir <- newdir
saveRDS(df$data, paste0(newdir, "/DF.rds"), compress = "xz")
for (i in seq_along(t2)) {
print(names(df$raster)[i])
writeRaster(
df$raster[[i]],
paste0(newdir, "/", t2[i], ".nc"),
overwrite = TRUE
)
print(
paste0(newdir, "/", t2[i], ".nc")
)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.