rutils/rwrf1.R

In ibarraespinosa/eixport: Export Emissions to Atmospheric Models

#!/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")
  )
   }