进行计算结果全为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

盐酸兔子

fcode 发表于 2014-5-12 22:59
Debug 一下，看 call 子程序的语句是否执行了？

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

fcode

Debug 一下，看 call 子程序的语句是否执行了？

盐酸兔子

fcode 发表于 2014-5-8 22:23
子程序里或许需要对 F 进行一些修改。

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

又修改了好多次，发现程序的子程序不运行，是什么原因呢？

盐酸兔子

fcode 发表于 2014-5-8 22:23
子程序里或许需要对 F 进行一些修改。

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

嗯 知道了 一定勤加练习 谢谢版主

fcode

子程序里或许需要对 F 进行一些修改。

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

自己要去多写，多练。遇到问题了再问。不要一上来直接就问，这样不好，大大的不好。

盐酸兔子

fcode 发表于 2014-5-8 22:09
一样，没有区别。
[mw_shl_code=fortran,true]
x=3.0

针对我的程序这样修改
call sub2（F0,G1,Fk,c1,F）
就可以实现子程序和主程序的传递了吧

fcode

一样，没有区别。

[Fortran] 纯文本查看 复制代码

x=3.0
call sub( x )
write(*,*) x

subroutine sub(y)
  y = 4.0
end subroutine sub

就会输出 4.0

盐酸兔子

fcode 发表于 2014-5-8 22:00
实参，虚参传递，是基本的函数调用。

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

但是要是想将子程序的返还到主程序呢，怎么传递

fcode

实参，虚参传递，是基本的函数调用。

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

subroutine 和 function 本质是同一个东西，只是形式不同而已。