<div dir="ltr"><div class="gmail_default" style="font-size:small">Chris,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Yes, the "all(ismissing(sst))" is a way to trap this special case.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I'm not sure what you are asking in your second question. When you subset "sst" with the given lat/lon range, then the lat/lon coordinate arrays should also be subsetted automatically.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I just tried to run your script again to verify that you indeed have coordinate arrays attached to your data, but it appears the file is no longer there.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">You can verify this yourself by doing a "printVarSummary" on "sst" after you read it off a file.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Here's a little sample script that you should be able to run as-is:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style=""><div class="gmail_default" style=""><font face="monospace, monospace"> a = addfile("$NCARG_ROOT/lib/ncarg/data/cdf/<a href="http://uv300.nc">uv300.nc</a>","r")</font></div><div class="gmail_default" style=""><font face="monospace, monospace"> u = a->U(1,:,:)</font></div><div class="gmail_default" style=""><font face="monospace, monospace"><br></font></div><div class="gmail_default" style=""><font face="monospace, monospace"> printVarSummary(u)</font></div><div class="gmail_default" style=""><font face="monospace, monospace"><br></font></div><div class="gmail_default" style=""><font face="monospace, monospace"> u_subset = a->U(1,{0:30},{-10:10})</font></div><div class="gmail_default" style=""><font face="monospace, monospace"> printVarSummary(u_subset)</font></div><div class="gmail_default" style=""><br></div><div class="gmail_default" style="">The output (trimmed) will look like this:</div><div class="gmail_default" style=""><br></div><div class="gmail_default" style=""><div class="gmail_default"><font face="monospace, monospace">Variable: u</font></div><div class="gmail_default"><font face="monospace, monospace">Type: float</font></div><div class="gmail_default"><font face="monospace, monospace">. . .</font></div><div class="gmail_default"><span style="font-family:monospace,monospace">Dimensions and sizes:</span><span class="" style="font-family:monospace,monospace;white-space:pre"> </span><span style="font-family:monospace,monospace">[lat | 64] x [lon | 128]</span><br></div><div class="gmail_default"><font face="monospace, monospace">Coordinates: </font></div><div class="gmail_default"><font face="monospace, monospace"> lat: [-87.8638..87.8638]</font></div><div class="gmail_default"><font face="monospace, monospace"> lon: [-180..177.1875]</font></div><div class="gmail_default"><font face="monospace, monospace">Number Of Attributes: 5</font></div><div class="gmail_default"><font face="monospace, monospace"> time :<span class="" style="white-space:pre"> </span>7</font></div><div class="gmail_default"><font face="monospace, monospace"> _FillValue :<span class="" style="white-space:pre"> </span>-999</font></div><div class="gmail_default"><font face="monospace, monospace"> long_name :<span class="" style="white-space:pre"> </span>Zonal Wind</font></div><div class="gmail_default"><font face="monospace, monospace"> short_name :<span class="" style="white-space:pre"> </span>U</font></div><div class="gmail_default"><font face="monospace, monospace"> units :<span class="" style="white-space:pre"> </span>m/s</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">Variable: u_subset</font></div><div class="gmail_default"><font face="monospace, monospace">Type: float</font></div><div class="gmail_default"><font face="monospace, monospace">. . .</font></div><div class="gmail_default"><font face="monospace, monospace">Dimensions and sizes:<span class="" style="white-space:pre"> </span>[lat | 11] x [lon | 7]</font></div><div class="gmail_default"><font face="monospace, monospace">Coordinates: </font></div><div class="gmail_default"><font face="monospace, monospace"> lat: [1.395307..29.30136]</font></div><div class="gmail_default"><font face="monospace, monospace"> lon: [-8.4375..8.4375]</font></div><div class="gmail_default"><font face="monospace, monospace">Number Of Attributes: 5</font></div><div class="gmail_default"><font face="monospace, monospace"> time :<span class="" style="white-space:pre"> </span>7</font></div><div class="gmail_default"><font face="monospace, monospace"> _FillValue :<span class="" style="white-space:pre"> </span>-999</font></div><div class="gmail_default"><font face="monospace, monospace"> long_name :<span class="" style="white-space:pre"> </span>Zonal Wind</font></div><div class="gmail_default"><font face="monospace, monospace"> short_name :<span class="" style="white-space:pre"> </span>U</font></div><div class="gmail_default"><font face="monospace, monospace"> units :<span class="" style="white-space:pre"> </span>m/s</font></div><div class="gmail_default"><br></div></div></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Note that for "u_subset", the lat/lon coordinates have the new subsetted range, based on the closest *less than* values to the lat/lon range I gave it.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_extra"><div class="gmail_quote">On Tue, May 17, 2016 at 5:54 PM, Herbster, Christopher G. <span dir="ltr"><<a href="mailto:herbstec@erau.edu" target="_blank">herbstec@erau.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">
<div lang="EN-US" link="blue" vlink="purple">
<div>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Mary,<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Thanks for checking what seems obvious in hindsight. Feeling a bit sheepish.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Can I trap this with an:<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">if(all(ismissing(sst))) then
<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"> … ; don’t do much<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">else<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">… ; make plots<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">end if<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Boy have I wasted some time thinking this had to do with the longitude being 74 – 434 degrees east. I had originally thought that I was somehow only getting
data over land, where SST would certainly be missing.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Also, is there a preferred way to have the meta data for the subset data have ranges that are indicative of the subset, rather than the whole domain?<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Thanks again for catching my oversight,<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Chris H.<u></u><u></u></span></p><span class="">
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">--
<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Dr. Christopher G. Herbster<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Associate Professor<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Director of Science and Technology<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">for the ERAU Weather Center<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Applied Aviation Sciences<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Embry-Riddle Aeronautical Univ.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">600 S. Clyde Morris Blvd.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Daytona Beach, FL 32114-3900<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"> <a href="tel:386.226.6444" value="+13862266444" target="_blank">386.226.6444</a> Office<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><a href="tel:386.226.6446" value="+13862266446" target="_blank">386.226.6446</a> Weather Center<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><a href="http://wx.erau.edu/" target="_blank">http://wx.erau.edu/</a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Schedule at: <a href="http://wx.erau.edu/faculty/herbster/Schedules/" target="_blank">http://wx.erau.edu/faculty/herbster/Schedules/</a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
</span><p class="MsoNormal"><b><span style="font-size:11pt;font-family:Calibri,sans-serif">From:</span></b><span style="font-size:11pt;font-family:Calibri,sans-serif"> Mary Haley [mailto:<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>]
<br>
<b>Sent:</b> Tuesday, May 17, 2016 5:58 PM<br>
<b>To:</b> Herbster, Christopher G. <<a href="mailto:herbstec@erau.edu" target="_blank">herbstec@erau.edu</a>><br>
<b>Cc:</b> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
<b>Subject:</b> Re: [ncl-talk] Problem reading data from NCEP DODS server<u></u><u></u></span></p><div><div class="h5">
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<div>
<p class="MsoNormal">Chris,<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">I think the issue is simply that the first time step of "sst" is all missing. <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">If you select a different time step, then it looks fine:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">; sst = f->sst(mytime,:,Slat:Nlat,Wlon:Elon) </span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'"> </span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">sst = f->sst(:,:,Slat:Nlat,Wlon:Elon) ; NOTE: READING ALL TIME STEPS</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">print("done reading sst")</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">print("dimsizes sst = " + str_join(dimsizes(sst),","))</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">do nt=0,dimsizes(sst(:,0,0,0))-1</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'"> printMinMax(sst(nt,:,:,:),0)</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">end do</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">This is the output:</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">(0) dimsizes sst = 65,1,265,481<u></u><u></u></span></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">nt = 0<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">sea_surface_temperature (c) : min=1.26765e+30 max=1.26765e+30<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">nt = 1<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">sea_surface_temperature (c) : min=21.6418 max=32.2542<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">nt = 2<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">sea_surface_temperature (c) : min=21.6017 max=32.3401<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">nt = 3<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">sea_surface_temperature (c) : min=21.5984 max=32.3429<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">nt = 4<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">sea_surface_temperature (c) : min=21.6628 max=32.1836<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">nt = 5<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">sea_surface_temperature (c) : min=22.0586 max=32.491<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">nt = 6<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">sea_surface_temperature (c) : min=22.2371 max=32.6994<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'">--Mary<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:'Courier New'"><u></u> <u></u></span></p>
</div>
</div>
</div>
</div>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">On Tue, May 17, 2016 at 10:28 AM, Herbster, Christopher G. <<a href="mailto:herbstec@erau.edu" target="_blank">herbstec@erau.edu</a>> wrote:<u></u><u></u></p>
<blockquote style="border-style:none none none solid;border-left-width:1pt;border-left-color:rgb(204,204,204);padding:0in 0in 0in 6pt;margin-left:4.8pt;margin-right:0in">
<div>
<div>
<p class="MsoNormal">Hi folks,<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">I’m trying to retrieve some data for a subset region from the Global Real Time Ocean Forecast System (RTOFS), and I can’t figure out why I get all missing values when I try to read
the data.<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">This seems to be the case for me even if I try to read this whole domain.<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">Here is some code:<u></u><u></u></p>
<p class="MsoNormal">------<u></u><u></u></p>
<p class="MsoNormal">; read-rtofs-dods.ncl
<u></u><u></u></p>
<p class="MsoNormal">load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/cd_string.ncl"<u></u><u></u></p>
<p class="MsoNormal">;
<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">Dev = "png"<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">;<u></u><u></u></p>
<p class="MsoNormal">myDate = 20160516<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">;NOMADS = "<a href="http://nomads.ncep.noaa.gov:9090/dods/rtofs/rtofs_global" target="_blank">http://nomads.ncep.noaa.gov:9090/dods/rtofs/rtofs_global</a>"+myDate+"/rtofs_glo_2ds_forecast_daily_prog"<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">NOMADS = "<a href="http://nomads.ncep.noaa.gov:9090/dods/rtofs/rtofs_global" target="_blank">http://nomads.ncep.noaa.gov:9090/dods/rtofs/rtofs_global</a>"+myDate+"/rtofs_glo_2ds_forecast_3hrly_prog"<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">filename = NOMADS ; = url
<u></u><u></u></p>
<p class="MsoNormal">print("Fetching URL: ")<u></u><u></u></p>
<p class="MsoNormal">print(" "+filename)<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal"> exists = isfilepresent(filename)<u></u><u></u></p>
<p class="MsoNormal"> if(.not.exists) then
<u></u><u></u></p>
<p class="MsoNormal"> print("OPeNDAP isfilepresent test unsuccessful.")<u></u><u></u></p>
<p class="MsoNormal"> print("Either file doesn't exist, or NCL does not have OPeNDAP capabilities on this system")<u></u><u></u></p>
<p class="MsoNormal"> exit<u></u><u></u></p>
<p class="MsoNormal"> else<u></u><u></u></p>
<p class="MsoNormal"> print("OPeNDAP isfilepresent test successful.")<u></u><u></u></p>
<p class="MsoNormal"> f = addfile(filename,"r")<u></u><u></u></p>
<p class="MsoNormal"> vars = getfilevarnames(f)<u></u><u></u></p>
<p class="MsoNormal"> print(vars) ; should be (in any order):
<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal"> if(.not.any(ismissing(vars))) then<u></u><u></u></p>
<p class="MsoNormal"> do i=0,dimsizes(vars)-1<u></u><u></u></p>
<p class="MsoNormal"> printFileVarSummary (f,vars(i))<u></u><u></u></p>
<p class="MsoNormal"> end do<u></u><u></u></p>
<p class="MsoNormal"> end if<u></u><u></u></p>
<p class="MsoNormal"> end if<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">; Define boundary of data<u></u><u></u></p>
<p class="MsoNormal">; Assign lat/lon boundaries to the data<u></u><u></u></p>
<p class="MsoNormal">; Longitudes between 0 and 74 E need to have 360 added to them<u></u><u></u></p>
<p class="MsoNormal">;<u></u><u></u></p>
<p class="MsoNormal">myNlat = 32.0<u></u><u></u></p>
<p class="MsoNormal">mySlat = 10.0<u></u><u></u></p>
<p class="MsoNormal">myWlon = 260.0<u></u><u></u></p>
<p class="MsoNormal">myElon = 300.0<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">;myNlat = 26.0<u></u><u></u></p>
<p class="MsoNormal">;mySlat = 21.0<u></u><u></u></p>
<p class="MsoNormal">;myWlon = 260.0<u></u><u></u></p>
<p class="MsoNormal">;myElon = 300.0<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">; Full domain
<u></u><u></u></p>
<p class="MsoNormal">;myNlat = 90.0<u></u><u></u></p>
<p class="MsoNormal">;mySlat = -90.0<u></u><u></u></p>
<p class="MsoNormal">;myWlon = 74.0<u></u><u></u></p>
<p class="MsoNormal">;myElon = 434.0<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">; Get lat/lon/time info first<u></u><u></u></p>
<p class="MsoNormal">lat = f->lat<u></u><u></u></p>
<p class="MsoNormal">print ("Done reading lat")<u></u><u></u></p>
<p class="MsoNormal">lon = f->lon<u></u><u></u></p>
<p class="MsoNormal">print ("Done reading lon")<u></u><u></u></p>
<p class="MsoNormal">time = f->time(0)<u></u><u></u></p>
<p class="MsoNormal">print ("Done reading time")<u></u><u></u></p>
<p class="MsoNormal">;<u></u><u></u></p>
<p class="MsoNormal">; Use just one time when "0"<u></u><u></u></p>
<p class="MsoNormal">mytime = 0<u></u><u></u></p>
<p class="MsoNormal">;mytime = ":"<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">Nlat = closest_val(myNlat, lat)<u></u><u></u></p>
<p class="MsoNormal">print ("Nlat = " + Nlat + " Lat = " + lat(Nlat) )<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">Slat = closest_val(mySlat, lat)<u></u><u></u></p>
<p class="MsoNormal">print ("Slat = " + Slat + " Lat = " + lat(Slat) )<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">Wlon = closest_val(myWlon, lon)<u></u><u></u></p>
<p class="MsoNormal">print("Wlon = " + Wlon + " Lon = " + lon(Wlon) )<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">Elon = closest_val(myElon, lon)<u></u><u></u></p>
<p class="MsoNormal">print("Elon = " + Elon + " Lon = " + lon(Elon) )<u></u><u></u></p>
<p class="MsoNormal">;<u></u><u></u></p>
<p class="MsoNormal">;exit<u></u><u></u></p>
<p class="MsoNormal">;<u></u><u></u></p>
<p class="MsoNormal">mylat = lat(Slat:Nlat)<u></u><u></u></p>
<p class="MsoNormal">mylon = lon(Wlon:Elon)<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">printVarSummary (mylat)<u></u><u></u></p>
<p class="MsoNormal">printVarSummary (mylon)<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">print (mylat)<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">;exit<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">delete(lat)<u></u><u></u></p>
<p class="MsoNormal">delete(lon)<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">lat = f->lat(Slat:Nlat)<u></u><u></u></p>
<p class="MsoNormal">lon = f->lon(Wlon:Elon)<u></u><u></u></p>
<p class="MsoNormal">print ("Done reading subset lat/lon")<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">printVarSummary(lat)<u></u><u></u></p>
<p class="MsoNormal">print (lat)<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">printVarSummary(lon)<u></u><u></u></p>
<p class="MsoNormal">print (lon)<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">;exit<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">sst = f->sst(mytime,:,Slat:Nlat,Wlon:Elon)<u></u><u></u></p>
<p class="MsoNormal">;sst = f->sst<u></u><u></u></p>
<p class="MsoNormal">print ("Done reading sst")<u></u><u></u></p>
<p class="MsoNormal">printVarSummary(sst)<u></u><u></u></p>
<p class="MsoNormal">print (sst)<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">; exit ; Uncomment to stop after sst data are read
<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">u_velocity = f->u_velocity(mytime,:,Slat:Nlat,Wlon:Elon)<u></u><u></u></p>
<p class="MsoNormal">;u_velocity = f->u_velocity<u></u><u></u></p>
<p class="MsoNormal">print ("Done reading U")<u></u><u></u></p>
<p class="MsoNormal">printVarSummary(u_velocity)<u></u><u></u></p>
<p class="MsoNormal">print (u_velocity)<u></u><u></u></p>
<p class="MsoNormal">v_velocity = f->v_velocity(mytime,:,Slat:Nlat,Wlon:Elon)<u></u><u></u></p>
<p class="MsoNormal">;v_velocity = f->v_velocity<u></u><u></u></p>
<p class="MsoNormal">print ("Done reading V")<u></u><u></u></p>
<p class="MsoNormal">printVarSummary(v_velocity)<u></u><u></u></p>
<p class="MsoNormal">print (v_velocity)<u></u><u></u></p>
<p class="MsoNormal">speed = sqrt((u_velocity*u_velocity) + (v_velocity*v_velocity))<u></u><u></u></p>
<p class="MsoNormal">print ("Done calculating Speed")<u></u><u></u></p>
<p class="MsoNormal">print (speed)<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">exit<u></u><u></u></p>
<p class="MsoNormal">;;;;;;<u></u><u></u></p>
<p class="MsoNormal">;;;;;; END of code snippet<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">------<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">If you run this code, I suggest a redirect into a file, as there are a lot of lines written to standard out.<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">What I don’t understand is that I can read a subset of the lat and lon data, but when I try to get any other data they all come up as missing/fill values.<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">Any help is GREATLY appreciated and I have really chased my tail trying to get past this.<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">Thanks,<u></u><u></u></p>
<p class="MsoNormal">Chris Herbster<u></u><u></u></p>
<p class="MsoNormal">--
<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">Dr. Christopher G. Herbster<u></u><u></u></p>
<p class="MsoNormal">Associate Professor<u></u><u></u></p>
<p class="MsoNormal">Director of Science and Technology<u></u><u></u></p>
<p class="MsoNormal">for the ERAU Weather Center<u></u><u></u></p>
<p class="MsoNormal">Applied Aviation Sciences<u></u><u></u></p>
<p class="MsoNormal">Embry-Riddle Aeronautical Univ.<u></u><u></u></p>
<p class="MsoNormal">600 S. Clyde Morris Blvd.<u></u><u></u></p>
<p class="MsoNormal">Daytona Beach, FL 32114-3900<u></u><u></u></p>
<p class="MsoNormal"> <a href="tel:386.226.6444" target="_blank">386.226.6444</a> Office<u></u><u></u></p>
<p class="MsoNormal"><a href="tel:386.226.6446" target="_blank">386.226.6446</a> Weather Center<u></u><u></u></p>
<p class="MsoNormal"><a href="http://wx.erau.edu/" target="_blank">http://wx.erau.edu/</a><u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">Schedule at:
<a href="http://wx.erau.edu/faculty/herbster/Schedules/" target="_blank">http://wx.erau.edu/faculty/herbster/Schedules/</a><u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
<p class="MsoNormal" style="margin-bottom:12pt"><br>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><u></u><u></u></p>
</blockquote>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div></div></div>
</div>
</blockquote></div><br></div></div>