returnlevel_map: Create return level maps for Austria

In SpatialModelsZAMG: Spatial Modeling of Extreme Events

Description

this function creates a map of Austria with return levels

Usage

returnlevel_map(covariables, rl, q, sd_or_swe,
                plottitle = NULL, save_name = NULL,
                save_dir = getwd(), printPlot = TRUE)

Arguments

covariables	a named matrix with the covariables. each row corresponds to one location, columns should include at least lon and lat
rl	a vector with the return level for every location
q	the return period – must be a number greater than 1
sd_or_swe	a character string; you can either chose snow depth ("sd"), snow water equivalent ("swe") or the quotient swe/sd ("quot")
plottitle	a character string defining the title of the plot. default is plottitle = "sd_or_swe return level map", depending on sd_or_swe. if no plottitle is wanted, use plottitle = ""
save_name	a character string defining the saving name of the map.
save_dir	a character string defining the directory for the map to be saved. default is the working directory
printPlot	logical value; if TRUE (default), the plot is printed

Value

a map of Austria with the return levels

See Also

returnlevels, returnlevel_plot, cond_returnlevel_map

Examples

# get covariables
lon = get(data(lon.at))
lon = as.vector(t(lon))
lat = get(data(lat.at))
lat = as.vector(t(lat))
alt = get(data(alt.at))
alt = as.vector(t(alt))

x = get(data(sample_grid_data))

mdday    = x$mdday
sd_mmax  = x$mmsd
swe_mmax = x$mmswe

# take only locations of the (Austrian) domain
lon = lon[which(!is.na(mdday))]
lat = lat[which(!is.na(mdday))]
alt = alt[which(!is.na(mdday))]
mdday    = mdday[which(!is.na(mdday))]
sd_mmax  = sd_mmax[which(!is.na(sd_mmax))]
swe_mmax = swe_mmax[which(!is.na(swe_mmax))]

# define matrix 'covariables'
covariables = cbind("lon" = lon, "lat" = lat, "alt" = alt,
                    "mdday" = mdday, "sd_mmax" = sd_mmax,
                    "swe_mmax" = swe_mmax)

# load function output from GEVparameters_from_models
sd_GEVparam  = get(data("sd_GEVparam"))
swe_GEVparam = get(data("swe_GEVparam"))

# calculate returnlevels
q = 100
sd_rl  = returnlevels(GEVparam = sd_GEVparam,  q = q)
swe_rl = returnlevels(GEVparam = swe_GEVparam, q = q)

# create return level maps
returnlevel_map(covariables = covariables,
                rl = sd_rl, q = q,
                sd_or_swe = "sd",
                printPlot = FALSE)
returnlevel_map(covariables = covariables,
                rl = swe_rl, q = q,
                sd_or_swe = "swe",
                printPlot = FALSE)

SpatialModelsZAMG documentation built on Nov. 11, 2019, 3 p.m.