Функция Fortran не возвращает значение основной программе

Question

Я вызвал функцию integ из другой функции rter3. Я напечатал integ значения из функции integ, и она печатает допустимые значения. Но когда я печатаю значения integ в функции rter3, она печатает 0.00. Я перепробовал все, что мог, чтобы понять, почему это так, но не смог понять. Пожалуйста, помогите мне здесь. Спасибо.

function rter3(z,c1,nx1,xarray,xlocation,angl3r,gamr,probr,aaa,da1)
    use data
    integer:: z,o,nx1,lc
    double precision:: c1(nx1),w1,w2,rter3,angl3r,gamr,probr(100),aaa,da1,rati2,xarray(nx1),xlocation,div1
    double precision:: variable
    dx=0.1
    rati2=20.0
    !write(*,*) z,o,nx1,xlocation,angl3,gam
    w1=div1(c1(z-1),c1(z),dx)
    w2=div1(c1(z),c1(z+1),dx)
    variable=integ(nx1,c1,xarray,xlocation,dx,angl3r,gamr,probr,aaa,da1,w1,w2)
    rter3=rati2*variable
    write(*,*) "rter3=", rter3, variable, rati2, xlocation, gam, angl3, da1
end function rter3

function integ(nxx,hfield,xfield,loc1,diff,angl3,gam,prob,afield,dprob,ab2,af2)
    implicit none
    integer :: l,nxx,dui,kk,indx,indt
    double precision:: du(101),prob(100),am1,a,ap1,angl3,gam,diff,dprob,integ,function12,afield(100),xfield(nxx)
    double precision:: loc1,loc2,hfield(nxx),am2,ap2,ab2,af2
    !write(*,*) 'test loop', loc1, diff, angl3, gam, dprob

    do kk=1,100
        do l=1,nxx
            if (0.95*xfield(l)>=loc1-afield(kk).and. 1.05*xfield(l)<loc1-afield(kk)) then
                dui=l
                !call exit()
            endif
        enddo
        !dui= loc1-kk
        !write(*,*) "Entered integ", dui
        !write(*,*) dui
        if (dui<=3) then
            dui=3
        endif
        if (dui>= nxx-2) then
            dui=nxx-2
        endif
        am1=hfield(dui-1)
        am2=hfield(dui-2)
        a=hfield(dui)
        ap2=hfield(dui+2)
        ap1=hfield(dui+1)
        du(kk)=abs(prob(kk)*(function12(am2,am1,ap1,ap2,diff,angl3,gam,ab2,af2))*dprob)
    enddo

    integ=sum(du)
    write(*,*) "integration value=", integ, nxx, loc1, diff, angl3, gam, dprob
end function integ

Итак, это отдельные функции для вашей справки. Пожалуйста, дайте мне знать, если вы чувствуете, что что-то не так. Весь код очень большой, но если вы чувствуете, что вам нужно его увидеть, я прикреплю его. Еще раз спасибо.

Answer 1

Как функция rter3 «знает» интерфейс функции Integ? Находятся ли оба в одном модуле? Если нет, возможно, rter3 не знает интерфейс и использует неявную типизацию. Поскольку вы не использовали «неявное ни одно» в rter3, это могло бы быть сделано, и в этом случае биты integ были бы интерпретированы как целое число, в то время как они должны интерпретироваться как двойная точность - значение будет отображаться неправильно. Мои предложения: всегда размещайте свои процедуры в модуле. Тогда вы в порядке с процедурами в том же модуле ... компилятор может проверить согласованность аргументов. Если вы вызываете подпрограмму или функцию из основной программы или другого модуля, «используйте» этот модуль для получения проверки аргумента. Используйте "неявное нет". Если вы забыли включить неявное none в свой исходный код, дополнительно используйте опцию компилятора, которая делает то же самое (например, -fimplicit-none для gfortran).