NCL计算Heat budget，出的图很奇怪，本人已死，请求大神帮助

燕灵飞雪 · 发表于 2015-12-14 13:48:10

燕灵飞雪发表于 2015-12-13 19:43
下面是画97年nonlinear heating的程序：

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

如果把乘的量级改一下，wata乘上100，即 nonlin=uata*1e6+vata*1e6+wata*100 ;据我输出数据看，乘之前uata和vata是1*10-7~10-9量级，wata是1*10-2~10-5量级。
图就变成了下面这样（值以及颜色没有人为设置，自动生成的）：

66kqkq · 发表于 2015-12-14 17:31:26

把计算差分的方法全部用ncl自带的函数替换试试呢

燕灵飞雪 · 发表于 2015-12-14 20:47:30

66kqkq 发表于 2015-12-14 17:31
把计算差分的方法全部用ncl自带的函数替换试试呢

你说的是center_finite_diff_n函数吗？之前试过的，但是这个函数会默认让缺测值也参与计算，不好用

风往北吹 · 发表于 2015-12-15 10:25:31

本帖最后由风往北吹于 2015-12-15 10:30 编辑

看错了~

井中月 · 发表于 2015-12-15 11:00:11

燕灵飞雪发表于 2015-12-14 20:47
你说的是center_finite_diff_n函数吗？之前试过的，但是这个函数会默认让缺测值也参与计算，不好用{:5_24 ...

用center_finite_diff_n算过，没发现问题啊

shaGuar · 发表于 2015-12-15 16:06:11

NCL直接用于计算，效率会不会不高？

66kqkq · 发表于 2015-12-15 16:21:17

center_finite_diff_n好像没有问题，是不是你的缺测值没有设置？
不是说这样会效率一些，而是这样可以避免一些可能会因为自己粗心造成的差分错误

燕灵飞雪 · 发表于 2015-12-15 17:21:52

井中月发表于 2015-12-15 11:00
用center_finite_diff_n算过，没发现问题啊

真的吗T_T 那你的资料有没有缺测值啊？算的时候有没有发现缺测值也参与计算了？

燕灵飞雪 · 发表于 2015-12-15 17:23:47

风往北吹发表于 2015-12-15 10:25
看错了~

啥。。？

燕灵飞雪 · 发表于 2015-12-15 17:24:22

shaGuar 发表于 2015-12-15 16:06
NCL直接用于计算，效率会不会不高？

效率是不太高。。计算一次起码二十分钟。。。

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