- 积分
- 6
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2020-11-28
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
如图,为什么等值线会出现交叉?
我的脚本如下:
a = addfile("/gpfsnewhome/mayxia/WANGH/Build_WRF/WRF-4.1.2/test/em_real/wrfout_d01_2016-01-16_00:00:00","r")
xlat = a->XLAT(0,:,:)
xlon = a->XLONG(0,:,:)
t = 0
time = wrf_user_getvar(a,"Times",t)
hgt = wrf_user_getvar(a,"geopt",t)
hgt = hgt/9.8
tc = wrf_user_getvar(a,"tc",t)
p = wrf_user_getvar(a,"pressure",t)
level = 500
hgt_interp = wrf_user_interp_level(hgt,p,level,True)
tc_interp = wrf_user_interp_level(tc,p,level,True)
hgt_interp@lon2d = xlon
hgt_interp@lat2d = xlat
tc_interp@lon2d = xlon
tc_interp@lat2d = xlat
wrf_smooth_2d(hgt_interp,2)
wrf_smooth_2d(tc_interp,2)
;=====================================================================================================================
wks = gsn_open_wks("pdf","polar")
gsn_define_colormap(wks,"GMT_globe")
;=====================================================================================================================
res1 = True
res1@gsnLeftString = ""
res1@gsnRightString = ""
res1@gsnAddCyclic = True
res1@gsnPolar = "NH"
res1@mpDataSetName = "/gpfsnewhome/mayxia/WANGH/heatwave/database/Earth..4"
res1@mpDataBaseVersion = "MediumRes"
res1@mpMinLatF = 50
res1@mpCenterLonF = 90
res1@mpOceanFillColor = 117
res1@mpLandFillColor = 255
res1@cnLinesOn = True
res1@pmTickMarkDisplayMode = "Always"
;=====================================================================================================================
plot = gsn_csm_contour_map_polar(wks,tc_interp,res1)
|
-
|