整层大气的水汽通量散度计算

lulu

删除删除

river

可以定义17个数据来代替lev啊，而这17个数据随Z变化。比如这样
lev.1=1000;lev.2=925;lev.3=850;lev.4=700;以此类推一直到lev.17=10
然后在你的循环里定义的时候就改成'define p=lev.nz'。至于nz应该是要不要加单引号你试试看吧，我只提供思路，具体操作验证还是你来做吧

lulu

river 发表于 2014-9-24 18:27
可以定义17个数据来代替lev啊，而这17个数据随Z变化。比如这样
lev.1=1000;lev.2=925;lev.3=850;lev.4=700 ...

谢谢，我试试哈！

lulu

river 发表于 2014-9-24 18:27
可以定义17个数据来代替lev啊，而这17个数据随Z变化。比如这样
lev.1=1000;lev.2=925;lev.3=850;lev.4=700 ...

试了一下，发现还是没法做出来。运行是每问题的，也出来了grd文件，但是我用ctl一个个图检验，发现还是只出现了1000hpa的图。试了一下画925hpa的图，出来的错误提示如图C:\Users\xxxcb4360\Desktop

lulu

river 发表于 2014-9-24 18:27
可以定义17个数据来代替lev啊，而这17个数据随Z变化。比如这样
lev.1=1000;lev.2=925;lev.3=850;lev.4=700 ...

river

lulu 发表于 2014-9-25 18:56

你的'set t 'nt都没放在循环里面啊

lulu

river 发表于 2014-9-25 22:01
你的'set t 'nt都没放在循环里面啊

啊啊，居然没 看到啊啊，谢谢版友，我改改看啊

lulu

river 发表于 2014-9-25 22:01
你的'set t 'nt都没放在循环里面啊

改了一下，发现还是显示一样的错误，我觉得错误的意思是说我的grd数据出来的层次实际上只有一层油数据，其他层次都是缺省，所以显示上面屏幕截图那个 ‘z=2 1’，但是在程序中我查找不出错误是哪里。。。

river

lulu 发表于 2014-9-26 10:54
改了一下，发现还是显示一样的错误，我觉得错误的意思是说我的grd数据出来的层次实际上只有一层油数据， ...

rhum有17层的资料吗？你用的2.5*2.5的资料，相对湿度只到300hPa吧，你看看。还有一个问题，你的时间循环应该是在所有变量的最外层，而不是像你现在这样在每个变量的最外层。
还有一个建议，最后计算水汽通量散度的那一步完全可以去掉，因为前面已经计算了qu qv，画图的时候直接使用hdivg出图就行，没必要重复计算。减少计算量，也减少grads运行时间。
还有你修改了gs，就麻烦你贴上来，不然就是你很肯定你修改过的是一定正确的了

lulu

river 发表于 2014-9-26 12:18
rhum有17层的资料吗？你用的2.5*2.5的资料，相对湿度只到300hPa吧，你看看。还有一个问题，你的时间循环 ...

谢谢建议！我接着改改