绘制夏季多年平均场的等值线图等值线不连续

wet510

绘制多年夏季平均场的水汽通量，结果出来的等值线图不连贯，断断续续的，不知是脚本问题还是数据问题，向大家请教，期待高手指导！！

1. 'reinit'
   
2. 'set grads off'
   
3. 'set grid off'
   
4. 'set timelab off'
   
5. 'sdfopen E:\nc_data_monthly\air.mon.mean.nc'
   
6. 'sdfopen E:\nc_data_monthly\rhum.mon.mean.nc'
   
7. 'sdfopen E:\nc_data_monthly\uwnd.mon.mean.nc'
   
8. 'sdfopen E:\nc_data_monthly\vwnd.mon.mean.nc'
   
9. 'sdfopen E:\nc_data_monthly\pres.sfc.mon.mean.nc'
   
10. 'set lon 115 145'
    
11. 'set lat 35 60'
    
12. 'set z 1 8'
    
13. 'set t 157 756'
    
14. *file='D:\NEdraw\vapor_flux\cv_ave\ave.dat'
    
15. *'set gxout fwrite'
    
16. *'set fwrite 'file
    
17. 'define prs=lev'
    
18. 'define es=6.112*exp(17.2693822*(air.1-273.16)/(air.1-29.65))'
    
19. 'define qs=0.622*es/(prs-0.378*es)'
    
20. 'define qv=rhum.2*qs/100'
    
21. 'define uq=qv*uwnd.3*1000'
    
22. 'define vq=qv*vwnd.4*1000'
    
23. 'define uqs=vint(pres.5(lev=0),uq,300)'
    
24. 'define vqs=vint(pres.5(lev=0),vq,300)'
    
25. 'set z 1'
    
26. 'define sq=mag(uqs,vqs)'
    
27. 'set t 6 8'
    
28. 'define sq1=ave(sq,t+156,t=756,12)'
    
29. 'modify sq1 seasonal'
    
30. 'define sq2=sum(sq1,t=162,t=164)'
    
31. 'd sq2'

复制代码

下面是出的图：

丑了点，见谅

lotina

13行的      13.'set t 157 756'   这是什么？
27行的      27.'set t 6 8'

不明白

wet510

lotina 发表于 2012-3-30 22:16

                               
登录/注册后可看大图

13行的      13.'set t 157 756'   这是什么？
27行的      27.'set t 6 8'

13行表示时间设定在1961年1月1日0时到2010年12月31日0时
27行表示设定t为6到8月，是为了后面求夏季多年平均

lotina

wet510 发表于 2012-3-31 10:53

                               
登录/注册后可看大图

13行表示时间设定在1961年1月1日0时到2010年12月31日0时
27行表示设定t为6到8月，是为了后面求夏季多年平 ...

前面是逐日资料，后边变为了逐月资料？
如果后面还是逐日资料，T应该设置为151  242
28.'define sq1=ave(sq,t+156,t=756,12)'  这个你是求什么的啊？  12应该改为365吧

lotina

wet510 发表于 2012-3-31 15:58

                               
登录/注册后可看大图

自己发现问题了，计算出来的值在一些经纬度的值为空，导致无法绘图。。。

PS，我还是不明白你的资料
嘿嘿

wet510

lotina 发表于 2012-3-31 18:37

                               
登录/注册后可看大图

PS，我还是不明白你的资料
嘿嘿

necp1949年至2010年的逐月再分析资料

wet510

lotina 发表于 2012-3-31 12:23

                               
登录/注册后可看大图

前面是逐日资料，后边变为了逐月资料？
如果后面还是逐日资料，T应该设置为151  242
28.'define sq1=av ...

打开的就是逐月资料，我前面说的1961年1月1日0时指的就是1月份，157是1961年1月，158就是1961年2月，依次类推756就是2010年12月，前面回复有误，请谅解。至于后面'set t 6 8'，因为我要求的是夏季的多年平均，所以时间上只要每年的6－8月，故设置 t 的变化范围为6－8，ave 中的12是时间间隔，t的实际变化范围是6 18 30 42.......直到截止时间t =756，这样设置后就是每年的6月份求平均，7、8月一样，不知解释清楚没？

wet510

自己发现问题了，计算出来的值在一些经纬度的值为空，导致无法绘图。。。

gszywfc

楼主你好：
我的程序运行这两句时出错“incompatible grids”
'define uqs=vint(pres.5(lev=0),uwnd.1*q/9.8,300)'
'define vqs=vint(pres.5(lev=0),vwnd.2*q/9.8,300)'
是pres.5高斯网格的问题吗

gszywfc

解决了，是我下错文件了