Reanalysis {SailoR}

Surface Wind from different Reanalyses

Description

Contains the information of both components of monthly averaged surface wind in January over the Northern Hemisphere as included in different state-of-the-art reanalyses during the period 2011-2018. ERA5 reanalysis was defined as the reference, and the remaining gridded products were included in the model matrix. Particularly, the original NCEP/NCAR first generation reanalysis, MERRA2, CFSv2 and ERA-Interim. All of them have been bilinearly interpolated to NCEP/NCAR original grid.

Usage

data("Reanalysis")

Format

A list with two matrices:

ref

a matrix including the data from reanalysis ERA5.

mod

a matrix including the data from other reanalyses.

Each of these data frames includes 5 variables:

V1

a vector with longitude data.

V2

a vector with latitude data.

V3

a factor defining if the data belong to the reanalysis CFSv2 (cfsv2), ERA-Interim (ei), MERRA2 (merra2) or NCEP/NCAR reanalysis (nnra).

V4

zonal component of surface wind (m/s).

V5

meridional component of surface wind (m/s).

Examples

# Load the data
data(Reanalysis)

# Parameters
Uxlim=c(-0.5,1)
Uylim=c(-1,0.5)
Uxlab<-"Ux (m/s)"
Uylab<-"Uy (m/s)"
plotmain<-"Reanalyses"
sfactor<-0.15

# Create ref, mod objects

weightWithLat=TRUE

if(weightWithLat){
  dlats=as.numeric(Reanalysis[["ref"]][,2])
  rlats=pi*dlats/180.
  lweights=sqrt(cos(rlats))

  ref<-data.frame(Reanalysis[["ref"]][,3],
                  lweights*as.numeric(Reanalysis[["ref"]][,4]),
                  lweights*as.numeric(Reanalysis[["ref"]][,5]))
}else{

  ref<-data.frame(Reanalysis[["ref"]][,3],
                  as.numeric(Reanalysis[["ref"]][,4]),
                  as.numeric(Reanalysis[["ref"]][,5]))

}
names(ref)<-c("mod","U","V")

mod<-data.frame(Reanalysis[["mod"]][,3],
                as.numeric(Reanalysis[["mod"]][,4]),
                as.numeric(Reanalysis[["mod"]][,5]))
names(mod)<-c("mod","U","V")
  
# Sailors

#--------------------------------------------------------
# Example 1: Figure 7 (left) from Sáenz et al. (2020)
#--------------------------------------------------------

SailoR.Plot(ref, mod, ColourList=NULL, sfactor, docenter=FALSE,
            Uxlim, Uylim, Uxlab, Uylab, plotmain, plotlegend=TRUE,
            plotRMSElegend=TRUE, plotscalelegend=TRUE,
            RMSE_legend_Rounding=1, RMSE_legend_units = " m/s",
            referenceName="ERA5")

#--------------------------------------------------------
# Example 2: Figure 7 (right) from Sáenz et al. (2020)
#--------------------------------------------------------

SailoR.Plot(ref, mod, ColourList=NULL, sfactor, docenter=TRUE,
            Uxlim, Uylim, Uxlab, Uylab, plotmain, plotlegend=TRUE,
            plotRMSElegend=TRUE, plotscalelegend=TRUE,
            RMSE_legend_Rounding=1, RMSE_legend_units = " m/s",
            referenceName="ERA5")