西太平洋暖池海温异常序列与500hpa位势高度场相关系数分布

勇敢的悲伤

海温异常与位势高度求相关的时候，是要把海温异常的距平值求出一个平均值，还是每年的距平值与位势高度相乘。这个好晕 啊，怎么编写呢？下面的是我自己写的，大家帮忙看看思路对不对
program main
implicit none
integer,parameter:: iix=145,iiy=73,mt=696,m=12,y=58
integer i,j,t,yr,mo,kk,k
real::a=0.0,avf,s,g,h
integer::anum=0
real atemp(y),sst(iix,iiy,mt),num(y),temp(iix,iiy,m,y),sst1(y),d(y),f(y),hgt(iix,iiy,mt),hgt1(iix,iiy,m,y),bnum(y),zhgt1(iix,iiy,y),b(iix,iiy,y),fm(iix,iiy),r(iix,iiy),zhgt(iix,iiy),c(iix,iiy),fz(iix,iiy),atemp1(y)
open(13,file='f:\1\hgt.grd',form='binary')
open(11,file='f:\1\sst2.grd',form='binary')
open(14,file='f:\1\temp1.grd',form='binary')
open(41,file='f:\1\xiangguan.grd',form='binary')
open(51,file='f:\1\xiangguan.txt')
do t=1,mt
read(11) ((sst(i,j,t),i=1,iix),j=1,iiy)
enddo
do yr=1,y
read(14) (atemp(yr))
enddo
do yr=1,y
if(atemp(yr)>1.0)then
anum=anum+1
endif
enddo
do yr=1,y
do mo=1,m
kk=(yr-1)*12+mo
do j=1,iiy
do i=1,iix
temp(i,j,mo,yr)=sst(i,j,kk)
enddo
enddo
enddo
enddo
do yr=1,y
num(yr)=0
do mo=6,8
do j=1,iiy
do i=1,iix
if(temp(i,j,mo,yr)<=-9.99e+08)then
num(yr)=num(yr)+1
endif
enddo
enddo
enddo
enddo
do yr=1,y
f(yr)=0
do mo=6,8
do j=1,iiy
do i=1,iix
if(temp(i,j,mo,yr)>-9.99e+08)then
f(yr)=f(yr)+temp(i,j,mo,yr)
endif
enddo
enddo
enddo
enddo
do yr=1,y
d(yr)=f(yr)/(iix*iiy*3.0-num(yr))
enddo
s=0.0
do yr=1,y
s=s+d(yr)
enddo
avf=s/y
do yr=1,y
sst1(yr)=(d(yr)-avf)
enddo
k=1
do yr=1,y
if(atemp(yr)>1.0)then
f(k)=sst1(yr)
k=k+1
endif
enddo
do k=1,anum
g=f(k)+g
enddo
g=g/anum
do t=1,mt
read(13) ((hgt(i,j,t),i=1,iix),j=1,iiy)
enddo
do yr=1,y
do mo=1,m
kk=(yr-1)*12+mo
do j=1,iiy
do i=1,iix
hgt1(i,j,mo,yr)=hgt(i,j,kk)
enddo
enddo
enddo
enddo
do j=1,iiy
do i=1,iix
do yr=1,y
bnum(yr)=0
do mo=6,8
if(hgt1(i,j,mo,yr)<=-9.9e+08)then
bnum(yr)=bnum(yr)+1
endif
enddo
enddo
enddo
enddo
do j=1,iiy
do i=1,iix
do yr=1,y
zhgt1(i,j,yr)=0
do mo=6,8
if(hgt1(i,j,mo,yr)>-9.9e+08)then
zhgt1(i,j,yr)=zhgt1(i,j,yr)+hgt1(i,j,mo,yr)
endif
enddo
enddo
enddo
enddo
do j=1,iiy
do i=1,iix
do yr=1,y
zhgt1(i,j,yr)=zhgt1(i,j,yr)/(3-bnum(yr))
enddo
enddo
enddo
do j=1,iiy
do i=1,iix
zhgt(i,j)=0.0
do yr=1,y
zhgt(i,j)=zhgt(i,j)+zhgt1(i,j,yr)
enddo
enddo
enddo
do j=1,iiy
do i=1,iix
zhgt(i,j)=zhgt(i,j)/y
enddo
enddo
do j=1,iiy
do i=1,iix
do yr=1,y
b(i,j,yr)=zhgt1(i,j,yr)-zhgt(i,j)
enddo
enddo
enddo
do j=1,iiy
do i=1,iix
fm(i,j)=0.0
do yr=1,y
fm(i,j)=fm(i,j)+(g*b(i,j,yr))
enddo
enddo
enddo

do j=1,iiy
do i=1,iix
fm(i,j)=fm(i,j)/anum
enddo
enddo
a=g**2
do j=1,iiy
do i=1,iix
do yr=1,y
c(i,j)=c(i,j)+b(i,j,yr)**2
enddo
enddo
enddo
do j=1,iiy
do i=1,iix
fz(i,j)=a*c(i,j)
enddo
enddo
do j=1,iiy
do i=1,iix
r(i,j)=fm(i,j)/sqrt(fz(i,j))
enddo
enddo
do j=1,iiy
do i=1,iix
write(41) r(i,j)
write(51,*) r(i,j)
enddo
enddo
end

mofangbao

看你这个程序会崩溃的。。。不过能自己努力先写已经不错啦，给你几个程序参考一下，弄明白算法就好办了，都是资料站的程序：
Smscorrelation 逐月序列求解相关系数fortran源程序
http://www.06climate.com/view/1119.html

smllcorrelation-逐月、逐季或其他异常序列不同月份之间的滞后超前nt月的相关系数fortran源程序
http://www.06climate.com/view/1118.html

sllyear-求逐月异常序列x(n,12)和y(n,12)(n是年)相同月份之间的滞后超前nt年的相关系数fortra
http://www.06climate.com/view/1117.html

Sllcorrelation-异常序列x(n)和y(n)之间的滞后超前的相关系数fortran源代码、源程序
http://www.06climate.com/view/1116.html

correlationmiss-求有缺省资料的两个序列x(n)和y(n)的相关系数r的fortran源代码
http://www.06climate.com/view/1114.html

lotina

你这程序太恐怖了！
你求他们两个异常场的相关就好啦！
海温异常场的每个格点跟相对应高度异常场的格点求相关

lotina

你这程序太恐怖了！
你求他们两个异常场的相关就好啦！
海温异常场的每个格点跟相对应高度异常场的格点求相关

勇敢的悲伤

lotina 发表于 2012-4-10 19:50

                               
登录/注册后可看大图

你这程序太恐怖了！
你求他们两个异常场的相关就好啦！
海温异常场的每个格点跟相对应高度异常场的格点求 ...

嘿嘿~~初次编写，不太会。。。但是这两个资料的分辨率不同，每个格点不对应的。而且还是58年的，这些混一块儿我就弄不清了。。。

lotina

那你要明白你主要目的啊。根据你的需要处理数据。如果是个点不对应，如果是我的话，就会考虑求区域平均，不过就是要根据你研究的对象来考虑了。

勇敢的悲伤

lotina 发表于 2012-4-10 22:59

                               
登录/注册后可看大图

那你要明白你主要目的啊。根据你的需要处理数据。如果是个点不对应，如果是我的话，就会考虑求区域平均，不 ...

谢谢，我再研究研究~~

勇敢的悲伤

mofangbao 发表于 2012-4-11 08:51

                               
登录/注册后可看大图

看你这个程序会崩溃的。。。不过能自己努力先写已经不错啦，给你几个程序参考一下，弄明白算法就好办了，都 ...

好的~谢谢~~

yunqicainiao

收集的好齐全

天道酬勤

程序怎么不能用。。。老出错