[ncl-talk] An error occurred reading times

[Wei Huang]

Grace,

NCL tells you that you have tried to access "times" beyond its boundary.
Add line:
printVarSummary(times)
before line 49 to make sure that "times" has dimension 137 or more.

Wei


================================================
1850 Table Mesa Dr.
Boulder, CO 80307
Phone: 303-497-8924

On Thu, Jan 8, 2015 at 1:01 AM, grace <313695096 at qq.com> wrote:

> Hello:
>  I want to add  lightning data on a plot,which is drawn with wrfout
> data,but it report errors as fellow:
> fatal:Subscript out of range, error in subscript #0
> fatal:An error occurred reading times
> fatal:["Execute.c":8567]:Execute: Error occurred at or near line 49 in
> file /home/Huanglei/wrf_Cloudgraup_add_L.ncl‍
> Here is my script,could you guys give me some tips?Any information will be
> appreciated.
> ;   Example script to produce plots for a WRF real-data run,
> ;   with the ARW coordinate dynamics option.
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
>
> begin
> ;
> ; The WRF ARW input file.
> ; This needs to have a ".nc" appended, so just do it.
>   a = addfile("/home/Huanglei/data/wrfoutd032013080521"+".nc","r")
>
>
>
> ; We generate plots, but what kind do we prefer?
>  ; type = "eps"
>  type = "pdf"
> ; type = "ps"
> ; type = "ncgm"
>   wks = gsn_open_wks(type,"20130806plt_Cloudgraup1020")
>
>   gsn_define_colormap(wks,"precip_11lev")
> ; Set some basic resources
>   res = True
>   res at MainTitle = "REAL-TIME WRF"
>  ; res at gsnDraw      =  False
>   ;res at gsnFrame     =  False
>
>   mpres  = True  ; Map resources
>   mpres at mpOutlineOn = False  ; Turn off map outlines
>   mpres at mpFillOn    = False  ; Turn off map fill
>   mpres at mpGridAndLimbOn = True
>   pltres = True ; Plot resources
>   pltres at PanelPlot  = True   ; Tells wrf_map_overlays not to remove
> overlays
>
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>
> ; What times and how many time steps are in the data set?
>   times = wrf_user_getvar(a,"times",-1)  ; get all times in the file
>   ntimes = dimsizes(times)         ; number of times in the file
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>
> ;  do it =1, ntimes-1,1        ; TIME LOOP
>    it = 136
>     print("Working on time: " + times(it) )
>     res at TimeLabel = times(it)   ; Set Valid time to use on plots
>    ; print(it + (/0,1,2,3,4,5,6,7,8,9,10,11,12,13,14/))
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
> ; First get the variables we will need
>     if(isfilevar(a,"QGRAUP"))
>   qi = wrf_user_getvar(a,"QGRAUP",it )
>       qi = qi*1000.
>       qi at units = "g/kg"
>     end if
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>
>     do level = 13,15,2      ; LOOP OVER LEVELS
>
>       display_level = level + 2
>       opts = res
>       opts at cnFillOn         = True
>       opts at gsnSpreadColors  = False
>       opts at ContourParameters       = (/ 1, 19, 2 /)
>       opts at PlotLevelID      = "Eta Level  " + display_level
>   opts at gsnDraw      =  False
>       opts at gsnFrame     =  False
>
>
>       if (isvar("qi"))
>         qis  = qi(level + (/0,1,2,3,4,5,6/),:,:)
>         qisum = dim_sum_n_Wrap(qis, 0)
>         contour = wrf_contour(a,wks,qisum,opts)
>         plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres)
>
>         delete(contour)
>       end if
> ;>============================================================<
> ;                      add  map
> ;>------------------------------------------------------------<
>
>
>   shp_name1    = "/home/Huanglei/map/China/diquJie_polyline.shp"
>
>   lnres                  = True
>   lnres at gsLineColor      = "gray25"
>   lnres at gsLineThicknessF = 0.5
>
>  id = gsn_add_shapefile_polylines(wks,plot,shp_name1,lnres)
>   shp_name2    = "/home/Huanglei/map/China/cnmap/cnhimap.shp"
>
>   prres=True
>   prres at gsLineThicknessF = 2.0
>   prres at gsLineColor = "black"
>   plotcn3 = gsn_add_shapefile_polylines(wks,plot,shp_name2,prres)
>
>    txres2  = True
>    txres2 at txFont  = 10
>    txres2 at txFontHeightF =0.01
>    txres2 at txFontColor = "Blue"
>    txdum1 =gsn_add_text(wks, plot, "Chengdu", 104.06,30.67, txres2)
>
>   draw(plot)       ; This will draw the map and the shapefile outlines.
>
>
>    delete(opts)
> ;>============================================================<
> ;                      add  lightning
> ;>------------------------------------------------------------<‍
> ascii_filename = "/home/Huanglei/data/20130806/utc11.txt"
>   seismic = asciiread(ascii_filename,(/183,3/),"float")
>
>     y = seismic(:,0)  ; Column 1 of file contains X values.
>     x = seismic(:,1)  ; Column 2 of file contains Y values.
>     z = seismic(:,2)  ; Column 3 of file contains Z values.
>   txres2  = True
>   txres2 at txFont  = 0.01
>   txres2 at txFontHeightF =0.01
>   txres2 at txFontColor = "Red"
>  idx = ind(z .gt. 0)
>  print(idx)
> if .not. all(ismissing(idx))
>     str = new(dimsizes(idx), "string")
>     str = "+"
>     txdum1 = gsn_add_text(wks, plot, str, x(idx),y(idx), txres2)
> end if
>
> txres2 at txFontColor = "Blue"
>
> idx := ind(z .lt. 0)
> if .not. all(ismissing(idx))
>     str := new(dimsizes(idx), "string")
>     str = "-"
>     txdum2 = gsn_add_text(wks, plot, str, x(idx),y(idx), txres2)
> end if
>    draw(plot)
>   frame(wks)
>     end do      ; END OF LEVEL LOOP
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>
> ;  end do        ; END OF TIME LOOP
>
> end
> ‍
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150108/3b48ea7e/attachment.html