Re: HDFEOS_v2.5 swath

Hi Chiara,

The geolocation fields supported by HDFEOS are only visible if you
read the file using the HDFEOS reader, which links
to the HDFEOS library, built on top of the HDF library.
So rather than opening the file as HDF you need to open it as HDFEOS.
That simply requires modifying your addfile call to:

f = addfile(diro+fname+".hdfeos", "r")

  -dave

On Sep 29, 2008, at 11:55 AM, Chiara Antoniazzi wrote:

> Hi,
> I'm trying to read data in hdfeos format (2B-GEOPROF-LIDAR: http://
> www.cloudsat.cira.colostate.edu/dataHome.php)
>
> If I open the file:
> f = addfile(diro+fname+".hdf", "r")
> printVarSummary(f)
> this is part of the output:
>
> granule_number : 9375
> ID_SITE : Cooperative Institute for Research in the Atmosphere
> ID_CENTER : CloudSat Data Processing Center
> ID_CENTER_URL : http://cloudsat.cira.colostate.edu
> ID_MACHINE : SKINKO
> ID_CREATED : Sun May 25 17:19:27 2008
> algorithm_name : 2B-GEOPROF-LIDAR
> algorithm_version : 1.1
> product_name : 2B-GEOPROF-LIDAR
> product_version : 003
> start_time : 20080201025505
> end_time : 20080201043358
> dimensions:
> nray_2B_GEOPROF_LIDAR = 37081
> nbin_2B_GEOPROF_LIDAR = 125
> ncloud_2B_GEOPROF_LIDAR = 5
>
> GROUP=GeoField
> OBJECT=GeoField_1
> GeoFieldName="Profile_time"
> DataType=DFNT_FLOAT32
> DimList=("nray")
> END_OBJECT=GeoField_1
> OBJECT=GeoField_2
> GeoFieldName="UTC_start"
> DataType=DFNT_FLOAT32
> DimList=("scalar")
> END_OBJECT=GeoField_2
> OBJECT=GeoField_3
> GeoFieldName="TAI_start"
> DataType=DFNT_FLOAT64
> DimList=("scalar")
> END_OBJECT=GeoField_3
> OBJECT=GeoField_4
> GeoFieldName="Latitude"
> DataType=DFNT_FLOAT32
> DimList=("nray")
> END_OBJECT=GeoField_4
> OBJECT=GeoField_5
> GeoFieldName="Longitude"
> DataType=DFNT_FLOAT32
> DimList=("nray")
> END_OBJECT=GeoField_5
> OBJECT=GeoField_6
> GeoFieldName="Height"
> DataType=DFNT_INT16
> DimList=("nray","nbin")
> END_OBJECT=GeoField_6
>
>
>
>
> variables:
> short Height ( nray_2B_GEOPROF_LIDAR, nbin_2B_GEOPROF_LIDAR )
> _FillValue : -9999
> factor : 1
> offset : 0
> long_name : Height of range bin in Reflectivity/Cloud
> Mask above reference surface (~ mean sea level).
> units : m
> valid_range : ( -5000, 30000 )
> missing : -9999
> missop : ==
> hdf_name : Height
>
> byte CloudFraction ( nray_2B_GEOPROF_LIDAR,
> nbin_2B_GEOPROF_LIDAR )
> _FillValue : 0xf7
> factor : 1
> offset : 0
> long_name : Cloud Fraction
> valid_range : ( 0x00, 0x64 )
> missing : 0xf7
> missop : ==
> hdf_name : CloudFraction
>
> But I cannot read time, latitude and longitude that I need.
> Is there any way to read these data correctly using NCL?
>
>
> Thanks.
> Chiara
>
>
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Sep 29 2008 - 12:15:54 MDT