From tmpwanadoo at wanadoo.fr  Thu Mar  2 19:58:42 2017
From: tmpwanadoo at wanadoo.fr (tmpwanadoo)
Date: Thu, 2 Mar 2017 19:58:42 +0100
Subject: [Rnomads-user] SubsetNOMADS
Message-ID: <73d1b021-378e-44c4-26e0-47469484b882@wanadoo.fr>

Hi Danny,

Would it be possible to have a script example showing how to use 
SubsetNOMADS to extract the forecast for a specific time when 
downloading a grib with several times ?

Below is an example of the many things a tried !

Best

S?bastien


 > #Get data from NOMADS real time server
 > tmp.data <- DODSGrab(latest.model, latest.model.run, variables, time, 
lon, lat, display.url=FALSE)#, verbose=TRUE)
 > names(tmp.data)
[1] "model.run.date" "forecast.date"  "variables" "levels"         
"ensembles"      "lon" "lat"            "value"
[9] "request.url"
 > summary(tmp.data)
                Length Class   Mode
model.run.date 129444 -none-  character
forecast.date  129444 POSIXct numeric
variables      129444 -none-  character
levels         129444 -none-  character
ensembles      129444 -none-  character
lon            129444 -none-  numeric
lat            129444 -none-  numeric
value          129444 -none-  numeric
request.url         1 -none-  character

 > unique(tmp.data$forecast.date)
[1] "2017-03-02 13:00:00 CET" "2017-03-02 16:00:00 CET" "2017-03-02 
19:00:00 CET" "2017-03-02 22:00:00 CET"

 > unique(tmp.data$variables)
[1] "ugrd10m"

 > unique(tmp.data$model.run.date)
[1] "20170302gfs_0p25_12z"

 > unique(tmp.data$request.url)
[1] 
"http://nomads.ncep.noaa.gov:9090/dods/gfs_0p25/gfs20170302/gfs_0p25_12z.ascii?ugrd10m[0:3][120:280][520:720]"

 > SubsetNOMADS(model.data=tmp.data, forecast.date=c("2017-03-01 
07:00:00 CET"))
$model.run.date
character(0)

$forecast.date
character(0)

$variables
character(0)

$levels
character(0)

$ensembles
character(0)

$lon
numeric(0)

$lat
numeric(0)

$value
numeric(0)

$request.url
[1] 
"http://nomads.ncep.noaa.gov:9090/dods/gfs_0p25/gfs20170302/gfs_0p25_12z.ascii?ugrd10m[0:3][120:280][520:720]"


---
L'absence de virus dans ce courrier ?lectronique a ?t? v?rifi?e par le logiciel antivirus Avast.
https://www.avast.com/antivirus


From danny.c.bowman at gmail.com  Fri Mar  3 04:59:57 2017
From: danny.c.bowman at gmail.com (Daniel Bowman)
Date: Thu, 2 Mar 2017 20:59:57 -0700
Subject: [Rnomads-user] SubsetNOMADS
In-Reply-To: <db1c215a09c7418e9f21c578bb97fba4@BL2PR03MB564.namprd03.prod.outlook.com>
References: <db1c215a09c7418e9f21c578bb97fba4@BL2PR03MB564.namprd03.prod.outlook.com>
Message-ID: <CAN=ZksNFdQyV_93+ES7d0AVUQrB-HU1GTbF01KLxjLWEpOd+og@mail.gmail.com>

Hi Sebastien,

Sorry for the delay in responding.  Yes, it looks like you've found a
problem! I've been able to reproduce it, and I see where it resides in
my code.  I will try and fix it this weekend if I have time.

For now, just set the time input to DODSGrab to one value.  That way,
you'll only get one model run.

Best,

Danny

On Thu, Mar 2, 2017 at 11:58 AM, tmpwanadoo <tmpwanadoo at wanadoo.fr> wrote:
> Hi Danny,
>
> Would it be possible to have a script example showing how to use
> SubsetNOMADS to extract the forecast for a specific time when
> downloading a grib with several times ?
>
> Below is an example of the many things a tried !
>
> Best
>
> S?bastien
>
>
>  > #Get data from NOMADS real time server
>  > tmp.data <- DODSGrab(latest.model, latest.model.run, variables, time,
> lon, lat, display.url=FALSE)#, verbose=TRUE)
>  > names(tmp.data)
> [1] "model.run.date" "forecast.date"  "variables" "levels"
> "ensembles"      "lon" "lat"            "value"
> [9] "request.url"
>  > summary(tmp.data)
>                 Length Class   Mode
> model.run.date 129444 -none-  character
> forecast.date  129444 POSIXct numeric
> variables      129444 -none-  character
> levels         129444 -none-  character
> ensembles      129444 -none-  character
> lon            129444 -none-  numeric
> lat            129444 -none-  numeric
> value          129444 -none-  numeric
> request.url         1 -none-  character
>
>  > unique(tmp.data$forecast.date)
> [1] "2017-03-02 13:00:00 CET" "2017-03-02 16:00:00 CET" "2017-03-02
> 19:00:00 CET" "2017-03-02 22:00:00 CET"
>
>  > unique(tmp.data$variables)
> [1] "ugrd10m"
>
>  > unique(tmp.data$model.run.date)
> [1] "20170302gfs_0p25_12z"
>
>  > unique(tmp.data$request.url)
> [1]
> "http://nomads.ncep.noaa.gov:9090/dods/gfs_0p25/gfs20170302/gfs_0p25_12z.ascii?ugrd10m[0:3][120:280][520:720]"
>
>  > SubsetNOMADS(model.data=tmp.data, forecast.date=c("2017-03-01
> 07:00:00 CET"))
> $model.run.date
> character(0)
>
> $forecast.date
> character(0)
>
> $variables
> character(0)
>
> $levels
> character(0)
>
> $ensembles
> character(0)
>
> $lon
> numeric(0)
>
> $lat
> numeric(0)
>
> $value
> numeric(0)
>
> $request.url
> [1]
> "http://nomads.ncep.noaa.gov:9090/dods/gfs_0p25/gfs20170302/gfs_0p25_12z.ascii?ugrd10m[0:3][120:280][520:720]"
>
>
> ---
> L'absence de virus dans ce courrier ?lectronique a ?t? v?rifi?e par le logiciel antivirus Avast.
> https://www.avast.com/antivirus
>
> _______________________________________________
> Rnomads-user mailing list
> Rnomads-user at lists.r-forge.r-project.org
> http://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rnomads-user