NCL Debugging and HDF-EOS Support for AMSR-E LAMAZ Projection?

From: Stewart Dickson <s-dickson_at_nyahnyahspammersnyahnyah>
Date: Fri Apr 08 2011 - 14:50:29 MDT

Hello, List,

Here's the short list of questions -- details below.

    Does NCL::addfile(".he2")->GridLon_, GridLat_ handle EOS (EASE)
    grids with projection type LAMAZ?
    Data fields of type DFNT_UINT8 ? The errors I see seem to relate
    to data type clashes, and I'm at a loss as to how to debug.

I am comparing a working example of visualizing HDF-EOS files in NCL:
      AMSR_E_L3_SeaIce6km_V11_20050118.hdf of EosZOO NSIDC Examples
                     
http://www.hdfeos.org/zoo/index_openNSIDC_Examples.php ->

AMSR_E_L3_SeaIce6km_V11_20050118_SI_06km_NH_89V_DAY.ncl
<http://www.hdfeos.org/zoo/NSIDC/AMSR_E_L3_SeaIce6km_V11_20050118_SI_06km_NH_89V_DAY.ncl>

       With HDF-EOS2 file:
ftp://ftp.hdfgroup.uiuc.edu/pub/outgoing/NASAHDF/AMSR_E_L3_5DaySnow_V09_20050126.hdf

                     [LAMAZ Projection] Listed further down on the same
page.

       ncl_filedump AMSR_E_L3_SeaIce6km_V11_20050118.hdf reports:

                     GridName="NpPolarGrid06km" ... Projection=GCTP_PS
                      DataFieldName="SI_06km_NH_89V_DAY"
                      DataType=DFNT_INT16 -- corresponds to NCL constant
Literal suffix 'h'

         ncl_filedump AMSR_E_L3_5DaySnow_V09_20050126.hdf reports:

                       GridName="Northern Hemisphere" ...
Projection=GCTP_LAMAZ
                       DataFieldName="SWE_NorthernPentad"
                       DataType=DFNT_UINT8 -- corresponds to NCL
constant Literal suffix 'B'

         My first attempt to visualize
AMSR_E_L3_5DaySnow_V09_20050126.hdf projection: 11 -- LAMAZ in NCL
begin
                     
eos_file=addfile("AMSR_E_L3_5DaySnow_V09_20050126.hdf.he2", "r")
                     
data=eos_file->SWE_NorthernPentad_Northern_Hemisphere(:,:)
                     data@lon2d=eos_file->GridLon_Northern_Hemisphere
                     data@lat2d=eos_file->GridLat_Northern_Hemisphere
                     data@_FillValue=0B
             Results in:
                         fatal:Type Mismatch: The type of missing value
could not be converted to type of variable (data)
             Changing:
                         data@_FillValue=0B
             To:
                         data@_FillValue=0b
             Results in:

fatal:NhlCvtStringToEnum: Unable to convert string "missing" to
requested type
                         warning:Error retrieving resource amJust from
args - Ignoring Arg
fatal:NhlCvtStringToEnum: Unable to convert string "missing" to
requested type
                         warning:Error retrieving resource amSide from
args - Ignoring Arg
... (Repeated 12 times)
            The remaining code consists of:
                
xwks=gsn_open_wks("pdf","AMSR_E_L3_5DaySnow_V09_20050126_SWE_NH")
res=True
res@... ;options...
gsn_define_colormap(xwks,"BlAqGrYeOrReVi200")
                plot=gsn_csm_contour_map(xwks,data,res) ; create plot

                delete(plot) ; cleaning up
                delete(xwks)
                delete(data)
                delete(res)
                delete(eos_file)
             end

Thanks,

-Stewart Dickson, HDF Group

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sat Apr 9 06:50:37 2011

This archive was generated by hypermail 2.1.8 : Tue Apr 19 2011 - 18:32:03 MDT