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

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Mon Apr 11 2011 - 16:50:01 MDT

NCL 6.0.0-beta has different "literal suffix" for the differnt types
EK ubyte has "B"

     0B

See:
http://www.ncl.ucar.edu/Document/Manuals/Ref_Manual/NclDataTypes.shtml

On 04/11/2011 04:41 PM, Stewart Dickson wrote:
> Thank you, Wei,
>
> Yes, I was previously using: ncl_ncarg-5.2.1.Linux_RedHat_i686_gcc412
> I have just now downloaded and installed:
> ncl_ncarg-6.0.0-beta.Linux_RedHat_i686_gcc412
> I find the results from the on-line example script:
> 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>
> To be the same as what I obtained previously (the control case) ;)
> I see the same results in ncl_ncarg-6.0.0-beta.Linux_RedHat_i686_gcc412
> as you report for
> ncl_filedump AMSR_E_L3_5DaySnow_V09_20050126.hdf
> and
> ncl_filedump AMSR_E_L3_5DaySnow_V09_20050126.hdf.he2
>
> I note the difference in data field type and _FillValue as reported by
> the NCL HDF4 interface vs. the NCL HDF-EOS interface:
> ncl_filedump AMSR_E_L3_5DaySnow_V09_20050126.hdf
> variables:
> ubyte SWE_NorthernPentad ( YDim_Northern_Hemisphere,
> XDim_Northern_Hemisphere )
> _FillValue : 255
> hdf_name : SWE_NorthernPentad
> ncl_filedump AMSR_E_L3_5DaySnow_V09_20050126.hdf.he2
> (signed)byte SWE_NorthernPentad_Northern_Hemisphere
> _FillValue :-1
>
> I am using:
> eos_file=addfile("AMSR_E_L3_5DaySnow_V09_20050126.hdf.he2", "r") ; The
> HDF-EOS2 Interface
> and
> data=eos_file->SWE_NorthernPentad_Northern_Hemisphere ; read data field(:,:)
>
> Do I also see correctly from ncl_filedump
> AMSR_E_L3_5DaySnow_V09_20050126.hdf.he2
> that
> data@lon2d=eos_file->GridLon_Northern_Hemisphere
> and
> data@lat2d=eos_file->GridLat_Northern_Hemisphere
> Should be well-defined?
>
> I am using the following:
> data@_FillValue=0b ; Making this double assignment allows the _FillValue
> of 0
> ; to be properly represented.
> data@_FillValue=-1b
> xwks=gsn_open_wks("pdf","AMSR_E_L3_5DaySnow_V09_20050126_SWE_NH")
> res=True ; plot mods desired
> res@... ; Options
> plot=gsn_csm_contour_map(xwks,data,res) ; create plot
>
> I am (still) seeing:
> 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)
> ERROR 2 IN MAPPOS/MDPPOS - ARGUMENTS ARE INCORRECT
> ERROR 1 IN MAPSTI - UNCLEARED PRIOR ERROR
> ERROR 1 IN MAPROJ - UNCLEARED PRIOR ERROR
> ERROR 2 IN MAPSET/MDPSET - UNCLEARED PRIOR ERROR
> fatal:MapSetTrans: error initializing map: MAPINT - UNCLEARED PRIOR ERROR
> warning:PDFWorkstationInitialize:Device X Coordinates invalid, defaulting
> warning:PDFWorkstationInitialize:Device Y Coordinates invalid, defaulting
>
> This time, in ncl_ncarg-6.0.0-beta.Linux_RedHat_i686_gcc412, NCL
> produces a PDF file with map axes, grid lines, decoration (continental
> outlines, etc.) but no contour plot data; whereas previously, in
> ncl_ncarg-5.2.1.Linux_RedHat_i686_gcc412, NCL produced no output file.
>
> Thanks again,
>
> -Stewart
>
> On 4/11/2011 11:01 AM, Wei Huang wrote:
>> Stewart,
>>
>> Which NCL version are you using (ncl -V)?
>> I'd suggest you try out NCL6.0.0beta,
>>
>> Here, use data as HDF(4), I have:
>> ncl_filedump AMSR_E_L3_5DaySnow_V09_20050126.hdf
>>
>> ...
>> dimensions:
>> YDim_Northern_Hemisphere = 721
>> XDim_Northern_Hemisphere = 721
>> YDim_Southern_Hemisphere = 721
>> XDim_Southern_Hemisphere = 721
>> variables:
>> ubyte SWE_NorthernPentad ( YDim_Northern_Hemisphere,
>> XDim_Northern_Hemisphere )
>> _FillValue : 255
>> hdf_name : SWE_NorthernPentad
>>
>> ubyte Flags_NorthernPentad ( YDim_Northern_Hemisphere,
>> XDim_Northern_Hemisphere )
>> _FillValue : 255
>> hdf_name : Flags_NorthernPentad
>>
>> ubyte SWE_SouthernPentad ( YDim_Southern_Hemisphere,
>> XDim_Southern_Hemisphere )
>> _FillValue : 255
>> hdf_name : SWE_SouthernPentad
>>
>> ubyte Flags_SouthernPentad ( YDim_Southern_Hemisphere,
>> XDim_Southern_Hemisphere )
>> hdf_name : Flags_SouthernPentad
>>
>>
>>
>> Use as HE2, we can see:
>> ncl_filedump AMSR_E_L3_5DaySnow_V09_20050126.hdf.he2
>>
>> ....
>>
>>
>> Variable: f
>> Type: file
>> filename:AMSR_E_L3_5DaySnow_V09_20050126.hdf
>> path:AMSR_E_L3_5DaySnow_V09_20050126.hdf
>> file global attributes:
>> dimensions:
>> YDim_Northern_Hemisphere = 721
>> XDim_Northern_Hemisphere = 721
>> YDim_Southern_Hemisphere = 721
>> XDim_Southern_Hemisphere = 721
>> variables:
>> double GridLat_Southern_Hemisphere ( YDim_Southern_Hemisphere,
>> XDim_Southern_Hemisphere )
>> projection :Lambert Azimuthal Equal Area
>> corners :( 1e+51, 1e+51, 1e+51, 1e+51 )
>> long_name :latitude
>> units :degrees_north
>>
>> double GridLon_Southern_Hemisphere ( YDim_Southern_Hemisphere,
>> XDim_Southern_Hemisphere )
>> projection :Lambert Azimuthal Equal Area
>> corners :( 1e+51, 1e+51, 1e+51, 1e+51 )
>> long_name :longitude
>> units :degrees_east
>>
>> byte Flags_SouthernPentad_Southern_Hemisphere (
>> YDim_Southern_Hemisphere, XDim_Southern_Hemisphere )
>> coordinates :GridLat_Southern_Hemisphere, GridLon_Southern_Hemisphere
>> hdfeos_name :Flags_SouthernPentad
>> projection :Lambert Azimuthal Equal Area
>> unsigned :True
>>
>> byte SWE_SouthernPentad_Southern_Hemisphere (
>> YDim_Southern_Hemisphere, XDim_Southern_Hemisphere )
>> coordinates :GridLat_Southern_Hemisphere, GridLon_Southern_Hemisphere
>> hdfeos_name :SWE_SouthernPentad
>> projection :Lambert Azimuthal Equal Area
>> unsigned :True
>> _FillValue :-1
>>
>> double GridLat_Northern_Hemisphere ( YDim_Northern_Hemisphere,
>> XDim_Northern_Hemisphere )
>> projection :Lambert Azimuthal Equal Area
>> corners :( 1e+51, 1e+51, 1e+51, 1e+51 )
>> long_name :latitude
>> units :degrees_north
>>
>> double GridLon_Northern_Hemisphere ( YDim_Northern_Hemisphere,
>> XDim_Northern_Hemisphere )
>> projection :Lambert Azimuthal Equal Area
>> corners :( 1e+51, 1e+51, 1e+51, 1e+51 )
>> long_name :longitude
>> units :degrees_east
>>
>> byte Flags_NorthernPentad_Northern_Hemisphere (
>> YDim_Northern_Hemisphere, XDim_Northern_Hemisphere )
>> coordinates :GridLat_Northern_Hemisphere, GridLon_Northern_Hemisphere
>> hdfeos_name :Flags_NorthernPentad
>> projection :Lambert Azimuthal Equal Area
>> unsigned :True
>> _FillValue :-1
>>
>> byte SWE_NorthernPentad_Northern_Hemisphere (
>> YDim_Northern_Hemisphere, XDim_Northern_Hemisphere )
>> coordinates :GridLat_Northern_Hemisphere, GridLon_Northern_Hemisphere
>> hdfeos_name :SWE_NorthernPentad
>> projection :Lambert Azimuthal Equal Area
>> unsigned :True
>> _FillValue :-1
>>
>>
>>
>> huangwei@ucar.edu <mailto:huangwei@ucar.edu>
>> VETS/CISL
>> National Center for Atmospheric Research
>> P.O. Box 3000 (1850 Table Mesa Dr.)
>> Boulder, CO 80307-3000 USA
>> (303) 497-8924
>>
> ...
>
>
>
> _______________________________________________
> 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 Tue Apr 12 08:50:12 2011

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