gridInfo: Read grid information from a NetCDF file

Description Usage Arguments Note Examples

View source: R/gridinfo.R

Description

Return a list containing information of a regular grid / domain

Usage

1
gridInfo(file = file.choose(), z = F, geo = F, verbose = T)

Arguments

file

file name/path to a wrfinput, wrfchemi or geog_em file

z

TRUE for read wrfinput vertical coordinades

geo

True for use geog_em files

verbose

display additional information

Note

just WRF-Chem is suported by now

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
grid_d1 <- gridInfo(paste(system.file("extdata", package = "EmissV"),"/wrfinput_d01",sep=""))

grid_d2 <- gridInfo(paste(system.file("extdata", package = "EmissV"),"/wrfinput_d02",sep=""))
grid_d3 <- gridInfo(paste(system.file("extdata", package = "EmissV"),"/wrfinput_d03",sep=""))
names(grid_d1)
# for plot the shapes
library(sp)
shape   <- raster::shapefile(paste0(system.file("extdata", package = "EmissV"),"/BR.shp"))
plot(shape,xlim = c(-55,-40),ylim = c(-30,-15), main="3 nested domains")
axis(1); axis(2); box(); grid()
lines(grid_d1$Box, col = "red")
text(grid_d1$xlim[2],grid_d1$Ylim[1],"d1",pos=4, offset = 0.5)
lines(grid_d2$Box, col = "red")
text(grid_d2$xlim[2],grid_d2$Ylim[1],"d2",pos=4, offset = 0.5)
lines(grid_d3$Box, col = "red")
text(grid_d3$xlim[1],grid_d3$Ylim[2],"d3",pos=2, offset = 0.0)

EmissV documentation built on May 2, 2019, 7:16 a.m.