NCL矢量图求指点

infocert

1. netcdf desert.UVEL\,VVEL.01 {
   
2. dimensions:
   
3.     nlat = 384 ;
   
4.     nlon = 320 ;
   
5.     time = UNLIMITED ; // (1 currently)
   
6.     z_t = 60 ;
   
7.     d2 = 2 ;
   
8. variables:
   
9.     double ULAT(nlat, nlon) ;
   
10.         ULAT:long_name = "array of u-grid latitudes" ;
    
11.         ULAT:units = "degrees_north" ;
    
12.         ULAT:_FillValue = 9.96920996838687e+36 ;
    
13.         ULAT:missing_value = 9.96920996838687e+36 ;
    
14.     double ULONG(nlat, nlon) ;
    
15.         ULONG:long_name = "array of u-grid longitudes" ;
    
16.         ULONG:units = "degrees_east" ;
    
17.         ULONG:_FillValue = 9.96920996838687e+36 ;
    
18.         ULONG:missing_value = 9.96920996838687e+36 ;
    
19.     float UVEL(time, z_t, nlat, nlon) ;
    
20.         UVEL:long_name = "Velocity in grid-x direction" ;
    
21.         UVEL:units = "centimeter/s" ;
    
22.         UVEL:coordinates = "ULONG ULAT z_t time" ;
    
23.         UVEL:grid_loc = "3221" ;
    
24.         UVEL:cell_methods = "time: mean" ;
    
25.         UVEL:_FillValue = 9.96921e+36f ;
    
26.         UVEL:missing_value = 9.96921e+36f ;
    
27.     float VVEL(time, z_t, nlat, nlon) ;
    
28.         VVEL:long_name = "Velocity in grid-y direction" ;
    
29.         VVEL:units = "centimeter/s" ;
    
30.         VVEL:coordinates = "ULONG ULAT z_t time" ;
    
31.         VVEL:grid_loc = "3221" ;
    
32.         VVEL:cell_methods = "time: mean" ;
    
33.         VVEL:_FillValue = 9.96921e+36f ;
    
34.         VVEL:missing_value = 9.96921e+36f ;
    
35.     double time(time) ;
    
36.         time:long_name = "time" ;
    
37.         time:units = "days since 0000-01-01 00:00:00" ;
    
38.         time:bounds = "time_bound" ;
    
39.         time:calendar = "noleap" ;
    
40.     double time_bound(time, d2) ;
    
41.         time_bound:long_name = "boundaries for time-averaging interval" ;
    
42.         time_bound:units = "days since 0000-01-01 00:00:00" ;
    
43.     float z_t(z_t) ;
    
44.         z_t:long_name = "depth from surface to midpoint of layer" ;
    
45.         z_t:units = "centimeters" ;
    
46.         z_t:positive = "down" ;
    
47.         z_t:valid_min = 500.f ;
    
48.         z_t:valid_max = 537500.f ;
    

复制代码

各位大神，小弟初学ncl，上面这种格式的文件怎么画（UVEL,VVEL）的矢量图，求各位大神指点，
风场的适量图我会处理，用风场的脚本套这个文件，总是出错，观察了下，问题应该出在nlat,nlon上，风场上的数据经纬都是一维的，应该用什么办法解决？多谢了！！！

兰溪之水

    nlat = 384 ;
    nlon = 320 ;

这是一维的？

infocert

兰溪之水 发表于 2013-6-21 08:24

                               
登录/注册后可看大图

nlat = 384 ;
    nlon = 320 ;

但我用风矢量图的脚本来套这个数据，总报错
这个数据跟风数据的区别在于ULAT和ULONG，在风数据里，这两个是一维的表示成lat(nlat),lon(nlon)，而这里确是二维的，我的脚本如下：

1. ncfile = "de-gr."+month+".nc"
   
2.   f     = addfile(ncfile, "r")
   
3.   ;data  = f->PS(0,:,:)
   
4.   u     = f->UVEL(0,0,:,:)
   
5.   v     = f->VVEL(0,0,:,:)
   
6. 
   
7. ;************************************************
   
8. ; create plot
   
9. ;************************************************
   
10.   wks_type = "png"
    
11.   wks_type@wkWidth = 1280
    
12.   wks_type@wkHeight = 1280
    
13.   wks = gsn_open_wks(wks_type,"uvvel-nh-"+month)
    
14.   gsn_define_colormap(wks, "BlWhRe")
    
15. 
    
16.   res                  = True                ; plot mods desired
    
17.   res@gsnMaximize      = True
    
18.   res@gsnAddCyclic     = True
    
19. 
    
20.   res@pmLabelBarHeightF=  0.100
    
21.   res@pmLabelBarWidthF =  0.60              ; default is 0.6
    
22. 
    
23.   mo = stringtoint(month)
    
24.   res@tiMainString     = "UVEL VVEL in "+monthArr(mo)
    
25.   res@gsnLeftString    = "desert-green"
    
26.   res@gsnRightString   = "units: "+u@units
    
27. 
    
28.   res@gsnPolar = "NH"
    
29. 
    
30.   res@vcRefMagnitudeF  = 1.0
    
31.   res@vcRefLengthF     = 0.045
    
32.   res@vcMinDistanceF   = 0.017
    
33. 
    
34.   res@vcRefMagnitudeF      = 5.0                 ; make vectors larger
    
35.   res@vcRefLengthF         = 0.2              ; ref vector length
    
36.   res@vcMonoLineArrowColor = True
    
37.   res@vcLineArrowColor      = "black"
    
38.   plot = gsn_vector_map(wks, u, v, res)

复制代码

兰溪之水

infocert 发表于 2013-6-21 10:09

                               
登录/注册后可看大图

但我用风矢量图的脚本来套这个数据，总报错
这个数据跟风数据的区别在于ULAT和ULONG，在风数据里，这两个 ...

你一直说错，那错误信息呢，大哥。。。

尽头的尽头

楼主把出错提示发上来

infocert

兰溪之水 发表于 2013-6-21 10:14

                               
登录/注册后可看大图

你一直说错，那错误信息呢，大哥。。。

呵呵，不好意思，错误如下

1. fatal:VectorPlotDraw: VVECTR - VECTOR NDC LENGTH TOO GREAT
   
2. fatal:VectorPlotDraw: error drawing vectors
   
3. fatal:VectorPlotDraw: draw error
   
4. fatal:PlotManagerDraw: error in plot draw
   
5. fatal:_NhlPlotManagerDraw: Draw error

复制代码

infocert

尽头的尽头 发表于 2013-6-21 10:22

                               
登录/注册后可看大图

楼主把出错提示发上来

不好意思，已经发了

尽头的尽头

infocert 发表于 2013-6-21 10:25

                               
登录/注册后可看大图

不好意思，已经发了

  res@vcRefMagnitudeF      = 5.0                 ; make vectors larger
  res@vcRefLengthF         = 0.2              ; ref vector length
这两句重复了，把下面的删掉试试，不行再根据数据改变这两句的数值

infocert

尽头的尽头 发表于 2013-6-21 10:33

                               
登录/注册后可看大图

res@vcRefMagnitudeF      = 5.0                 ; make vectors larger
  res@vcRefLengthF         ...

谢谢，我试试，但我觉得应该不是这两句的问题，这个脚本之前画过一个风适量的图，是没问题的，就算res的属性设置重复，也只会产生一个warning，不会产生error的

兰溪之水

ncfile = "de-gr."+month+".nc"
  f     = addfile(ncfile, "r")
  u     = f->UVEL(0,0,:,:)
  v     = f->VVEL(0,0,:,:)

;************************************************
; create plot
;************************************************
  wks_type = "png"
  wks_type@wkWidth = 1280
  wks_type@wkHeight = 1280
  wks = gsn_open_wks(wks_type,"uvvel-nh-"+month)
  gsn_define_colormap(wks, "BlWhRe")

  res                  = True                ; plot mods desired
  plot = gsn_vector_map(wks, u, v, res)



看看这样会报什么错~