整层水汽通量gs

拖延症

撒花点赞

清梅浅影

支持一个！

Yangtzeriver

楼主应该给一些介绍啊，便于大家学习和掌握

Yangtzeriver

等有了金钱再来下载看看！

银酱赛高

楼主，我参照你的GS来做整层水汽积分，GRADS报了vint函数出错"too few levels in default file".
我的思路是，每个要素（qs，ps，u，v）都输出为一个文件，再读取文件的数据来生成整层水汽通量积分的图。每个文件都能正常读取并画图，但是使用vint函数就报错，请楼主指点一下。

这是我画图的gs：

1. 'reinit'
   
2. 'open e:\data\gs\hanchao\ps.ctl'
   
3. 'open e:\data\gs\hanchao\qs.ctl'
   
4. 'open e:\data\gs\hanchao\u.ctl'
   
5. 'open e:\data\gs\hanchao\v.ctl'
   
6. 
   
7. 'set grads off'
   
8. 'set grid off'
   
9. 
   
10. 'set lat 0 60'
    
11. 'set lon 80 150'
    
12. 'set z 1'
    
13. 'set t 1'
    
14. 'set mpdraw on'
    
15. 'set mpdset cnworld'
    
16. 'set map 15 1 3'
    
17. 'set poli on'
    
18. 
    
19. 'define a=vint(ps.1/100,u.3*qs.2/9.,300)'
    
20. 'define b=vint(ps.1/100,v.4*qs.2/9.8,300)'
    
21. 'define div=hdivg(a,b)*100000'
    
22. 
    
23. 'set gxout contour'
    
24. 'set lon 90 140'
    
25. 'set lat 15 30'
    
26. 'd smth9(smth9(div))'
    
27. 'd div'
    
28. 
    
29. 'printim e:\data\pic\hanchao\vint_VF.jpg white'
    
30. ;

复制代码

这是文件的ctl：

1. dset e:\data\hanchao\ps.dat
   
2. undef 9.999E+20
   
3. ydef  61  linear 0.000000  1
   
4. xdef  71  linear 80.000000 1.000000
   
5. tdef 1 linear 00Z22JAN2016 6hr
   
6. zdef 1 levels 100000
   
7. vars 1
   
8. ps 1 99 ps
   
9. endvars

复制代码

dset e:\data\hanchao\qs.dat
undef 9.999E+20
ydef  61  linear 0.000000  1
xdef  71  linear 80.000000 1.000000
tdef 1 linear 00Z22JAN2016 6hr
zdef 26 levels 100000 97500 95000 92500 90000 85000 80000 75000 70000 65000 60000 55000 50000 45000 40000 35000 30000 25000 20000 15000 10000 7000 5000 3000 2000 1000
vars 1
qs 26 99 qs
endvars

dset e:\data\hanchao\u.dat
undef 9.999E+20
ydef  61  linear 0.000000  1
xdef  71  linear 80.000000 1.000000
tdef 1 linear 00Z22JAN2016 6hr
zdef 26 levels 100000 97500 95000 92500 90000 85000 80000 75000 70000 65000 60000 55000 50000 45000 40000 35000 30000 25000 20000 15000 10000 7000 5000 3000 2000 1000
vars 1
u 26 99 uwnd
endvars

dset e:\data\hanchao\v.dat
undef 9.999E+20
ydef  61  linear 0.000000  1
xdef  71  linear 80.000000 1.000000
tdef 1 linear 00Z22JAN2016 6hr
zdef 26 levels 100000 97500 95000 92500 90000 85000 80000 75000 70000 65000 60000 55000 50000 45000 40000 35000 30000 25000 20000 15000 10000 7000 5000 3000 2000 1000
vars 1
v 26 99 vwnd
endvars

银酱赛高

补充一下，我猜想可能是lev设置得不对，所以设了几次z都是同样的错误提示

╭☆DJ乱儛シ

谢谢楼主！

TZ一叶舟

谢谢楼主，抱走程序参考了

xuehongyang

现在疯狂滴在查找相关资料，都怪自己平时不努力，加油，共勉！

xuehongyang

气象家园里好多厉害的资料，但是更厉害的是楼主