有关求相关的程序求调试

碎云飘摇

！一个副高指数时间与全球冬季海温场进行相关的程序求调试（资料是60年的，求前28年相关）

1. 
   
2. 
   
3. program xiangguan                                                            
   
4. parameter(nt=60,nt1=28,in=180,jn=89,nk=1)
   
5. real an(in,jn,nk,nt),dt(nt),ir(nt),rr(in,jn,nk)
   
6. open(2,file='e:\av05.grd',form='binary')
   
7. do k=1,nt
   
8. read(2)dt(k)
   
9. enddo
   
10. 
    
11. open(1,file='e:\sst\sst.win.grd',form='binary')
    
12. do k=1,nt
    
13. do ik=1,nk
    
14. read(1)((an(i,j,ik,k),i=1,in),j=1,jn)
    
15. enddo
    
16. enddo
    
17. 
    
18. do ik=1,nk
    
19. do i=1,in
    
20. do j=1,jn
    
21.     do k=1,nt1
    
22.     ir(k)=an(i,j,ik,k)
    
23.     enddo
    
24.     call correlation(nt1,ir,dt,r)
    
25.     rr(i,j,ik)=r
    
26. enddo
    
27. enddo
    
28. enddo
    
29. open(3,file='e:\corr.grd',form='binary')
    
30. do ik=1,nk
    
31. do j=1,jn
    
32. do i=1,in
    
33. do k=1,nt1
    
34. if(an(i,j,ik,k).eq.32767.)then
    
35. rr(i,j,ik)=32767.
    
36. endif
    
37. enddo
    
38. enddo
    
39. enddo
    
40. enddo
    
41. do ik=1,nk
    
42. do j=1,jn
    
43. do i=1,in
    
44. write(*,*)rr(i,j,ik)
    
45. write(3)rr(i,j,ik)
    
46. enddo
    
47. enddo
    
48. enddo
    
49. close(3)
    
50. end
    
51. 
    
52. 
    
53. subroutine correlation(n,x,y,r)                                                                                                                                                !
    
54. real x(n),y(n)                                 
    
55. ave1=0.0; ave2=0.0; Var1=0.0; Var2=0.0                                   
    
56. do i=1,n                                                                 
    
57.    ave1=ave1+x(i)/real(n)                                               
    
58.    ave2=ave2+y(i)/real(n)                                               
    
59. enddo                                                                    
    
60. do i=1,n                                                                 
    
61.    Var1=Var1+(x(i)-ave1)**2                                             
    
62.    Var2=Var2+(y(i)-ave2)**2                                             
    
63. enddo                                                                    
    
64. tmp=0.0                                                                  
    
65. do i=1,n                                                                 
    
66.    tmp=tmp+(x(i)-ave1)*(y(i)-ave2)                                       
    
67. enddo                                                                    
    
68. r=tmp/sqrt(Var1*Var2)                                                   
    
69. end  
    
70. 
    
71. 
    

</font><br /> <br /> program xiangguan                                                            <br /> parameter(nt=60,nt1=28,in=180,jn=89,nk=1)<br /> real an(in,jn,nk,nt),dt(nt),ir(nt),rr(in,jn,nk)<br /> open(2,file='e:\av05.grd',form='binary')<br /> do k=1,nt<br /> read(2)dt(k)<br /> enddo<br /> <br /> open(1,file='e:\sst\sst.win.grd',form='binary')<br /> do k=1,nt<br /> do ik=1,nk<br /> read(1)((an(i,j,ik,k),i=1,in),j=1,jn)<br /> enddo<br /> enddo<br /> <br /> do ik=1,nk<br /> do i=1,in<br /> do j=1,jn<br />     do k=1,nt1<br />     ir(k)=an(i,j,ik,k)<br />     enddo<br />     call correlation(nt1,ir,dt,r) <br />     rr(i,j,ik)=r<br /> enddo<br /> enddo<br /> enddo<br /> open(3,file='e:\corr.grd',form='binary')<br /> do ik=1,nk<br /> do j=1,jn<br /> do i=1,in<br /> do k=1,nt1<br /> if(an(i,j,ik,k).eq.32767.)then<br /> rr(i,j,ik)=32767.<br /> endif<br /> enddo<br /> enddo<br /> enddo<br /> enddo<br /> do ik=1,nk<br /> do j=1,jn<br /> do i=1,in<br /> write(*,*)rr(i,j,ik)<br /> write(3)rr(i,j,ik)<br /> enddo<br /> enddo<br /> enddo<br /> close(3)<br /> end<br /> <br /> <br /> subroutine correlation(n,x,y,r)                                                                                                                                                !<br /> real x(n),y(n)                                  <br /> ave1=0.0; ave2=0.0; Var1=0.0; Var2=0.0                                   <br /> do i=1,n                                                                 <br />    ave1=ave1+x(i)/real(n)                                               <br />    ave2=ave2+y(i)/real(n)                                               <br /> enddo                                                                    <br /> do i=1,n                                                                 <br />    Var1=Var1+(x(i)-ave1)**2                                              <br />    Var2=Var2+(y(i)-ave2)**2                                              <br /> enddo                                                                    <br /> tmp=0.0                                                                  <br /> do i=1,n                                                                 <br />    tmp=tmp+(x(i)-ave1)*(y(i)-ave2)                                       <br /> enddo                                                                    <br /> r=tmp/sqrt(Var1*Var2)                                                    <br /> end  <br /> <br /> <br /> <font style=\"font-size:12px\">

言深深

数组长度貌似不匹配，你在call进去数组的时候，数组长度nt=60
但是子程序里面的数组长度变成了nt1=28

碎云飘摇

谢谢，仔细一看昨晚一点半发的哎，很感动
现已经调试好了，长度写乱了，循环也写乱了
居然是高手，群里面没认出，深表抱歉