- 积分
- 161
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-11-13
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 清木77 于 2022-7-7 18:12 编辑
resm@mpGeophysicalLineColor = -1 ;关闭自然边界
resm@pmTickMarkDisplayMode = "always" ;坐标标签添加度符号
res@gsnAddCyclic = True ;添加循环点,否则在0度或180度出现白条
***************************************************
;三个copy的区别
copy_VarMeta(from, to) ;复制包括经纬度、缺测值在内的所有属性
copy_VarCoords(from, to) ;仅复制各维度属性,如时间和经纬度
copy_VarAtts(from, to) ;仅复制除维度属性以外的信息
***************************************************
;平滑
(1)先使用参数设置修改
res@cnSmoothingOn = True
res@cnSmoothingDistanceF = 0.0001
res@cnSmoothingTensionF = -2.5
(2)函数选择
X = smth9_Wrap(x, 0.50, q, False) ;q=-0.25轻平滑,q=0.25重平滑
(3)最下策的选择
wrf_smooth_2d(x, num) ;num为平滑的次数*******************************
;区域平均
(1)纬度余弦为权重
lat = f->lat
rad = 4.0*atan(1.0)/180.0
clat = cos(lat*rad)
uAve_clat = wgt_areaave(u, clat, 1.0, 1)
(2)格点面积为权重
lat = f->lat
lon = f->lon
jlat = dimsizes(lat)
rad = 4.0*atan(1.0)/180.0
re = 6371220.0
rr = re * rad
dlon = abs(lon(2) - lon(1)) * rr
dx = dlon * cos(lat*rad)
dy = new(jlat, typeof(dx))
dy(0) = abs(lat(2) - lat(1)) * rr
dy(1:jlat - 2) = abs(lat(2:jlat - 1) - lat(jlat - 2)) * rr
area = dx*dy
uAve_area = wgt_areaave(u, area, 1.0, 1)
(3)高斯格点权重
gwgt = f->gwgt
uAve_gwgt = wgt_areaave(u, gwgt, 1.0, 1)
|
|