Re: Converting axis coordinates

From: Adam Phillips <asphilli_at_nyahnyahspammersnyahnyah>
Date: Wed Sep 14 2011 - 15:56:26 MDT

Hi Gerry,
Does your model1Ssh array have associated latitude/longitude coordinate
variables? You can see if it does by using printVarSummary:

You should see something like this:

Variable: model1Ssh
             time: [17067072..17618592]
             lat: [90..-90]
             lon: [ 0..357.5]
(Last 3 lines show the coordinate variables.)

If coordinate variables aren't present, you need to assign them. Looking
at your script though, I see you are doing this:
model1Ssh=100 * model1->sea_surface_height_above_sea_level
This will result in only the data being passed from the right-hand side
of the equal sign to the left-hand side, and not the coordinate
variables associated with sea_surface_height_above_sea_level array.

I'd recommend changing this:
model1Ssh=100 * model1->sea_surface_height_above_sea_level
to this:
model1Ssh = model1Ssh*100

I would venture to guess that that will fix your coordinate variables,
and therefore the lats/lons shown on your plot...

As far as the land mask issue: As you are plotting SSH there is no data
over land, so those areas are getting color filled white. You need to
tell NCL to color fill the map after the contour field gets color
filled. Try setting res@mpFillDrawOrder = "PostDraw", and set
res@mpLandFillColor = "gray70", as you are adding gray70 to your
colormap, and it's best to exactly specify the color your want NCL to use.

Let the group know if that doesn't solve everything..

On 09/14/2011 03:34 PM, Gerald Creager wrote:
-------- Original Message --------
Subject: Converting axis coordinates
Date: Wed, 14 Sep 2011 16:34:13 -0500
From: Gerald Creager <>
Organization: AATLT, Texas A&M University
To: <>

I'm creating an image of a model output (attached), and the X and Y
coordinates are currently in model indices. I'd like to have them
represented as Lon/Lat.

Also, I tried to create a gray landmask but it didn't seem to work. If
someone's got a suggestion for that, I'd appreciate it. The ncl script
is also attached.

Thanks, gerry

Gerry Creager --
Texas Mesonet -- AATLT, Texas A&M University
Cell: 979.229.5301 Office: 979.458.4020 FAX: 979.862.3983
Office: 1700 Research Parkway Ste 160, TAMU, College Station, TX 77843
Adam Phillips                      
NCAR/Climate and Global Dynamics Division       (303) 497-1726
P.O. Box 3000				
Boulder, CO 80307-3000
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
Received on Wed Sep 14 15:56:32 2011

This archive was generated by hypermail 2.1.8 : Fri Sep 16 2011 - 11:24:25 MDT