| 
 
	积分3058贡献 精华在线时间 小时注册时间2019-1-19最后登录1970-1-1 
 | 
 
| 
本帖最后由 supersupergirl 于 2019-8-29 14:40 编辑
x
登录后查看更多精彩内容~您需要 登录 才可以下载或查看,没有帐号?立即注册 
  
 如图所示,这是绘制风场的散度和涡度图 用填色的方法就都有数值 可是用绘制等值线的方法就只显示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
 
 求问原因!!!
 
 | 
 
  |