- 积分
 - 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  
 |   
 
 
 
 |