HDFEOS_v2.5 swath

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
Received on Mon Sep 29 2008 - 11:55:59 MDT