The Weather Research and Forecasting (WRF) Model is a numerical weather prediction (NWP) system. NWP refers to the simulation and prediction of the atmosphere with a computer model, and WRF is a set of software for this.
meteoForecast
is a new R package that implements functions to download data from the Meteogalicia and OpenMeteo NWP-WRF services using the NetCDF Subset Service.
Read this introduction for additional information and examples.
Ricardo Bessa dijo:
Oscar, another interesting add-on for the function could be a selection of the spatial resolution. In the WRF case it would be interesting to choose between 36, 12 and 4 km data.
Oscar Perpiñán Lamigueiro dijo:
I agree. Could you please open an issue here https://github.com/oscarperpinan/meteoForecast/issues? Thanks!
Ricardo Bessa dijo:
Mauricio, the GFS model is already used by the WRF from Meteogalicia. It provides the initial conditions for the mesoscale model (the WRF model). The ECMEF is a paid service. Open Weather Map I think that it does not provide high resolution forecasts.
Oscar: some aspects can be improved in a future version. For instance, when getting data for an historical period some files for some days are missing and it returns an error.
Oscar Perpiñán Lamigueiro dijo:
Hi,
You are right about GFS and ECMWF. I didn’t know OpenWeatherMap. It’s like forecast.io, isn’t it?
Anyhow, I have just published a wiki with a collection of data sources I have been collating during last year. Some of them could be included in
meteoForecast
.Regarding the error, what function are you using?
getRasterDays
andgetPointDays
are coded to manage missing days. I would appreciate if you could open an issue. It helps me to track the bugs.Mauricio Zambrano-Bigiarini dijo:
Hi Ricardo,
Even if GFS is already used by the WRF from Meteogalicia, when I try to get a forecast for south America I get an error:
library(meteoForecast)
wrfDays2 <- getRasterDays('cft',
start = '2014-01-01',
end = '2014-01-05',
box = c(-72, 70, -35, 33))
Error in R_nc_open: NetCDF: Unknown file format
Error in open.ncdf(filename) :
Error in open.ncdf trying to open file /tmp/RtmpBzM3Ie/cft_20140105_00.nc
Error in x[[1]] : subscript out of bounds
Oscar Perpiñán Lamigueiro dijo:
Hi Mauricio,
GFS is used to initialize the WRF model, but the output of the WRF Meteogalicia system is restricted to a smaller area.
By the way, I have to improve the error message.
Ricardo Bessa dijo:
Great package! When I run the following code «tempK <- getPoint(c(0, 40), vars = 'temp')" I get the following error message: Error in getPoint(c(0, 40), vars = "temp") :
Data not found. Check the date and variables name.
URL: http://mandeo.meteogalicia.es/thredds/ncss/grid/wrf_2d_12km/fmrc/files/20140709/wrf_arw_det_history_d02_20140709_0000.nc4?var=temp&point=true&longitude=0&latitude=40
It seemsn that the last file is not availlable in the server, but it is.
Oscar Perpiñán Lamigueiro dijo:
It seems that there is a problem at the server today: «HTTP Status 500 – Internal Server Error». It works without problems if you choose a different day.
I will check it again later.
Thanks
Oscar.
Pingback: meteoForecast, a package to obtain NWP-WRF fore...
Mauricio Zambrano-Bigiarini dijo:
Thanks Oscar for such a nice and useful package !
In the future, would your consider to add some global forecast services (not only for Europe, as Meteogalicia and OpenMeteo) ?.
Oscar Perpiñán Lamigueiro dijo:
Yes, sure! Do you have any suggestion?
Mauricio Zambrano-Bigiarini dijo:
1) Global Forecast System (GFS) :
http://www.ncdc.noaa.gov/data-access/model-data/model-datasets/global-forcast-system-gfs
2) Open Weather Map:
http://openweathermap.org/Maps?zoom=10&lat=51&lon=7&layers=B0TTTFFT
3) ECMWF (but I think this is a paid product ):
http://www.ecmwf.int/en/forecasts/datasets/dataset-i-i-atmospheric-fields-high-resolution-forecast
If you want some collaboration, I’m very willing to help 🙂
Oscar Perpiñán Lamigueiro dijo:
Thanks!. I will check them as soon as possible, and I will contact you then.