Question

我已经制作了线性插值代码，但是当我尝试在主程序或其他子程序中输出该值时，它显示0.0000的值，但是当我在子程序中运行程序时，它获得正确的值。

写（，）'Interpolverdi ='，interpolverdi获取正确的值但当我运行写（，）'主interpolverdi'时，interpolverdi我得到的值为0.0000

我的输入文件工作正常。

我做错了什么？如何从子程序中获取特定值，然后在另一个子程序中使用该值？

感谢我能得到的所有帮助：）

!main.f90

program mainn

use inputt
use interpoll
!use mongo

 implicit none

 integer :: n, itmax,k,bc
 real(8) :: epsilon, Re, tmax,dt,h,fux,interpolverdi
 real(8), dimension (:), allocatable :: oo,u,v,p
 integer, dimension (:), allocatable :: nn

 !real(8), dimension(:,:), allocatable :: Uvel, Vvel, pres, psi


 call lesinput(n,itmax,k,epsilon,Re,tmax,dt,oo,nn,bc)



 call interpol(nn,oo,n)

 write(*,*)'main interpolverdi',interpolverdi


 end program mainn

  !interpoll.f90
 module interpoll

 implicit none

 contains


 subroutine interpol(nn,oo,n)
  integer :: nn(:), n,s,i
  real(8) :: oo(:),interpolverdi




  s = size(nn)

  do i=1,s-1

  if (n==nn(i)) then
   interpolverdi= oo(i)
  exit


 else if ((n>nn(i)) .and. (n<nn(i+1))) then
 interpolverdi= oo(i)+((oo(i+1)-oo(i))*(n-nn(i)))/(nn(i+1)-nn(i))
 exit

 else
 interpolverdi=nn(s)

 end if

end do

write(*,*) ' Interpolverdi=', interpolverdi
write(*,*)'n=',n

end subroutine interpol

 end module interpoll

Answer 1

您有两个不同的变量名为interpolverdi。一个在主程序中定义，一个在子程序中定义。

当您在主程序和子程序中执行write(*,*) interpolverdi时，它每次都会打印一个不同的变量。

如果您想直接从内部子程序更改主程序变量，请从子程序中删除interpolverdi的声明。

然而，更好的方法是在子例程中添加另一个参数

 subroutine interpol(nn,oo,n,interpolverdi)
  integer :: nn(:), n,s,i
  real(8) :: oo(:),interpolverdi

end subroutine interpol

我建议将参数或主程序变量重命名为其他名称，以便它们具有不同的名称。

BTW real(8)非常丑陋且不完全便携。如果您正在学习Fortran，请不要学习糟糕习惯。 Fortran 90 kind parameter Confusing double precision real in Fortran

从子程序中获取变量

1 个答案: