fortranプログラムの中央値の求め方

Question

fortran90を用いてa(i-1),a(i),a(i+1)の中央値を求める方法について教えてください。

Answer 1

３個のデータの中央値でいいなら二回、最小値を取れば求まることになる。
以下では mask を使って不用なデータを排除している。

Fortran90 ではなく、Fortran95 の機能を多少使っているが。

program Hello
    implicit none
    real :: a(10), ans
    logical :: mask(size(a))
    integer :: i, k
    
    call random_number(a)
    
    do i = 2, 9
        print *, 'i=', i, a(i-1:i+1)
        mask = .false.
        mask(i-1:i+1) = .true.
        k = minloc(a, dim = 1, mask = mask)
        mask(k) = .false.
        ans = minval(a, mask)
        print *, ans
        print *
    end do       
end program Hello