Re: map of Standardized Precipitation Index (SPI)

Hi Maria,

It fixes it, but I got another error:

atal:NclOneDValGetClosestIndex: finish coordinate index out of range, can't=
 continue
fatal:Could not obtain coordinate indexes, unable to perform subscript
fatal:["Execute.c":8128]:Execute: Error occurred at or near line 71 in file=
 rainfall_SPI.ncl


===

res@gsnCenterString = yyyymm({200007})
    plot(0) = gsn_csm_contour_map_ce(wks,spi({200007},:,:), res)

==

How to get my coordinate indexes?


________________________________
De : Maria Gehne - NOAA Affiliate <maria.gehne@noaa.gov>
Envoyé : vendredi 13 juin 2014 09:32
À : mamadoulamine.mbaye@ucad.edu.sn
Cc : ncl-talk@ucar.edu USERS
Objet : Re: [ncl-talk] map of Standardized Precipitation Index (SPI)

Hi Mamadou,

if you look at the documentation for yyyymm_time (http://www.ncl.ucar.edu/D=
ocument/Functions/Contributed/yyyymm_time.shtml)
you can see that the first and second argument have to be integers. Your yr=
Strt and yrLast are not integers. You can try

year = yyyymm/100
yrStrt = toint(year(0))
yrLast = toint(year(ntim-1))

and see if that fixes your problem.

Maria



On Fri, Jun 13, 2014 at 9:53 AM, mamadoulamine.mbaye@ucad.edu.sn<mailto:mam=
adoulamine.mbaye@ucad.edu.sn> <mamadoulamine.mbaye@ucad.edu.sn<mailto:mamad=
oulamine.mbaye@ucad.edu.sn>> wrote:

Hi,

I am trying to compute and plot the map of Standardized Precipitation In=
dex (SPI), but I got the following error:

---------------------------------------------------------------------------=
------

Variable: prc
Type: double
Total Size: 2918400 bytes
            364800 values
Number of Dimensions: 3
Dimensions and sizes: [time | 1140] x [lat | 20] x [lon | 16]
Coordinates:
            time: [20060131.99930555..21001231.99930555]
            lat: [18.75..9.25]
            lon: [-13.75..-6.25]
Number Of Attributes: 2
  _FillValue : -8.999999873090293e+33
  missing_value : -8.999999873090293e+33
(0) min=0 max=33.5802
fatal:Argument type mismatch on argument (0) of (yyyymm_time) can not coerc=
e
fatal:["Execute.c":8128]:Execute: Error occurred at or near line 33 in file=
 rainfall_SPI.ncl

==========================
==========================
===================

here is the script :

  prc = flt2dble(f1->var1)
  pmsg = prc@_FillValue ; convenience

  printVarSummary(prc)
  printMinMax(prc,0)

  runlen = (/12,24/)
  nrun = dimsizes(runlen)
;*********************************
; plot parameters
;*********************************
  dimprc = dimsizes(prc)
  ntim = dimprc(0)
  nlat = dimprc(1)
  mlon = dimprc(2)

  yyyymm = f1->time
  yyyymm = yyyymm/100

  year = yyyymm/100
  yrStrt = year(0)
  yrLast = year(ntim-1)
  nyear = yrLast-yrStrt+1
  yyyymm = yyyymm_time(yrStrt, yrLast, "integer")
  yrfrac = (/ yyyymm_to_yyyyfrac(yyyymm, 0.0) /)

  prc&time = (/ yyyymm /)

;=

  wks = gsn_open_wks ("x11","spi")
  gsn_define_colormap(wks,"BlueWhiteOrangeRed")
  gsn_reverse_colormap(wks) ; Reverse the color map.

  res = True
  res@gsnDraw = False ; don't draw
  res@gsnFrame = False ; don't advance frame
  res@gsnSpreadColors = True

  res@cnFillOn = True ; color Fill
  res@cnFillMode = "RasterFill" ; Raster Mode
  res@cnLinesOn = False ; Turn off contour lines
  res@cnLevelSelectionMode = "ManualLevels" ; set manual contour levels
  res@cnMinLevelValF = -3.0 ; set min contour level
  res@cnMaxLevelValF = 3.0 ; set max contour level
  res@cnLevelSpacingF = 0.5 ; set contour spacing
  res@lbLabelBarOn = False ; turn off individual cb's

  resP = True ; panel resources
  resP@txString = "SPI: 2006-2100 (BC_REMO_RCP45)" ; add center s=
tring
  resP@gsnPanelLabelBar= True ; add common colorbar

  plot = new ( 2, "graphic")

  do nr=0,nrun-1
    spi = dim_spi_n(prc, runlen(nr), False, 0)

    spi@long_name = "SPI"
    spi@units = "run="+runlen(nr)
    copy_VarCoords(prc, spi)

    res@gsnCenterString = yyyymm({200007})
    plot(0) = gsn_csm_contour_map_ce(wks,spi({200007},:,:), res)

    res@gsnCenterString = yyyymm({201012})
    plot(1) = gsn_csm_contour_map_ce(wks,spi({200612},:,:), res)

    gsn_panel(wks,plot,(/2,1/),resP) ; now draw as one plot
  end do
end
=======

I was trying to adapt the script from NCL scripts.

Thanks in advance

Mamadou

Regards


_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk


Received on Fri Jun 13 06:31:14 2014