关于画风场图的问题求助~画出来的图很怪。。。

海棠aroma

数据是从NCEP上下的1948-2014年的月均u和v的.nc文件

程序如下：
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/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
;************************************************
begin
;************************************************
; read in netCDF file
;************************************************
   ufname = "uwnd_mon_mean.nc"  ;get uwind file name
   vfname = "vwnd_mon_mean.nc"  ;get vwind file name

  udata  = addfile(ufname,"r")  ; Open uwind file
  vdata  = addfile(vfname,"r")  ; Open vwind file

  utmp = udata->uwnd           ; get uwnd data
  vtmp = vdata->vwnd            ; get uwnd data
  level = udata->level
  time  = udata->time
  ;time1 =  ut_calendar(time,-1)
  ttime = dimsizes(time)
   lat = udata->lat
   lon = udata->lon

   nlat=dimsizes(lat)
   nlon=dimsizes(lon)

avr_monthu = new((/nlat,nlon/),float)
  avr_monthu!0 = "lat"
  avr_monthu&lat = utmp&lat
  avr_monthu&lat@units = "degrees_north"
  avr_monthu!1= "lon"
  avr_monthu&lon = utmp&lon
  avr_monthu&lon@units = "degrees_east"
  avr_monthv = avr_monthu
   nmod=12
;*************************************
do nlevel = 0,16
   do ntime =0,nmod-1


    u = utmp(:,nlevel,:,:)
    v = vtmp(:,nlevel,:,:)

  avr_monthu = dim_avg_n(u(ntime:ttime-1:nmod,:,:),0)
  avr_monthv = dim_avg_n(v(ntime:ttime-1:nmod,:,:),0)
   printVarSummary(avr_monthv)




;************************************************
; create plot
;************************************************
  filename="uv_mon_1980-2013_ltm_"+level(nlevel)+"_"+ntime

  wks = gsn_open_wks("eps",filename)                ; open a eps file

  res                         = True             ; plot mods desired
  res@gsnFrame                = False            ; so we can draw time stamp
  res@vcRefAnnoOrthogonalPosF = -1.0             ; move ref vector up
  res@vcRefMagnitudeF         = 10.0             ; define vector ref mag
  res@vcRefLengthF            = 0.045            ; define length of vec ref
  res@vcGlyphStyle            = "CurlyVector"    ; turn on curly vectors
  res@vcMinDistanceF          = 0.017

   res@vpWidthF   = 0.8                          ; make map bigger
   res@vpHeightF  = 0.8
   res@mpMaxLatF  = 90                         ; select subregion
   res@mpMinLatF  = -90
   res@mpMinLonF  = 0
   res@mpMaxLonF  = 357.5

   res@gsnLeftString  = "uv at "+level(nlevel)+"hPa "+ntime


  vector = gsn_csm_vector_map_ce(wks, avr_monthu, avr_monthv,res)


;************************************************
; add time stamp to plot
;************************************************

  frame (wks)  
       cmd="convert -geometry 1000x1000 -density 1000  -trim "+filename+".eps "+filename+".png"


                system(cmd)
    end do
end do
end

得出的图是下面这样子。。。







刚开始提示错误是(0) check_for_y_lat_coord: Warning: Data either does not contain a valid latitude coordinate array or doesnt contain one at all (0) check_for_lon_coord: Warning: Data either does not contain a valid longitude coordinate array or doesnt contain one at all
后来我在前面加了数组的属性，他不报这个错了，但是画出来的图还是错的，这是为啥？？
我感觉是因为我求平均求错了？？

飏逸雅

我也遇到过这种问题。在读入变量后，貌似要用short2flt把变量转换一下，有的nc数据好像是要这么做吧。。。

lleoiu

你读出来的U,V是短整形的，在计算的时候要转成浮点数组来计算。
u     = short2flt(u)

longlivehj

飏逸雅 发表于 2014-11-17 09:38
我也遇到过这种问题。在读入变量后，貌似要用short2flt把变量转换一下，有的nc数据好像是要这么做吧。。。

学习了

海棠无香aroma

lleoiu 发表于 2014-11-17 11:52
你读出来的U,V是短整形的，在计算的时候要转成浮点数组来计算。
u     = short2flt(u)

我short2flt之后 就成这样了。。。。。

飞云侠客

很明显嘛，你的两个分量是相等的，用了两个U或两个V，而不是UV，你的数据有问题

longlivehj

U、V分量相等了，楼主仔细检查！

海棠无香aroma

飞云侠客 发表于 2014-11-17 19:08
很明显嘛，你的两个分量是相等的，用了两个U或两个V，而不是UV，你的数据有问题

T.T谢谢提醒。。。我看了两天就没看到这里我写错了。。。T.T谢谢谢谢

海棠无香aroma

longlivehj 发表于 2014-11-17 19:09
U、V分量相等了，楼主仔细检查！

十分感谢- -！！！我太大意了！！！！谢谢！

tSSkill

不知道为什么，现在下这个资料只有三维的了，没有高度维度了，很奇怪