Re: blank plot

Hi,
Thanks for your reply and I tried. Only one more blank plot appears in phf
file and no other things change.

The lat/lon in this file is 4060*2780. I'm not sure if it matters.

Thanks

Lizzie

> Hi,
>
> Add the following after plot=gsn_csm_contour_map_polar(xwks,data,res):
>
> draw(plot)
> frame(xwks)
>
> -------------------------------------------------------------
> Brandon Fisel
>
>
> On Wed, Sep 15, 2010 at 13:36, <ziying-li@hdfgroup.org> wrote:
>
>> Hi,
>> This is a follow up case for below example:
>> http://www.ncl.ucar.edu/Support/talk_archives/2010/1904.html
>>
>> I successfully read out lat and lon from that ASCII file. The size
>> matches
>> and the value looks reasonable.
>>
>> But the final pdf plot I got is blank. You can see title, map and color
>> bar well except data.
>>
>> NCL code:
>> ------------------------------------------------------------
>> load "$NCARG_ROOT/lib/ncarg/nclex/gsun/gsn_code.ncl"
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>>
>> begin
>> eos_file=addfile("MYD02HKM.A2010031.0035.005.2010031183706.hdf.he2",
>> "r")
>>
>> data_unscaled = eos_file->EV_500_RefSB_MODIS_SWATH_Type_L1B(0,:,:) ; get
>> the clues from print(eos_file)
>> data_unscaled@_FillValue = 65535h
>> data = 3.4915833E-5 * ( data_unscaled - (-0.0) )
>>
>> nlon=4060
>> nlat=2708
>>
>>
>> lat=asciiread("lat_MYD02HKM.A2010031.0035.005.2010031183706.output",(/nlon,nlat/),"float")
>>
>> lon=asciiread("lon_MYD02HKM.A2010031.0035.005.2010031183706.output",(/nlon,nlat/),"float")
>>
>> data@lat2d=lat
>> data@lon2d=lon
>>
>>
>> xwks=gsn_open_wks("pdf","MYD02HKM.A2010031.0035.005.2010031183706_EV_500_RefSB.ncl")
>>
>> setvalues NhlGetWorkspaceObjectId() ; make maximum filesize larger
>> "wsMaximumSize" : 200000000
>> end setvalues
>>
>> res=True
>> res@cnFillOn=True
>> res@gsnMaximize=True
>> res@gsnPaperOrientation = "portrait"
>> res@cnLinesOn=False
>> res@cnLineLabelsOn = False
>> res@gsnSpreadColors=True
>> res@cnFillMode="RasterFill"
>> res@lbOrientation="vertical"
>> res@cnMissingValFillPattern = 0
>> res@cnMissingValFillColor=0
>> res@lbLabelAutoStride = True
>> res@sfXArray=lon
>> res@sfYArray=lat
>> res@gsnAddCyclic=False
>>
>> gsn_define_colormap(xwks,"BlAqGrYeOrReVi200")
>> res@tiMainString = "MYD02HKM.A2010031.0035.005.2010031183706.hdf"
>> res@gsnLeftString = "Field name: EV_500_RefSB"
>> res@gsnRightString ="Units: none"
>>
>> ;res@mpLimitMode = "LatLon"
>> ;res@mpMinLatF = min(data@lat2d) ;
>> ;res@mpMaxLatF = max(data@lat2d) ;
>> ;res@mpMinLonF = min(data@lon2d) ;
>> ;res@mpMaxLonF = max(data@lon2d) ;
>>
>> res@gsnPolar = "NH" ; specify the hemisphere
>> plot=gsn_csm_contour_map_polar(xwks,data,res)
>>
>>
>> delete(plot)
>> delete(xwks)
>> delete(data)
>> delete(res)
>> delete(eos_file)
>>
>> end
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Sep 15 11:54:42 2010