请问wmbarbmap函数该怎么对经纬度进行定义

爱地理不爱物理

我画出来的风场和实际位置对不上，[img][/img]

代码如下：

1. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
   
2. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
   
3. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
   
4. load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"
   
5. 
   
6. begin
   
7. 
   
8. 
   
9. ;-----------------读取数据-----------------------------
   
10. fpath = "D:/research/20220629/surf/wind-2914.txt"
    
11. 
    
12. 
    
13. nrow = numAsciiRow(fpath)
    
14.     data=asciiread(fpath, (/102,4/), "float")
    
15.     lat=data(:,2)           ;¶áèëÎ3¶è
    
16.     lon=data(:,3)           ;¶áèë¾-¶è
    
17.     d  =data(:,0)
    
18.     s  =data(:,1)
    
19. 
    
20. 
    
21. 
    
22. 
    
23. wks=gsn_open_wks("eps","d:/research/20220629/surf/wind-2914")
    
24.     gsn_define_colormap(wks,"BlueRed")      ;
    
25.     res=True
    
26.     res@gsnDraw                 = False
    
27.     res@gsnFrame                = False
    
28.     res@gsnMaximize=True
    
29.     res@mpFillOn                = True
    
30.     res@mpDataSetName="Earth..4"       ;ÖD1úμØí¼°üo¬ÔúÕa¸ö½DEarth..4μÄμØí¼¿a
    
31.     res@mpDataBaseVersion = "MediumRes";íøéÏμÄÄǸö′úÂëàïûóDÕa¾ä￡¬o|ÎòÕÛìúáËoþÃ2ÅÃ÷°×]
    
32.    
    
33.    
    
34.    
    
35.    
    
36.     res@mpAreaMaskingOn         = True
    
37.     res@mpMaskAreaSpecifiers    = (/"Zhejiang"/)
    
38.     res@mpFillAreaSpecifiers = (/"Zhejiang"/)
    
39.     res@mpOutlineSpecifiers=(/"Zhejiang"/)
    
40.     res@mpOutlineBoundarySets = "AllBoundaries"
    
41.    
    
42.     ;
    
43.     ;res@mpOutlineOn=False  ;»-±ß½ç
    
44.     ;res@mpMaskAreaSpecifiers    = (/"China","Taiwan","Disputed area between India and China","India:Arunachal Pradesh"/)
    
45.     res@mpGeophysicalLineThicknessF=2.0 ;Õaá½DDêÇÎaá˼ó′ֱ߽çoí1ú½çÏß
    
46.     res@mpNationalLineThicknessF=2.0
    
47.     res@mpMinLatF=27
    
48.     res@mpMaxLatF=31.5
    
49.     res@mpMinLonF=118
    
50.     res@mpMaxLonF=122.5
    
51. 
    
52.    
    
53.    
    
54.     res@mpLandFillColor="white"
    
55.     res@mpInlandWaterFillColor="white"
    
56.     res@mpProvincialLineColor      = "black"
    
57.     res@mpProvincialLineThicknessF = 2.
    
58.     res@mpShapeMode= "FreeAspect"
    
59.     res@mpFillDrawOrder="PreDraw"  ;post
    
60.    
    
61. 
    
62.    
    
63.     res@tmXBOn = True
    
64.      res@tmXTOn = False
    
65.      res@tmYLOn = True
    
66.      res@tmYROn = False
    
67.      res@tmXBPrecision=10
    
68.      res@tmYLPrecision=10
    
69.      res@tmXBMode = "Explicit"
    
70.      res@tmXBValues = (/118,120,122.5/)
    
71.      res@tmXBLabels = (/"118~S~o~N~E","120~S~o~N~E","122.5~S~o~N~E"/)
    
72.      res@tmYLMode = "Explicit"
    
73.      res@tmYLValues = (/27,29,31.5/)
    
74.      res@tmYLLabels = (/"27~S~o~N~N","29~S~o~N~N","31.5~S~o~N~N"/)
    
75.      res@tmXBLabelFontHeightF =0.015
    
76.      res@tmYLLabelFontHeightF =0.015
    
77.      
    
78.     map=gsn_csm_map(wks,res)
    
79. 
    
80. uv=wind_component(s/10,d,0);将风场分解为U;v
    
81. ;u = uv(0,:)*1.94386 ;转单位
    
82. u = uv(0,:)*2.5 ;转单位
    
83. ;v = uv(1,:)*1.94386
    
84. v = uv(1,:)*2.5
    
85. wmsetp("col", 2)   
    
86. wmsetp("wbs", .04)
    
87. wmsetp("WDF", 1)
    
88. wmbarbmap(wks,lat,lon,u,v)  ;wmbarbmap(wks,lat, lon,u, v)和wmbarb(wks,lon, lat,u, v)
    
89. 
    
90. shp3="D:/research/2004Hagupit/gis/浙江省/浙江省.shp"  
    
91. ;shp3="d:/research/2004Hagupit/rain-hour/jinhua_polygon.shp"  
    
92. lnres3        = True     
    
93. lnres3@gsLineColor      ="black"
    
94. lnres3@gsLineThicknessF = 1.5      
    
95. shp_plot3     =gsn_add_shapefile_polylines(wks,map,shp3,lnres3)  
    
96. ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;河流
    
97. ;shp4="/home/fc/ncl/lib/ncarg/nclscripts/chinamap/hyd1_4l.shp"  
    
98. ; lnres4        = True
    
99. ;  lnres4@gsLineColor      = "blue"
    
100.   ;lnres4@gsLineThicknessF = 1.0   
     
101. ; shp_plot4     =gsn_add_shapefile_polylines(wks,plot,shp4,lnres4)  
     
102. maximize_output(wks,False)      ;这句话不写，add shapefiles 就没意义！
     
103. ;-----------------------主要部分------------------------------------------------
     
104. draw(map)
     
105. frame(wks)
     
106. 
     
107. 
     
108. end

复制代码

爱地理不爱物理

爱地理不爱物理

已解决，draw(map)要在wmbarbmap之前