errror loading ERA5 data
Hi Nuria,
here is the code which return the error message I mentioned last week:
library(ncdf)
setwd("/scratch/Earth/lcaron/WINDSTORMS")
library(s2dverification)
library(ncdf)
library(easyNCDF)
library(ggplot2)
library(pracma)
outputdir='/scratch/Earth/lcaron/WINDSTORMS/'
figpath='/home/Earth/lcaron/Desktop/Windstorms/NewFig/'
datapath='/home/Earth/lcaron/Desktop/Windstorms/data'
start_year=1979
end_year=2010
years<-seq(start_year,end_year,1)
numyears<-length(years)
yrs=c(start_year,end_year)
latbnd<-c(-10,70)
lonbnd<-c(120,360)
#latbnd<-c(10,70)
#lonbnd<-c(-80,0)
mm <-'08' # First month
months<-c('August','September','October','November','December','January','February','March')
seasons<- c('ASO',"SON","OND","NDJ","DJF","JFM")
numseas<-length(seasons)
#
Mystart <- as.Date(paste0(yrs[1],mm,'01'), "%Y%m%d")
#
Myend <- as.Date(paste0(yrs[2],mm,'01'), "%Y%m%d")
#
# Build sequence of dates to load forecast for (date of initialization)
Mydateseq <- format(seq(Mystart, Myend, by = "year"), "%Y%m%d")
var="g500"
obsdata <- c('era5')
#obsdata <- c('jra55')
# Load the data. Averaging over the region of interest is also performed here.
tosdata <- Load(var,
exp = NULL, # experiment
obs = obsdata, # obs/reanalysis used for comparison
sdates = Mydateseq, # series of forecast to load
leadtimemin = 1, leadtimemax = 8, #
lonmin = lonbnd[1], lonmax = lonbnd[2], # max/min longitude
latmin = latbnd[1], latmax = latbnd[2], # max/min latitude
storefreq = "monthly", #sampleperiod = 1,
nmemberobs = 1, output = "lonlat",
configfile = '/home/Earth/lcaron/git-stuff/s2dverification/inst/config/BSC.conf')
If I use JRA55 instead of ERA5, it works fine. Any insight as to what the problem is would be appreciated...