<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:13px"><div class="" style="">Hi NCL User,</div><div class="" style="">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</div><div class="" style=""><br class="" style=""></div><div class="" style="color: rgb(0, 0, 0); font-size: 13px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal;">Error are highlighted in the script below</div><div class="" style="color: rgb(0, 0, 0); font-size: 13px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal;"><br></div><div class="" style="color: rgb(0, 0, 0); background-color: transparent; font-style: normal;"><span
class="">*************</span><span style="font-family: Menlo; font-size: 11px;" class="">fatal:Subscript out of range, error in subscript #0</span></div>
<div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class="">fatal:An error occurred reading pr9Area***********</div><div class="" style="color: rgb(0, 0, 0); font-size: 13px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal;">Script for GPCP sub</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class="">j =addfile ("/Users/ahaque/Desktop/CMIP5/historical/Daily/historical/jjas_hist_50y/jjas_hist_50y/GPCP_1DD_v1.2_199610-201312_sa.nc", "r")</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> pr9=j->PREC;(6301x16X21;x21;{lat1:lat2},{lon1:lon2}) ; (time,lat,lon) (0,1,2)</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> time9 = i->time</div><div style="margin-top: 0px;
margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class="">;Convert kg to mm</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> pr9 = pr9*1.</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> pr9@units = "mm/day"</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> printVarSummary(pr9)</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> printMinMax(pr9,1)</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> clat9 = latRegWgt(pr9&lat, typeof(pr9), 0) ; create lat weights</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> pr9Area =
wgt_areaave_Wrap(pr9, clat9, 1, 0) ; weighted area avg</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> ymdStrt9 = 199610 ; start yyyymmdd</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> ymdLast9 = 201312 ;last</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> yrStrt = ymdStrt9/10000</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> yrLast = ymdLast9/10000</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> nhar = 4 ; no of fourier comp</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> var =
"pr"</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> TIME9 = cd_calendar(time9,0) ; type of float</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> year9 = floattointeger( TIME9(:,0) )</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> month9 = floattointeger( TIME9(:,1) )</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> day9 = floattointeger( TIME9(:,2) )</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> ddd9 = day_of_year(year9, month9, day9)</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> yyyyddd9
= year9*1000 + ddd9</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> nday9 = 365</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> pday9 = new ( nday9, typeof(pr9), getFillValue(pr9))</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> do nd9 =0,nday9-2</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> iday9 = ind(ddd9.eq.(nd9+1))</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> <span style="background-color: rgb(253, 239, 43);">pday9(nd9) = avg(pr9Area(iday9))</span></div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;"
class=""> delete(iday9) ; may change next ireration</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> end do</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> pday9(nday9-1)= 0.5*(pday9(0)+pday9(nday9-2))</div><div style="background-color: transparent;" class="">
</div><div style="margin-top: 0px; margin-bottom: 0px; font-size: 11px; font-family: Menlo;" class=""> pr9AvgRegion = avg(pr9Area)</div><div class="" style=""><br class="" style=""></div></div></body></html>