Number of subscripts do not match number of dimensions of variable

Dear NCL users,

I got an error during the panel plot: fatal:Number of subscripts do not
match number of dimensions of variable,(3) Subscripts used, (4) Subscripts
expected

fatal:Execute: Error occurred at or near line 16 in file

Any help will be appreciated about fixing the problem

Thanks in advance

wei

*********************************

= = = = = = =

= = = = = = =

= = = = = =

  times = wrf_user_list_times (a); get times in the file

  ntimes = dimsizes(times) ; number of times in the file

  FirstTime = res

  mdims = getfilevardimsizes(a,"P") ; get some dimension sizes for the file

  nd = dimsizes(mdims)

  plots = new (4 , graphic )

; Get the data

  z = wrf_user_getvar(a, "z",-1) ; grid point height

  p = wrf_user_getvar(a,"pressure",-1)

; Interpolate to pressure

  if (FirstTime) then

  zz = z/100.

  dims = dimsizes(zz)

  zmax = -999.

  do imax = 1,dims(0)-1

  *if ( .not.ismissing(zz(imax,1,1)) .and. zz(imax,1,1) .gt. zmax ) then*

  zmax = zz(imax,1,1)

  end if

  end do

  imax = floattoint(zmax)

  zmax = imax/10.

  print("Y labels set to: 0 - " + zmax + " km")

  FirstTime = False

end if

; Interpolate data vertically

  mpv_plane = wrf_user_intrp3d(mpv,p,"v",plane,0.,opts)

  eth_plane = wrf_user_intrp3d(eth,p,"v",plane,0.,opts)

 ;

  gmpv_plane = wrf_user_intrp3d(gmpv,p,"v",plane,0.,opts)

  ths_plane = wrf_user_intrp3d(ths,p,"v",plane,0.,opts)

  ; Options for XY Plots

  opts_xy = res

  opts_xy@tiYAxisString = "Pressure (hPa)"

  opts_xy@AspectRatio = 0.75

  opts_xy@cnMissingValPerimOn = True

  opts_xy@cnMissingValFillColor = 0

  opts_xy@cnMissingValFillPattern = 11

  opts_xy@tmYLMode = "Explicit"

  opts_xy@tmYLValues = fspan(0,100,10) ; over the 0-100 span of
the Y label

  opts_xy@tmYLLabels = fspan(1000,100,10); Corresponding 11
labels

  opts_xy@PlotOrientation = ths_plane@Orientation

; Plotting options for gmpv

  opts_gmpv = opts_xy

  opts_gmpv@ContourParameters = (/ - 4., 5., 1. /)

  opts_gmpv@cnFillOn = True

  opts_gmpv@lbLabelBarOn = False

  opts_gmpv@gsnSpreadColors = True

; Plotting options for Temperature

  opts_eth = opts_xy

  opts_eth@ContourParameters = (/ 270., 450., 4. /)

  opts_eth@cnInfoLabelOn = False ; Turn off individual label bars so we
can

; Plotting options for Equivalent Potential Temperature

  opts_ths = opts_xy

  opts_ths@ContourParameters = (/ 220., 330., 4. /)

; Get the contour info for the mpv and theta_e

  opts_xy@PlotOrientation = eth_plane@Orientation

  contour_eth = wrf_contour(a,wks,eth_plane(7,:,:),opts_eth)

  contour_mpv = wrf_contour(a,wks,mpv_plane(7,:,:),opts_gmpv)

  contour_eth1 = wrf_contour(a,wks,eth_plane(8,:,:),opts_eth)

  contour_mpv1 = wrf_contour(a,wks,mpv_plane(8,:,:),opts_gmpv)

; Get the contour info for the gmpv and ths

  opts_xy@PlotOrientation = ths_plane@Orientation

  contour_ths = wrf_contour(a,wks,ths_plane(7,:,:),opts_ths)

  contour_gmpv = wrf_contour(a,wks,gmpv_plane(7,:,:),opts_gmpv)

  contour_ths1 = wrf_contour(a,wks,ths_plane(8,:,:),opts_ths)

  contour_gmpv1 = wrf_contour(a,wks,gmpv_plane(8,:,:),opts_gmpv)

; MAKE PLOTS

  plots(0) = wrf_overlays(a,wks,(/contour_mpv,contour_eth/),pltres)

  plots(1) = wrf_overlays(a,wks,(/contour_gmpv,contour_ths/),pltres)

  plots(2) = wrf_overlays(a,wks,(/contour_mpv1,contour_eth1/),pltres)

  plots(3) = wrf_overlays(a,wks,(/contour_gmpv1,contour_ths1/),pltres)

; Delete options and fields, so we don't have carry over

  delete(opts_eth)

  delete(opts_ths)

  delete(opts_gmpv)

  delete(eth_plane)

  delete(mpv_plane)

  delete(ths_plane)

  delete(gmpv_plane)

 ; create panel

;************************************************

 ;Panel the WRF plots

  resP = True ; modify the panel plot

  resP@gsnFrame = False ; don't advance panel plot

  resP@gsnPanelLabelBar = True ; add common colorbar

  resP@gsnPanelScalePlotIndex = 1

  resP@gsnPanelBottom = 0.2 ; shrink panel plot by setting bottom
edge of plot

  resP@gsnPanelTop = 0.9

  gsn_panel(wks,(/plots/),(/2,2/),resP) ; now draw as one plot

  frame(wks)

  end

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Feb 4 09:37:57 2011