爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 9916|回复: 20

[grib] 求每个月平均

[复制链接]

新浪微博达人勋

发表于 2013-11-5 09:57:16 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
在论坛上找了好久都没找到,昨天自己尝试着写了一个,运行成功啦,看有没有人能用的到哒。
dset j:/1022/era40/dayave.grb
title j:/1022/era40/dayave.grb
undef -9.99e08
ydef 81 linear 15.000000 0.5
xdef 91 linear 60.000000 0.500000
tdef 16436 linear 18Z01sep1957 1dy
zdef 1 linear 1 1
vars 1
no2Tsfc  0 167,1,0  ** surface 2 metre temperature K
ENDVARS

这是ctl。


'reinit'
'set gxout fwrite'
'set fwrite j:\1022\erainterim\tem\monave.grb'

'open j:\1022\erainterim\tem\dayave.ctl'
'set x 1 91'
'set y 1 81'
'set z 1'
t0 = 1
m=0
while (t0 <= 12631)
    'set t 't0
    'q time'

    time1 = subwrd(result,3)
    say time1
    yearmon=substr(time1,6,7)
    say yearmon

    t1=t0
    t2=t1+40
    n=0
    while (t1< t2)

        'set t 't1
        'q time'
        time2 = subwrd(result,3)
    yearmon2=substr(time2,6,7)
    say yearmon
if (yearmon2=yearmon) ; n=n+1;t1=t1+1;endif;
if (yearmon2!=yearmon);break;endif
endwhile

  t5=t0+n-1
  say n
  'set t 't0
say t0
say t5
'display ave(no2Tsfc,t='t0',t='t5')'




  t0=t0+n
  m=m+1
          endwhile
  say m
          'disable fwrite'
          ;

这是.gs。

评分

参与人数 2金钱 +25 贡献 +7 收起 理由
言深深 + 10 + 2 赞一个!
mofangbao + 15 + 5 鼓励原创分享

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-5 10:13:42 | 显示全部楼层
哇哦 谢谢楼主慷慨分享~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

0
早起挑战累计收入
发表于 2013-11-5 12:38:43 | 显示全部楼层
原创分享必须支持!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

0
早起挑战累计收入
发表于 2013-11-5 12:38:51 | 显示全部楼层
原创分享必须支持!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-5 16:18:00 | 显示全部楼层
{:eb502:}{:eb502:}
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-5 17:06:25 | 显示全部楼层
谢谢搂住分享
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-8 15:17:31 | 显示全部楼层
不错不错   很好啊!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-4-3 09:23:02 | 显示全部楼层
正好是我要用的,楼主好人
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-9-21 20:38:57 | 显示全部楼层
有没有用于nc文件的呢?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-11-6 07:50:35 | 显示全部楼层
谢谢,很不错
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表