Fortran Coder

标题: 进行计算结果全为0 [打印本页]
作者: 盐酸兔子    时间: 2014-5-6 20:53
标题: 进行计算结果全为0
[Fortran] 纯文本查看 复制代码
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

QQ截图20140506205104.png (22.96 KB, 下载次数: 455)

计算结果

计算结果

beijing.zip

225 Bytes, 下载次数: 1

资料
作者: vvt    时间: 2014-5-6 21:12
没有输入文件 beijing.grd，无法替你调试。

学会 Debug，你会自己解决所有的计算结果不符合预期的问题。
请看本站文章：http://www.fcode.cn/guide-44-1.html
作者: 盐酸兔子    时间: 2014-5-6 22:04
vvt 发表于 2014-5-6 21:12
没有输入文件 beijing.grd，无法替你调试。

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

一调试就跳出了，我又加了数据，修改了好多遍了，还是没找到原因，帮忙修改下吧 谢谢
作者: fcode    时间: 2014-5-6 22:13
真是....说一句做一句。

你让别人调试，至少要让别人能按照你的意图跑一遍啊。一开始，不提供输入文件。然后，又不提供输入值。G1,c1 应该输入多少？
作者: 盐酸兔子    时间: 2014-5-6 22:18
fcode 发表于 2014-5-6 22:13
真是....说一句做一句。

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

sorry，我疏忽了，以后注意G1=0.2 C=30000
G2=0.4 C=15000
作者: fcode    时间: 2014-5-6 22:20
Fd(i,j)=Fa(i,j)-Fb(i,j)
write(*,*)Fd

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

Fa=F 这算么？？ F 也没有计算啊。
作者: 盐酸兔子    时间: 2014-5-6 22:44
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
盐酸兔子 发表于 2014-5-6 22:44
想将第一次计算的F数组赋给Fa，第二次计算的赋给Fb。Fa=F不可以直接将数组赋值过去么 ...

可以直接赋值数组。问题是，你的 F 也没有值啊。你算 F 了么？

请注意，子程序里的 F，主程序不能直接使用。需要用实参虚参传递，或 common 共享，或 module 共享。
作者: aliouying    时间: 2014-5-7 12:21
fcode 发表于 2014-5-6 22:46
可以直接赋值数组。问题是，你的 F 也没有值啊。你算 F 了么？

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

开个贴讲下如何发帖求助吧
作者: fcode    时间: 2014-5-7 19:32
aliouying 发表于 2014-5-7 12:21
开个贴讲下如何发帖求助吧

这个版块头部有的，但是没人看
作者: 盐酸兔子    时间: 2014-5-8 21:46
fcode 发表于 2014-5-6 22:46
可以直接赋值数组。问题是，你的 F 也没有值啊。你算 F 了么？

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

我在两个函数里使用common 但是总出现Error: The dimensions of this array have been defined more than once. 具体使用common怎么修改呢。谢谢版主
作者: fcode    时间: 2014-5-8 21:49
绝大多数错误，只要自己翻译一下错误提示就可以了。

错误就是，数组的维度定义了超过一次。一般数组只定义一次既可。

另外，不建议用 common，强烈建议放弃使用它
作者: 盐酸兔子    时间: 2014-5-8 21:56
fcode 发表于 2014-5-8 21:49
绝大多数错误，只要自己翻译一下错误提示就可以了。

错误就是，数组的维度定义了超过一次。一般数组只定义 ...

使用实参虚参传递呢，怎么修改，麻烦给我大概说说，是使用function么实行传递么
作者: fcode    时间: 2014-5-8 22:00
实参，虚参传递，是基本的函数调用。

就像你的 sub1 一样，把 FF 和 c 传递过去。

subroutine 和 function 本质是同一个东西，只是形式不同而已。
作者: 盐酸兔子    时间: 2014-5-8 22:07
fcode 发表于 2014-5-8 22:00
实参，虚参传递，是基本的函数调用。

就像你的 sub1 一样，把 FF 和 c 传递过去。

但是要是想将子程序的返还到主程序呢，怎么传递
作者: fcode    时间: 2014-5-8 22:09
一样，没有区别。
[Fortran] 纯文本查看 复制代码


x=3.0

call sub( x )

write(*,*) x



subroutine sub(y)

  y = 4.0

end subroutine sub

就会输出 4.0
作者: 盐酸兔子    时间: 2014-5-8 22:22
fcode 发表于 2014-5-8 22:09
一样，没有区别。
[mw_shl_code=fortran,true]
x=3.0

针对我的程序这样修改
call sub2（F0,G1,Fk,c1,F）
就可以实现子程序和主程序的传递了吧
作者: fcode    时间: 2014-5-8 22:23
子程序里或许需要对 F 进行一些修改。

实参虚参传递，是一个最基本的使用，你必须学会。（实际上，sub1里面你已经在使用了，为什么换一个sub2地方就不懂了呢？）

自己要去多写，多练。遇到问题了再问。不要一上来直接就问，这样不好，大大的不好。
作者: 盐酸兔子    时间: 2014-5-8 22:28
fcode 发表于 2014-5-8 22:23
子程序里或许需要对 F 进行一些修改。

实参虚参传递，是一个最基本的使用，你必须学会。（实际上，sub1里 ...

嗯 知道了 一定勤加练习 谢谢版主
作者: 盐酸兔子    时间: 2014-5-12 22:51
fcode 发表于 2014-5-8 22:23
子程序里或许需要对 F 进行一些修改。

实参虚参传递，是一个最基本的使用，你必须学会。（实际上，sub1里 ...

又修改了好多次，发现程序的子程序不运行，是什么原因呢？
作者: fcode    时间: 2014-5-12 22:59
Debug 一下，看 call 子程序的语句是否执行了？
作者: 盐酸兔子    时间: 2014-5-13 20:14
fcode 发表于 2014-5-12 22:59
Debug 一下，看 call 子程序的语句是否执行了？

执行了 我找到自己的错误了 我再努力修改下 谢谢您



欢迎光临 Fortran Coder (http://bbs.fcode.cn/) Powered by Discuz! X3.2