[ncl-talk] SPI/SPI DROUGHT CLASSIFICATION

Dennis Shea shea at ucar.edu
Mon Apr 15 20:34:04 MDT 2019


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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190415/8120c401/attachment.html>


More information about the ncl-talk mailing list