Re: visualize cloudsat data

Yin,

For starters, you want to use "gsn_csm_contour" and not "gsn_contour".
The "gsn_contour" routine will not look for coordinate arrays;
"gsn_csm_contour" will.

Dennis asked offline for you to send the file. He'll take look at
it.

In the future, however, it would be helpful if you provide us with
more information about your data. For example, the output from
statements like:

   print(in)
   ...
   printVarSummary(yyy)
   ...
   dnames=getfilevardims(in,"ReceivedEchoPowers_1B_CPR")
   print(dnames)

would be very useful. Once we have your dataset, we can do this
ourselves.

--Mary

On Thu, 16 Jul 2009, yinjinfang wrote:

>
> Hi,dear;
> I want to visualize cloudsat data with ncl 5.1.0. I am confusing when I run the following script.
> It shows a wrong plot. I think the problem arises from the cootdinates setting (default).
> I attempt to set the cootdinates, but fail. Because the cloudsat data is a two-demension datasat (a profile), and I can not obtain the latitude and/or longitude varibles for
> x-cootdinate. How can I get the reasonable cootdinates?
> Thanks in advance!
>
>
> -----------------------------------------------------------
>
> in=addfile("./data/2009168171301_16694_CS_1B-CPR_GRANULE_P_R04_E02.hdf.hdfeos","r")
> printVarSummary(in)
> dnames=getfilevardims(in,"ReceivedEchoPowers_1B_CPR")
> print(dnames)
> xwks=gsn_open_wks("X11","yjf")
> yyy=in->ReceivedEchoPowers_1B_CPR
> lon=in->Longitude_1B_CPR
> ;resources = True
> ;resources_at_cnFilleOn = True
> ;resources_at_sfXCStartV= lon
> plot=gsn_contour(xwks,yyy(0:1196,:),resources)
> end
>
>
> -------------------------------------------------------------
> Yin
> 2009-07-15
> 2009-07-15
>
>
>
> yinjinfang
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jul 15 2009 - 22:34:10 MDT