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

To: Wei Huang <huangwei@ucar.edu>

OK,

Thanks to Wei and Dennis.

I am doing the following:
eos_file=addfile("AMSR_E_L3_5DaySnow_V09_20050126.hdf.he2", "r")
bdata=eos_file->SWE_NorthernPentad_Northern_Hemisphere
fdata = tofloat(bdata)
fdata@lon2d=eos_file->GridLon_Northern_Hemisphere
fdata@lat2d=eos_file->GridLat_Northern_Hemisphere

xwks=gsn_open_wks("pdf","AMSR_E_L3_5DaySnow_V09_20050126_SWE_NH")

res=True
res@cnFillOn=True
res@gsnSpreadColors=True
res@cnFillMode="RasterFill"
res@cnMissingValFillPattern = 0 ; missing value pattern is set to
"SolidFill"
res@cnMissingValFillColor=0; white color for missing values

res@mpProjection = "LambertEqualArea" ; LAMAZ
res@mpCenterLatF = 90.0; == North Pole
res@mpCenterLonF = 0.0
res@mpMinLatF = min(fdata@lat2d)

plot=gsn_csm_contour_map(xwks,fdata,res)

Using the res@mpProjection = "LambertEqualArea" options, I no longer see
any errors,
But the contour plot is blank (all white), but for the continental
outlines and color bar.

Suggestions as to how to coerce gsn_csm_contour_map into finding contours?

Thanks,

-Stewart

On 4/12/2011 3:14 PM, Wei Huang wrote:
> Stewart,
>
> You cannot use this udata to plot functions.
>
> You may need to do something like:
>
> fdata = tofloat(udata)
>
> then:
>
> plot=gsn_csm_contour_map(xwks,fdata,res)
>
> Wei
>
> 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
Received on Fri Apr 15 00:33:01 2011