求每个月平均

ric_bob · 发表于 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。

生活在别处 · 发表于 2013-11-5 10:13:42

哇哦谢谢楼主慷慨分享~

mofangbao · 发表于 2013-11-5 12:38:43

原创分享必须支持！

mofangbao · 发表于 2013-11-5 12:38:51

原创分享必须支持！

zhouling · 发表于 2013-11-5 16:18:00

{:eb502:}{:eb502:}

jeffggg · 发表于 2013-11-5 17:06:25

谢谢搂住分享

Akbar · 发表于 2013-11-8 15:17:31

不错不错很好啊！

hikingbaby · 发表于 2014-4-3 09:23:02

正好是我要用的，楼主好人

愚遥 · 发表于 2014-9-21 20:38:57

有没有用于nc文件的呢？

qiaoyt · 发表于 2014-11-6 07:50:35

谢谢，很不错

		自动登录	找回密码
密码			立即注册

[grib] 求每个月平均

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

评分