Plot a tower footprint

Description

Plots the output of a Kormann and Meixner footprint model.

Usage

1
2
plotFootprint(footprint, main.title = NULL, key.title =
NULL, color.palette = NULL, ...)

Arguments

footprint

The footprint output produced by either Calculate or Average function

main.title

The title of the plot. Default is "Footprint"

key.title

The title over the key bar. Default is "Probability"

color.palette

The color palette to use.

...

Further graphical arguments for the filled.contour() function

Author(s)

Georgios Xenakis

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
## Close any previously open graphic devices
graphics.off()

## Load data
data(fluxes)

## Clean and de-spike the fluxes
fluxes=cleanFluxes(fluxes,sdCor=TRUE,sdTimes=3,timesList=3,distCor=TRUE,
thresholdList=list(H=c(-100,1000),LE=c(-100,1000)))

## Use only non-NA data
fluxes_noNA<-fluxes[which(!is.na(fluxes$L)),]

## Displacement height
d=17.42 ## Displacement height

## Calculate input variables for the footprint model
fluxes_noNA$z=(fluxes_noNA$X.z.d..L*fluxes_noNA$L)+d
fluxes_noNA$zol=fluxes_noNA$z/fluxes_noNA$L

## Calculate the footprint for day 853
footprint=Calculate(fetch=500,height=33,grid=200,fluxes_noNA$wind_speed[853],
fluxes_noNA$wind_dir[853],fluxes_noNA$u.[853],fluxes_noNA$zol[853],
sqrt(fluxes_noNA$u_var[853]))

## Plot the footprint
plotFootprint(footprint)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.