Re: Converting axis coordinates

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

You should see something like this:

Variable: model1Ssh
<snip>
Coordinates:
             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=model1->sea_surface_height_above_sea_level
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..
Adam

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 <gerry.creager@tamu.edu>
Reply-To: gerry.creager@tamu.edu
Organization: AATLT, Texas A&M University
To: ncl-talk@ucar.edu <ncl-talk@ucar.edu>

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 -- gerry.creager@tamu.edu
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                                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
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk