amber (version 1.0.3)

plotEnsembleHovmoeller: Hovmoeller Diagram for model ensemble

Description

This function plots Hovmoeller diagrams of monthly climatological mean values and biases computed by scores.grid.time for multiple ensemble members

Usage

plotEnsembleHovmoeller(mod.path.list = mod.path.list,
  modelIDs = modelIDs, myVariables = myVariables, myBin = 20,
  gridCellWidth = 2, my.ylim = c(-100, 100), plot.width = 8.4,
  plot.height = 5, outputDir = FALSE)

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').

myBin

An integer number that defines the latitudinal range used for computing the zonal mean. For instance, a value of 10 implies that a zonal mean is computed for every 10 degrees latitude.

gridCellWidth

A number that is used as a factor to adjust the width of grid cells, e.g. 1.

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.

Value

Figures in PDF format.

Examples

# NOT RUN {
# }
# NOT RUN {
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 <- 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-MODIS', 'GPP-GOSIF')

plotEnsembleHovmoeller(mod.path.list = mod.path.list,
modelIDs = modelIDs, myVariables = myVariables, myBin = 20, gridCellWidth = 2,
my.ylim = c(-100, 100), plot.width = 8.4, plot.height = 5.0)
# }
# NOT RUN {
 #donttest
# }