<div dir="ltr"><br><div>dear ncl talk</div><div><br></div><div>I need help if there is anyone can help me.</div><div><br></div><div>I was using grid data of my country with NetCDF format. but I cannot understand what the error is if anyone can help me.this is my code and the error it gives me is.</div><div><br></div><div><div><div>begin</div><div>;--------------------------------------------------------------------------</div><div>f =   "Sweetheart.nc"</div><div>;f = "<a href="http://mon.precip.nc">mon.precip.nc</a>"</div><div>fin  = addfile(f, "r")</div><div>timet = fin->time</div><div>ymd = round(timet, 3)</div><div><br></div><div>;ymd = cd_calendar(timet,2)</div><div>iTim =  ind (ymd.ge.19830101.and.ymd.lt.20100101)</div><div>precpi = fin->precip(iTim,{4:8},{34:39})</div><div>;==============================================================================</div><div>timem =precpi&time</div><div>rain = rmMonAnnCycTLL(precpi)</div><div>prec_jjas1 = rain(2::12,:,:)</div><div>rf1= runave_Wrap(rain(lat|:,lon|:,time|:),4,0)</div><div>arr_jjas1= (/rf1(time|2::12,lat|:,lon|:)/)</div><div><br></div><div>arr_jjas1!0 = "time"</div><div>arr_jjas1!1  = "lat"</div><div>arr_jjas1!2 = "lon"</div><div>arr_jjas1&time = prec_jjas1&time</div><div>arr_jjas1&lat = prec_jjas1&lat</div><div>arr_jjas1&lon = prec_jjas1&lon</div><div>printVarSummary(arr_jjas1)</div><div><br></div><div>yene = dim_avg_Wrap (arr_jjas1(time|:,lat|:,lon|:))</div><div> yene2 = dim_avg_Wrap(yene)</div><div> stand = dim_standardize_Wrap(yene2,1)</div><div>;printVarSummary(stand)</div><div>;*************************************************************************</div><div> fil = "Sweetheart.nc"</div><div> ;fil = "<a href="http://mon.precip.nc">mon.precip.nc</a>"</div><div>  fin1 =addfile(fil,"r")</div><div>  timep = fin1->time</div><div>ymd1 = round(timet, 3)</div><div><br></div><div>  ;ymd1 = cd_calendar(timep,2)</div><div>  iTim = ind(ymd1.ge.19830101.and.ymd1.lt.20100101)</div><div>  precip = fin1->precip(iTim,:,:)</div><div>;***************************************************************************</div><div>mare =rmMonAnnCycTLL(precip)</div><div>time = mare&time</div><div>prec_jjas = mare(3::12,:,:)</div><div>rf= runave_Wrap(mare(lat|:,lon|:,time|:),3,0)</div><div>arr_jjas= (/rf(time|3::12,lat|:,lon|:)/)</div><div><br></div><div>arr_jjas!0 = "time"</div><div>arr_jjas!1  = "lat"</div><div>arr_jjas!2 = "lon"</div><div>arr_jjas&time = prec_jjas&time</div><div>arr_jjas&lat = prec_jjas&lat</div><div>arr_jjas&lon = prec_jjas&lon</div><div>printVarSummary(arr_jjas)</div><div>;===========calculate corelation==========================================</div><div>ccr = escorc(stand,arr_jjas(lat|:,lon|:,time|:))</div><div>copy_VarCoords(prec_jjas(0,:,:),ccr)</div><div>printVarSummary(ccr)</div><div>printMinMax(ccr,0)</div><div>;========================================================================</div><div>siglvl =0.05</div><div>test = rtest(ccr(lat|:,lon|:),34,0)</div><div>test = where(ismissing(ccr),test@_FillValue,test)</div><div><br></div><div>test!0 =  "lat"</div><div>test!1 ="lon"</div><div>test&lat = ccr&lat</div><div>test&lon =ccr&lon</div><div>printVarSummary(test)</div><div>;==============================================================================</div><div>wks = gsn_open_wks("x11","Correlation")</div><div>gsn_define_colormap(wks,"precip_diff_12lev")</div><div>res = True</div><div>res@cnFillOn     = True</div><div>;res@cnLineOn =   True</div><div>res@cnLevelSelectionMode = "ManualLevels"</div><div>res@mpCenterLonF =180</div><div>res@cnMinLevelValF       = -0.6  </div></div></div><div><div>res@cnMaxLevelValF       =  0.6</div><div>res@cnLineLabelsOn = False</div><div><br></div><div>res@cnLevelSpacingF      = 0.1</div><div>;res@gsnSpreedColors = True</div><div>res@mpPerimOn =False</div><div>res@tiMainString         = " (a) Corr (SERI ,MAM  precipitation and 850-hpa winds) "</div><div>res@mpFillOn  = False</div><div>res@gsnDraw =False</div><div>res@cnLinesOn            = False</div><div>res1 =True</div><div>res1@gsnDraw =False</div><div>res1@gsnFrame     =  False</div><div><br></div><div>res1@cnLevelSelectionMode = "ManualLevels"</div><div>res1@cnMinLevelValF       = 0.00       ; min level</div><div>res1@cnMaxLevelValF       =  1.05</div><div>res1@cnLevelSpacingF      = 0.05</div><div>res1@cnInfoLabelOn = False</div><div>res1@cnLinesOn = False</div><div>res1@cnLineLabelsOn = False</div><div> ;res@mpMinLatF            = 3     ; range to zoom in on</div><div><br></div><div>       ;res@mpMaxLatF            = 15</div><div>       ;res@mpMinLonF             = 33</div><div>       ;res@mpMaxLonF            = 48</div><div>       ;res@mpOutlineBoundarySets = "National"</div><div>       ;res@mpAreaMaskingOn =True</div><div><br></div><div>res1@cnLineLabelOn  = False</div><div>res1@cnFillScaleF = 0.6 ;add extra density</div><div>res1@gsnDraw      =  False                   ; do not draw the plot</div><div>res1@gsnFrame     =  False                   ; do not advance the frame</div><div>;res@tiMainString =   "Correlation between NIO  and rainfall 9month led "</div><div>plot0    = gsn_csm_contour_map(wks,ccr,res)</div><div><br></div><div>opt =True</div><div>opt@gsnShadeFillType = "pattern"</div><div>opt@gsnShadeLow = 17</div><div>opt@gsnShadeDotSizeF = 0.6</div><div>plot1    = gsn_csm_contour(wks,test,res1)</div><div>plot1    = gsn_contour_shade(plot1,0.05,0.75,opt)</div><div>;=========================================================================================</div><div>ECMWF   = "/home/cyoo/scratch/reanl/ecmwf/monthly/"</div><div>fileU   = systemfunc ("ls "+ECMWF+"uwnd/u.mon.*.nc")</div><div>F    =  addfiles (fileU, "r")</div><div>timeu = F[:]->time</div><div>ymds = cd_calendar(timeu,2)</div><div>iTim1 =  ind (ymds.ge.19790101.and.ymds.lt.20120101)</div><div>uwind  = F[:]->u (iTim1,{850},:,:)</div><div>;======================================================================</div><div>Uanom =rmMonAnnCycTLL(uwind)</div><div>;T = Uanom&time</div><div> T = uwind&time</div><div><br></div><div>;ym1 =cd_calendar(T,0)</div><div> ;month1 = tointeger(ym1(:,1))</div><div> ;imon1 = ind(month1.ge.02.and.month1.le.05)</div><div><br></div><div>U1 = Uanom(2::12,:,:)</div><div>ub1= runave_Wrap(Uanom(lat|:,lon|:,time|:),4,0)</div><div>uu21= (/ub1(time|2::12,lat|:,lon|:)/)</div><div>uu21!0 = "time"</div><div>uu21!1  = "lat"</div><div>uu21!2 = "lon"</div><div>uu21&time = U1&time</div><div>uu21&lat = U1&lat</div><div>uu21&lon = U1&lon</div><div>;printVarSummary(uu21)</div><div><br></div><div>U = dim_avg_Wrap(uu21(lat|:,lon|:,time|:))</div><div><br></div><div>;=====================================================================</div><div>fileV   = systemfunc ("ls "+ECMWF+"vwnd/v.mon.*.nc")</div><div>fv      =  addfiles (fileV, "r")</div><div>time2 = fv[:]->time</div><div>ymd2 = cd_calendar(time2,2)</div><div>iTim2 =  ind (ymd2.ge.19790101.and.ymd2.lt.20120101)</div><div>vwind =  fv[:]->v (iTim2,{850},:,:)</div></div><div><div>Vanom =rmMonAnnCycTLL(vwind)</div><div>;T1 = Vanom&time</div><div>;T1 = vwind&time</div><div>V1 = Vanom(2::12,:,:)</div><div>vb1= runave_Wrap(Vanom(lat|:,lon|:,time|:),4,0)</div><div>vu21= (/vb1(time|2::12,lat|:,lon|:)/)</div><div>vu21!0 = "time"</div><div>vu21!1  = "lat"</div><div>vu21!2 = "lon"</div><div>vu21&time = V1&time</div><div>vu21&lat = V1&lat</div><div>vu21&lon = V1&lon</div><div>V = dim_avg_Wrap(vu21(lat|:,lon|:,time|:))</div><div><br></div><div><br></div><div>  res2           = True</div><div>  res2@vcMinFracLengthF     = 0.3</div><div>  res2@vcRefMagnitudeF       = 20               ; define vector ref mag</div><div>  res2@vcRefLengthF          =   0.04      ; define length of vec ref</div><div>res2@vcGlyphStyle               = "LineArrow"</div><div>res2@vcMinDistanceF        = 0.02           ; thin out vectors</div><div>  ;res2@vcGlyphStyle          = "CurlyVector"        ; turn on curley vectors</div><div>  res2@gsnDraw               = False                ; don't draw yet</div><div>  res2@gsnFrame              = False                ; don't advance frame yet</div><div>  res2@vcLineArrowThicknessF      = 1.5</div><div>  res2@gsnLeftString         = ""                   ; turn off strings</div><div>  res2@gsnRightString        = ""</div><div><br></div><div>plotV = gsn_csm_vector(wks,U,V,res2)</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div>overlay(plot0,plot1)                       ; overlay the U-wind plot on the rainfall plot</div><div>overlay(plot0,plotV)</div><div>draw(plot0)                                  ; draw the rainfall plot (with the U-wind plot overlaid)</div><div>frame(wks)</div><div>end</div><div><br></div><div><div>fatal:["NclFile.c":2099]:Subscript out of range, error in subscript #0</div><div>fatal:["Execute.c":8575]:Execute: Error occurred at or near line 22 in file stan</div><div>dardizerain.ncl</div></div><div><br></div><div><br></div></div></div>