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