盐酸兔子 发表于 2014-5-6 20:53:10

进行计算结果全为0

program barnes
    implicit none
    external sub1,sub2
    integer,parameter::xx=4,yy=4
    real Fa(xx,yy),Fb(xx,yy),Fd(xx,yy),Fk(xx,yy),F0(xx,yy),F(xx,yy)
    real G1,G2
    integer c1,c2,i,j
    write(*,*)'please write G1,c1'
    read *,G1,c1
open(3,file='d:/barns/beijing.grd',form='binary')
read(3) ((Fk(i,j),j=1,yy),i=1,xx)
    call sub1(Fk,c1)
    call sub2(F0,G1,Fk,c1)
Fa=F
write(*,*)'please write G2,c2'
read*,G2,c2
call sub1(Fk,c2)
    call sub2(F0,G1,Fk,c2)
Fb=F
    do i=1,xx
       do j=1,yy
      Fd(i,j)=Fa(i,j)-Fb(i,j)
      write(*,*)Fd
   enddo
    enddo
    write(*,*)            
    OPEN(11,FILE='d:/barns/out2.grd',form='binary')
    write(11) ((Fd(i,j),j=1,yy),i=1,yy)
    close(11)
    close(8)
end

subroutine sub1(Ff,c)!计算初始值
implicit none
integer i,j,m,n,c,xx,yy
real F0(xx,yy),Ff(xx,yy),F2,F3,F4,Wk,Rk
common xx,yy
   do i=1,xx
   do j=1,yy
      do m=1,xx
         do n=1,yy
            Rk=sqrt((i-m)**2.0+(j-n)**2.0)
            Wk=exp(-(Rk**2.0)/(4.0*c))!权重函数计算
            F2=Wk*Ff(m,n)    !分子计算
            F3=F2+F3
            F4=F4+Wk
         enddo
      enddo
   F0(i,j)=F3/F4
   enddo
enddo
endsubroutine

subroutine sub2(F0,G,Ff,c)!计算F值
implicit none
integer i,j,m,n,c,xx,yy
real F0(4,4),Ff(4,4),F(4,4),Fs(4,4),D,F2,F3,F4,Rk,Wk,G
common xx,yy
do i=1,xx
   do j=1,yy
      do m=1,xx
         do n=1,yy
            D=Ff(m,n)-F0(m,n)
            Rk=sqrt((i-m)**2.0+(j-n)**2.0)
            Wk=exp(-(Rk**2.0)/(4.0*c*G))!权重函数计算
            F2=Wk*D    !分子计算
            F3=F2+F3
            F4=F4+Wk
         enddo
      enddo
   Fs(i,j)=F3/F4
   F(i,j)=Fs(i,j)+F(i,j)
   enddo
enddo
endsubroutine
G1=0.2 C=30000

G2=0.4 C=15000

vvt 发表于 2014-5-6 21:12:41

没有输入文件 beijing.grd，无法替你调试。

学会 Debug，你会自己解决所有的计算结果不符合预期的问题。
请看本站文章：http://www.fcode.cn/guide-44-1.html

盐酸兔子 发表于 2014-5-6 22:04:44

vvt 发表于 2014-5-6 21:12
没有输入文件 beijing.grd，无法替你调试。

学会 Debug，你会自己解决所有的计算结果不符合预期的问题。


一调试就跳出了，我又加了数据，修改了好多遍了，还是没找到原因，帮忙修改下吧 谢谢

fcode 发表于 2014-5-6 22:13:53

真是....说一句做一句。

你让别人调试，至少要让别人能按照你的意图跑一遍啊。一开始，不提供输入文件。然后，又不提供输入值。G1,c1 应该输入多少？

盐酸兔子 发表于 2014-5-6 22:18:36

fcode 发表于 2014-5-6 22:13
真是....说一句做一句。

你让别人调试，至少要让别人能按照你的意图跑一遍啊。一开始，不提供输入文件。然 ...

sorry，我疏忽了，以后注意G1=0.2 C=30000
G2=0.4 C=15000

fcode 发表于 2014-5-6 22:20:29

Fd(i,j)=Fa(i,j)-Fb(i,j)
write(*,*)Fd

Fd 是 Fa 减去 Fb，但是 Fa 和 Fb 都没有计算啊。

Fa=F 这算么？？ F 也没有计算啊。

盐酸兔子 发表于 2014-5-6 22:44:24

fcode 发表于 2014-5-6 22:20
Fd(i,j)=Fa(i,j)-Fb(i,j)
write(*,*)Fd



想将第一次计算的F数组赋给Fa，第二次计算的赋给Fb。Fa=F不可以直接将数组赋值过去么

fcode 发表于 2014-5-6 22:46:39

盐酸兔子 发表于 2014-5-6 22:44
想将第一次计算的F数组赋给Fa，第二次计算的赋给Fb。Fa=F不可以直接将数组赋值过去么 ...
可以直接赋值数组。问题是，你的 F 也没有值啊。你算 F 了么？

请注意，子程序里的 F，主程序不能直接使用。需要用实参虚参传递，或 common 共享，或 module 共享。

aliouying 发表于 2014-5-7 12:21:18

fcode 发表于 2014-5-6 22:46
可以直接赋值数组。问题是，你的 F 也没有值啊。你算 F 了么？

请注意，子程序里的 F，主程序不能直接使 ...

开个贴讲下如何发帖求助吧

fcode 发表于 2014-5-7 19:32:30

aliouying 发表于 2014-5-7 12:21
开个贴讲下如何发帖求助吧

这个版块头部有的，但是没人看
页: [1] 2 3
查看完整版本: 进行计算结果全为0