GrADS中斜温图的绘制（Skew-T）

圣诞2011 · 发表于 2012-6-13 21:59:33

本帖最后由圣诞2011 于 2012-6-13 22:03 编辑

传说楼主想问你一个问题：为什么在求风向的时候，选择那么多的if-endif呢？实际上只是执行了ugrdprs>=0 ,vgrdprs>=0这个if-endif。虽然ugrdprs 和 vgrdprs 是数据文件中的变量，但是在if-endif中的这两个量却是作为描述性变量出现的，而且它们没有初始化，所以它们的值等于变量名本身，而在grads里面字符变量总是大于数字的，所以四个if-endfi语句其实仅仅只是执行了其中一个而已。刚好比较凑巧的是其实不管纬向、经向风是大于0或者是小于0，风向都可以用'define dir=180+atan2(ugrdprs,vgrdprs)*180/pi'求出来，楼主其他3个求风向的表达式其实是不正确的。原因在于在grads里面函数atan2(a,b)与我们之前在高中教材中求反正弦函数的结果是不同的。
atan2 of (1.73205080756888,1) = 60.0000
atan2 of (1.73205080756888,-1) = 120.0000
atan2 of (-1.73205080756888,-1) = -120.0000
atan2 of (-1.73205080756888,1) = -60.0000
以上是我本人的一些见解，可能有不对的地方，还希望大家指正！

超台十八 · 发表于 2012-6-13 22:04:30

ncl 吧... grads 做这些复杂的东西真的惨了... grads 除了互动操作模式, 真的没多大好处了.

冬日情愫+ · 发表于 2012-6-13 22:22:14

问个问题：就是如果有一个仅在z维变化的一维变量（x,y,t固定），如果我想改变这个变量在某个层次(z=××)中的值，该怎么办呢？

超台十八 · 发表于 2012-6-13 23:15:34

冬日情愫+ 发表于 2012-6-13 22:22

登录/注册后可看大图

问个问题：就是如果有一个仅在z维变化的一维变量（x,y,t固定），如果我想改变这个变量在某个层次(z=××)中 ...

fortran 改呀... grads 不清楚可不可以.

冬日情愫+ · 发表于 2012-6-13 23:51:44

当然可以用fortran。但是不知道可不可以在grads里面实现。比如在楼主中涉及到风向的确定，如果找不到一个所有条件下都适用的计算式的时候，要采用楼主原始的想法那样（通过判断某一层中变量满足某种条件是定义不同的变量），不知道该怎么实现。

rosewei · 发表于 2012-7-5 11:51:21

飞雪残剑发表于 2011-12-2 23:04

登录/注册后可看大图

我这次弄了，还是不行，你看看我的程序问题在哪里？谢啦！程序见附件！我用的资料是NCEP fnl格式的！

哈哈，下载了

shuihu88888 · 发表于 2012-7-18 16:07:49

威武！用下试试~~

土豆笑公主 · 发表于 2012-8-2 16:21:18

正好需要。受教了，收获很大，顶

zhf_ahqx · 发表于 2012-7-31 16:58:37

运行过程中出现"Data Request Warning: Request is completely outside file limits"，请问怎么解决？

alexmitsui · 发表于 2012-8-12 21:01:56

感谢分享！十分感谢！

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

[分享资料] GrADS中斜温图的绘制（Skew-T）