[ncl-talk] SPI/SPI DROUGHT CLASSIFICATION

Oluwaseun Ilori wilseun at gmail.com
Tue Apr 16 09:56:08 MDT 2019


Thanks so much and I've attached my data ncdump -h .nc output here but
somehow lost at the original dimension of the code "spi(nrun,time,lat,lon)"
that is leading to an error.

Thanks so much again.

On Tue, Apr 16, 2019 at 3:34 AM Dennis Shea <shea at ucar.edu> wrote:

> When you pose a 'help' question, you should *always* provide the output
> from:
>
> %> ncdump -h SPI,nc
>
> or, equivalently,
>
> %> ncl_filedump SPI.nc
>
> or
>
>   f = addfile(....")
>   spi = f->SPI
>   printVarSummary(spi)
>   pritMinMax(spi,0)
>
> otherwise, people are guessing.
> =================
>    nrun  = 3                     ; 3, 6, 12
>    nmos = 12                  ; number of months/year
>    NYRS  = 90
>    NTIM =NYRS*nmos   ; TOTAL # of months
>
>    NYR30 = 30
>    NTM30 = NYR30*nmos
>    NSG30 = 3       ; # of 30-year segments
>
> Original dimensions:
>
>    spi(nrun,time)              ; station or location
> or
>    spi(nrun,time,lat,lon)    ; grid
>
> ;-----
>    KNT = new((/nrun,NSG30/),"integer")
> or
>    KNT = new((/nrun,NSG30,nlat,mlon/),"integer")
>
>    spiCrit = -1.0   ;
>    ntStrt  = 0
>    ntLast = NTM30 -1        ; 360 months [30 years];  -1  is for NCL index
>
>    do ns=0,NSG30-1
>         KNT(:, ns) = dim_num_n( spi(:,*ntStrt:ntLast*).le.spiCrit, 1)
> or
>        KNT(:, ns,:,:) = dim_num_n( spi(:,*ntStrt:ntLast,*:,:).le.spiCrit,
> 1)
>        ntStrt = ntStrt + NSG30
>         ntLast = ntLast + NSG30
>     end do
>
>   printVarSummary(KNT)
>   printMinMax(KNT,0)
>   print("=====")
>
>
>
> On Mon, Apr 15, 2019 at 4:35 AM Oluwaseun Ilori <wilseun at gmail.com> wrote:
>
>> Thanks, I'll check them out.
>>
>> On Mon, Apr 15, 2019 at 2:39 AM Md. Jalal Uddin <dmjalal90 at gmail.com>
>> wrote:
>>
>>> You can look at two functions that might help you.
>>> count_unique_values
>>>
>>> https://www.ncl.ucar.edu/Document/Functions/Built-in/count_unique_values.shtml
>>>
>>>
>>> where
>>>  https://www.ncl.ucar.edu/Document/Functions/Built-in/where.shtml
>>>
>>> Good luck!
>>>
>>> On Mon, Apr 15, 2019 at 2:35 AM Oluwaseun Ilori <wilseun at gmail.com>
>>> wrote:
>>>
>>>> Good day,
>>>> Please I need help to perform drought classification base on McKee et.,
>>>> 1993. I'm having SPI.nc file that contain 3, 6, and 12months time scales
>>>> for the period of 90 years (2011-2100). I want to use spi < - 1.0 to
>>>> determine drought frequency for the first 30 years, second 30 years and the
>>>> third 30 years and save it as nc file. I've been looking for any function
>>>> that I could use to perform this task but none yet. Pleas could someone put
>>>> me through cos I'm not so good in this but I love the output of ncl.
>>>> Thanks in advance
>>>> _______________________________________________
>>>> ncl-talk mailing list
>>>> ncl-talk at ucar.edu
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>>
>>>
>>> --
>>> *Md. Jalal Uddin*
>>> MSc in Applied Meteorology (English Language)
>>> Nanjing University of Information, Science and Technology, China
>>> B.Sc. in Disaster Management (Hons.)
>>> Patuakhali Science and Technology University, Bangladesh.
>>> Cell: +8613260859092, +8801792052662
>>> Web: www.dmjalal90.weebly.com
>>> Facebook: jalal.hossen.39
>>> LinkedIn: https://bd.linkedin.com/in/md-jalal-uddin-80a026b0
>>> Twitter: dmjalal90
>>> Skype: dmjalal90
>>>
>>
>>
>> --
>> *OLUWASEUN ILORI*
>> *B.Tech(Hons) in Meteorology*
>> *M.Tech in Meteorology and climate science*
>>
>> *Phone: +2347031000687 +2348069154213*
>> *Federal University of Technology Akure. Nigeria.*
>> *Facebook: *http://facebook.com/hilorywilson
>> <https://www.facebook.com/hilorywilson>
>> linkedln: https://www.linkedin.com/in/wilson-ilori-447952108
>> Twitter: HiloryW
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>

-- 
*OLUWASEUN ILORI*
*B.Tech(Hons) in Meteorology*
*M.Tech in Meteorology and climate science*

*Phone: +2347031000687 +2348069154213*
*Federal University of Technology Akure. Nigeria.*
*Facebook: *http://facebook.com/hilorywilson
<https://www.facebook.com/hilorywilson>
linkedln: https://www.linkedin.com/in/wilson-ilori-447952108
Twitter: HiloryW
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190416/1f3d6d31/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screenshot from 2019-04-16 15-06-21.png
Type: image/png
Size: 155653 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190416/1f3d6d31/attachment.png>


More information about the ncl-talk mailing list