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]]
______________________________________________
[email protected] mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
