Description Usage Arguments Value Examples
View source: R/plotEnsembleMean.R
This function plots ensemble mean, minimum, and maximum values of a statistical metric computed by scores.grid.time and scores.grid.notime.
1 2 3 4 5 | plotEnsembleMean(long.name, metric, mod.path.list, modelIDs, myVariables,
shp.filename = system.file("extdata/ne_110m_land/ne_110m_land.shp",
package = "amber"), my.xlim = c(-180, 180), my.ylim = c(-60, 85),
plot.width = 5, plot.height = 7, outputDir = FALSE,
subcaption = "")
|
long.name |
A string that gives the full name of the variable, e.g. 'Gross primary productivity' |
metric |
A string that specifies what statistical metric should be plotted. This includes for instance 'bias', 'crmse', 'phase', 'iav', 'bias-score', 'rmse-score', 'phase-score', and 'iav-score'. |
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'). |
shp.filename |
A string that gives the coastline shapefile |
my.xlim |
An R object that gives the longitude range that you wish to plot, e.g. c(-180, 180) |
my.ylim |
An R object that gives the longitude range that you wish to plot, e.g. c(-90, 90) |
plot.width |
Number that gives the plot width, e.g. 8 |
plot.height |
Number that gives the plot height, e.g. 4 |
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. |
subcaption |
A string that defines the subcaption of the figure, e.g. '(a)'. |
Figures in PDF format.
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 28 29 30 | 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)
long.name <- 'Gross Primary Productivity'
metric <- 'mod-mean'
mod01.path <- paste(system.file('extdata', package = 'amber'), 'model01', sep = '/')
mod02.path <- paste(system.file('extdata', package = 'amber'), 'model02', sep = '/')
mod.path.list <- list(mod01.path, mod02.path)
modelIDs <- c('CLASSIC.CRUJRAv2', 'CLASSIC.GSWP3W5E5')
myVariables <- c('GPP-GOSIF', 'GPP-MODIS')
plotEnsembleMean(long.name, metric, mod.path.list, modelIDs, myVariables,
plot.width = 5, plot.height = 5.5)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.