PC correlation disabled

ashleywu

fortran程序运行时出现incrementally linked image- PC correlation disabled 这是什么错误？谁知道怎么改吗？

andrewsoong

矩阵设置成超过100x100？？？

mofangbao

楼主请看 提问的智慧 这篇帖子

兰溪之水

sqrt里面小于0 了吧。。。

river

好歹贴个程序出来看看啊

ashleywu

附上程序。
****相关*****文件1为单纯的一列数据,文件2为空间点,每个点都有同样长时间序列的数据*********
      parameter(nt=34,nx=144,ny=73)   
!nt::time series, nx:X方向格点数,ny:Y方向格点数
      real f(nt),ff(nx,ny,nt)    !分别存放两个原始数据
      real Y(nt)  
      open(1,file='sjxs.prewinter.grd',form='binary')    !!!!!!!!!
do i=1,nt
read(1) f(i)
print*,f(i)
enddo     
close(1)          !读入文件1
open(2,file='chi-summer200.grd',
     $form='binary') !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
do k=1,nt
do j=1,ny
do i=1,nx
read(2) ff(i,j,k)
! print*,ff(i,j,k)
enddo
enddo
enddo     
close(2)          !读入文件2
open(5,file='r.grd',form='binary') !结果输出
      open(6,file='t.grd',form='binary')
      
      do 50 jj=1,ny
do 40 ii=1,nx
do t=1,nt
Y(t)=ff(ii,jj,t)
enddo
      call correlation(f,Y,nt,r,tt)
      write(5) r
      write(6) tt
40    continue
50    continue
      close(5)
close(6)
end
ccc计算相关系数  CCCCCCCC
      subroutine correlation(X,Y,nt,RHO,TTT)
      real X(nt),Y(nt)
real*8 SX,SY,SXX,SYY,SXY,FN
SX=0.0
SY=0.0
SXX=0.0
SYY=0.0
SXY=0.0
FN=0.0
      do 20 i=1,nt
if(Y(i).ne.-9.99E+33)then
FN=FN+1
SX=SX+X(i)
SY=SY+Y(i)
SXX=SXX+X(i)*X(i)
SYY=SYY+Y(i)*Y(i)
SXY=SXY+X(i)*Y(i)
endif
20    continue
      if(FN.ne.0)then
RHO=(SXY-SX*SY/FN)/sqrt((SXX-SX*SX/FN)
     &  *(SYY-SY*SY/FN))
else
RHO=-9.99E+33
endif
      if(RHO.ne.-9.99E+33)then
TTT=RHO*sqrt(30.0)/sqrt(1-RHO*RHO)
else
TTT=-9.99E+33
endif
      end

ashleywu

mofangbao 发表于 2013-11-1 12:43

                               
登录/注册后可看大图

楼主请看 提问的智慧 这篇帖子

不用每次发帖都说别人不会问问题吧，每个人都是从新人开始的啊。诚心诚意请求帮助的。至于这么打击别人的积极性吗！@AIRES

lqouc

ashleywu 发表于 2013-11-1 18:10

                               
登录/注册后可看大图

不用每次发帖都说别人不会问问题吧，每个人都是从新人开始的啊。诚心诚意请求帮助的。至于这么打击别人的 ...

既然知道这个‘提问的智慧’为什么不尽量照着做呢，看来也不是一次说你不会问问题了，为什么只知道抱怨不知道改善自己呢。不是想打击你的提问积极性，是想帮你提高回答者的积极性，让你看你自己的这个问题你怎么回答。
报错很明显了就是你的sqrt出现了问题，一般是由于开方的数小于0造成的。
至于后面的就不用想了，前面搞对了再说吧，
这些基本的报错你在论坛都可以搜得到，有些帖子专门总结了常见的报错及解决方法。
多搜索你会发现根本用不着提问。
还有，让你把程序贴上来就是直接贴代码，不是上传附件，除了版主和管理员都要扣除积分，想想帮了你还要付出积分，谁愿意？这次我就先给你贴出来了。

lqouc

ashleywu 发表于 2013-11-1 18:07

                               
登录/注册后可看大图

附上程序。
****相关*****文件1为单纯的一列数据,文件2为空间点,每个点都有同样长时间序列的数据********* ...

你的程序，我没有数据，不知道具体哪一步出的错，但是我觉得是你的公式不对，不知道是不是你自己编的，那个sqrt里面应该是两个变量偏差的平方，但是你的表达看起来是先平方再做差。
你可以在论坛搜索一下，应该还有fortran的相关系数程序，对比一下。

ashleywu

lqouc 发表于 2013-11-1 18:31

                               
登录/注册后可看大图

你的程序，我没有数据，不知道具体哪一步出的错，但是我觉得是你的公式不对，不知道是不是你自己编的，那 ...

谢谢你的回复，我之前不知道下附件回答问题要扣分，我以为直接贴出来有的注释什么的语句看起来会让人头晕。平方这个我用过算短一点时间段的相关，是没有报错的。这次是把nt改成了34，原来是32，我在想是不是数组大小问题。ps:每次提问都在改进中，但是直接说你连问题都不会提，让人没有勇气提问了。