[ncl-talk] Fatal error

[azizul haque]

Hi NCL User,
I am facing big problem with GPCP daily data. I doing daily annual cycle with other CMIP5 models this work well. However, GPCP data gives following error message that I could not solve.Your kind help are appreciated

Error are highlighted in the script below

*************fatal:Subscript out of range, error in subscript #0
fatal:An error occurred reading pr9Area***********
Script for GPCP sub
j =addfile ("/Users/ahaque/Desktop/CMIP5/historical/Daily/historical/jjas_hist_50y/jjas_hist_50y/GPCP_1DD_v1.2_199610-201312_sa.nc", "r")
   pr9=j->PREC;(6301x16X21;x21;{lat1:lat2},{lon1:lon2})   ; (time,lat,lon)  (0,1,2)
 time9     = i->time
;Convert kg to mm
  pr9        = pr9*1.
  pr9 at units  = "mm/day"
  printVarSummary(pr9)
  printMinMax(pr9,1)
  clat9      = latRegWgt(pr9&lat, typeof(pr9), 0) ; create lat weights
  pr9Area    = wgt_areaave_Wrap(pr9, clat9, 1, 0)  ; weighted area avg
  ymdStrt9   = 199610 ; start yyyymmdd
  ymdLast9   = 201312      ;last
  yrStrt     = ymdStrt9/10000
  yrLast     = ymdLast9/10000
  nhar       = 4  ; no of fourier comp
  var        = "pr"
  TIME9      = cd_calendar(time9,0)   ; type of float
  year9      = floattointeger( TIME9(:,0) )
  month9     = floattointeger( TIME9(:,1) )
  day9       = floattointeger( TIME9(:,2) )
  ddd9       = day_of_year(year9, month9, day9)
  yyyyddd9   = year9*1000 + ddd9
  nday9      = 365
  pday9      = new ( nday9, typeof(pr9), getFillValue(pr9))
  do nd9      =0,nday9-2
     iday9   = ind(ddd9.eq.(nd9+1))
     pday9(nd9) = avg(pr9Area(iday9))
     delete(iday9)                       ; may change next ireration
    end do
   pday9(nday9-1)= 0.5*(pday9(0)+pday9(nday9-2))
   pr9AvgRegion  = avg(pr9Area)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140929/a815b3a1/attachment.html