爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 304|回复: 8

[作图] 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)
111.png

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2024-9-11 18:47:49 | 显示全部楼层
本科毕设做过相关的,后来也是搞不出来,老师用Fortran帮忙弄了,楼主如果有好办法可以把代码贴出来
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2024-9-11 19:51:43 | 显示全部楼层
P-E-R 是什么量
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2024-9-12 08:52:49 | 显示全部楼层

就是delta_R 是降水 蒸散发 径流的差
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2024-9-12 11:10:19 | 显示全部楼层
哈哈镜 发表于 2024-9-12 08:52
就是delta_R 是降水 蒸散发 径流的差

是地面上的量吗
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2024-9-12 11:13:26 | 显示全部楼层

是的 都是三维数据  查不出是啥原因
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 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....
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 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是同一个东西 是降水 所以能这么做吗?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2024-9-12 16:22:55 | 显示全部楼层
大神,我也不会
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表