Description Usage Arguments Value Examples
View source: R/plotEnsembleHovmoeller.R
This function plots Hovmoeller diagrams of monthly climatological mean values and biases computed by scores.grid.time for multiple ensemble members
1 2 3 4 | 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)
|
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. |
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 | 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)
#donttest
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.