<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif">Grace,</div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">NCL tells you that you have tried to access "times" beyond its boundary.</div><div class="gmail_default" style="font-family:verdana,sans-serif">Add line:</div><div class="gmail_default" style="font-family:verdana,sans-serif">printVarSummary(times)</div><div class="gmail_default" style="font-family:verdana,sans-serif">before line 49 to make sure that "times" has dimension 137 or more.</div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">Wei</div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><font face="tahoma, sans-serif">================================================</font><div><font face="tahoma, sans-serif">1850 Table Mesa Dr.</font></div><div><font face="tahoma, sans-serif">Boulder, CO 80307</font></div><div><font face="tahoma, sans-serif">Phone: 303-497-8924</font></div></div></div></div>
<br><div class="gmail_quote">On Thu, Jan 8, 2015 at 1:01 AM, grace <span dir="ltr"><<a href="mailto:313695096@qq.com" target="_blank">313695096@qq.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div>Hello:</div><div> I want to add  lightning data on a plot,which is drawn with wrfout data,but it report errors as fellow:</div><div>fatal:Subscript out of range, error in subscript #0<br>fatal:An error occurred
reading times<br>fatal:["Execute.c":8567]:Execute: Error occurred at or near
line 49 in file /home/Huanglei/wrf_Cloudgraup_add_L.ncl<span style="line-height:0px">â€</span><br></div><div>Here is my script,could you guys give me some tips?Any information will be appreciated.</div><div>;  Example script to produce plots for a WRF real-data run,</div><div>;  with the ARW coordinate dynamics option.</div><div><br></div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" </div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" </div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"</div><div><br></div><div>begin</div><div>;</div><div>; The WRF ARW input file.  </div><div>; This needs to have a ".nc" appended, so just do it.</div><div> a = addfile("/home/Huanglei/data/wrfoutd032013080521"+".nc","r")       </div><div><br></div><div><br></div><div>; We generate plots, but what kind do we prefer?</div><div> ; type = "eps"</div><div> type = "pdf"</div><div>; type = "ps"</div><div>; type = "ncgm"</div><div> wks = gsn_open_wks(type,"20130806plt_Cloudgraup1020")</div><div><br></div><div> gsn_define_colormap(wks,"precip_11lev")</div><div>; Set some basic resources</div><div> res = True</div><div> res@MainTitle = "REAL-TIME WRF"</div><div> ; res@gsnDraw    =  False          </div><div> ;res@gsnFrame   =  False</div><div><br></div><div> mpres  = True  ; Map resources</div><div> mpres@mpOutlineOn = False  ; Turn off map outlines</div><div> mpres@mpFillOn   = False  ; Turn off map fill</div><div> mpres@mpGridAndLimbOn = True</div><div> pltres = True ; Plot resources</div><div> pltres@PanelPlot  = True  ; Tells wrf_map_overlays not to remove overlays</div><div><br></div><div><br></div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div><br></div><div>; What times and how many time steps are in the data set?</div><div> times = wrf_user_getvar(a,"times",-1)  ; get all times in the file</div><div> ntimes = dimsizes(times)     ; number of times in the file</div><div><br></div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div><br></div><div>;  do it =1, ntimes-1,1     ; TIME LOOP</div><div>  it = 136</div><div>  print("Working on time: " + times(it) )</div><div>  res@TimeLabel = times(it)  ; Set Valid time to use on plots</div><div>  ; print(it + (/0,1,2,3,4,5,6,7,8,9,10,11,12,13,14/))</div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div>; First get the variables we will need     </div><div>  if(isfilevar(a,"QGRAUP"))</div><div><span style="white-space:pre-wrap"> </span>  qi = wrf_user_getvar(a,"QGRAUP",it )</div><div>   qi = qi*1000.</div><div>   qi@units = "g/kg"   </div><div>  end if</div><div><br></div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div><br></div><div>  do level = 13,15,2    ; LOOP OVER LEVELS</div><div><br></div><div>   display_level = level + 2</div><div>   opts = res</div><div>   opts@cnFillOn     = True</div><div>   opts@gsnSpreadColors  = False</div><div>   opts@ContourParameters    = (/ 1, 19, 2 /)</div><div>   opts@PlotLevelID    = "Eta Level  " + display_level</div><div><span style="white-space:pre-wrap"> </span>  opts@gsnDraw    =  False          </div><div>   opts@gsnFrame   =  False</div><div><br></div><div><br></div><div>   if (isvar("qi"))</div><div>    qis  = qi(level + (/0,1,2,3,4,5,6/),:,:)</div><div>    qisum = dim_sum_n_Wrap(qis, 0)</div><div>    contour = wrf_contour(a,wks,qisum,opts)</div><div>    plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres)</div><div><br></div><div>    delete(contour)</div><div>   end if</div><div>;>============================================================<</div><div>;            add  map</div><div>;>------------------------------------------------------------<</div><div><br></div><div>   </div><div> shp_name1   = "/home/Huanglei/map/China/diquJie_polyline.shp"</div><div><br></div><div> lnres          = True</div><div> lnres@gsLineColor    = "gray25"</div><div> lnres@gsLineThicknessF = 0.5   </div><div><br></div><div> id = gsn_add_shapefile_polylines(wks,plot,shp_name1,lnres)</div><div> shp_name2   = "/home/Huanglei/map/China/cnmap/cnhimap.shp"</div><div><br></div><div> prres=True</div><div> prres@gsLineThicknessF = 2.0     </div><div> prres@gsLineColor = "black"</div><div> plotcn3 = gsn_add_shapefile_polylines(wks,plot,shp_name2,prres)</div><div><br></div><div>  txres2  = True</div><div>  txres2@txFont  = 10</div><div>  txres2@txFontHeightF =0.01</div><div>  txres2@txFontColor = "Blue"</div><div>  txdum1 =gsn_add_text(wks, plot, "Chengdu", 104.06,30.67, txres2)</div><div><br></div><div> draw(plot)    ; This will draw the map and the shapefile outlines.</div><div><br></div><div><br></div><div>  delete(opts)</div><div><div style="line-height:21px">;>============================================================<</div><div style="line-height:21px">;            add  lightning</div><div style="line-height:21px">;>------------------------------------------------------------<<span style="line-height:0px">â€</span></div></div><div>ascii_filename = "/home/Huanglei/data/20130806/utc11.txt"</div><div> seismic = asciiread(ascii_filename,(/183,3/),"float")   </div><div>  </div><div>  y = seismic(:,0)  ; Column 1 of file contains X values.</div><div>  x = seismic(:,1)  ; Column 2 of file contains Y values.</div><div>  z = seismic(:,2)  ; Column 3 of file contains Z values.</div><div> txres2  = True</div><div> txres2@txFont  = 0.01</div><div> txres2@txFontHeightF =0.01</div><div> txres2@txFontColor = "Red"</div><div> idx = ind(z .gt. 0)</div><div> print(idx)</div><div>if .not. all(ismissing(idx))</div><div>  str = new(dimsizes(idx), "string")</div><div>  str = "+"</div><div>  txdum1 = gsn_add_text(wks, plot, str, x(idx),y(idx), txres2)</div><div>end if</div><div><br></div><div>txres2@txFontColor = "Blue"</div><div><br></div><div>idx := ind(z .lt. 0)</div><div>if .not. all(ismissing(idx))</div><div>  str := new(dimsizes(idx), "string")</div><div>  str = "-"   </div><div>  txdum2 = gsn_add_text(wks, plot, str, x(idx),y(idx), txres2)</div><div>end if</div><div>  draw(plot) </div><div> frame(wks)</div><div>  end do    ; END OF LEVEL LOOP</div><div><br></div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div><br></div><div>;  end do     ; END OF TIME LOOP    </div><div><br></div><div>end</div></div><div><span style="line-height:0px">â€</span><br></div><br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>