有人知道calculate_monthly_values的用法吗

xiahou90

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"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
begin
f=addfile("/date/rday/pday.nc","r")
n=f->pday
time=n&time
time@units="days since 1961-1-1 00:00:0.0"
xMonthAvg=new((/18628,72,128/),"byte")
xMonthAvg = calculate_monthly_values(n, "avg", 0, False)
xMonthSum = calculate_monthly_values(n, "sum", 0, False)
xMonthMin = calculate_monthly_values(n, "min", 0, False)
xMonthMax = calculate_monthly_values(n, "max", 0, False)
print(xMonthAvg)
end


报错
fatal:Assignment type mismatch, right hand side can't be coerced to type of left hand side
fatal:["Execute.c":8567]:Execute: Error occurred at or near line 12762 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl

fatal:["Execute.c":8567]:Execute: Error occurred at or near line 19 in file test.ncl

longlivehj

把xMonthAvg=new…那行去掉。

navy.blue

calculate_monthly_values不支持字节型运算，要么把xMonthAvg的定义删掉，要么变换类型。
还有，calculate_monthly_values函数只有最左边维是时间维的变量才能正常运行，如果还出错，就只能去查看错误提示里contributed.ncl里面的具体行了。

xiahou90

定义是我后来加的 还是出错 类型我从double到byte都试过 还是出错
话说提示Assignment type mismatch, right hand side can't be coerced to type of left hand side 不就是变量类型不匹配吗
最左边一维确实是时间 不会错

longlivehj

xiahou90 发表于 2014-6-10 16:10
定义是我后来加的 还是出错 类型我从double到byte都试过 还是出错
话说提示Assignment type mismatch, rig ...

定义删了也不行？
ncl一般不用显式定义变量。其实，我觉得你定义也有问题，函数计算出来的是月值，18628个月大概有1500多年！所以，感觉18626可能是n的时间维大小。

xiahou90

18628是n的日值，我这里要计算月值，改成732也是报错

xiahou90

问题已解决

lideqin

能说下怎么解决的吗?

xiahou90

查过$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl
改了变量类型，另外多算了一天。这两个错误

¤拉普兰德¤

xiahou90 发表于 2014-6-17 15:42
查过$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl
改了变量类型，另外多算了一天。这两个错误

你好，请问你最后把变量类型改成什么了，我用typeof(var)还是不行，谢谢