爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 6023|回复: 0

为什么使用vector画出的放大1000倍后的图不一样

[复制链接]
发表于 2022-4-7 15:13:49 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
本帖最后由 樱花阿喵 于 2022-4-7 15:15 编辑

我从wrfout文件里面读取变量名为DTAUX3D和DTAUY3D的拖曳力,并插值到500hPa ,使用gsn_csm_vector画出图后数值太小,想要扩大1000倍,结果两个图完全不一样,但是填色底图是一样的,有没有大神可以指点一下(上面是扩大1000倍后的结果) 图片1.png

图片2.png

附上我的部分脚本,被注释掉的是放大1000倍的代码
dragu = f1[:]->DTAUX3D
dragv = f1[:]->DTAUY3D


p1 = wrf_user_getvar(f1, "pressure", -1)

;;;interp to 500hPa
dragu_plane = wrf_user_intrp3d(dragu,p1,"h",500,0.0,False);;;time ,lat,lon
dragv_plane = wrf_user_intrp3d(dragv,p1,"h",500,0.0,False)


delete(dragu)
delete(dragv)
delete(p1)

;;;插值到ERA5网格点上
lat2d = f1[0]->XLAT(0,:,:)
lon2d = f1[0]->XLONG(0,:,:)
hgt = f1[0]->HGT(0,:,:)

lat = fspan(25, 46, 85)
lon = fspan(62, 107, 181)
dugrid = rcm2rgrid(lat2d, lon2d, dragu_plane, lat, lon, 0)
dvgrid = rcm2rgrid(lat2d, lon2d, dragv_plane, lat, lon, 0)
hgtgrid = rcm2rgrid(lat2d, lon2d, hgt, lat, lon, 0)

;;;求时间平均
du_avgt = dim_avg_n(dugrid, 0)
dv_avgt = dim_avg_n(dvgrid, 0)
dw_avgt = sqrt(du_avgt^2+dv_avgt^2)
;du_avgt = du_avgt*1000
;dv_avgt = dv_avgt*1000
;dw_avgt = dw_avgt*1000

;;;赋值坐标变量
du_avgt!0 = "lat"
lat@units = "degrees_north"
du_avgt&lat = lat
du_avgt!1 = "lon"
lon@units = "degrees_east"
du_avgt&lon = lon


copy_VarCoords(du_avgt, dv_avgt)
copy_VarCoords(du_avgt, dw_avgt)
copy_VarCoords(du_avgt, hgtgrid)

密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表