求助！时间维连续的问题

liusun2010

之前有听说grads计算对时间要求特别严格，但是由于我要算的数据很大，用fortran计算不太合适，所以想问一下。我有55年天的数据，高度场，每年365天，2.29的数据已经舍去，已经做了31天的滑动平均。
我要计算5、6、7、8月的月平均，必须要用天的资料转成月平均，在单独分别输出55年5、6、7、8月的数据。
如果令a=不为1的常数，计算d ave（h,t=a,t=a+31），d ave外面是a=a+365的循环，这样计算可以吗？

dzxconan

应该d ave(h,t='a',t='a'+31)ba

liusun2010

dzxconan 发表于 2013-4-27 20:25

                               
登录/注册后可看大图

应该d ave(h,t='a',t='a'+31)ba

谢谢，我先试试

lqouc

真不知道谁跟你说fortran不适合算大数据，有点常识都知道fortran就是为了算大数据的好吧。楼上的方法可以，或者ave(h,t='a',t='a+31')'这样应该也行。

liusun2010

lqouc 发表于 2013-4-27 21:26

                               
登录/注册后可看大图

真不知道谁跟你说fortran不适合算大数据，有点常识都知道fortran就是为了算大数据的好吧。楼上的方法可以， ...

我知道是可以算，但会很久，之前我要做1825天的滑动平均，frotran算了14h+都没算出来，grads算了7h左右吧，就出来了

lqouc

liusun2010 发表于 2013-4-27 21:36

                               
登录/注册后可看大图

我知道是可以算，但会很久，之前我要做1825天的滑动平均，frotran算了14h+都没算出来，grads算了7h左右吧 ...

那是程序编的有问题，绝对可以改善的。
而且想不到你做那个的意义，1825点滑动，算完之后自由度就小很多，如果想滤波直接拿年平均数据好了。

liusun2010

'reinit'
'open f:\data\hgt\31t.ctl'
'set gxout fwrite'
'set fwrite f:\data\hgt\mon_5.dat'

'set x 1 144'
'set y 1 73'
'set z 1'

a=121
while (a<=1000)
say a
'd ave(hgt,t='a',t='a'+31)'
a=a+365
endwhile

'disable fwrite'
'reinit'
只画了3年的，又用fortran画了相同的3年，但画出的图有些不太一样。。。

liusun2010

'reinit'
'open f:\data\hgt\31t.ctl'
'set gxout fwrite'
'set fwrite f:\data\hgt\mon_5.dat'

'set x 1 144'
'set y 1 73'
'set z 1'

a=121
while (a<=1000)
say a
'd ave(hgt,t='a',t='a'+31)'
a=a+365
endwhile

'disable fwrite'
'reinit'
只画了3年的，又用fortran画了相同的3年，但画出的图有些不太一样。。。

liusun2010

'reinit'
'open f:\data\hgt\31t.ctl'
'set gxout fwrite'
'set fwrite f:\data\hgt\mon_5.dat'

'set x 1 144'
'set y 1 73'
'set z 1'

a=121
while (a<=1000)
say a
'd ave(hgt,t='a',t='a'+31)'
a=a+365
endwhile

'disable fwrite'
'reinit'
只画了3年的，又用fortran画了相同的3年，但画出的图有些不太一样。。。