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