<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 id="_editor_bookmark_start_1" style="display: none; 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 class="Apple-tab-span" style="white-space:pre"> </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 class="Apple-tab-span" style="white-space:pre"> </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 id="_editor_bookmark_start_2" style="display: none; 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 id="_editor_bookmark_start_0" style="display: none; line-height: 0px;"></span><br></div>