[gradsusr] Plotting of Freezing level

Tue Sep 24 13:25:33 EDT 2013

Use the fndlvl <http://www.iges.org/grads/gadoc/gradfuncfndlvl.html>function.

Jeff Duda

On Tue, Sep 24, 2013 at 3:44 AM, Geeta Geeta <geeta124 at hotmail.com> wrote:

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

-- 
Jeff Duda
Graduate research assistant
University of Oklahoma School of Meteorology
Center for Analysis and Prediction of Storms
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130924/237b2936/attachment-0001.html 