[ncl-talk] Plotting spatial map of SPI

Dawit Abebe dawit.kan at gmail.com
Mon Apr 24 07:16:04 MDT 2017


Thanks Mary! I have been going through your Webinar presentation on NCL and Shapefiles of couple of years back. I’m trying to build my script following that. I may come back if I need help after testing the script. 

Thank you so much

Dawit 


> On Apr 24, 2017, at 3:47 PM, Mary Haley <haley at ucar.edu> wrote:
> 
> Dawit,
> 
> The example you pointed to at:
> 
> https://www.dkrz.de/Nutzerportal-en/doku/vis/sw/ncl/dkrz-ncl-shapefile-german-federal-states-example <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
>  <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 <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 <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 <mailto: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 <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 <mailto: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 <mailto: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 <mailto: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/b5000b48/attachment.html 


More information about the ncl-talk mailing list