Plot zonal mean values of model and reference data in a single figure
This function plots zonal mean values of model and reference data, computed from zonalMeanStats, for a selection of variables.
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 = "")
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)'. |
A figure that shows zonal mean values of model and reference data.
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')Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.