[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