Simple plot from wrf emission file

Source: R/wrf_plot.R

Create a quick plot from wrf emission file

wrf_plot(
  file = file.choose(),
  name = NA,
  time = 1,
  nivel = 1,
  barra = T,
  lbarra = 0.2,
  col = cptcity::cpt(n = 20, rev = T),
  map = NULL,
  skip = FALSE,
  verbose = TRUE,
  ...
)

Arguments

file

emission file name
name

pollutant name
time

time from emission file
nivel

level from the emission file
barra

barblot if TRUE
lbarra

length of barplot
col

color vector
map

function call to plot map lines, points and annotation (experimental)
skip

logical, skip plot of constant valuess
verbose

if TRUE print some information
...

Arguments to be passed to plot methods

Note

If the file contains levels (kemit>1), and one frame (auxinput5_interval_m = 1) time with control the level which will be ploted

In case of an error related to plot.new() margins lbarra must be adjusted

See also

Lights, to_wrf and wrf_create

Author

Daniel Schuch

Examples

{

dir.create(file.path(tempdir(), "EMISS"))
wrf_create(wrfinput_dir = system.file("extdata", package = "eixport"),
          wrfchemi_dir = file.path(tempdir(), "EMISS"))

# get the name of created file
files <- list.files(path = file.path(tempdir(), "EMISS"),
                   pattern = "wrfchemi",
                   full.names = TRUE)

# load end write some data in this emission file
data(Lights)
to_wrf(Lights, files[1], total = 1521983, name = "E_CO")

wrf_plot(files[1], "E_CO")
}

#> Warning: '/tmp/Rtmp0unmTl/EMISS' already exists
#> creating emission for domain 1 ...
#> output file: /tmp/Rtmp0unmTl/EMISS/wrfchemi_d01_2011-08-01_00:00:00 
#> writing emissions: E_CO weight 1 on file /tmp/Rtmp0unmTl/EMISS/wrfchemi_d01_2011-08-01_00:00:00 
#> /tmp/Rtmp0unmTl/EMISS/wrfchemi_d01_2011-08-01_00:00:00
#> E_CO:
#> Max value: 2481.08618164062, Min value: 0