新手求FORTRAN怎么处理日平均

mayuki

新手FORTRAN苦手。已经处理了很多天，但是一直没有成功过。

-------------------------------------------------------------------------------------------


2012002        2012-1-2        0        18.8
                1        17.3
                2        21.1
                3        26.8
                4        30.6
                5        35.5
                6        42
                7        49.9
                8        64.1
                9        81.9
                10        70.9
                11        32.4
                12        27.9
                13        27.5
                14        30.9
                15        27
                16        30.8
                17        35.4
                18        37.7
                19        44.2
                20        33.6
                21        36
                22        50.4
                23        37.3
2012003        2012-1-3        0        26.7
                1        20.9
                2        25.7
                3        29
                4        33.9
                5        31.1
                6        43.4
                7        44.8
                8        46
                9        58.5
                10        83.5
                11        78.2
                12        116
                13        114.1
                14        86.8
                15        58.4
                16        53.2
                17        30.9
                18        20.6
                19        30.5
                20        42.5
                21        48
                22        46.3
                23        45.7
---------------------------------------------------

具体数据如下，一年的每天24小时的值，想要求的是这一年的日平均。

avday =0.
      do i =1,24
           read(1.*) a     !从文件中读入一个小时的读数，也可以是已有数组a()           
           avday=avday+a
      end do
      avday=avday/24.
      print*, “日平均",avday

用过这个，定义过数组，但是程序运行不了。然后数据文件中有很多缺测项，不知道该怎么写。
拜托拜托，希望有人能教教。

言深深

给你写一天的处理程序，具体需要针对你自己的数据，实现读取
x(24)代表的是24个时次的资料，自己想办法read进去，你的程序是read不进去的。下面的程序只会read第一列，自己修改一下。
s是对24个时次求和，ic是计数器。注意求和和计数器都是在非缺测情况下进行的，缺测标记为999，属于浮点型。
avag是计算平均。

do i=1,24
read(1,*)x(i)
enddo
s=0;ic=0
do i=1,24
if(x(i).ne.999.)then
s=s+x(i)
ic=ic+1
endif
enddo
avag=s/ic

iken

这个数据读进来就简单了哇。。用sum这个函数试试

mayuki

言深深 发表于 2014-8-31 13:29
给你写一天的处理程序，具体需要针对你自己的数据，实现读取
x(24)代表的是24个时次的资料，自己想办法rea ...

问个弱智的问题，我read数据是之read数据列还是把时间列一起读进去，我是要求的一年365天的小时平均，就是按小时值求平均的，所以最后输出的应该是24个值
谢谢回复

mayuki

iken 发表于 2014-8-31 23:05
这个数据读进来就简单了哇。。用sum这个函数试试

我把时间和浓度一起读进来然后分小时用sum再除以366?

iken

mayuki 发表于 2014-9-1 16:57
我把时间和浓度一起读进来然后分小时用sum再除以366?

日均是当天24小时的平均吧，有366个值，你说的sum/366为全年的日变化，有24个值？

lanxueren121

新手{:eb303:}{:eb303:}

mayuki

iken 发表于 2014-9-1 19:18
日均是当天24小时的平均吧，有366个值，你说的sum/366为全年的日变化，有24个值？

对 日变化，我用Matlab做出来了 ╮(╯▽╰)╭

lu19320

学习了，谢谢分享～！

悠然一梦123

对此我也很惆怅