- 积分
- 4199
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-2-11
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2016-1-22 14:48:57
|
显示全部楼层
program main
implicit none
integer,parameter::n=13149
integer i,d
real r1(n),r2(n),r(n),nor(n),avg,sum,std,max
open(10,file='f:\study\mjo\RMM1RMM2(1979-2015).txt')
read(10,*) (r1(i),r2(i),i=1,n)
print*,r2(1)
do i=1,n
r(i)=sqrt(r1(i)**2+r2(i)**2)
end do
sum=0
do i=1,n
sum=sum+r(i)
end do
avg=sum/n
sum=0
do i=1,n
sum=sum+(r(i)-avg)**2
end do
std=sqrt(sum/n)
do i=1,n
nor(i)=(r(i)-avg)/std
end do
open(20,file='f:\study\mjo\nor.grd',form='binary')
write(20) (nor(i),i=1,n)
open(30,file='f:\study\mjo\nor.dat')
write(30,*) (nor(i),i=1,n)
!select the date of maxinum MJO intensity
max=nor(1)
d=1
do i=2,n
if(nor(i)>=max) then
max=nor(i)
d=i
end if
end do
print*,d,max
close(10)
close(20)
close(30)
end
|
|