求每个月平均

ric_bob

在论坛上找了好久都没找到，昨天自己尝试着写了一个，运行成功啦，看有没有人能用的到哒。
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。

生活在别处

哇哦 谢谢楼主慷慨分享~

mofangbao

原创分享必须支持！

mofangbao

原创分享必须支持！

zhouling

{:eb502:}{:eb502:}

jeffggg

谢谢搂住分享

Akbar

不错不错   很好啊！

hikingbaby

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

愚遥

有没有用于nc文件的呢？

qiaoyt

谢谢，很不错