Re: x-axis tick mark labels in contour plot (Adam Phillips)

Hi Adam,

Thank you very much for your help. It is very much appreciated.

Starry

________________________________________
From: ncl-talk-bounces@ucar.edu [ncl-talk-bounces@ucar.edu] On Behalf Of ncl-talk-request@ucar.edu [ncl-talk-request@ucar.edu]
Sent: Thursday, January 26, 2012 1:00 PM
To: ncl-talk@ucar.edu
Subject: ncl-talk Digest, Vol 98, Issue 37

Send ncl-talk mailing list submissions to
        ncl-talk@ucar.edu

To subscribe or unsubscribe via the World Wide Web, visit
        http://mailman.ucar.edu/mailman/listinfo/ncl-talk
or, via email, send a message with subject or body 'help' to
        ncl-talk-request@ucar.edu

You can reach the person managing the list at
        ncl-talk-owner@ucar.edu

When replying, please edit your Subject line so it is more specific
than "Re: Contents of ncl-talk digest..."

Today's Topics:

   1. Re: GRIB2 (Mary Haley)
   2. Re: Satellite cloud map (Mary Haley)
   3. Re: x-axis tick mark labels in contour plot (Adam Phillips)

----------------------------------------------------------------------

Message: 1
Date: Thu, 26 Jan 2012 10:33:12 -0700
From: Mary Haley <haley@ucar.edu>
Subject: Re: GRIB2
To: "Chartier, Alex" <alex.chartier@metoffice.gov.uk>
Cc: ncl-talk@ucar.edu
Message-ID: <24D98413-2DA6-4ABF-9079-CB66F62CAF58@ucar.edu>
Content-Type: text/plain; charset=us-ascii

Hi Alex,

I don't know if anybody contacted you offline, but can you provide your GRIB2 file?

ftp ftp.cgd.ucar.edu
anonymous
Use your email address for the password
cd incoming
put ...
quit

I will need to know the exact name of the file in order to retrieve it.

--Mary

On Jan 24, 2012, at 10:03 AM, Chartier, Alex wrote:

> Hi,
>
> I recently installed NCL in the hope of turning some GRIB2 files into
> netCDF, but I came up against the following error:
> warning: NclGRIB2: codetable file
> "/data/local/frch/NCL/lib/ncarg/grib2_codetables/ncep/4/4.2.4.2.table"
> not a valid GRIB2 code table.
>
> There is no table '4.2.4.2.table', although the directory exists and has
> other tables in it. Where do I go to get the GRIB2 code table?
> The demo worked fine when I installed NCL.
>
> Thanks,
>
> Alex Chartier
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

------------------------------

Message: 2
Date: Thu, 26 Jan 2012 10:36:28 -0700
From: Mary Haley <haley@ucar.edu>
Subject: Re: Satellite cloud map
To: Michael Notaro <mnotaro@wisc.edu>
Cc: ncl-talk NCL <ncl-talk@ucar.edu>
Message-ID: <E4EF4AEC-473A-44F1-AB98-24471D7DEF80@ucar.edu>
Content-Type: text/plain; charset=us-ascii

Hi Mike

It depends. If you are using NCL V6.0.0 or earlier, and this image has more than 256 colors, then no, you wouldn't be able to reproduce it.

In our next release of NCL, V6.1.0, we hope to have support for larger color maps, and something like this could be possible.

Do you have the data for this image, and could you provide it? We can try it out here on our 6.1.0 test version.

--Mary

On Jan 25, 2012, at 3:54 PM, Michael Notaro wrote:

> Is it possible in NCL to produce a satellite cloud
> map similar to the attached plot?
>
> Thanks, Mike
>
> <cloud.png>_______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

------------------------------

Message: 3
Date: Thu, 26 Jan 2012 11:33:50 -0700
From: Adam Phillips <asphilli@ucar.edu>
Subject: Re: x-axis tick mark labels in contour plot
To: ncl-talk@ucar.edu
Message-ID: <4F219C8E.8020306@ucar.edu>
Content-Type: text/plain; charset="iso-8859-1"

Hi Starry,
Your data file didn't come through so I filled your data array with
random values and ran your script.

In the future when writing to ncl-talk, please mention any errors you
are getting from NCL. In your case you are getting the following error
when you try to set the tmXBLabels:

warning:_NhlCreateSplineCoordApprox: Attempt to create spline
approximation for X axis failed: consider adjusting trXTensionF value
^Mwarning:IrTransInitialize: error creating spline approximation for
trXCoordPoints; defaulting to linear

What this error means is that NCL is having trouble figuring out how to
handle your x-axis coordinate values. So I started to focus on your
input latitudes.

(This turns out not to be the cause of your issue, but you were setting
this:
   lat@units="degrees North"
   lon@units="degrees East"
You were right to set the units attribute of your lat/lon arrays, but
NCL expects them to be one of a few possibilities, the most common one
being:
   lat@units="degrees_north"
   lon@units="degrees_east"

(See this page for more information:
http://www.ncl.ucar.edu/Document/Graphics/error_msg.shtml)

I changed the units attribute and I started getting this warning:
get_lat_values: Warning: Your latitude values do not fall between -90
and 90 inclusive.
You will not get 'nice' latitude labels.

This of course made me really zero in on your latitudes..)

So I printed out your latitudes and longitudes to see what they looked like:
15.0068N 76.9888E
15.0068N 76.9888E
15.0068N 76.9888E
15.0068N 76.9888E
15.0068N 76.9888E
15.0068N 76.9888E
15.0068N 76.9888E
15.0068N 76.9888E
15.0068N 76.9888E
15.0068N 76.9888E
15.0068N 76.9888E
15.0068N 76.9888E
15.0068N 76.9888E
15.0068N 76.9888E
15.0068N 76.9888E
15.0516N 76.9788E
15.0516N 76.9788E
15.0516N 76.9788E
15.0516N 76.9788E
15.0516N 76.9788E
15.0516N 76.9788E
15.0516N 76.9788E
15.0516N 76.9788E

You're plotting latitude vs. height, so the longitudes aren't the issue.
However, as you can see, you have multiple instances of the same
latitude, all with different data values. When you tell NCL to draw a
tickmark at 15.0068N, it doesn't know which point to draw them at, the
1st instance or the 15th instance of 15.0068N... That is definitely the
cause of your X-axis labeling problems, and points to a possible issue
with your data. I'm not sure if the data is the same for each instance
of the same latitude value as I don't have your data. You can get an
idea of that by inserting:
print(lat+"N "+lon+"E "+vfm_swap(0,:))
Good luck,
Adam

On 01/26/2012 08:44 AM, Manoharan, Vani Starry wrote:
> Hi Adam,
>
> I tried runnning the program by not specifying the XBValues/XBLabels. But it generates numbers on the x-axis. I am attaching my program and the data files used in the program with this mail and the following mail (as the files are quite big).
>
> Thank you.
> Starry
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

--
______________________________________________________________
Adam Phillips                                asphilli@ucar.edu
NCAR/Climate and Global Dynamics Division       (303) 497-1726
P.O. Box 3000
Boulder, CO 80307-3000    http://www.cgd.ucar.edu/cas/asphilli
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20120126/7b761e6b/attachment.html
------------------------------
_______________________________________________
ncl-talk mailing list
ncl-talk@ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
End of ncl-talk Digest, Vol 98, Issue 37
****************************************
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk