Re: reading HDF data

It looks like Examples 4 and 5 at

http://www.ncl.ucar.edu/Applications/eosdis.shtml

Note ... your file is hdfeos

HDFEOSVersion : HDFEOS_V2.8

Add a .hdfeos extension in the addfile function. No need to change the
actual file name:

f = addfile("2007014182825_03807_CS_ZZ_GRANULE.hdf.hdfeos", "r")
^^^^^^

In any event .... I don't see that the file has lat and lon values!!!!!

Tip: when plotting a (37088,125) gris use Raster mode.

D

Cory Wolff wrote:
> I am working with some data that is in HDF format and am not sure how to
> read all of the information out of the file. The output from printing
> this file in ncl (v4.3) is attached. I am able to pull out the
> Reflectivity data because it is a variable, but I also need other
> information such as "Time", "Latitude", "Longitude", and "TopBinHeight",
> which appear to be in the file. I can't seem to get at that information
> at this point. Any help is appreciated.
>
> Cory
>
>
> ------------------------------------------------------------------------
>
> ncl 39> print(f)
> Variable: f (file variable)
>
> filename: 2007014182825_03807_CS_ZZ_GRANULE
> path: 2007014182825_03807_CS_ZZ_GRANULE.hdf
> file global attributes:
> HDFEOSVersion : HDFEOS_V2.8
> StructMetadata_0 : ????athStructure
> GROUP=SWATH_1
> SwathName="CloudSat - REFLECTIVITY 01/14/2007"
> GROUP=Dimension
> OBJECT=Dimension_1
> DimensionName="scalar"
> Size=1
> END_OBJECT=Dimension_1
> OBJECT=Dimension_2
> DimensionName="nbin"
> Size=125
> END_OBJECT=Dimension_2
> OBJECT=Dimension_3
> DimensionName="nray"
> Size=37088
> END_OBJECT=Dimension_3
> END_GROUP=Dimension
> GROUP=DimensionMap
> END_GROUP=DimensionMap
> GROUP=IndexDimensionMap
> END_GROUP=IndexDimensionMap
> GROUP=GeoField
> OBJECT=GeoField_1
> GeoFieldName="Time"
> DataType=DFNT_FLOAT32
> DimList=("nray")
> END_OBJECT=GeoField_1
> OBJECT=GeoField_2
> GeoFieldName="Latitude"
> DataType=DFNT_FLOAT32
> DimList=("nray")
> END_OBJECT=GeoField_2
> OBJECT=GeoField_3
> GeoFieldName="Longitude"
> DataType=DFNT_FLOAT32
> DimList=("nray")
> END_OBJECT=GeoField_3
> OBJECT=GeoField_4
> GeoFieldName="DEM"
> DataType=DFNT_FLOAT32
> DimList=("nray")
> END_OBJECT=GeoField_4
> OBJECT=GeoField_5
> GeoFieldName="LandSeaFlag"
> DataType=DFNT_UINT8
> DimList=("nray")
> END_OBJECT=GeoField_5
> END_GROUP=GeoField
> GROUP=DataField
> OBJECT=DataField_1
> DataFieldName="BinSize"
> DataType=DFNT_FLOAT32
> DimList=("scalar")
> END_OBJECT=DataField_1
> OBJECT=DataField_2
> DataFieldName="SceneMap"
> DataType=DFNT_UINT8
> DimList=("nray")
> END_OBJECT=DataField_2
> OBJECT=DataField_3
> DataFieldName="TopBinHeight"
> DataType=DFNT_FLOAT32
> DimList=("nray")
> END_OBJECT=DataField_3
> OBJECT=DataField_4
> DataFieldName="Reflectivity"
> DataType=DFNT_FLOAT32
> DimList=("nray","nbin")
> END_OBJECT=DataField_4
> OBJECT=DataField_5
> DataFieldName="Zsfc"
> DataType=DFNT_FLOAT32
> DimList=("nray")
> END_OBJECT=DataField_5
> END_GROUP=DataField
> GROUP=MergedFields
> END_GROUP=MergedFields
> END_GROUP=SWATH_1
> END_GROUP=SwathStructure
> GROUP=GridStructure
> END_GROUP=GridStructure
> GROUP=PointStructure
> END_GROUP=PointStructure
> END
>
> dimensions:
> nray_CloudSat_REFLECTIVITY_01_14_2007 = 37088 // unlimited
> nbin_CloudSat_REFLECTIVITY_01_14_2007 = 125 // unlimited
> variables:
> float Reflectivity ( nray_CloudSat_REFLECTIVITY_01_14_2007, nbin_CloudSat_REFLECTIVITY_01_14_2007 )
> hdf_name : Reflectivity
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>

-- 
======================================================
Dennis J. Shea                  tel: 303-497-1361    |
P.O. Box 3000                   fax: 303-497-1333    |
Climate Analysis Section                             |
Climate & Global Dynamics Div.                       |
National Center for Atmospheric Research             |
Boulder, CO  80307                                   |
USA                        email: shea 'at' ucar.edu |
======================================================
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk