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

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Wed Apr 13 2011 - 15:51:59 MDT

Is the LAMAZ the Lambert Eqiual Area???

res@mpProjection = "LambertEqualArea"

===

On 04/08/2011 02:50 PM, Stewart Dickson wrote:
> 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

-- 
======================================================
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
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Apr 14 07:52:09 2011

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