[ncl-talk] Plotting spatial map of SPI

Mary Haley haley at ucar.edu
Mon Apr 24 06:47:53 MDT 2017


Dawit,

The example you pointed to at:

https://www.dkrz.de/Nutzerportal-en/doku/vis/sw/
ncl/dkrz-ncl-shapefile-german-federal-states-example

shows you exactly how to add shapefile outlines to your plot.

You can also see this example on our own website at:

http://www.ncl.ucar.edu/Applications/NUG.shtml

Look for "NUG_shapefile_plot.ncl".

This example depends on the RANGS data base and the DEU_adm1.shp shapefile,
but you can remove the requirement for the RANGS database by setting
mpDataBaseVersion to "MediumRes" instead of "HighRes"

To get the DEU_adm1.shp shapefile or the adm shapefiles for any country, go
to:

http://www.gadm.org/country

Use the pull down list to select the desired country and make sure the
"File format" is set to "shapefile".

Once you click "OK", you should see a new page with a "Download" button.
Click on this, and it will start downloading the adm files.

You can then use gsn_add_shapefile_polyline to add the shapefile outlines
to your plot.

There are many shapefile examples on our shapefiles page:

http://www.ncl.ucar.edu/Applications/shapefiles.shtml

Good luck,

--Mary


On Fri, Apr 21, 2017 at 2:11 AM, Dawit Abebe <dawit.kan at gmail.com> wrote:

> Dear Dennis,
>
> Thank you so much. I’m aware of the area size, it is only 36 pixels of 0.5
> by 0.5 degree grid. I’m just wondering if it is possible to add the admin
> shape file of the area as shown in examples provided: https://www.dkrz.de/
> Nutzerportal-en/doku/vis/sw/ncl/dkrz-ncl-shapefile-german-
> federal-states-example
>
> If this is possible, I will try this just to visualise the the area in
> terms of administrative boundaries overlaid with the geographic coordinate
> system?
>
> Thanks for you so much again for your patiences with me and going too far
> supporting me.
>
> Dawit
>
>
>
>
>
> On Apr 21, 2017, at 2:01 AM, Dennis Shea <shea at ucar.edu> wrote:
>
> See attached:
>
> you have a very small are.
>
>
>
> On Wed, Apr 19, 2017 at 11:18 AM, Dawit Abebe <dawit.kan at gmail.com> wrote:
>
>> Dear Dennis
>>
>> Thank you so much!
>>
>> I tried with corrected commands but still have a problem with the command:
>>
>>           res at gsnCenterString  = spi&time({192112}), returning the
>> following error message:
>>
>> fatal:Coordinate indexing can not be used with coordinate variables
>> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 100 in
>> file spi12_plotmap_script.ncl
>>
>> What did I miss?
>>
>> Dawit
>>
>>
>>
>> On Apr 19, 2017, at 6:25 PM, Dennis Shea <shea at ucar.edu> wrote:
>>
>> latS = min(prc&lat)
>>   latN = ma
>>
>>
>>
> <dawit.ncl><spi.000001.png><spi.000002.png>
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170424/fca69c2a/attachment-0001.html 


More information about the ncl-talk mailing list