ncl如何画wrfout的垂直风场图

material1022

使用gsn_csm_pres_hgt函数，结果出错，求大神帮帮忙
(0)        gsn_csm_pres_hgt: Fatal: The first dimension of the input data must
(0)        have a coordinate variable called 'lev.'
(0)        Cannot create plot.
fatal:Illegal right-hand side type for assignment
fatal:["Execute.c":8128]:Execute: Error occurred at or near line 12473 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl


ncl脚本如下：

1. ;================================================;
   
2. ;  gsn_contour_2.ncl
   
3. ;================================================;
   
4. ; Concepts illustrated:
   
5. ;   - Drawing a basic contour plot using gsn_contour
   
6. ;   - Spanning the full color map for contour fill
   
7. ;   - Turning on color fill for a contour plot
   
8. ;   - Turning off contour line labels
   
9. ;   - Turning off the contour informational label
   
10. ;   - Changing the size/shape of a contour plot using viewport resources
    
11. ;   - Increasing the size of tickmark labels
    
12. ;
    
13. ;================================================;
    
14. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
    
15. load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
    
16. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
    
17. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
    
18. ; ================================================;
    
19. begin
    
20. ;=================================================;
    
21. ; open file and read in data
    
22. ;=================================================;
    
23. cfilename="wrfout_ncl.cfg"
    
24. cfils    =asciiread(cfilename,-1,"string")
    
25. 
    
26. files = systemfunc("ls "+cfils(0)+":00:00") + ".nc"
    
27. f = addfiles(files,"r")
    
28. nfil=dimsizes(files)
    
29. times  = wrf_user_list_times(f)
    
30. U10=f[:]->U10(:,:,:)
    
31. V10=f[:]->V10(:,:,:)
    
32. U=wrf_user_getvar(f,"ua",-1)
    
33. V=wrf_user_getvar(f,"va",-1)
    
34. W=wrf_user_getvar(f,"wa",-1)
    
35. Tc=wrf_user_getvar(f,"tc",-1)
    
36. nxy=dimsizes(U)
    
37. nt=nxy(0)
    
38. nz=nxy(1)
    
39. ny=nxy(2)
    
40. nx=nxy(3)
    
41. pa=wrf_user_getvar(f,"pressure",-1) ;full mode pressure(hpa)
    
42. 
    
43. 
    
44. pressure_levels=(/1000.,900.,850.,750.,700.,600.,500.,300.,200./)
    
45. nleves=dimsizes(pressure_levels)
    
46. 
    
47. VV=new((/nt,nleves,ny,nx/),float)
    
48. UU=new((/nt,nleves,ny,nx/),float)
    
49. WW=new((/nt,nleves,ny,nx/),float)
    
50. TTc=new((/nt,nleves,ny,nx/),float)
    
51. 
    
52. 
    
53.   VV=wrf_user_intrp3d(V,pa,"h",pressure_levels,0.0,False)
    
54.   UU=wrf_user_intrp3d(U,pa,"h",pressure_levels,0.0,False)
    
55.   WW=wrf_user_intrp3d(W,pa,"h",pressure_levels,0.0,False)
    
56.   TTc=wrf_user_intrp3d(Tc,pa,"h",pressure_levels,0.0,False)
    
57. 
    
58. wAve=avg(WW(0,:,:,10))
    
59. uAve=avg(UU(0,:,:,10))
    
60. scale=fabs(uAve/wAve)
    
61. wscale=WW*scale
    
62. copy_VarCoords(WW,wscale)
    
63. ;=================================================;
    
64. ; PLOT 2
    
65. ;=================================================;
    
66. 
    
67. 
    
68.   
    
69. 
    
70. 
    
71. 
    
72.   
    
73.   res                     = True
    
74. 
    
75. ;---This resource not needed in V6.1.0
    
76.   res@gsnSpreadColors     = True         ; use full colormap
    
77. res@vcRefAnnoOn         = True      ; turns off the ref vector
    
78.   res@vcRefMagnitudeF  = 5.0                       ; add a reference vector
    
79.   res@vcRefLengthF     = 0.045                     ; what the ref length is
    
80.     res@vcGlyphStyle     = "CurlyVector"             ; turn on curley vectors
    
81. 
    
82.   
    
83. do ifle=0,0
    
84.   wks  = gsn_open_wks("png",cfils(1)+"/gsi_UW_"+times(ifle))
    
85.   gsn_define_colormap(wks,"gui_default")
    
86.   
    
87.   
    
88.   plot=gsn_csm_pres_hgt_vector(wks,TTc(ifle,:,:,10),UU(ifle,:,:,10),wscale(0,:,:,10),res)
    
89.   end do
    
90.   
    
91. 
    
92.   
    
93. 
    
94. 
    
95.   end
    

复制代码

藕茭白空心菜

楼主解决了吗

梁静茹93

我看到的一个问题是，wrf_user_intrp3d使用错误，别的就不清楚了，你可以看下

花小菜

请问楼主解决了吗，我跟你遇到了一样的问题，不知道该怎么解决

subtropical

画wrfout数据的图不是用gsn_csm_pres_hgt_vector吧，你可以去找下专门wrf画图的命令

yangjunhua

一样的问题，搞定了吗，楼主

recite

最开始就不能按照层次插值 而且不能用gsn开头的函数画WRFOUT的剖面
你可以看看 http://www2.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/NCL_examples.php