grads作图用 river.gs 作图结果黄河位置与经纬度不符？求解决

adnap

sq850

grads作图如上：
目的：9月3日—29日850 hPa平均水汽通量(100—130°E20—45°N)：
这是作图的gs程序

1. ***'reinit'
   
2. 'sdfopen e:\xu\ff\air.nc'
   
3. 'sdfopen e:\xu\ff\hgt.nc'
   
4. 'sdfopen e:\xu\ff\rhum.nc'
   
5. 'sdfopen e:\xu\ff\uwnd.nc'
   
6. 'sdfopen e:\xu\ff\vwnd.nc'
   
7. 'enable print e:\xu\ff\sq850.gmf'
   
8. 'set parea 0.5 10.5 1 8'
   
9. 'set grads off'
   
10. 'set grid off'
    
11. 'set mpdraw off'
    
12. 'run river.gs'
    
13. 'set lon 70 140'
    
14. 'set lat 20 70'
    
15. 'set xlint 5'
    
16. 'set ylint 5'
    
17. 'set lev 850'
    
18. 'set t 1'
    
19. 
    
20. 'define g=9.8'
    
21. 'define es850=(6.112*exp((17.67*(air.1-273.15))/(air.1-29.65)))'
    
22. 'define qs850=(0.62197*es850/(850-0.378*es850))'
    
23. 'define q850=qs850*rhum.3/100'
    
24. 'define sq850=mag(q850*uwnd.4/g,q850*vwnd.5/g)*1000'
    
25. 'set gxout vector'
    
26. 'define uuu=ave(uwnd.4/g,t=3,t=29)'
    
27. 'define vvv=ave(vwnd.5/g,t=3,t=29)'
    
28. 'set arrscl 0.5 0.7'
    
29. 'set arrlab off'
    
30. 'd uuu;vvv'
    
31. *'d ave(sq850,t=3,t=29)'
    
32. 'print'
    
33. 'disable print'
    
34. 'reinit'
    
35. ;
    

***'reinit'<br /> 'sdfopen e:\xu\ff\air.nc'<br /> 'sdfopen e:\xu\ff\hgt.nc'<br /> 'sdfopen e:\xu\ff\rhum.nc'<br /> 'sdfopen e:\xu\ff\uwnd.nc'<br /> 'sdfopen e:\xu\ff\vwnd.nc'<br /> 'enable print e:\xu\ff\sq850.gmf'<br /> 'set parea 0.5 10.5 1 8'<br /> 'set grads off'<br /> 'set grid off'<br /> 'set mpdraw off'<br /> 'run river.gs'<br /> 'set lon 70 140'<br /> 'set lat 20 70'<br /> 'set xlint 5'<br /> 'set ylint 5'<br /> 'set lev 850'<br /> 'set t 1'<br /> <br /> 'define g=9.8'<br /> 'define es850=(6.112*exp((17.67*(air.1-273.15))/(air.1-29.65)))'<br /> 'define qs850=(0.62197*es850/(850-0.378*es850))'<br /> 'define q850=qs850*rhum.3/100'<br /> 'define sq850=mag(q850*uwnd.4/g,q850*vwnd.5/g)*1000'<br /> 'set gxout vector'<br /> 'define uuu=ave(uwnd.4/g,t=3,t=29)'<br /> 'define vvv=ave(vwnd.5/g,t=3,t=29)'<br /> 'set arrscl 0.5 0.7'<br /> 'set arrlab off'<br /> 'd uuu;vvv'<br /> *'d ave(sq850,t=3,t=29)'<br /> 'print'<br /> 'disable print'<br /> 'reinit'<br /> ;<br />

麻烦高手帮我看看是哪里出错了？

adnap

river.gs文件来源于”辣椒酱“同学
http://bbs.06climate.com/forum.p ... p;extra=&page=2

mofangbao

我在你帖子里面的cnriver不能用吗？

adnap

mofangbao 发表于 2012-6-19 13:36

                               
登录/注册后可看大图

我在你帖子里面的cnriver不能用吗？

可以用，但是作图（下图）有中国行政边界（PS：老师要求图中只出现黄河和长江标志），不知该如何解决？

MeteoInfo

adnap 发表于 2012-6-19 13:44

                               
登录/注册后可看大图

可以用，但是作图（下图）有中国行政边界（PS：老师要求图中只出现黄河和长江标志），不知该如何解决？ ...

用MeteoInfo输出一个只包含长江黄河的GrADS地图就可以了。

adnap

MeteoInfo 发表于 2012-6-19 13:48

                               
登录/注册后可看大图

用MeteoInfo输出一个只包含长江黄河的GrADS地图就可以了。

手头这个活要的比较紧，楼上手上有没有“只包含长江黄河的适合GrADS作图的底图文件”，万分感谢........
MeteoInfo还没来得及学呢

MeteoInfo

adnap 发表于 2012-6-19 13:54

                               
登录/注册后可看大图

手头这个活要的比较紧，楼上手上有没有“只包含长江黄河的适合GrADS作图的底图文件”，万分感谢....... ...

现做一个吧：
river_china (918 Bytes, 下载次数: 52)

2012-6-19 14:00 上传

点击文件名下载附件

mofangbao

用MeteoInfo导出地图很容易，点几下鼠标而已，楼主应该自己去试试的，你先去试试吧，我做了一个，其实通过软件自带的帮助手册就能做出来了，论坛上也有教程：http://bbs.06climate.com/forum.php?mod=viewthread&tid=893，方法和这个帖子一样，只是最后选择输出GrADS Map File就行了

mofangbao

MeteoInfo 发表于 2012-6-19 14:01

                               
登录/注册后可看大图

现做一个吧：

额。。。你比我快了

adnap

多谢8楼和9楼，刚已经下载了MeteoInfo软件（一直听说MeteoInfo很强大，算是领略到了），准备开始学习