[R] rNOMADS Package

[Roy Mendelssohn - NOAA Federal]

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 using 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 using 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 using 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 using 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.