# Re: plotting climate divisions

From: Mary Haley <haley_at_nyahnyahspammersnyahnyah>
Date: Tue, 25 Sep 2007 11:19:07 -0600 (MDT)

Wendi,

I *think* there's a way to do this. I need to get some information
from our developers first, and then I'll try to provide
a PyNGL example.

--Mary

On Tue, 25 Sep 2007, Wendi Madsen wrote:

> Hi Mary
>
> That helps, but how do I plot data to each of the divisions. I have
> data that looks like this:
> 1 1 8.110001
> 1 2 7.546665
> 1 3 7.916664
> 1 4 7.1066665
> 1 5 6.209999
> 1 6 6.8200073
> 1 7 6.6500015
> 1 8 5.813339
> 2 1 -0.49666976
> 2 2 2.4066658
> 2 3 1.3466682
>
> where the first column is the state number, 1 = Alabama, and 2 =
> Arizona. The second column is the division number and the third column
> is the temperature anomoly. Is there some way to use this data or do I
> need the lat and lon for each division in order to plot it?? The
> polyg_2 script uses the lat and lon for the x and y in the call to
> gsn_polygon which I assume does the data plotting.
>
>
> Wendi
>
> Mary Haley wrote:
>>
>> Hi Wendi,
>>
>> We have climate divisions built into NCL and PyNGL. The polyg_2.ncl
>> example you are looking at was written before this feature was added
>> (I will update the example to reflect this).
>>
>> Anyway, I added a new example in the NCL "maps onlys" applications
>> page to show how to draw the built-in climate divisions:
>>
>> http://www.ncl.ucar.edu/Applications/maponly.shtml#ex15
>>
>> To get these divisions in a PyNGL plot, I think all you need are
>> the following resources:
>>
>> res.mpDataSetName = "Earth..3"
>> res.mpDataBaseVersion = "Ncarg4_1"
>> res.mpOutlineBoundarySets = "AllBoundaries"
>>
>> Let me know if you need something more than this.
>>
>> --Mary
>>
>>
>> On Mon, 24 Sep 2007, Wendi Madsen wrote:
>>
>>> Hi
>>>
>>> I am trying to use polyg_2.ncl as an example to plot climate divisions
>>> for the US temp anomolies. I have a text file that has 3 columns which
>>> are state number, 1 to 48, division number, and temp anolomies values.
>>> The example uses lat and lon for the x and y values for the gsn_polygon
>>> call. Is there any way to use the state number/division number or do I
>>> need the lat and lon for each division??
>>>
>>>
>>> Wendi
>>> --
>>> Systems Analyst
>>> Physical Sciences Division
>>> NOAA/OAR/Earth System Research Laboratory
>>> R/PSD0
>>> Boulder, Co 80305
>>> voice 303 497-6869
>>> fax: 303 497-6020
>>>
>>> _______________________________________________
>>> pyngl-talk mailing list
>>> pyngl-talk_at_ucar.edu
>>> http://mailman.ucar.edu/mailman/listinfo/pyngl-talk
>>>
>>
>
> --
> Systems Analyst
> Physical Sciences Division
> NOAA/OAR/Earth System Research Laboratory
> R/PSD0
> Boulder, Co 80305
> voice 303 497-6869
> fax: 303 497-6020