- 积分
- 484
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-5-6
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin
data = addfile("air.mon.mean.nc","r")
T = data->air(336:696,:,:,:)
Tavgtime = dim_avg_Wrap(T(level|:,lat|:,lon|:,time|:))
Tavgtimelon = dim_avg_Wrap(Tavgtime(level|:,lat|:,lon|:))
Tavgtimelon!0 = "lev"
Tavgtimelon!1 = "lat"
result = new(73,integer)
lat=fspan(-90,90,73)
lat@units="degree_north"
result!0 = "lat"
result&lat = lat
do ilat = 0,72
do ilev = 1,15
if ((Tavgtimelon(ilev,ilat) < Tavgtimelon(ilev-1,ilat)) .and.(Tavgtimelon(ilev,ilat) < Tavgtimelon(ilev+1,ilat))) then
result(ilat) = ilev
end if
end do
end do
wks = gsn_open_wks("pdf","tropopouse")
res = True
plot = gsn_csm_xy(wks,result&lat,result,res)
end
出来的图是一条水平直线,麻烦熟练工们帮忙看看问题所在 |
|