amber (version 1.0.3)

zonalMeanStats: Zonal mean plots of AMBER results (bias, bias scores, etc)

Description

This function computes zonal mean values of model and reference data and the zonal mean bias, centralized root-mean-squre error, phase, inter-annual variability, and corresponding scores. The computation is based on the NetCDF files produced by scores.grid.time.

Usage

zonalMeanStats(inputDir, outputDir = FALSE)

Arguments

inputDir

A string that gives the location of NetCDF files produced by scores.grid.time and scores.grid.notime, e.g. '/home/project/study'.

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

A list with two tables. The first table gives the zonal mean values of the mean, bias, standard deviation, centralized root-mean-squre error, month of seasonal peak, phase, inter-annual variability, and corresponding scores for each variable and globally gridded reference data set. The second table gives the physical units of each variable in LaTeX notation (e.g. 'W m$^-2$'). Both tables are written to two text files (zonalMeanStats and zonalMeanStatsUnits) if the user specifies an output directory.

Examples

# 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)

inputDir <- paste(system.file('extdata', package = 'amber'), 'zonalMeanStats', sep = '/')
zonalMeanStats(inputDir, outputDir = FALSE)
# }