When you start rNOMADS it says:
> Welcome to rNOMADS 2.4.2 "Pandaemonium Fortress"!
> Questions? Follow @rNOMADS_r on Twitter or send a message to rnomads-user
at lists.r-forge.r-project.org
>
Likely to get much more knowledgeable answers there.
-Roy
> On Jul 28, 2020, at 1:45 PM, Philip <herd_dog at cox.net> wrote:
>
> 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]]
>
> ______________________________________________
> R-help at r-project.org 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.
**********************
"The contents of this message do not reflect any position of the U.S.
Government or NOAA."
**********************
Roy Mendelssohn
Supervisory Operations Research Analyst
NOAA/NMFS
Environmental Research Division
Southwest Fisheries Science Center
***Note new street address***
110 McAllister Way
Santa Cruz, CA 95060
Phone: (831)-420-3666
Fax: (831) 420-3980
e-mail: Roy.Mendelssohn at noaa.gov www: https://www.pfeg.noaa.gov/
"Old age and treachery will overcome youth and skill."
"From those who have been given much, much will be expected"
"the arc of the moral universe is long, but it bends toward justice"
-MLK Jr.