[R] rNOMADS Package
Philip
herd_dog @end|ng |rom cox@net
Tue Jul 28 22:45:57 CEST 2020
Don’t know if anyone out there can help me since the rNOMADS package is so specialized and it requires a download of wgrib2 software from NOAA – National Oceanographic and Atmospheric Administration.
tion.
Is there a way to grab multiple forecasts with the ReadGrib function (page 33 of the documentation). I’m trying to get the rapid update forecast (RAP) on the 13 km scale
RAPOut <- CrawlModels(abbrev = "rap", depth = 1)
RAPParameters <- ParseModelPage(RAPOut[1])
RAPPred <- RAPParameters$pred[grep("t12z", RAPParameters$pred)]#12 Zulu is 5 am in Phoenix
RAPPred2 <- RAPPred[grep(c("awp130"), RAPPred)] #picks out 13 km forecasts
RAPPred3 <- RAPPred2[grep(c("pgrbf00|pgrbf01"),RAPPred2)] #Forecast for 12 Z and 13 Z
levels <- c("800 mb")
variables <- c("TMP")
RAPInfo <- GribGrab(RAPOut[1], RAPPred3, levels, variables) #Everything works perfectly to this point. The two requested forecasts are listed below. As you can see [[1]] is for 12 Z and [[2]] is for 13 Z – red in 18 point type.
RAPInfo
[[1]]
[[1]]$file.name
[1] "C:\\Users\\Owner\\Documents\\Ballooning\\WeatherBriefing/rap.t12z.awp130pgrbf00.grib2.grb"
[[1]]$url
[1] "https://nomads.ncep.noaa.gov/cgi-bin/filter_rap.pl?file=rap.t12z.awp130pgrbf00.grib2&lev_800_mb=on&var_TMP=on&dir=%2Frap.20200728"
[[2]]
[[2]]$file.name
[1] "C:\\Users\\Owner\\Documents\\Ballooning\\WeatherBriefing/rap.t12z.awp130pgrbf01.grib2.grb"
[[2]]$url
[1] "https://nomads.ncep.noaa.gov/cgi-bin/filter_rap.pl?file=rap.t12z.awp130pgrbf01.grib2&lev_800_mb=on&var_TMP=on&dir=%2Frap.20200728"
RAPData <- ReadGrib(RAPInfo[[2]]$file.name, levels, variables,
domain=c(-113.20,-112.78,33.70,33.40), #an area west of Phoenix
domain.type="latlon",
file.type="grib2")
I can grab the 12 Z forecast with [[1]] in the first line above or the 13 Z forecast if I enter [[2]]. But is there a way to get both at t he same time?
Thanks,
Philip Heinrich
[[alternative HTML version deleted]]
More information about the R-help
mailing list