[ncl-talk] printing spi
Dawit Abebe
dawit.kan at gmail.com
Mon Apr 10 07:01:30 MDT 2017
Dear Guido,
What a shame of myself!
Thanks Guido. I did the correction and went through smoothly but the last command failed to return the graph (4 one side and 3 on the other side in one plot). What do you that I missed to complete this?
Thanks again
…………………………………...
ncl 109> do k=0,klen-1
ncl 110> res at xyLineColors = xyLineColors(k)
ncl 111> res at gsnCenterString = "len="+klen(k)
ncl 112> plt(k) = gsn_csm_xy (wks, yrfrac, spi(k,:), res)
ncl 113> end do
ncl 114> resP = True
ncl 115> resP at txString = "SPI: Southern Ethiopia 1901-2014"
ncl 116> resP at gsnMaximize = True
ncl 117> gsn_panel (wks,plt,(/4,2/), resP)
………………….
> On Apr 10, 2017, at 3:39 PM, Guido Cioni <guidocioni at gmail.com> wrote:
>
> There is a typo in your code. It should be xyLineColors and not xyLineColours.
> See here https://www.ncl.ucar.edu/Document/Graphics/Resources/xy.shtml <https://www.ncl.ucar.edu/Document/Graphics/Resources/xy.shtml>.
>
> (Look closely at the warning/error messages that you are receiving: most of the times they can help you a lot in identifying the problem)
>
> Hope that helps
>
> Guido Cioni
> http://guidocioni.altervista <http://guidocioni.altervista/>.org
>
>> On 10 Apr 2017, at 14:26, Dawit Abebe <dawit.kan at gmail.com <mailto:dawit.kan at gmail.com>> wrote:
>>
>> Dear Dennis and all
>>
>> I’m sorry for bothering you all but tried to do xy plot of spi following example 1but failed. You hinted in your email below that I should change the spi array but could not figure out how to do thi. I tried the following but could not able to continue after some step. I have gone through the documentations and examples but failed to fix my problem.
>>
>> ………………………………
>>
>> Variable: spi
>> Type: double
>> Total Size: 3151872 bytes
>> 393984 values
>> Number of Dimensions: 4
>> Dimensions and sizes: [len | 8] x [time | 1368] x [lat | 4] x [lon | 9]
>> Coordinates:
>> len: [1..36]
>> time: [ 380..41987]
>> lat: [3.75..5.25]
>> lon: [36.75..40.75]
>> Number Of Attributes: 2
>> long_name : SPI
>> _FillValue : -999
>>
>> Plot Parameter
>> ncl 55> yyyymm = f->time
>> ncl 56> yyyymm = yyyymm/100yrStrt = 1901
>> ncl 57> yrLast = 2014
>> ncl 58> nyear = yrLast-yrStrt+1
>> ncl 59> yyyymm = yyyymm_time(yrStrt, yrLast, "integer")
>> ncl 60> yrfrac = (/yyyymm_to_yyyyfrac(yyyymm, 0.0)/)
>> ncl 61>
>> ncl 62> wks = gsn_open_wks ("ps","spi")
>> ncl 63> res = True
>> ncl 64> res at gsnDraw=False
>> ncl 65> res at gsnFrame=False
>> ncl 66> res at vpHeightF = 0.4
>> ncl 67> res at vpWidthF = 0.8
>> ncl 68> res at vpXF = 0.1
>> ncl 69>
>> ncl 70> res at trYMinF = -3.0
>> ncl 71> res at trYMaxF = 3.0
>> ncl 72> res at gsnYRefLine = 0.0
>> ncl 73> res at xyMonoDashPattern = True
>> ncl 74> res at xyLineThicknessF = 1
>> ncl 75>
>> ncl 76> plt = new (klen, "graphic")
>> ncl 77> xyLineColors = (/"black","red","blue","green","brown","magenta","grape"/)
>> ncl 78> do k=0,klen-1
>> ncl 79> res at xyLineColours = xyLineColours(k)
>> ncl 80> res at gsnCenterString = "len="+klen(k)
>> ncl 81> plt(k) = gsn_csm_xy (wks, yrfrac, spi(k,:), res)
>> ncl 82> end do
>> fatal:Undefined identifier: (xyLineColours) is undefined, can't continue
>> fatal:["Execute.c":8638]:Execute: Error occurred at or near line 79
>> ……………………………
>>
>> Thanks for your help
>>
>> Dawit
>>
>>
>>
>>
>>> On Apr 6, 2017, at 9:00 PM, Dennis Shea <shea at ucar.edu <mailto:shea at ucar.edu>> wrote:
>>>
>>> I do not understand the question.
>>>
>>> https://www.ncl.ucar.edu/Applications/ <https://www.ncl.ucar.edu/Applications/>
>>>
>>> Under Data Analysis: Click 'Standardized Precipitation Index'
>>>
>>> See: Example 1
>>>
>>> You need to make the same array changes as in the print.
>>>
>>> EG: spi(0,:,)==> spi(0,:,nl,ml)
>>>
>>>
>>> On Thu, Apr 6, 2017 at 1:22 AM, Dawit Abebe <dawit.kan at gmail.com <mailto:dawit.kan at gmail.com>> wrote:
>>> Hello All,
>>>
>>> Thanks to Dennis, I am able to appreciate the power of NCL computing SPI. My next step is to produce
>>> (1) a time serious XY plot for seasonal (SPI-3 ending in May and ending in November), annual (SPI-12) and multi annual (SPI-24 and SPI-36 months) and
>>> (2) Shape map of selected annual SPI with extremely sever drought events to examine the spatial characteristics of the drought.
>>>
>>> I’m not asking for help before I try by my self (thanks Dennis for the push to do this). However, I need a quick comment how I should proceed from the last step producing the SPI output, which was printing:
>>>
>>> print(yyyymm+sprintf("%8.2f", prc(:,nl,ml)) \
>>> +sprintf("%8.2f", spi(0,:,nl,ml))+sprintf("%8.2f", spi(1,:,nl,ml)) \
>>> +sprintf("%8.2f", spi(2,:,nl,ml))+sprintf("%8.2f", spi(3,:,nl,ml)) \
>>> +sprintf("%8.2f", spi(4,:,nl,ml))+sprintf("%8.2f", spi(5,:,nl,ml)) \
>>> +sprintf("%8.2f", spi(6,:,nl,ml)))
>>>
>>> Do I need to go through this and then start with plot parameters, or I can continue with plot parameters without the print command?
>>>
>>> Thank you so much
>>> Dawit
>>>
>>>
>>>
>>>> On Apr 5, 2017, at 4:43 PM, Dennis Shea <shea at ucar.edu <mailto:shea at ucar.edu>> wrote:
>>>>
>>>> The error mesages is telling you the exact problem:
>>>> fatal:Plus: Number of dimensions do not match, can't continue
>>>>
>>>> spi has 4-dimensions:
>>>> [len | 7] x [time | 1368] x [lat | 4] x [lon | 9]
>>>>
>>>> prc has 3-dimensions
>>>> [time | 1368] x [lat | 4] x [lon | 9]
>>>>
>>>> yyyymm has 1-dimension
>>>> [time | 1368]
>>>>
>>>> ==============
>>>> The print statement you used was copied from:
>>>> https://www.ncl.ucar.edu/Document/Functions/Built-in/dim_spi_n.shtml <https://www.ncl.ucar.edu/Document/Functions/Built-in/dim_spi_n.shtml>
>>>> Example 2
>>>> That was working upon a single time series prc[*]
>>>>
>>>>
>>>>
>>>> You have a 4-D array. *You* must make the appropriate changed to handle your situation.
>>>> The eamples are not generic to all dimensions.
>>>>
>>>> nlat=4
>>>> mlon=9
>>>>
>>>> do nl=0,nlat-1
>>>> do ml=0,mlon-1
>>>> print("------------------------------------------------")
>>>>
>>>> print(yyyymm+sprintf("%8.2f", prc(:,nl,ml)) \
>>>> +sprintf("%8.2f", spi(0,:,nl,ml))+sprintf("%8.2f", spi(1,:,nl,ml)) \
>>>> +sprintf("%8.2f", spi(2,:,nl,ml))+sprintf("%8.2f", spi(3,:,nl,ml)) \
>>>> +sprintf("%8.2f", spi(4,:,nl,ml))+sprintf("%8.2f", spi(5,:,nl,ml)) \
>>>> +sprintf("%8.2f", spi(6,:,nl,ml)))
>>>>
>>>> end do
>>>> en do
>>>>
>>>> On Wed, Apr 5, 2017 at 2:20 AM, Dawit Abebe <dawit.kan at gmail.com <mailto:dawit.kan at gmail.com>> wrote:
>>>> Hello all,
>>>>
>>>> Kindly help me to fix my problem with printing SPI. I pasted below first the return of printVarSummary(spi), which seems okay. However, my attempt to print the output (text in blue) did not work. What was wrong?
>>>>
>>>> Greatly appreciate your help.
>>>>
>>>> Dawit
>>>>
>>>> …………………………………….
>>>>
>>>> copy_VarCoords(prc,spi(0,:,:,:))
>>>> spi at long_name = "SPI"
>>>> spi!0 = "len"
>>>> spi&len = len
>>>> printVarSummary(spi)
>>>>
>>>>
>>>> Variable: spi
>>>> Type: double
>>>> Total Size: 2757888 bytes
>>>> 344736 values
>>>> Number of Dimensions: 4
>>>> Dimensions and sizes: [len | 7] x [time | 1368] x [lat | 4] x [lon | 9]
>>>> Coordinates:
>>>> len: [1..36]
>>>> time: [ 380..41987]
>>>> lat: [3.75..5.25]
>>>> lon: [36.75..40.75]
>>>> Number Of Attributes: 2
>>>> long_name : SPI
>>>> _FillValue : -999
>>>>
>>>>
>>>> yyyymm = cd_calendar(f->time, -1)
>>>>
>>>> print(yyyymm+sprintf("%8.2f", prc) \
>>>> +sprintf("%8.2f", spi(0,:))+sprintf("%8.2f", spi(1,:)) \
>>>> +sprintf("%8.2f", spi(2,:))+sprintf("%8.2f", spi(3,:)) \
>>>> +sprintf("%8.2f", spi(4,:))+sprintf("%8.2f", spi(5,:)) \
>>>> +sprintf("%8.2f", spi(6,:)))
>>>> fatal:Plus: Number of dimensions do not match, can't continue
>>>> fatal:["Execute.c":8638]:Execute: Error occurred at or near line 37
>>>>
>>>>
>>>>
>>>> _______________________________________________
>>>> ncl-talk mailing list
>>>> ncl-talk at ucar.edu <mailto:ncl-talk at ucar.edu>
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk <http://mailman.ucar.edu/mailman/listinfo/ncl-talk>
>>>>
>>>>
>>>
>>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu <mailto:ncl-talk at ucar.edu>
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170410/2d3c8896/attachment.html
More information about the ncl-talk
mailing list