<br>Hi Wesley,<br><br>Apologies for not being more clear. The example below would be a way of using lats4d on reanalysis data behind a GDS.<br><br>lats4d -dods -i <a href="http://nomads.ncdc.noaa.gov:9091/dods/NCEP_GFS_ANALYSIS/analysis">http://nomads.ncdc.noaa.gov:9091/dods/NCEP_GFS_ANALYSIS/analysis</a><br>
_complete -o test1 -format grads_grib -lat 15 50 -lon -100 -70 -time 00Z01SEP2004 18Z01SEP2004 -vars tmpsfc tmp2m<br><br>The result is a grib file for 6 hourly subsetted temperatures.<br><br>It&#39;s possible that this works for me because of the local grads installation. The q config from my workstation is appended below.<br>
<br>The NARR grid throws a curve, some things may work, for example I&#39;m not sure if regrid works on NARR. I put the following data set into the same command as above:<br><br><a href="http://nomads.ncdc.noaa.gov:9091/dods/NCEP_NARR_DAILY/200409/200409/narr-a_221_200409dd_hh00_000">http://nomads.ncdc.noaa.gov:9091/dods/NCEP_NARR_DAILY/200409/200409/narr-a_221_200409dd_hh00_000</a><br>
<br>I don&#39;t recall encountering grib2 data yet, so I cannot comment on that. We like serving up HDF data which is likely the anti-fun compared to grib2.<br><br>Mike<br><br><br>Config: OpenGrADS-1.9.8-m3 32-bit little-endian readline sdf/xdf hdf-sds netcdf lats athena printim<br>
<br>Grid Analysis and Display System (GrADS) Version OpenGrADS-1.9.8-m3<br>Copyright (c) 1988-2005 by Brian Doty and IGES<br>Center for Ocean-Land-Atmosphere Studies (COLA)<br>Institute for Global Environment and Society (IGES)<br>
This program is distributed WITHOUT ANY WARRANTY<br>See file COPYRIGHT for more information.<br><br>Built Fri Jul&nbsp; 6 00:20:28 EDT 2007 for i686-pc-linux-gnu<br><br>This version of GrADS has been configured with the following options:<br>
<br>&nbsp;&nbsp; o This is a 32-bit LITTLE ENDIAN machine version.<br>&nbsp;&nbsp; o Command line editing (readline) ENABLED.<br>&nbsp;&nbsp; o CIRES/CDC (<a href="http://www.cdc.noaa.gov">http://www.cdc.noaa.gov</a>) SDF/XDF interface ENABLED.<br>&nbsp;&nbsp;&nbsp;&nbsp; Use sdfopen or xdfopen to read both NetCDF and HDF-SDS files.<br>
&nbsp;&nbsp; o DTYPE netcdf and DTYPE hdfsds are ENABLED.<br>&nbsp;&nbsp; o OPeNDAP (a.k.a. DODS) gridded data interface DISABLED.<br>&nbsp;&nbsp; o OPeNDAP (a.k.a. DODS) station data interface DISABLED.<br>&nbsp;&nbsp; o PCMDI (<a href="http://www-pcmdi.llnl.gov">http://www-pcmdi.llnl.gov</a>) LATS interface ENABLED.<br>
&nbsp;&nbsp;&nbsp;&nbsp; This version is configured to write GRIB and HDF-SDS files.<br>&nbsp;&nbsp; o DAO (<a href="http://dao.gsfc.nasa.gov">http://dao.gsfc.nasa.gov</a>) Athena Widget GUI ENABLED.<br>&nbsp;&nbsp; o NRL/DAO/PCMDI XA or ImageMagick Image Output DISABLED.<br>
&nbsp;&nbsp; o printim command for direct png/gif output ENABLED.<br>&nbsp;&nbsp;&nbsp;&nbsp; (via the GD Library -- <a href="http://www.boutell.com/gd">http://www.boutell.com/gd</a>)<br><br>For additional information please consult <a href="http://grads.iges.org/grads/">http://grads.iges.org/grads/</a><br>
<br><br><br><div class="gmail_quote">On Jan 30, 2008 11:10 AM, Wesley Ebisuzaki &lt;<a href="mailto:Wesley.Ebisuzaki@noaa.gov">Wesley.Ebisuzaki@noaa.gov</a>&gt; wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Mike,<br><br> &nbsp; &nbsp; &nbsp;I hate to disappoint you but the NOMADS sites (NCEP, NCDC) don&#39;t<br>use lats4d for the distribution of NARR, global reanalyses and<br>other model products. &nbsp;With the recent release of GrADS v2.0a, all the<br>
building blocks are in place for NOMADS to support grib-2.<br>For many users, grib-2 will be more fun than grib-1.<br><font color="#888888"><br> &nbsp; &nbsp;Wesley Ebisuzaki<br></font><div><div></div><div class="Wj3C7c"><br><br>Mike Bosilovich wrote:<br>
&gt; Thanks for the prompt reply Brian,<br>&gt;<br>&gt; I realize and greatly appreciate the amount of work you, Jennifer and<br>&gt; others have done, with very little support. Grads is a great tool,<br>&gt; hands down. And I tell that to anyone who will listen.<br>
&gt;<br>&gt; I assumed that you all, as developers, would like feedback on what<br>&gt; users need to do. Lats4d fills a very important role in my work, and I<br>&gt; expect others as well. I wanted to pass on a point of view. Certainly,<br>
&gt; I&#39;m not trying to tell you what you must do, only providing feedback<br>&gt; on the capability that I find useful.<br>&gt;<br>&gt; Now, as for what Jennifer said, there was a post that says lats will<br>&gt; not be supported. That was the impetus for my post. If I<br>
&gt; misunderstood, and there will be lats4d capability in some way, I<br>&gt; apologize for over reacting.<br>&gt;<br>&gt; In any event, I also need to clarify my statement on the new NCEP<br>&gt; reanalysis. They will have 1000TB of Grib1 format data at NCDC (this<br>
&gt; is still changing, as I understand). In Grib2, it comes down to 500Tb<br>&gt; or so. For what that&#39;s worth, it&#39;s still immense. I got the sense that<br>&gt; they are planning on grib1, because of the lack of utilities with<br>
&gt; grib2 capability.<br>&gt;<br>&gt; Again, this is merely my opinion as a user, and I hoped the opinion<br>&gt; and the information would be useful to the developers.<br>&gt;<br>&gt; Sincerely,<br>&gt;<br>&gt; Mike Bosilovich<br>
&gt;<br>&gt; On Jan 29, 2008 10:02 PM, Brian Doty &lt;<a href="mailto:doty@cola.iges.org">doty@cola.iges.org</a><br></div></div><div><div></div><div class="Wj3C7c">&gt; &lt;mailto:<a href="mailto:doty@cola.iges.org">doty@cola.iges.org</a>&gt;&gt; wrote:<br>
&gt;<br>&gt; &nbsp; &nbsp; Mike, the new Version 2 of GrADS is being released in alpha version<br>&gt; &nbsp; &nbsp; and incomplete so that people would have access to the new GRIB2<br>&gt; &nbsp; &nbsp; support. &nbsp; With no funding for many years and key personnel working<br>
&gt; &nbsp; &nbsp; limited part time hours we do not have the resources to support old<br>&gt; &nbsp; &nbsp; code and also move forward to provide essential new functionality for<br>&gt; &nbsp; &nbsp; present and future research needs. &nbsp;GrADS development first started<br>
&gt; &nbsp; &nbsp; in the late 1980s and it has taken us several years of work to clean<br>&gt; &nbsp; &nbsp; up the 15 years of accumulated mods and patches on top of the<br>&gt; &nbsp; &nbsp; original code base in order to come out with a new code base so we<br>
&gt; &nbsp; &nbsp; could move forward. &nbsp; I have made every effort over the years to<br>&gt; &nbsp; &nbsp; ensure upward compatibility and maintain stability. &nbsp;There is no way<br>&gt; &nbsp; &nbsp; that a 20 year old code base can be patched and hacked for another 20<br>
&gt; &nbsp; &nbsp; years or even 5 years and maintain stability and still support the<br>&gt; &nbsp; &nbsp; 100TB data sets you talk about and the other projects probably in the<br>&gt; &nbsp; &nbsp; pipeline. &nbsp; As new versions of data formats appear -- grib2, hdf5,<br>
&gt; &nbsp; &nbsp; new netcdf on the way, etc --- &nbsp;we have no choice. &nbsp;We cannot keep<br>&gt; &nbsp; &nbsp; using old libraries and maintain them ourselves. &nbsp;We cannot patch and<br>&gt; &nbsp; &nbsp; hack old code that we didn&#39;t write to begin with and which in some<br>
&gt; &nbsp; &nbsp; cases has restrictive copyrights.<br>&gt;<br>&gt; &nbsp; &nbsp; As I say, version 2 as it stands now is incomplete. &nbsp; Key areas of<br>&gt; &nbsp; &nbsp; functionality which need to be re-implemented have not yet been<br>&gt; &nbsp; &nbsp; done. &nbsp;Other major parts of grads have not yet been upgraded, such as<br>
&gt; &nbsp; &nbsp; the graphics package. &nbsp; As Jennifer said in her original email, the<br>&gt; &nbsp; &nbsp; output of various data formats is high on the list of things yet to<br>&gt; &nbsp; &nbsp; be done... &nbsp; Brian<br>&gt;<br>&gt; &nbsp; &nbsp; On Jan 29, 2008, at 4:34 PM, Mike Bosilovich wrote:<br>
&gt;<br>&gt; &nbsp; &nbsp; &gt;<br>&gt; &nbsp; &nbsp; &gt; I have to admit an ignorance to some of the inner workings of Grads<br>&gt; &nbsp; &nbsp; &gt; development. But I am quite surprised to hear that lats4d is not<br>&gt; &nbsp; &nbsp; &gt; presently supported in Grads 2.<br>
&gt; &nbsp; &nbsp; &gt;<br>&gt; &nbsp; &nbsp; &gt; I have been using Grads for 10 years now (I can still recall the<br>&gt; &nbsp; &nbsp; &gt; relief after years of ncargf77 programming :-) Of course, I still<br>&gt; &nbsp; &nbsp; &gt; write code when appropriate, but grads and lats4d have been the<br>
&gt; &nbsp; &nbsp; &gt; main tools in my work. Lats4d fills a critical void. By pointing it<br>&gt; &nbsp; &nbsp; &gt; at any grads formatted file (netcdf, hdf, binary or grib), I can<br>&gt; &nbsp; &nbsp; &gt; reformat the data to what ever a colleague may need. I can also,<br>
&gt; &nbsp; &nbsp; &gt; easily, no, effortlessly subset variables space or time. Coupled<br>&gt; &nbsp; &nbsp; &gt; with a call to regrid, this subsetting utility is beyond compare.<br>&gt; &nbsp; &nbsp; &gt; In a c shell, it can rip through huge data files with simplicity. I<br>
&gt; &nbsp; &nbsp; &gt; hope this does not sound like exaggeration, lats4d is the strongest<br>&gt; &nbsp; &nbsp; &gt; data tool I have used.<br>&gt; &nbsp; &nbsp; &gt;<br>&gt; &nbsp; &nbsp; &gt; In our office, we have just begun production of a new reanalysis<br>
&gt; &nbsp; &nbsp; &gt; data product. While it will take some time to complete, we are<br>&gt; &nbsp; &nbsp; &gt; beginning to develop examples on how users can access and analyze<br>&gt; &nbsp; &nbsp; &gt; the data. Comparison with other reanalyses is obvious and will be<br>
&gt; &nbsp; &nbsp; &gt; in high demand. To accomplish that, they will need to regrid our<br>&gt; &nbsp; &nbsp; &gt; data sets to the existing coarser reanalysis data sets. Or, they<br>&gt; &nbsp; &nbsp; &gt; may need to change the format out of our native HDF. The easiest<br>
&gt; &nbsp; &nbsp; &gt; way to explain to others how to do this is with lats4d.<br>&gt; &nbsp; &nbsp; &gt;<br>&gt; &nbsp; &nbsp; &gt; There will be 100Tb of reanalysis data available through a GDS.<br>&gt; &nbsp; &nbsp; &gt; With Lats4d and gradsdods (or gradsdap), this would be much more<br>
&gt; &nbsp; &nbsp; &gt; accessible. Users, with some examples, will access the data through<br>&gt; &nbsp; &nbsp; &gt; online capabilities, rather than bulk downloading of the native HDF<br>&gt; &nbsp; &nbsp; &gt; files (a throttle may be needed for access if that is the preferred<br>
&gt; &nbsp; &nbsp; &gt; by users). In addition, a plan is being prepared to develop a DVD,<br>&gt; &nbsp; &nbsp; &gt; similar to the NCEP reanalyses CD and NARR DVD. Personally, I would<br>&gt; &nbsp; &nbsp; &gt; like to see a flavor of grads and lats4d included therein<br>
&gt; &nbsp; &nbsp; &gt; (admittedly, I have not gotten to discuss this with the grads<br>&gt; &nbsp; &nbsp; &gt; developers yet).<br>&gt; &nbsp; &nbsp; &gt;<br>&gt; &nbsp; &nbsp; &gt; It doesn&#39;t stop there. At AMS last week, NCEP and NCDC held a town<br>&gt; &nbsp; &nbsp; &gt; hall meeting to discuss their plans for the next reanalyses. They<br>
&gt; &nbsp; &nbsp; &gt; expect to have nearly 1000Tb of data from three different<br>&gt; &nbsp; &nbsp; &gt; reanalyses, and their production has also started. I don&#39;t see<br>&gt; &nbsp; &nbsp; &gt; storage as a barrier, but bandwidth is. Too many users making too<br>
&gt; &nbsp; &nbsp; &gt; many big requests will limit accessibility. Lats4d access to their<br>&gt; &nbsp; &nbsp; &gt; GDS will become an important function.<br>&gt; &nbsp; &nbsp; &gt;<br>&gt; &nbsp; &nbsp; &gt; Again, I have to admit I do not know the extent of the issues here.<br>
&gt; &nbsp; &nbsp; &gt; By necessity, I will have to use versions of grads that include<br>&gt; &nbsp; &nbsp; &gt; lats4d not only in my work, but as I show other how to use our<br>&gt; &nbsp; &nbsp; &gt; data. I felt the need to speak up and I hope that these issues will<br>
&gt; &nbsp; &nbsp; &gt; be considered.<br>&gt; &nbsp; &nbsp; &gt;<br>&gt; &nbsp; &nbsp; &gt; A sincere grads and lats4d user,<br>&gt; &nbsp; &nbsp; &gt;<br>&gt; &nbsp; &nbsp; &gt; Mike Bosilovich<br>&gt; &nbsp; &nbsp; &gt;<br>&gt;<br>&gt;<br></div></div></blockquote></div><br>