plotFootprint: Plot a tower footprint

In FREddyPro: Post-Processing EddyPro Full Output File

Description

Plots the output of a Kormann and Meixner footprint model.

Usage

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

## 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)

Example output

FREddyPro documentation built on May 29, 2017, 7:22 p.m.