Issues about ''retrieve''

[QuocPhamBao]

Hello
I used to implement this to extract some data and it is successful:
r1 <- retrieve('slp.1966.nc', lev = 500, it = as.Date(c('1966-01-08',
'1966-01-17')), lon=c(-180,180),lat=c(35,90))

However, I recently used the same for the data in NCEP: https://downloads.psl.noaa.gov/Datasets/ncep.reanalysis/Dailies/surface/
and it is not working.
This is the error notification:
"Warning : Calendar attribute has not been found in the meta data and will be set automatically."
Error in grep(x, seq(length(idim))) : invalid 'pattern' argument

I still can use it for the old data that I download from NCEP last year (July), but I am unable to use new data from the NCEP.
I attached two files (an old file and a new file which were downloaded from NCEP).
I do not know if there is any update about this package?

Regards

[abdelkaderm]

Hi

I will have a look and come back to you.

Bests,

[abdelkaderm]

We have recently made some changes that affected the command line you are using and we will try to fix it. In the meantime, you can try the following two lines separately that should work.

x <- retrieve('Downloads/slp.1966.nc', param = 'slp', lev = 500, lon=c(-180,180),lat=c(35,90))
x <- subset(x,it =as.Date(c('1966-01-08','1966-01-17')))

[QuocPhamBao]

Thank you ! It is very useful.

[abdelkaderm]

You are welcome 👍