爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 3247|回复: 3

[源代码] 求夏季1996年逐日平均

[复制链接]
发表于 2017-8-8 15:22:40 | 显示全部楼层 |阅读模式

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

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

x
原txt是有156各站点,1996年的6-8月逐日资料,现在相求每日的多站点平均,自己编了以下,程序,但是提示数组越界,原数据一共是14352行。想咨询下大神,有什么方法可以修改下源程序吗,比较急太感谢了program main
implicit none
integer,parameter::n=14352                  
!integer state(n),year(n),id(state_num)
integer state(n),year(n),mon(n),R(n),day(n)
real  var(n),sum
integer  i,j,k,m,it,iy,a

!**********************************************
open(30,file="sum2.txt")
do i=1,n
read(30,*)state(i),year(i),mon(i),day(i),R(i)
print*,R(i)
end do

open(50,file='sum4.txt')
        do it=1,92
          sum=0
                a=it
                do iy=92,14260,92        !155个92行,每隔92行计算一次,          
                    if(R(a)/=32766)then       
                    sum=sum+R(a)
                        a=a+iy
                  end if
                 end do
                sum=sum/156/10
          write(3)sum
    enddo
close(50)
end


密码修改失败请联系微信:mofangbao
发表于 2017-8-8 22:39:31 | 显示全部楼层
参考一下,假设需要求平均数据是R,共156*92个数据

RR=reshape(R,(/92,156/))
mean=sum(R,2)
mean=mean/156.0

mean为一维数组,共92个元素
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2017-8-29 14:49:30 | 显示全部楼层
゛绍〤z、 发表于 2017-8-8 22:39
参考一下,假设需要求平均数据是R,共156*92个数据

RR=reshape(R,(/92,156/))

还是不太明白这个函数的意义
密码修改失败请联系微信:mofangbao
发表于 2017-9-5 20:39:18 | 显示全部楼层
自由行 发表于 2017-8-29 14:49
还是不太明白这个函数的意义

将数组R重新排列为92行*156列的数组
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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