extended_jc: extended_jc

Description

Gets daily classification of weather types at every grid-point over the map or selected area of interest according to the Jenkison and Collison scheme.

Usage

extended_jc(mslp, loni, lati, times, gale = FALSE, num_cores = 2)

Arguments

mslp

3-Dimensional array ([loni,lati,time]) with mean sea level pressure in Pa.

loni

Array with longitude values.

lati

Array with latitude values.

times

Array with the dates used.

gale

Logial. If TRUE, the function returns also Gale days.

num_cores

Number of cores (2 by default).

Value

A list with two objects:

A list of eleven matrix of daily frequencies of weather types ("wtypes"). Each matrix is a 3D array [loni,lati,times] and it refers to each weather type (N,NE,E,SE,S,SW,W,A,C and U).
A list of six matrix of daily frequency of airflow indices ("indices"). Each matrix is a 3D array [loni,lati,times] and it refers to each airflow index (W,S,TF,ZW,ZS,Z and D).
A list with the central points for which the classification is applied.

References

Otero, N., Sillmann, J. & Butler, T. Assessment of an extended version of the Jenkinson<U+2013>Collison classification on CMIP5 models over Europe Climate Dynamics. https://doi.org/10.1007/s00382-017-3705-y

Examples

Run this code

# NOT RUN {
# Load data
data(press)
# Get coordinates
longitudes <- press$loni
latitudes  <- press$lati
times      <- press$dates

# Example when the classification is restricted to an area
# Select longitudes and latitudes within the European domain: -10W,40E, 40N,70N
ilon  <- which(longitudes>(-10)&longitudes<40)
loni  <- longitudes[ilon]
ilat  <- which(latitudes>40&latitudes<70)
lati  <- latitudes[ilat]
cwtEU <- extended_jc(press$msl[ilon,ilat,], loni, lati, times, gale=FALSE, num_cores=2)
 
# }
# NOT RUN {
# Not run
# This is a long running example
# Get the classification for the whole map, all longitudes and latitudes
cwtGlobal <- extended_jc(press$msl, longitudes, latitudes, times, gale=FALSE, num_cores=2)
# }
# NOT RUN {
# }