NCL中央差分求偏导，delta_R和pr维度均为（time lat lon) 出图不对，求大神指导

哈哈镜 · 发表于 2024-9-11 17:45:50

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

delta_R = (/pr - et - ro/)
dPERdP = center_finite_diff_n(delta_R, pr, False, 0, 0)
copy_VarMeta(pr, dPERdP)
R_P = dim_avg_n_Wrap(dPERdP, 0)

1733217252 · 发表于 2024-9-11 18:47:49

本科毕设做过相关的，后来也是搞不出来，老师用Fortran帮忙弄了，楼主如果有好办法可以把代码贴出来

内马尔 · 发表于 2024-9-11 19:51:43

P-E-R 是什么量

哈哈镜 · 发表于 2024-9-12 08:52:49

内马尔发表于 2024-9-11 19:51
P-E-R 是什么量

就是delta_R 是降水蒸散发径流的差

内马尔 · 发表于 2024-9-12 11:10:19

哈哈镜发表于 2024-9-12 08:52
就是delta_R 是降水蒸散发径流的差

是地面上的量吗

哈哈镜 · 发表于 2024-9-12 11:13:26

内马尔发表于 2024-9-12 11:10
是地面上的量吗

是的都是三维数据查不出是啥原因

内马尔 · 发表于 2024-9-12 13:15:19

本帖最后由内马尔于 2024-9-13 10:49 编辑

哈哈镜发表于 2024-9-12 11:13
是的都是三维数据查不出是啥原因

delta_R = (/pr - et - ro/)
dPERdP = center_finite_diff_n(delta_R, pr, False, 0, 0)
这两行儿的 pr 不是一个东西吧
说实话，不太明白你要算什么，但是提一点，ncl的那个函数只能在均匀坐标上进行运算，也就是等差数列，1，3，5，7....

哈哈镜 · 发表于 2024-9-12 14:40:45

内马尔发表于 2024-9-12 13:15
delta_R = (/pr - et - ro/)
dPERdP = center_finite_diff_n(delta_R, pr, False, 0, 0)
这两行儿 ...

目的是求降水对delta_R的贡献，所以先求偏delta_R/偏Pr
这两行的pr是同一个东西是降水所以能这么做吗?

浅色记忆 · 发表于 2024-9-12 16:22:55

大神，我也不会

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

[作图] NCL中央差分求偏导，delta_R和pr维度均为（time lat lon) 出图不对，求大神指导

登录后查看更多精彩内容~