- 积分
- 114
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-4-8
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2021-2-24 10:36:20
|
显示全部楼层
9楼的程序中就有关于地形叠加的命令
begin
lon = fspan(0,360,4320) ;gsn_xy
setfileoption("nc","MissingToFillValue",False)
stdlat=31 ;纬向剖面
elev= read_elev_data("E:/ncl/ETOPO5.DAT") ;地形数据
;printVarSummary(elev)
elev0 = 1013.25*(1-elev*0.0065/288.15)^5.25145 ;½«µØÐÎÊý¾ÝµÄµ¥Î»£ºmת»¯ÎªhPa
copy_VarCoords(elev,elev0)
;printVarSummary(elev0)
elev2 = elev0({stdlat},:)
;printVarSummary(elev2)
elev2@_FillValue= -9.96921e+36
elev3 = where(elev2.lt.5000,elev2@_FillValue,elev2)
copy_VarCoords(elev2,elev3)
;printVarSummary(elev3)
res.....................................(所画变量的描述语句)
xyres =True ;实现地形阴影
xyres@gsnMaximize = True
xyres@trYMaxF = 4000
xyres@gsnYRefLine = 4000.
xyres@gsnBelowYRefLineColor = "black"
xyres@gsnDraw = False
xyres@gsnFrame = False
plot变量
plot_hgt = gsn_csm_xy(wks,lon,elev2,xyres) ;plot阴影
overlay(plot(0),plot_hgt) ;叠加变量和阴影,注意overlay的顺序
draw plot
|
|