grads读取纬向风数据

太阳雨

请大家帮帮忙
我从ncep的逐日资料uwnd.1999.nc中读取t=163到t=178这段时间平均纬向风，这样写gs对不对
'reinit'
'sdfopen d:\99\uwnd.1999.nc'
'set fwrite d:\99\Umean.grd'
'set gxout fwrite'
   iz=1
   while(iz<=6)
     'set z 'iz''
     'define u1=ave(uwnd,t=163,t=178,1)'
     'd u1'
      iz=iz+1
   endwhile
'disable fwrite'
;

river

感觉问题不大，画图出来看看吧。如果只是需要图，不需要对数据进行输出的。你可以在gs里直接出图就行啊

1. 'reinit'
   
2. 'sdfopen d:\99\uwnd.1999.nc'
   
3. 'enable print d:\99\umean.gmf'
   
4. 'set grads off'
   
5. 'set grid off'
   
6. iz=1
   
7. while(iz<=6)
   
8. 'set z 'iz''
   
9. 'set t 1'
   
10. *****计算过程****
    
11. 'define u1=ave(uwnd,t=163,t=178,1)'
    
12. ****display之前可以控制图形的输出，加一些设置什么的*****
    
13. 'd u1'
    
14. 'draw title Umean on Z'iz''
    
15. 'print'
    
16. 'c'
    
17. iz=iz+1
    
18. endwhile
    
19. 'disable print'
    
20. ;

复制代码

可以拿这个出的图和你那个比较一下，看看是不是一样的

太阳雨

自己顶一下

做个霸气的木头

对与不对，楼主可以通过写ctl画图，自己验证~

太阳雨

做个霸气的木头 发表于 2012-10-2 14:32

                               
登录/注册后可看大图

对与不对，楼主可以通过写ctl画图，自己验证~

嗯，图是画过了，但还是不知道图对不对 不知道图是不是就应该这个样子

做个霸气的木头

太阳雨 发表于 2012-10-3 12:24

                               
登录/注册后可看大图

嗯，图是画过了，但还是不知道图对不对 不知道图是不是就应该这个样子

图对不对要看跟你的要求是否一致，你的程序是每次循环，求一个平面的时间平均，循环结束得到的数据是六层数据，时间点只有一个，是不是你的要求，自己可以对应了。图对不对，你可以寻找相关文章看人家的的图
。再者你的gs有个小瑕疵，在define前应该加上‘set t 1’和经纬度的定义，时间不定义有时会出错的。

太阳雨

做个霸气的木头 发表于 2012-10-3 13:42

                               
登录/注册后可看大图

图对不对要看跟你的要求是否一致，你的程序是每次循环，求一个平面的时间平均，循环结束得到的数据是六层 ...

嗯，我就是想求一个平面的时间平均，循环求六个面的，我把set t 1 加上试试，谢谢提醒

太阳雨

river 发表于 2012-10-5 09:24

                               
登录/注册后可看大图

感觉问题不大，画图出来看看吧。如果只是需要图，不需要对数据进行输出的。你可以在gs里直接出图就行啊可以 ...

谢谢您这么细心的帮助，我试过了，画出来是一样啊，谢谢

river

太阳雨 发表于 2012-10-5 19:03

                               
登录/注册后可看大图

谢谢您这么细心的帮助，我试过了，画出来是一样啊，谢谢

不客气。那就没问题了，以后自己多做做这些也就熟悉了

hello呀喵

正好用到，太棒了~~