<div dir="ltr">Hi, thanks, <div><br></div><div>My problem is more on the output side. I do addfiles, and then I'm reading in the variables with NCL, like this: </div><div><br></div>u = f[:]->UGRD_P0_L103_GLL0(:,{latNGFS:latSGFS},{lonWGFS:lonEGFS})<br>v = f[:]->VGRD_P0_L103_GLL0(:,{latNGFS:latSGFS},{lonWGFS:lonEGFS})<br>prate = f[:]->PRATE_P0_L1_GLL0(:,{latNGFS:latSGFS},{lonWGFS:lonEGFS})<br>t = f[:]->TMP_P0_L103_GLL0(:,{latNGFS:latSGFS},{lonWGFS:lonEGFS})<div><br></div><div>Currently, I'm passing this data to a Fortran wrapper, because it's always seemed easier than writing data in NCL. I loop through time, lat, lon to create concatenated strings for each variable (u,v,prate,t). Then, I loop through lat, and lon to write 1 line to a file for each latitude and longitude. So, the output lines will look something like this:</div>-77.5,0,"15.36;15.25,14.5...","354;350;349,...","273.3;274.5;275.8...","0;0;0;...","999;1002;1002;..."<div><br></div><div>The first 2 values represent latitude and longitude, and then the concatenated strings for wind speed, wind direction (converting u and v), prate, and t. There would be exactly latitude x longitude lines in the output file. So, again, just wondering if there is a more efficient way to write this output file.</div><div><div><br></div><div><br></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Sep 22, 2019 at 12:09 PM Dennis Shea <<a href="mailto:shea@ucar.edu">shea@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>As noted by Rashed:</div><div><br></div><div><a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/addfiles.shtml" target="_blank"><b>addfiles</b></a> also:<a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/setfileoption.shtml" target="_blank"><b> setfileoption</b></a></div><div>These work with netCDF, GRIB and HDF</div><div><br></div><div>Say the GFS files begin with: <b>gfs </b>and are GRIB [grib1 or grib2]<br></div><div> <br><div> <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/setfileoption.shtml" target="_blank"><strong>setfileoption</strong></a>("<b>grb</b>","SingleElementDimensions", (/"Initial_time","Forecast_time"/))</div></div><div> diri = "..."</div><div> fili = <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/systemfunc.shtml" target="_blank"><b>systemfunc</b></a>("cd "+diri+"<b> ; </b>ls gfs*)</div><div> fgfs = <b>addfiles</b>(diri+fili,"r")</div><div> x = fgfs->X</div><div> printVarSummary(x) ; ( initial_time0_hours, <b>forecast_time0,</b> lat, lon)</div><div><br></div><div>These are concatenated initial times and forecast times<br></div><div>or, if no file extension on the source files, add<br></div><div> <b> f</b>ili = fili + ".grb"</div><div> fgfs = <b>addfiles</b>(diri+fili,"r")</div><div><br></div><div>Good Luck</div><div>D<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Sep 22, 2019 at 9:39 AM Micah Sklut via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div dir="auto">Hi Rashed,</div></div><div dir="auto"><br></div><div dir="auto">I’m using addfiles to add all of the files, but I want to take the 384 hours of Data and make one array that has the hourly data concatenated for each latitude and longitude. </div><div dir="auto"><br></div><div dir="auto">As far as I can see, loops are the only way to go. I just need a lot of memory to pull it off.</div><div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Sep 22, 2019 at 8:02 AM Rashed Mahmood <<a href="mailto:rashidcomsis@gmail.com" target="_blank">rashidcomsis@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi Micah,</div><div>Did you try using the addfiles function in NCL:</div><div><br></div><div><a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/addfiles.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Functions/Built-in/addfiles.shtml</a></div><div><br></div><div>For example:</div><div>
<pre> fils = <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/systemfunc.shtml" target="_blank"><strong>systemfunc</strong></a> ("ls *.nc") ; file paths
f = <strong>addfiles</strong> (fils, "r") </pre>
</div><div><br></div><div>You would need to make sure that "ls *.nc" results in the correct order with respect to time and day.</div><div><br></div><div>Cheers,</div><div>Rashed<br></div><div><br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"></div><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Sep 21, 2019 at 1:53 PM Micah Sklut via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br></div></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"></blockquote></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi, <div><br></div><div>I have a task to take GFS data and concatenate variables across all the forecast hours, for each grid point. </div><div><br clear="all"><div>For example, taking a temperature variable for each grid point, for all latitude and longitudes, and creating a string that represents the values for each forecast hour. </div><div>So, if we have for latitude y, and longitude x, there would be a string value like "70.5,71.5,71.5,72.0,...nHours". The end product will be a line for each grid point for the variables i"m looking at that will be written to a file to be imported into a database. </div><div><br></div><div>Creating loops through all hours, latitudes, and longitudes will get the job done, but is expensive and was looking to see if there were any NCL functions that would help here. </div><div><br></div><div>Thank you. </div><div><br></div></div></div></blockquote></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
_______________________________________________<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" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>
</blockquote></div></div>-- <br><div dir="ltr">Micah Sklut<br><br></div>
_______________________________________________<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" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature">Micah Sklut<br><br></div>