Skip to content

Instantly share code, notes, and snippets.

@ibarraespinosa

ibarraespinosa/readwrf.R

Last active August 13, 2020 13:43
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save ibarraespinosa/fec92f4cea0194572824e4bbe76a67c3 to your computer and use it in GitHub Desktop.
Save ibarraespinosa/fec92f4cea0194572824e4bbe76a67c3 to your computer and use it in GitHub Desktop.
Download ZIP
Raw
readwrf.R
#!/usr/bin/env Rscript
args = commandArgs(trailingOnly=TRUE)
library(wrftools)
data(cetesb)
cetesb <- cetesb[!is.na(cetesb$Station), "Station"]
t2 = c("T2", "o3", "co", "no", "no2", "U10", "V10")
df <- xtractor(atmos = args[1],
vars = t2,
points = cetesb,
stations = cetesb$Station,
return_list = TRUE)
saveRDS(df, args[2], compress = "xz")
file.copy(args[2], paste0("../../../git_anfavea_2/rds/", args[2]), overwrite = T)
@ibarraespinosa
Copy link
Author

Rscript rwrfpetro.R ../em_real_aer_2016/wrfout_d01_2016-04-06_00:00:00 "2016.rds"

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment