Re: Graph doesn't show the contours

Hi Kemal,

The main problem is that you are not using your masking variable 'lulc' correctly. You have a comment that states that "ocean=0,land=1,sea_ice=2", but if you actually look at the data in the variable you will see that it ranges from 1-24, and if you contour the mask variable itself you can see that apparently the ocean is all equal to 16 while the values over land are all the other numbers. So you must change the mask calls to reflect this:

land_only = mask(tmp,(lulc.ne.16),True)
ocean_only = mask(tmp,(lulc.eq.16),True)

The other less serious problem is that you tried to use "BoxCenters" for the LabelBar. But the contour levels that are the input to the LabelBar represent the dividing lines between the boxes. Therefore using "BoxCenters" gives a misleading impression of the data and leads to there not being enough labels (hence the last label is set to "Label_16" rather than an actual value).
 
I am attaching an edited version of your slp.ncl that gives a reasonable plot of the land-only data.
 -dave

On Oct 19, 2010, at 5:56 PM, Kemal Gurer wrote:

> Hello,
>
> I am trying to generate contour plots for sea level pressure and surface pressure using ncl, but the plots are almost empty. While I have been successful with temperature data that I used before, I am not successful with these two data. I placed the data and the scripts to:
>
> ftp://eos.arb.ca.gov/pub/outgoing/kgurer/ncl_error.tar.gz
>
> I would appreciate any help you may give. Thank you.
>
> Kemal.
> _______________________________________________
> 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 Oct 20 13:19:27 2010