Skip to content

National Weather Service Data

8 messages · Philip, stephen sefick, Thomas Adams +2 more

Philip
# 
Is anyone out there familiar with rNOMADS?  It is a package to get into National Weather Service forecasting data with R?

I'm not sure the Weather Service software named wgrib2 loaded correctly because some of the stuff won't run and I can't make much sense out of some of the output.

Thanks.
stephen sefick
# 
I am unfamiliar with Rnomads. Could you provide a minimal reproducable
example? You are more likely to receive help this way.
On Thu, Jul 2, 2020, 18:06 Philip <herd_dog at cox.net> wrote:

            

      
      
      
    

        

      
      
    

        

  
  


  


      

    

    

      
      

        


  

        

      Thomas Adams
    

    

      
      #
    

  


  

      
Hi Philip!

I'm a little familiar with rNOMADS... I tried following the example for
'ArchiveGribGrab' using a more recent date

#An example for the Global Forecast System
#Get data for January 1 2014
#Temperature at 2 m above ground
#3 hour prediction
# using GRIB
abbrev <- "gfsanl"
model.date <- 20200601
model.run <- 00
preds <- 3

I got this result...
URL 'https://nomads.ncdc.noaa.gov/data/gfsanl/202006/20200601/': status was
'Couldn't resolve host name'

Even this part of the URL was not found: https://nomads.ncdc.noaa.gov -- so
there is documentation problems.

I did get the first example, using contour to work fine. You probably
should contact Daniel C. Bowman <danny.c.bowman at gmail.com> directly with
issues, because the problems you are seeing may be resolvable only by him...

Best,
Tom

        
On Thu, Jul 2, 2020 at 6:06 PM Philip <herd_dog at cox.net> wrote:

        

            

      
      
      
    

        

      
      
    

        

  
    

      
      
    

  


  


      

    

    
3 days later

    

      
      

        


  

        

      Philip
    

    

      
      #
    

  


  

      
I am trying to access National Weather Service forecasting data through the rNOMADS package.  I?m not sure if the Weather Service software ? grib2 ? loaded correctly.  Second, some of the examples in the rNOMADS documentation seem to run correctly but I?m not sure what the output means.  Any anvise would be greatly appreciated.

1 - I tried to load the wgrib2 software from instructions in the following website:

    https://bovineaerospace.wordpress.com/2015/04/26/how-to-install-rnomads-with-grib-file-support-on-windows/

2 ? the instructions say that if it loaded correctly I should get a laundry list similar to what is below from the command: >system(?wgrib2?).



The list I get looks different.  Below is the first 20 or so entries.  How can I check to see if the wgrib2 loaded correctly?

wgrib2 v0.1.9.9 9/2013 Wesley Ebisuzaki, Reinoud Bokhorst, Jaakko Hyv??tti, Dusan Jovic, Kristian Nilssen, Karl Pfeiffer, Pablo Romero, Manfred Schwarb, Arlindo da Silva, Niklas Sondell, Sergey Varlamov
-0xSec           inv  X      Hex dump of section X (0..8)
-MM              inv         reference time MM
-N_ens           inv         number of ensemble members
-RT              inv         type of reference Time
-S               inv         simple inventory with minutes and seconds (subject to change)
-Sec0            inv         contents of section0
-Sec3            inv         contents of section 3 (Grid Definition Section)
-Sec4            inv         Sec 4 values (Product definition section)
-Sec5            inv         Sec 5 values (Data representation section)
-Sec6            inv         show bit-map section
-Sec_len         inv         length of various grib sections
-T               inv         reference time YYYYMMDDHHMMSS
-V               inv         diagnostic output
-VT              inv         verf time = reference_time + forecast_time (YYYYMMDDHHMMSS)
-YY              inv         reference time YYYY

3 ? As I said, some of the documentation examples work and for some I get error messages.  Below is an example of one that seemed to work but I don?t understand the output.

#GribInfo - page 20
urlsOut <- CrawlModels(abbrev="gfs_0p50",depth=2)
ModelParameters <- ParseModelPage(urlsOut[2])#[1] is most recent model
MyPred <- ModelParameters$pred[grep("06$",ModelParameters$pred)]
   Levels <- c("2_m_above_ground","800_mb")
   Variables <- c("TMP","RH")
   GribInfo <- GribGrab(urlsOut[2],MyPred,Levels,Variables)
GribInv <- GribInfo(GribInfo[[1]]$file.name,"grib2")

The command GribInv$inventory returns:

$inventory
[1] "1:0:d=2020070606:TMP:800 mb:6 hour fcst:"                "2:148450:d=2020070606:RH:800 mb:6 hour fcst:"           
[3] "3:414132:d=2020070606:TMP:2 m above ground:6 hour fcst:" "4:571266:d=2020070606:RH:2 m above ground:6 hour fcst:" 

This is supposed to be temperature and relative humidity 2 meters above the ground and at 800 milibars for 2020 ? July 6 ? at ZULU time 0600.  But I have no idea what the numbers 414132 ? second line ? mean.

Any advise would be greatly appreciated.

Philip Heinrich

From: stephen sefick 
Sent: Thursday, July 2, 2020 3:20 PM
To: Philip 
Cc: r-help 
Subject: Re: [R] National Weather Service Data

I am unfamiliar with Rnomads. Could you provide a minimal reproducable example? You are more likely to receive help this way.

        
On Thu, Jul 2, 2020, 18:06 Philip <herd_dog at cox.net> wrote:

        
Is anyone out there familiar with rNOMADS?  It is a package to get into National Weather Service forecasting data with R?

  I'm not sure the Weather Service software named wgrib2 loaded correctly because some of the stuff won't run and I can't make much sense out of some of the output.

  Thanks.

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

  
  


  


      

    

    

      
      

        


  

        

      Rasmus Liland
    

    

      
      #
    

  


  

      
On 2020-07-06 09:29 -0700, Philip wrote:

            

      
      
      
    

        
Dear Philip,

This page[1] says it's the byte location.  
?byte location of the start of the message 
starting from 0? 

Best,
Rasmus

[1] https://www.cpc.ncep.noaa.gov/products/wesley/wgrib2/default_inv.html

-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 833 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20200706/64019fec/attachment.sig>

  
  


  


      

    

    

      
      

        


  

        

      Roy Mendelssohn - NOAA Federal
    

    

      
      #
    

  


  

      
Hi Philip:

Results look correct to me.  This might help you:

https://www.cpc.ncep.noaa.gov/products/wesley/wgrib2/default_inv.html

-Roy

            

      
      
      
    

        

      
      
    

        
**********************
"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.

  
  


  


      

    

    

      
      

        


  

        

      Philip
    

    

      
      #
    

  


  

      
Thanks for getting back to me.  It is good to know that I am on the right 
track.

-----Original Message----- 
From: Rasmus Liland
Sent: Monday, July 6, 2020 10:42 AM
To: r-help
Subject: Re: [R] National Weather Service Data

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

  
  


  


      

    

    

      
      

        


  

        

      Rasmus Liland
    

    

      
      #
    

  


  

      
On 2020-07-06 17:41 -0700, Philip wrote:

            

      
      
      
    

        
Oh, it's always such a pleasure to be of help 
to someone in need :-)

-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 833 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20200707/e466908f/attachment.sig>