- 积分
- 2130
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-9-10
- 最后登录
- 1970-1-1
![[] 粉丝数: 微博数: 新浪微博达人勋](source/plugin/sina_login/img/light.png)
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
将wrfout数据插值后标量场没有数据,求各位大神赐教!
warning:ContourPlotInitialize: no valid values in scalar field; ContourPlot not possible:[errno=1101]
脚本:
begin
DATADir = "./"
a=addfile("ctl.nc","r")
u1 = wrf_user_getvar(a,"U10",19)
v1 = wrf_user_getvar(a,"V10",19)
lon2d = wrf_user_getvar(a,"XLONG",19)
lat2d = wrf_user_getvar(a,"XLAT",19)
w1 = sqrt(u1*u1+v1*v1)
copy_VarMeta(u1,w1)
olon = fspan(119.125,122.125,13)
olat = fspan(31.625,35.125,15)
data1 = new((/13,15/),"float")
olon!0 = "lon"
olon@long_name = "lon"
olon@units = "degrees-east"
olon&lon = olon
olat!0 = "lat"
olat@long_name = "lat"
olat@units = "degrees-north"
olat&lat = olat
wa = rcm2rgrid_Wrap(lon2d,lat2d,w1,olat,olon,0)
delete(u1)
delete(v1)
printVarSummary(wa)
contour_w = gsn_csm_contour_map(wks,wa,res)
end
|
|