plot dim_avg

Dear all
I try to average data only for daytime and plot it out.
After I used dim_avg, then dimesion changed from 3 to 2 dimension and when
I used
gsm_contour_map to plot it out.

Please help
Mark

    filesn1 = systemfunc("ls " + dirs + "2006*LDASOUT*.nc")
       f = addfiles(filesn, "r") ; change to filesn1 if 2006 )
       x = f[:]->HFX(:,:,:) ; (Time,
south_north, west_east)
       x1 = x(5:10:19,:,:)
    x1@units = " (a) (Default)"
; print(x)
        dirk = "/usr/rmt_share/tgdata/hrldas/workspace/jam2/india_gemout/"
        filesk = systemfunc("ls " + dirk + "200511*LDASOUT*.nc")
        filesk1 = systemfunc("ls " + dirs + "2006*LDASOUT*.nc")

       g = addfiles(filesk, "r")
       y = g[:]->HFX(:,:,:)
       y1 = y(5:10:19,:,:)

 y1@units = " (b) (GEM)"
        xavg = dim_avg_Wrap(x1)
        yavg = dim_avg_Wrap(y1)
; print(x1)
print(yavg)
      x@lat2d = lat2d
      x@lon2d = lon2d
      y@lat2d = lat2d
      y@lon2d = lon2d

  res@mpFillOn = False ; turn off map fill
  res@mpOutlineDrawOrder = "PostDraw" ; draw continental outline
last
  res@mpOutlineBoundarySets = "National" ; state boundaries

  res@tfDoNDCOverlay = False ; True for 'native' grid
  res@gsnAddCyclic = False ; data are not cyclic
; res@lbOrientation = "vertical"

;************************************************
; allocate array for 3 plots
;************************************************
  plts = new (2,"graphic")

;************************************************
; Tell NCL not to draw or advance frame for individual plots
;************************************************
  res@gsnDraw = False ; (a) do not draw
  res@gsnFrame = False ; (b) do not advance 'frame'

  plts(0) = gsn_csm_contour_map(wks,xavg(:,:),res)
  plts(1) = gsn_csm_contour_map(wks,yavg(:,:),res)
;************************************************

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jun 11 11:43:34 2012