<table cellspacing="0" cellpadding="0" border="0"><tr><td valign="top"><p dir="ltr"><font size ="2"><font size ="2">Thank you! Much appreciated. </font></font></p>
<p dir="ltr"><font size ="2"><font size ="2"><a href="https://overview.mail.yahoo.com/mobile/?.src=Android">Sent from Yahoo Mail on Android</a></font></font></p>
<table cellspacing="0" cellpadding="0" border="0"> <tbody> <tr> <td valign="top"> <div style="font-family:Roboto, sans-serif;color:#7e7d80;">From:"Bowman, Daniel" <haksaeng@live.unc.edu><br/>Date:Wed, Jan 7, 2015 at 4:38 PM<br/>Subject:Re: [Rnomads-user] ModelGrid function with Dods data<br/><br/></div>
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri, Arial, Helvetica, sans-serif;">
<p>Hi Adam,</p>
<p><br clear="none">
</p>
<p>Yes, ModelGrid is the culprit. First of all, GFS is in "Lat Lon" and the resolution is 0.5. So your code should read:</p>
<p><br clear="none">
</p>
<p>GFS.modelgrid <- ModelGrid( GFS.Model.Data, grid.type = "latlon", resolution = c(0.5, 0.5),levels = c("surface","2_m_above_ground"), variables)</p>
<p><br clear="none">
</p>
<p>However, since you pulled multiple forecasts from DODS, the above won't work unless you separate out each forecast. This is a problem I ran into recently but was too lazy to solve, but now you've given me a reason to put together an example. I'll get it
to you by the end of the day tomorrow.</p>
<p><br clear="none">
</p>
<p>Danny<br clear="none">
</p>
<p><br clear="none">
</p>
<div id="Signature">
<div style="font-family:Calibri, Arial, Helvetica, sans-serif;margin:0;">
Daniel C. Bowman<br clear="none">
DoctGFS.modelgrid <- ModelGrid(model.data = GFS.Model.Data, grid.type = "cartesian", resolution = c(-140,40),levels = c("surface","2_m_above_ground"), variables, model.domain = NULL, cartesian.nodes = (list(x=0.5,y=0.5))oral Candidate in Geophysics<br clear="none">
Mitchell 315<br clear="none">
Department of Geological Sciences<br clear="none">
The University of North Carolina at Chapel Hill<br clear="none">
phone: 575-418-8555<br clear="none">
email: daniel.bowman@unc.edu<br clear="none">
web:http://geosci.unc.edu/page/daniel-c-bowman<br clear="none">
<br clear="none">
</div>
</div>
<div class="yQTDBase yqt4731101039" id="yqt36340"><div style="color:rgb(33, 33, 33);">
<hr tabindex="-1" style="display:inline-block;width:98%;">
<div dir="ltr" id="divRplyFwdMsg"><font style="font-size:11pt;" color="#000000" face="Calibri, sans-serif"><b>From:</b> rnomads-user-bounces@lists.r-forge.r-project.org <rnomads-user-bounces@lists.r-forge.r-project.org> on behalf of Adam Simkowski <akssimkowski@sbcglobal.net><br clear="none">
<b>Sent:</b> Wednesday, January 7, 2015 4:12 PM<br clear="none">
<b>To:</b> rnomads-user@lists.r-forge.r-project.org<br clear="none">
<b>Subject:</b> [Rnomads-user] ModelGrid function with Dods data</font>
<div> </div>
</div>
<div>
<div style="color:#000;background-color:#fff;font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:12px;">
<div>Hi all, </div>
<div> </div>
<div>I am attempting to build an array of GFS model data, but I am running into some trouble when it comes to defining my Cartesian nodes in the ModelGrib function.
</div>
<div> </div>
<div>I'll walk you through my code before that point... I pull the current day 0z GFS_hd run and then specify a given lat lon area along with certain variables to reduce the size of the data that I am reading in. From there, I run the DODSGrab function to pull
in 36 hours worth of forecast data (gfs_hd) which works fine. </div>
<div> </div>
<div>From there, my troubles begin with the ModelGrid function. I believe I have all the arguments correct other than Cartesian.nodes, but I am skeptical on the resolution... not completely sure if I did that correctly. When it comes to the Cartesian.nodes
function, I'm not exactly sure what to put in the argument. I understand that it has to be a list with an x and y element, but I can't find any literature on the gfs_hd grid in kilometers. Any thoughts on how to move forward here?</div>
<div> </div>
<div>Thank you. </div>
<div> </div>
<div> </div>
<div>Code is listed below. I am running on a windows 7 machine so pulling Grib data is not an option for me, unfortunately.
</div>
<div> </div>
<div>#############################################3</div>
<div> <br clear="none">
library("rNOMADS")<br clear="none">
library("RCurl")<br clear="none">
<br clear="none">
GFS.list <- GetDODSDates(abbrev="gfs_hd", archive=F, request.sleep=1) #returns list of available dates with GFS data in a numerical list<br clear="none">
<br clear="none">
<br clear="none">
current.GFS.url <- tail(GFS.list$url, 1)<br clear="none">
GFS.model.runs <- GetDODSModelRuns(current.GFS.url)<br clear="none">
<br clear="none">
GFS.0z <- (head(GFS.model.runs$model.run,1))<br clear="none">
<br clear="none">
#run when needed <br clear="none">
#GFS.info.variables <- GetDODSModelRunInfo(current.GFS.url, GFS.0z)<br clear="none">
<br clear="none">
lon <- c(480, 620) # USA Longitude<br clear="none">
lat <- c(240, 280) # USA Latitude<br clear="none">
<br clear="none">
variables <- c("tmax2m","tmin2m", "weasdsfc", "ugrdmwl", "vgrdmwl")<br clear="none">
<br clear="none">
GFS.Model.Data <- DODSGrab(current.GFS.url, GFS.0z, variables, time = c(0,36), lon, lat)</div>
<div style="color:rgb(0,0,0);font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:12px;font-style:normal;background-color:transparent;">
</div>
<div>GFS.modelgrid <- ModelGrid(model.data = GFS.Model.Data, grid.type = "cartesian", resolution = c(-140,40),levels = c("surface","2_m_above_ground"), variables, model.domain = NULL, cartesian.nodes = (list(x=0.5,y=0.5)))</div>
<div> </div>
<div>########################################################################</div>
<div> </div>
<div>Adam </div>
</div>
</div>
</div></div>
</div>
</td> </tr> </tbody> </table></td></tr></table>