- 积分
- 2397
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-9-10
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
我按照ncl官网的说明使用latRegWgt函数求权重,但未得到结果。
我的ncl脚本:
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
f = addfile ("~/LOVECLIM/STORAGE/LOVECLIM1.3/control/repo/output/atmos/atmym001950_001.nc","r")
ts=f->ts
lat=f->lat
wgt = latRegWgt(lat, "float", 0)
printVarSummary(wgt)
print(lat+" "+wgt)
print("sum(wgt)="+sum(wgt))
end
输出:
(0) latRegWgt: Expecting equally spaced latitudes
Variable: wgt
Type: float
Total Size: 128 bytes
32 values
Number of Dimensions: 1
Dimensions and sizes: [32]
Coordinates:
Number Of Attributes: 1
_FillValue : 1e+20
(0) -85.7606 1e+20
(1) -80.2688 1e+20
(2) -74.7445 1e+20
(3) -69.213 1e+20
(4) -63.6786 1e+20
(5) -58.143 1e+20
(6) -52.6065 1e+20
(7) -47.0696 1e+20
(8) -41.5325 1e+20
(9) -35.9951 1e+20
(10) -30.4576 1e+20
(11) -24.9199 1e+20
(12) -19.3822 1e+20
(13) -13.8445 1e+20
(14) -8.3067 1e+20
(15) -2.7689 1e+20
(16) 2.7689 1e+20
(17) 8.3067 1e+20
(18) 13.8445 1e+20
(19) 19.3822 1e+20
(20) 24.9199 1e+20
(21) 30.4576 1e+20
(22) 35.9951 1e+20
(23) 41.5325 1e+20
(24) 47.0696 1e+20
(25) 52.6065 1e+20
(26) 58.143 1e+20
(27) 63.6786 1e+20
(28) 69.213 1e+20
(29) 74.7445 1e+20
(30) 80.2688 1e+20
(31) 85.7606 1e+20
(0) sum(wgt)=1e+20
请问,为什么我使用wgt没有计算出权重呢?是哪里出现问题了呢?
|
|