<div dir="ltr"><div>Use the <a href="http://www.iges.org/grads/gadoc/gradfuncfndlvl.html">fndlvl</a> function.<br><br></div>Jeff Duda<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Sep 24, 2013 at 3:44 AM, Geeta Geeta <span dir="ltr">&lt;<a href="mailto:geeta124@hotmail.com" target="_blank">geeta124@hotmail.com</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">


<div><div dir="ltr">Hi all, <br>I want to know how can plot the freezing level using Grads. <br>I have the WRF output which I am post processing using ARWpost. <br>Any hints???<br><br>geeta<br> <br><div>&gt; From: <a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a><br>
&gt; Subject: gradsusr Digest, Vol 43, Issue 47<br>&gt; To: <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>&gt; Date: Thu, 19 Sep 2013 18:51:45 -0400<br>&gt; <br>&gt; Send gradsusr mailing list submissions to<br>
&gt;         <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>&gt; <br>&gt; To subscribe or unsubscribe via the World Wide Web, visit<br>&gt;         <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt; or, via email, send a message with subject or body &#39;help&#39; to<br>&gt;         <a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a><br>&gt; <br>&gt; You can reach the person managing the list at<br>
&gt;         <a href="mailto:gradsusr-owner@gradsusr.org" target="_blank">gradsusr-owner@gradsusr.org</a><br>&gt; <br>&gt; When replying, please edit your Subject line so it is more specific<br>&gt; than &quot;Re: Contents of gradsusr digest...&quot;<br>
&gt; <br>&gt; <br>&gt; Today&#39;s Topics:<br>&gt; <br>&gt;    1. Re: GrADS Control File for Array with Spectral Band        Dimension<br>&gt;       (Jennifer Adams)<br>&gt;    2. Re: GrADS Control File for Array with Spectral Band        Dimension<br>
&gt;       (James E. Johnson)<br>&gt;    3. Problem Solved - Difference Files (Phil Klotzbach)<br>&gt; <br>&gt; <br>&gt; ----------------------------------------------------------------------<br>&gt; <br>&gt; Message: 1<br>
&gt; Date: Thu, 19 Sep 2013 16:55:58 -0400<br>&gt; From: Jennifer Adams &lt;<a href="mailto:jma@cola.iges.org" target="_blank">jma@cola.iges.org</a>&gt;<br>&gt; Subject: Re: [gradsusr] GrADS Control File for Array with Spectral<br>
&gt;         Band        Dimension<br>&gt; To: GrADS Users Forum &lt;<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>&gt;<br>&gt; Message-ID: &lt;<a href="mailto:FEDA1F66-6D49-470D-A8E0-F78EE3E928CA@cola.iges.org" target="_blank">FEDA1F66-6D49-470D-A8E0-F78EE3E928CA@cola.iges.org</a>&gt;<br>
&gt; Content-Type: text/plain; charset=&quot;us-ascii&quot;<br>&gt; <br>&gt; First of all, thank you for providing all the information and data I needed to reproduce the error! The failure of method 1 (using array indices in the units field of the variable declaration instead of x,y,z,t,e) is indeed a bug that is unique to the hdf5 data type. The critical bit of information in the error message from the HDF5 library was this: &quot;src and dest data spaces have different sizes.&quot; In method 1, the variables are described to GrADS as being only 2-dimensional (i.e., not Z-varying) even though they really have 3 varying dimensions, and the library trips on this deception. The netcdf library is not as fussy. I will put a fix in version 2.1, but for now you will have to use method 2. <br>
&gt; --Jennifer<br>&gt; <br>&gt; <br>&gt; On Sep 19, 2013, at 3:29 PM, James E. Johnson wrote:<br>&gt; <br>&gt; &gt; I&#39;m using grads-2.0.2-bin-CentOS5.8-x86_64.tar.gz<br>&gt; &gt; <br>&gt; &gt; I have an HDF-EOS5 file with some data arrays that are 3 dimensional, with the <br>
&gt; &gt; third dimension representing wavelength/spectral band. I created a control <br>&gt; &gt; file for this data following the description at <br>&gt; &gt; <a href="http://www.iges.org/grads/gadoc/gadocindex.html" target="_blank">http://www.iges.org/grads/gadoc/gadocindex.html</a> (section VAR subsection <br>
&gt; &gt; units). It says there are two ways to implement this, first with an array <br>&gt; &gt; index in the units dimension list, the second with dimension values specified <br>&gt; &gt; in the Z axis. The second method works, however the first method keeps giving <br>
&gt; &gt; me errors:<br>&gt; &gt; <br>&gt; &gt; HDF5-DIAG: Error detected in HDF5 (1.8.7) thread 0:<br>&gt; &gt;  #000: H5Dio.c line 174 in H5Dread(): can&#39;t read data<br>&gt; &gt;    major: Dataset<br>&gt; &gt;    minor: Read failed<br>
&gt; &gt;  #001: H5Dio.c line 348 in H5D_read(): src and dest data spaces have <br>&gt; &gt; different sizes<br>&gt; &gt;    major: Invalid arguments to routine<br>&gt; &gt;    minor: Bad value<br>&gt; &gt; H5Dread failed for short data<br>
&gt; &gt; Data Request Error:  Error for variable &#39;aaod342&#39;<br>&gt; &gt;  Error ocurred at column 1<br>&gt; &gt; DISPLAY error:  Invalid expression<br>&gt; &gt;  Expression = aaod342<br>&gt; &gt; <br>&gt; &gt; <br>
&gt; &gt; <br>&gt; &gt; The sample data file is:<br>&gt; &gt; <a href="ftp://acdisc.gsfc.nasa.gov/data/s4pa/Aura_OMI_Level3/OMAEROe.003/2013/OMI-" target="_blank">ftp://acdisc.gsfc.nasa.gov/data/s4pa/Aura_OMI_Level3/OMAEROe.003/2013/OMI-</a><br>
&gt; &gt; Aura_L3-OMAEROe_2013m0601_v003-2013m0603t015851.he5<br>&gt; &gt; <br>&gt; &gt; <br>&gt; &gt; First method (array index) control file:<br>&gt; &gt; DSET OMI-Aura_L3-OMAEROe_2013m0601_v003-2013m0603t015851.he5<br>
&gt; &gt; DTYPE hdf5_grid<br>&gt; &gt; TITLE OMI/Aura Multi-wavelength Aerosol Optical Depth and Single Scattering <br>&gt; &gt; Albedo Daily L3 Global 0.25x0.25 deg Lat/Lon Grid for 1 June 2013<br>&gt; &gt; UNDEF -32767 _FillValue<br>
&gt; &gt; XDEF 1440 linear -180 0.25<br>&gt; &gt; YDEF 720 linear -90 0.25<br>&gt; &gt; ZDEF 1 levels 1<br>&gt; &gt; TDEF 1 linear 1jun2013 1dy<br>&gt; &gt; VARS 5<br>&gt; &gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/AbsorbingAerosolOpticalThicknessMW=&gt;aaod342  <br>
&gt; &gt; 0  0,y,x  Absorbing Aerosol Optical Thickness at 342.5 nm<br>&gt; &gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/AbsorbingAerosolOpticalThicknessMW=&gt;aaod388  <br>&gt; &gt; 0  1,y,x  Absorbing Aerosol Optical Thickness at 388.0 nm<br>
&gt; &gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/AbsorbingAerosolOpticalThicknessMW=&gt;aaod442  <br>&gt; &gt; 0  2,y,x  Absorbing Aerosol Optical Thickness at 442.0 nm<br>&gt; &gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/AbsorbingAerosolOpticalThicknessMW=&gt;aaod463  <br>
&gt; &gt; 0  3,y,x  Absorbing Aerosol Optical Thickness at 463.0 nm<br>&gt; &gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/AbsorbingAerosolOpticalThicknessMW=&gt;aaod483  <br>&gt; &gt; 0  4,y,x  Absorbing Aerosol Optical Thickness at 483.5 nm<br>
&gt; &gt; ENDVARS<br>&gt; &gt; <br>&gt; &gt; <br>&gt; &gt; Second method (Z dimension) control file:<br>&gt; &gt; DSET OMI-Aura_L3-OMAEROe_2013m0601_v003-2013m0603t015851.he5<br>&gt; &gt; DTYPE hdf5_grid<br>&gt; &gt; TITLE OMI/Aura Multi-wavelength Aerosol Optical Depth and Single Scattering <br>
&gt; &gt; Albedo Daily L3 Global 0.25x0.25 deg Lat/Lon Grid for 1 June 2013<br>&gt; &gt; UNDEF -32767 _FillValue<br>&gt; &gt; XDEF 1440 linear -180 0.25<br>&gt; &gt; YDEF 720 linear -90 0.25<br>&gt; &gt; ZDEF 5 levels 342.5 388.0 442.0 463.0 483.5<br>
&gt; &gt; TDEF 1 linear 1jun2013 1dy<br>&gt; &gt; VARS 1<br>&gt; &gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/AbsorbingAerosolOpticalThicknessMW=&gt;aaod  <br>&gt; &gt; 5  z,y,x  Absorbing Aerosol Optical Thickness<br>
&gt; &gt; ENDVARS<br>&gt; &gt; <br>&gt; &gt; <br>&gt; &gt; How do I get method 1 (array index) to work? Is there something wrong with <br>&gt; &gt; ZDEF? Or is this a GrADS bug?<br>&gt; &gt; <br>&gt; &gt; ------------------------------------------------------------------------------<br>
&gt; &gt; | James E. Johnson                   | address:                              |<br>&gt; &gt; | ADNET Systems, Inc.                |  NASA Goddard Space Flight Center     |<br>&gt; &gt; | e-mail: <a href="mailto:james.johnson@nasa.gov" target="_blank">james.johnson@nasa.gov</a>     |  Data and Information Services Center |<br>
&gt; &gt; | phone:  301-614-5121               |  Code 610.2, Bldg 32, Room S130G      |<br>&gt; &gt; | fax:    301-614-5268               |  Greenbelt, MD  20771                 |<br>&gt; &gt; ------------------------------------------------------------------------------<br>
&gt; &gt; _______________________________________________<br>&gt; &gt; gradsusr mailing list<br>&gt; &gt; <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>&gt; &gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt; <br>&gt; --<br>&gt; Jennifer M. Adams<br>&gt; Center for Ocean-Land-Atmosphere Studies (COLA)<br>&gt; 111 Research Hall, Mail Stop 2B3<br>&gt; George Mason University<br>&gt; 4400 University Drive<br>&gt; Fairfax, VA 22030 <br>
&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; -------------- next part --------------<br>&gt; An HTML attachment was scrubbed...<br>&gt; URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20130919/406de0ac/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20130919/406de0ac/attachment-0001.html</a> <br>
&gt; <br>&gt; ------------------------------<br>&gt; <br>&gt; Message: 2<br>&gt; Date: Thu, 19 Sep 2013 18:11:44 -0400<br>&gt; From: &quot;James E. Johnson&quot; &lt;<a href="mailto:James.Johnson@nasa.gov" target="_blank">James.Johnson@nasa.gov</a>&gt;<br>
&gt; Subject: Re: [gradsusr] GrADS Control File for Array with Spectral<br>&gt;         Band        Dimension<br>&gt; To: &lt;<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>&gt;<br>&gt; Message-ID: &lt;<a href="mailto:201309191811.44416.James.Johnson@nasa.gov" target="_blank">201309191811.44416.James.Johnson@nasa.gov</a>&gt;<br>
&gt; Content-Type: text/plain; charset=&quot;us-ascii&quot;<br>&gt; <br>&gt; On Thursday 19 September 2013 04:55:58 pm Jennifer Adams wrote:<br>&gt; &gt; First of all, thank you for providing all the information and data I needed<br>
&gt; &gt;  to reproduce the error! The failure of method 1 (using array indices in<br>&gt; &gt;  the units field of the variable declaration instead of x,y,z,t,e) is<br>&gt; &gt;  indeed a bug that is unique to the hdf5 data type. The critical bit of<br>
&gt; &gt;  information in the error message from the HDF5 library was this: &quot;src and<br>&gt; &gt;  dest data spaces have different sizes.&quot; In method 1, the variables are<br>&gt; &gt;  described to GrADS as being only 2-dimensional (i.e., not Z-varying) even<br>
&gt; &gt;  though they really have 3 varying dimensions, and the library trips on<br>&gt; &gt;  this deception. The netcdf library is not as fussy. I will put a fix in<br>&gt; &gt;  version 2.1, but for now you will have to use method 2. --Jennifer<br>
&gt; <br>&gt; Thanks Jennifer,<br>&gt; <br>&gt; I think there may be another bug. I have a more extensive control file with <br>&gt; all the variables from the file. I can successfully display data from the <br>&gt; first variable I read, but subsequent variables complain:<br>
&gt; <br>&gt; HDF5-DIAG: Error detected in HDF5 (1.8.7) thread 0:<br>&gt;   #000: H5Shyper.c line 6593 in H5Sselect_hyperslab(): not a data space<br>&gt;     major: Invalid arguments to routine<br>&gt;     minor: Inappropriate type<br>
&gt; HDF5 Error: unable to select dataspace hyperslab<br>&gt; Data Request Error:  Error for variable &#39;uvai&#39;<br>&gt;   Error ocurred at column 1<br>&gt; DISPLAY error:  Invalid expression<br>&gt;   Expression = uvai<br>
&gt; <br>&gt; Here is the complete config file for this product:<br>&gt; DSET OMI-Aura_L3-OMAEROe_2013m0601_v003-2013m0603t015851.he5<br>&gt; DTYPE hdf5_grid<br>&gt; TITLE OMI/Aura Multi-wavelength Aerosol Optical Depth and Single Scattering <br>
&gt; Albedo Daily L3 Global 0.25x0.25 deg Lat/Lon Grid<br>&gt; UNDEF -32767 _FillValue<br>&gt; UNPACK ScaleFactor Offset<br>&gt; XDEF 1440 linear -180 0.25<br>&gt; YDEF 720 linear -90 0.25<br>&gt; ZDEF 5 levels 342.5 388.0 442.0 463.0 483.5<br>
&gt; TDEF 5000 linear 1oct2004 1dy<br>&gt; VARS 15<br>&gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/AbsorbingAerosolOpticalThicknessMW=&gt;aaod  <br>&gt; 5  z,y,x  Absorbing Aerosol Optical Thickness<br>&gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/AerosolModelMW=&gt;amod 0 y,x  <br>
&gt; Aerosol Model<br>&gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/AerosolOpticalThicknessMW=&gt;aod  <br>&gt; 5  z,y,x  Aerosol Optical Thickness<br>&gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/AerosolOpticalThicknessPassedThresholdMean=&gt;aodptm  <br>
&gt; 0  y,x  Aerosol Optical Thickness Passed Threshold Mean<br>&gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/AerosolOpticalThicknessPassedThresholdStd=&gt;aodpts  <br>&gt; 0  y,x  Aerosol Optical Thickness Passed Threshold Standard Deviation<br>
&gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/Latitude=&gt;lat  0  y,x  Latitude<br>&gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/Longitude=&gt;lon  0  y,x  <br>&gt; Longitude<br>&gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/SingleScatteringAlbedoMW=&gt;ssa  5  <br>
&gt; z,y,x  Single Scattering Albedo<br>&gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/SingleScatteringAlbedoPassedThresholdMean=&gt;ssaptm  <br>&gt; 0  y,x  Single Scattering Albedo Passed Threshold Mean<br>&gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/SingleScatteringAlbedoPassedThresholdStd=&gt;ssapts  <br>
&gt; 0  y,x  Single Scattering Albedo Passed Threshold Standard Deviation<br>&gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/SolarZenithAngle=&gt;sza  0  y,x  <br>&gt; Solar Zenith Angle<br>&gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/TerrainReflectivity=&gt;tref  0  <br>
&gt; y,x  Terrain Reflectivity<br>&gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/UVAerosolIndex=&gt;uvai  0  y,x  UV <br>&gt; Aerosol Index<br>&gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/VISAerosolIndex=&gt;visai  0  y,x  <br>
&gt; VIS Aerosol Index<br>&gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/ViewingZenithAngle=&gt;vza  0  y,x <br>&gt; Viewing Zenith Angle<br>&gt; ENDVARS<br>&gt; <br>&gt; <br>&gt; Does the optional CACHESIZE attribute needed to be set?<br>
&gt; <br>&gt; James Johnson<br>&gt; <br>&gt; &gt; On Sep 19, 2013, at 3:29 PM, James E. Johnson wrote:<br>&gt; &gt; &gt; I&#39;m using grads-2.0.2-bin-CentOS5.8-x86_64.tar.gz<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; I have an HDF-EOS5 file with some data arrays that are 3 dimensional,<br>
&gt; &gt; &gt; with the third dimension representing wavelength/spectral band. I created<br>&gt; &gt; &gt; a control file for this data following the description at<br>&gt; &gt; &gt; <a href="http://www.iges.org/grads/gadoc/gadocindex.html" target="_blank">http://www.iges.org/grads/gadoc/gadocindex.html</a> (section VAR subsection<br>
&gt; &gt; &gt; units). It says there are two ways to implement this, first with an array<br>&gt; &gt; &gt; index in the units dimension list, the second with dimension values<br>&gt; &gt; &gt; specified in the Z axis. The second method works, however the first<br>
&gt; &gt; &gt; method keeps giving me errors:<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; HDF5-DIAG: Error detected in HDF5 (1.8.7) thread 0:<br>&gt; &gt; &gt;  #000: H5Dio.c line 174 in H5Dread(): can&#39;t read data<br>&gt; &gt; &gt;    major: Dataset<br>
&gt; &gt; &gt;    minor: Read failed<br>&gt; &gt; &gt;  #001: H5Dio.c line 348 in H5D_read(): src and dest data spaces have<br>&gt; &gt; &gt; different sizes<br>&gt; &gt; &gt;    major: Invalid arguments to routine<br>&gt; &gt; &gt;    minor: Bad value<br>
&gt; &gt; &gt; H5Dread failed for short data<br>&gt; &gt; &gt; Data Request Error:  Error for variable &#39;aaod342&#39;<br>&gt; &gt; &gt;  Error ocurred at column 1<br>&gt; &gt; &gt; DISPLAY error:  Invalid expression<br>
&gt; &gt; &gt;  Expression = aaod342<br>&gt; &gt; &gt;<br>&gt; &gt; &gt;<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; The sample data file is:<br>&gt; &gt; &gt; <a href="ftp://acdisc.gsfc.nasa.gov/data/s4pa/Aura_OMI_Level3/OMAEROe.003/2013/OMI" target="_blank">ftp://acdisc.gsfc.nasa.gov/data/s4pa/Aura_OMI_Level3/OMAEROe.003/2013/OMI</a><br>
&gt; &gt; &gt;- Aura_L3-OMAEROe_2013m0601_v003-2013m0603t015851.he5<br>&gt; &gt; &gt;<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; First method (array index) control file:<br>&gt; &gt; &gt; DSET OMI-Aura_L3-OMAEROe_2013m0601_v003-2013m0603t015851.he5<br>
&gt; &gt; &gt; DTYPE hdf5_grid<br>&gt; &gt; &gt; TITLE OMI/Aura Multi-wavelength Aerosol Optical Depth and Single<br>&gt; &gt; &gt; Scattering Albedo Daily L3 Global 0.25x0.25 deg Lat/Lon Grid for 1 June<br>&gt; &gt; &gt; 2013 UNDEF -32767 _FillValue<br>
&gt; &gt; &gt; XDEF 1440 linear -180 0.25<br>&gt; &gt; &gt; YDEF 720 linear -90 0.25<br>&gt; &gt; &gt; ZDEF 1 levels 1<br>&gt; &gt; &gt; TDEF 1 linear 1jun2013 1dy<br>&gt; &gt; &gt; VARS 5<br>&gt; &gt; &gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/AbsorbingAerosolOpticalThic<br>
&gt; &gt; &gt;knessMW=&gt;aaod342 0  0,y,x  Absorbing Aerosol Optical Thickness at 342.5 nm<br>&gt; &gt; &gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/AbsorbingAerosolOpticalThic<br>&gt; &gt; &gt;knessMW=&gt;aaod388 0  1,y,x  Absorbing Aerosol Optical Thickness at 388.0 nm<br>
&gt; &gt; &gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/AbsorbingAerosolOpticalThic<br>&gt; &gt; &gt;knessMW=&gt;aaod442 0  2,y,x  Absorbing Aerosol Optical Thickness at 442.0 nm<br>&gt; &gt; &gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/AbsorbingAerosolOpticalThic<br>
&gt; &gt; &gt;knessMW=&gt;aaod463 0  3,y,x  Absorbing Aerosol Optical Thickness at 463.0 nm<br>&gt; &gt; &gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/AbsorbingAerosolOpticalThic<br>&gt; &gt; &gt;knessMW=&gt;aaod483 0  4,y,x  Absorbing Aerosol Optical Thickness at 483.5 nm<br>
&gt; &gt; &gt; ENDVARS<br>&gt; &gt; &gt;<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; Second method (Z dimension) control file:<br>&gt; &gt; &gt; DSET OMI-Aura_L3-OMAEROe_2013m0601_v003-2013m0603t015851.he5<br>&gt; &gt; &gt; DTYPE hdf5_grid<br>
&gt; &gt; &gt; TITLE OMI/Aura Multi-wavelength Aerosol Optical Depth and Single<br>&gt; &gt; &gt; Scattering Albedo Daily L3 Global 0.25x0.25 deg Lat/Lon Grid for 1 June<br>&gt; &gt; &gt; 2013 UNDEF -32767 _FillValue<br>&gt; &gt; &gt; XDEF 1440 linear -180 0.25<br>
&gt; &gt; &gt; YDEF 720 linear -90 0.25<br>&gt; &gt; &gt; ZDEF 5 levels 342.5 388.0 442.0 463.0 483.5<br>&gt; &gt; &gt; TDEF 1 linear 1jun2013 1dy<br>&gt; &gt; &gt; VARS 1<br>&gt; &gt; &gt; /HDFEOS/GRIDS/ColumnAmountAerosol/Data~Fields/AbsorbingAerosolOpticalThic<br>
&gt; &gt; &gt;knessMW=&gt;aaod 5  z,y,x  Absorbing Aerosol Optical Thickness<br>&gt; &gt; &gt; ENDVARS<br>&gt; &gt; &gt;<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; How do I get method 1 (array index) to work? Is there something wrong<br>
&gt; &gt; &gt; with ZDEF? Or is this a GrADS bug?<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; -------------------------------------------------------------------------<br>&gt; &gt; &gt;-----<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; | James E. Johnson                   | address:                          <br>
&gt; &gt; &gt; |    | ADNET Systems, Inc.                |  NASA Goddard Space Flight<br>&gt; &gt; &gt; | Center     | e-mail: <a href="mailto:james.johnson@nasa.gov" target="_blank">james.johnson@nasa.gov</a>     |  Data and Information<br>
&gt; &gt; &gt; | Services Center | phone:  301-614-5121               |  Code 610.2,<br>&gt; &gt; &gt; | Bldg 32, Room S130G      | fax:    301-614-5268               | <br>&gt; &gt; &gt; | Greenbelt, MD  20771                 |<br>
&gt; &gt; &gt;<br>&gt; &gt; &gt; -------------------------------------------------------------------------<br>&gt; &gt; &gt;----- _______________________________________________<br>&gt; &gt; &gt; gradsusr mailing list<br>
&gt; &gt; &gt; <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>&gt; &gt; &gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt; &gt; <br>&gt; &gt; --<br>&gt; &gt; Jennifer M. Adams<br>&gt; &gt; Center for Ocean-Land-Atmosphere Studies (COLA)<br>&gt; &gt; 111 Research Hall, Mail Stop 2B3<br>&gt; &gt; George Mason University<br>&gt; &gt; 4400 University Drive<br>
&gt; &gt; Fairfax, VA 22030<br>&gt; &gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; ------------------------------<br>&gt; <br>&gt; Message: 3<br>&gt; Date: Thu, 19 Sep 2013 15:51:39 -0700<br>&gt; From: Phil Klotzbach &lt;<a href="mailto:philk@atmos.colostate.edu" target="_blank">philk@atmos.colostate.edu</a>&gt;<br>
&gt; Subject: [gradsusr] Problem Solved - Difference Files<br>&gt; To: <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>&gt; Message-ID: &lt;<a href="mailto:523B7FFB.3070307@atmos.colostate.edu" target="_blank">523B7FFB.3070307@atmos.colostate.edu</a>&gt;<br>
&gt; Content-Type: text/plain; charset=&quot;iso-8859-1&quot;<br>&gt; <br>&gt; Hi everyone,<br>&gt; <br>&gt; Thanks for all of your very helpful comments.  I solved the problem by <br>&gt; concatenating the two files together using ncrcat (<a href="http://uwnd.2013.nc" target="_blank">uwnd.2013.nc</a> and <br>
&gt; <a href="http://uwnd.ltm.nc" target="_blank">uwnd.ltm.nc</a>).  Then, I set the dimensions and did the following difference:<br>&gt; <br>&gt; define uwnd2 = aave(uwnd-(uwnd(t+258), lon=-70, lon=-20, lat=0, lat=30)<br>
&gt; <br>&gt; <br>&gt; I did 258 days because that&#39;s how many days were available for <br>&gt; computation in 2013 so far.<br>&gt; <br>&gt; <br>&gt; Thanks again for the comments.  This list is an invaluable resource!<br>
&gt; <br>&gt; Phil<br>&gt; <br>&gt; <br>&gt; -- <br>&gt; ----------------------------------------------------<br>&gt; Phil Klotzbach, Ph.D.<br>&gt; Research Scientist<br>&gt; Department of Atmospheric Science<br>&gt; Colorado State University<br>
&gt; Email: <a href="mailto:philk@atmos.colostate.edu" target="_blank">philk@atmos.colostate.edu</a><br>&gt; ----------------------------------------------------<br>&gt; <br>&gt; -------------- next part --------------<br>
&gt; An HTML attachment was scrubbed...<br>&gt; URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20130919/94c48199/attachment.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20130919/94c48199/attachment.html</a> <br>
&gt; <br>&gt; ------------------------------<br>&gt; <br>&gt; _______________________________________________<br>&gt; gradsusr mailing list<br>&gt; <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>&gt; <br>&gt; <br>&gt; End of gradsusr Digest, Vol 43, Issue 47<br>&gt; ****************************************<br>
</div>                                               </div></div>
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms<br>
</div>