- 积分
- 952
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2020-2-28
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 樱花阿喵 于 2022-4-7 15:15 编辑
我从wrfout文件里面读取变量名为DTAUX3D和DTAUY3D的拖曳力,并插值到500hPa ,使用gsn_csm_vector画出图后数值太小,想要扩大1000倍,结果两个图完全不一样,但是填色底图是一样的,有没有大神可以指点一下(上面是扩大1000倍后的结果)
附上我的部分脚本,被注释掉的是放大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)
|
|