plotZonalMeans: Plot zonal mean values of model and reference data in a...

In amber: Automated Model Benchmarking R Package

Description

This function plots zonal mean values of model and reference data, computed from zonalMeanStats, for a selection of variables.

Usage

plotZonalMeans(mod.path.list, modelIDs, myVariables, metric = "mean",
  lat.range = c(-50, 80), plot.width = 6, plot.height = 10,
  outputDir = FALSE, myFilename = "zonalMeans.pdf",
  legendLocation = "topright", plotRef = TRUE, subcaption = "")

Arguments

mod.path.list	A List of directories where AMBER output is stored for different model runs, e.g. list(mod01.path, mod02.path, mod03.path)
modelIDs	An R object with the different model run IDs, e.g. c('CLASSIC.CRUJRAv2', 'CLASSIC.GSWP3W5E5', 'CLASSIC.CRUNCEP')
myVariables	An R object with the variable names of interest, e.g. c('GPP.FluxCom', 'RECO.FluxCom').
metric	Specify for what statistical metric you wish to plot zonal means. Current options are 'mean', 'bias', 'crmse', 'phase', 'iav', 'sd', 'bias.score', 'rmse.score', 'phase.score', and 'iav.score'.
lat.range	Latitudinal range of ticks and labels on the horizontal axis, e.g. c(-50, 80).
plot.width	A number that gives the desired plot width, e.g. 6.
plot.height	A number that gives the desired plot height, e.g. 10.
outputDir	A string that gives the output directory, e.g. '/home/project/study'. The output will only be written if the user specifies an output directory.
myFilename	A string that gives name of the Figure, e.g. 'zonalMeans.pdf'.
legendLocation	A string that specifies the location of the model ID legend, e.g. 'topright'.
plotRef	Logical. If FALSE, reference data is not plotted. Default is TRUE.
subcaption	A string that defines the subcaption of the figure, e.g. '(a)'.

Value

A figure that shows zonal mean values of model and reference data.

Examples

library(amber)
library(classInt)
library(doParallel)
library(foreach)
library(Hmisc)
library(latex2exp)
library(ncdf4)
library(parallel)
library(raster)
library(rgdal)
library(rgeos)
library(scico)
library(sp)
library(stats)
library(utils)
library(viridis)
library(xtable)

mod01.path <- system.file('extdata/SIMod01', package = 'amber')
mod02.path <- system.file('extdata/SIMod02', package = 'amber')
mod.path.list <- list(mod01.path, mod02.path)
modelIDs <- c('CLASSIC.CRUJRAv2', 'CLASSIC.GSWP3W5E5')
myVariables <- c('GPP.MODIS', 'LAI.AVHRR')
plotZonalMeans(mod.path.list, modelIDs, myVariables, metric = 'mean')

amber documentation built on Aug. 28, 2020, 5:08 p.m.