map of Standardized Precipitation Index (SPI)

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


Received on Fri Jun 13 03:53:47 2014