- 积分
- 2913
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2019-1-19
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 supersupergirl 于 2019-8-29 14:40 编辑
如图所示,这是绘制风场的散度和涡度图 用填色的方法就都有数值 可是用绘制等值线的方法就只显示0 求问各位大神是什么原因啊 还有我画的散度和涡度为什么数值这么大呢 以下是我的脚本 求各位帮忙看下 拜托
xlat = wrf_user_getvar(a, "XLAT",0)
xlon = wrf_user_getvar(a, "XLONG",0)
ter = wrf_user_getvar(a, "HGT",0)
lat = xlat(:,0)
lon = xlon(0,:)
u = wrf_user_getvar(a,"ua",it) ; 3D U at mass points
v = wrf_user_getvar(a,"va",it) ; 3D V at mass points
;计算风场散度和涡度
u@_FillValue = 9.96921e+36
scale=1.e04
dv = new ( dimsizes(u), typeof(u), u@_FillValue )
dv@long_name = "divengence"
dv@units = "scaled"
dv = uv2dv_cfd (u,v,lat,lon, 2) * scale ;散度
copy_VarCoords(u, dv)
vr = new ( dimsizes(u), typeof(u), u@_FillValue )
vr@long_name = "vorticity"
vr@units = "scaled"
vr = uv2vr_cfd (u,v,lat,lon, 2) * scale ;涡度
copy_VarCoords(u, vr)
printVarSummary(dv)
printVarSummary(vr)
计算出来散度和涡度的最大最小值为:
dv : min=-26.9542 max=18.9527
vr : min=-19.3438 max=19.2699
感觉数值有些不合理 过于大了 请问是什么原因呢
以下是绘制等值线的脚本:
; Plotting options for dv
opts_dv = opts_xy
opts_dv@cnInfoLabelZone = 1
opts_dv@cnInfoLabelSide = "Top"
opts_dv@cnInfoLabelPerimOn = True
opts_dv@cnInfoLabelOrthogonalPosF = -0.00005
opts_dv@ContourParameters = (/ 10. /)
; Plotting options for vr
opts_vr = opts_xy
opts_vr@ContourParameters = (/ 5 /)
opts_vr@cnFillOn = True
opts_vr@cnLineColor="Red"
opts_vr@cnLineThicknessF = 2.0
opts_vr@gsnContourPosLineDashPattern=1
opts_vr@cnInfoLabelOrthogonalPosF=0.00
opts_vr@cnLineDashPattern=1
求问原因!!!
|
-
|