<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<style type="text/css" style="display:none"><!--P{margin-top:0;margin-bottom:0;} p
{margin-top:0;
margin-bottom:0}--></style>
</head>
<body dir="ltr" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p><br>
</p>
<p><br>
</p>
<div id="Signature">
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
Daniel C. Bowman<br>
Doctoral Candidate in Geophysics<br>
Mitchell 315<br>
Department of Geological Sciences<br>
The University of North Carolina at Chapel Hill<br>
phone: 575-418-8555<br>
email: daniel.bowman@unc.edu<br>
web:http://geosci.unc.edu/page/daniel-c-bowman<br>
<br>
</div>
</div>
<div style="font-size:12pt; color:#000000; background-color:#FFFFFF; font-family:Calibri,Arial,Helvetica,sans-serif" dir="ltr">
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" color="#000000" face="Calibri, sans-serif"><b>From:</b> Daniel Bowman<br>
<b>Sent:</b> Monday, September 8, 2014 8:31 PM<br>
<b>To:</b> Adam Simkowski<br>
<b>Subject:</b> RE: [Rnomads-user] issues on retreveing specific location data in rNOMADS</font>
<div> </div>
</div>
<div>
<p>Hi Adam,</p>
<p><br>
</p>
<p>The issue here is the North American Mesoscale (NAM) model. The package vignettes were written specifically to interface with the Global Forecast System (GFS) model. The GFS is discretized into 0.5 x 0.5 degree cells, thus there are 720 longitude values
and 180 latitude values. That is the origin of the following lines in the vignette:</p>
<p><br>
</p>
<p>lons <- seq(0, 359.5, by = 0.5)<br>
lats <- seq(-90, 90, by = 0.5)<br>
</p>
<p><br>
</p>
<p>So...we need to figure out how the NAM model is set up before we can figure out how to get the data you want.</p>
<p><br>
</p>
<p>When I run </p>
<p><br>
</p>
<p>GetDODSModelRunInfo(model.url, NAM.model.run)</p>
<p><br>
</p>
<p>I find:</p>
<p><br>
</p>
<p>"Longitude: -153.00000000000°E to -49.30000000000°E (1038 points, avg. res. 0.1°)" <br>
[4] "Latitude: 12.00000000000°N to 61.90000000000°N (500 points, avg. res. 0.1°)"
</p>
<p><br>
</p>
<p>It is important to make sure each point is in fact 0.1 degrees apart (some models have variable grid spacings)...so I have to download the *entire model* and check:</p>
<p><br>
</p>
<p>everything <- DODSGrab(model.url, NAM.model.run, "dswrfsfc", time, c(0, 1037), c(0, 499))</p>
<p><br>
</p>
<p>diff(everything$lat)</p>
<p>diff(everything$lon)</p>
<p><br>
</p>
<p>Well, looks like we are lucky and everything is in fact 0.1 degrees apart. So let's redefine our indices:</p>
<p><br>
</p>
<p>lat <- 36.424789 #lat at Murphy Solar<br>
<br>
lon <- -77.062891 #lon at Murphy Solar<br>
<br>
#Get nearest model nodes<br>
lons <- seq(-153, -49.3, by = 0.1)<br>
lats <- seq(12, 61.9, by = 0.1)<br>
lon.diff <- abs(lon - lons)<br>
lat.diff <- abs(lat - lats)<br>
model.lon.ind <- which(lon.diff == min(lon.diff))<br>
model.lat.ind <- which(lat.diff == min(lat.diff))<br>
</p>
<p><br>
</p>
<p>Okay, now let's try getting the data for that specific point.</p>
<p><br>
</p>
<p>time <- c(0,0) #Model status at initialization<br>
lon.inds <- c(model.lon.ind - 2, model.lon.ind + 2)<br>
lat.inds <- c(model.lat.ind - 2, model.lat.ind + 2)<br>
<br>
</p>
<p>solar.data <- DODSGrab(model.url, NAM.model.run, "dswrfsfc", time, lon.inds, lat.inds)# pull surface downward short wave flux [w/m^2] from the model
</p>
<p><br>
</p>
<p>According to the coordinates in solar.data, we've captured the point we want. Anyway, the take home message is that each *model* is unique, and you have to investigate each one individually.</p>
<p><br>
</p>
<p>Hope this helps!<br>
<br>
Danny<br>
<br>
</p>
<p><br>
</p>
<div id="Signature">
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
Daniel C. Bowman<br>
Doctoral Candidate in Geophysics<br>
Mitchell 315<br>
Department of Geological Sciences<br>
The University of North Carolina at Chapel Hill<br>
phone: 575-418-8555<br>
email: daniel.bowman@unc.edu<br>
web:http://geosci.unc.edu/page/daniel-c-bowman<br>
<br>
</div>
</div>
<div style="color:rgb(33,33,33)">
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><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>
<b>Sent:</b> Monday, September 8, 2014 3:31 PM<br>
<b>To:</b> rnomads-user@lists.r-forge.r-project.org<br>
<b>Subject:</b> [Rnomads-user] issues on retreveing specific location data in rNOMADS</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:12pt">
<div>Hi all, </div>
<div> </div>
<div><font face="Times New Roman"></font></div>
<div style="margin:0in 0in 0pt"><span style="color:rgb(33,33,33)"><font face="Calibri">I have been attempting to pull “dswrfsfc” (incoming solar radiation) at a specific latitude and longitude from the historical NAM data set available through rNOMADS. I followed
the example shown in various parts of the package notes where data is pulled for Chapel Hill, NC. I am attempting to pull data for a location not too far from Chapel Hill (Murfreesboro, NC). However, when I attempt to pull either solar radiation or 2 meter
temperature data for a given hour, I get 24 rows of data with a longitude far away from my desired location and a latitude somewhat nearby. I was under the impression that I should only get one row of data closest to my desired location since I’m only querying
for one specific hour. I am not completely sure I understand the code to query data closest to a desired location given a latitude and longitude... I have included this code in the attached items with my specific lat and lon ( file name: rNOMADS lat-lon code).
Could someone please explain code here a little more in-depth than is shown in on the package vignette? </font></span><font face="Times New Roman"></font></div>
<span style="color:rgb(33,33,33)"></span>
<div><font face="Times New Roman"></font> </div>
<div style="margin:0in 0in 0pt"><span style="color:rgb(33,33,33)"><font face="Calibri">I have also attached my script that I wrote to retrieve point specific model data for Murfreesboro (Murphy solar) and an example of the 24 row data frame that it creates.
The script was written as a function where the input is a specific number assigned to a url from the historical NAM list. An example input would be 1473 which is the url for NAM runs on September 1<sup><font size="2">st</font></sup>, 2014. I am hoping to fix
the code so that it only returns one row of data instead of 24. </font></span></div>
<div><font face="Times New Roman"></font></div>
<span style="color:rgb(33,33,33)">
<div style="margin:0in 0in 0pt"><font face="Calibri">I have also written a script that pulls more than one model at a time using a for loop, but this lat lon issue seems to be inhibiting me from creating a proper final output data frame with good data. Once
I figure out the lat lon issue here, I will certainly share this code for all to use!</font></div>
<div style="margin:0in 0in 0pt; color:rgb(33,33,33); font-family:Calibri; font-size:16px; font-style:normal; background-color:transparent">
<font face="Calibri"></font> </div>
<div style="margin:0in 0in 0pt; color:rgb(33,33,33); font-family:Calibri; font-size:16px; font-style:normal; background-color:transparent">
</div>
<div style="margin:0in 0in 0pt; color:rgb(33,33,33); font-family:Calibri; font-size:16px; font-style:normal; background-color:transparent">
Thank you, </div>
<div style="margin:0in 0in 0pt; color:rgb(33,33,33); font-family:Calibri; font-size:16px; font-style:normal; background-color:transparent">
</div>
<div style="margin:0in 0in 0pt; color:rgb(33,33,33); font-family:Calibri; font-size:16px; font-style:normal; background-color:transparent">
Adam Simkowski </div>
<div style="margin:0in 0in 0pt; color:rgb(33,33,33); font-family:Calibri; font-size:16px; font-style:normal; background-color:transparent">
</div>
</span>
<div style="margin:0in 0in 0pt; color:rgb(33,33,33); font-family:Calibri; font-size:16px; font-style:normal; background-color:transparent">
</div>
<div><font face="Times New Roman"></font> </div>
<div style="margin:0in 0in 0pt"><span style="color:rgb(33,33,33)"><font face="Calibri"> </font></span></div>
<div><font face="Times New Roman"></font></div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>