Re: Converting axis coordinates

From: Gerald Creager <gerry.creager_at_nyahnyahspammersnyahnyah>
Date: Thu Sep 15 2011 - 12:43:50 MDT

Hit Send too fast. That did, in fact, cure the problem. Thanks for the
gentle reminder.

Now, an aesthetics problem: The plot looks "squished" horizontally. Is
there a suggested aspect ratio I could set?

On 09/15/2011 08:16 AM, Gerald Creager wrote:
> Argh... No, Adam, it doesn't:
>
> printVarSummary(model1Ssh)
>
> Variable: model1Ssh
> Type: float
> Total Size: 37677960 bytes
> 9419490 values
> Number of Dimensions: 3
> Dimensions and sizes: [90] x [261] x [401]
> Coordinates:
> Number Of Attributes: 1
> _FillValue : -9999
>
>
> On 09/14/2011 04:56 PM, Adam Phillips wrote:
>> 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
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Sep 15 12:43:57 2011

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