ncl画图报错，求助

xzx199913 · 发表于 2022-3-1 13:05:01

有大佬可以告诉我这是哪出问题了吗，谢谢

你是最美七月天 · 发表于 2022-3-1 13:25:20

month to annual没有提供单位属性

xzx199913 · 发表于 2022-3-1 13:32:07

你是最美七月天发表于 2022-3-1 13:25
month to annual没有提供单位属性

大佬，我应该怎么改啊

你是最美七月天 · 发表于 2022-3-1 13:43:00

xzx199913 发表于 2022-3-1 13:32
大佬，我应该怎么改啊

脚本和资料发我，我改改发你？或者私信我？

xzx199913 · 发表于 2022-3-1 13:58:02

你是最美七月天发表于 2022-3-1 13:43
脚本和资料发我，我改改发你？或者私信我？

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

begin

finame = "./et_m_vic_2000_2019.nc4"
fi = addfile(finame,"r")

time = fi->time
time(date) = cd_calendar(time,0)
year = time(:,0)
idx_2000_2019 = ind(year.ge.2000.and.year.le.2019)
ts_ann = month_to_annual(ts(idx_2000_2019,:,:),1)

rc = regCoef_n(ispan(2000,2019,1),ts_ann,0,0)
rc = rc*10
rc@long_name = "Linear Trend"
rc@units = ts_ann@units+"/decade"
copy_VarCoords(ts_ann(0,:,:),rc)

df = new((/nlat,nlon/),"integer")
tval = new((/nlat,nlon/),"float")
b = new((/nlat,nlon/),"float")
df = onedtond(rc@nptxy-2,(/nlat,nlon/))
tval = onedtond(rc@tval,(/nlat,nlon/))
b = 0.5
prob = 1.-betainc(df/(df+tval^2),df/2.0,b)
copy_VarCoords(ts_ann(0,:,:),prob)

rc@mpDateSetName = "Earth..4"

plot = gsn_csm_contour_map(wks, ts_ann, rc)
end

你是最美七月天 · 发表于 2022-3-1 14:04:03

xzx199913 发表于 2022-3-1 13:58
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/c ...

资料发我一份

xzx199913 · 发表于 2022-3-1 14:05:32

jiaoben.ncl (987 Bytes, 下载次数: 0)

xzx199913 · 发表于 2022-3-1 14:07:48

et_m_vic_2000_2019.rar (27.66 MB, 下载次数: 12)

xzx199913 · 发表于 2022-3-2 10:58:47

大佬，我的脚本还有拯救的机会吗

coliane · 发表于 2022-3-3 16:27:09

因为你的数据文件里的time并没有用cd calendar可以读取的格式去写、我看你后面就是想用year这个变量，完全可以自己创建一个year这个数组供后面使用。
date=yyyymm_time(2000, 2019, "integer")
year=toint(str_get_cols(sprinti("%i", date), 0, 3))
用上面这两行替换掉
原来的time = fi->time
time(date) = cd_calendar(time,0)
year = time(:,0)这些

		自动登录	找回密码
密码			立即注册

[其他] ncl画图报错，求助