subprogramas - funciones intrinsecas fortran

¿Cómo pasar una función con múltiples argumentos a una subrutina que espera una función con un solo argumento?

Question

Tengo una subrutina (ejemplo mínimo)

subroutine treatfunction(f,input,output)
   external, real::f
   real, intent(in):: input
   real, intent(out):: output

   output = f(input) + f(1.0)  ! i.e. f has only one argument
end subroutine

y una función con dos argumentos

real function fun(x,a)
   real,intent(in)::x,a

Ahora, por un tiempo determinado, en tiempo de ejecución, quiero pasarme la fun para treatfunction . Entonces, idealmente, me gustaría llamar algo así como

call treatfunction(fun(:,a=a0), input=myinput, output=myoutput)

¿Cuál es la forma más elegante de hacer esto con las características de Fortran2003 gfortran-5 ?

Por supuesto, podría insertar un argumento ficticio opcional a en la función treatfunction y llamar a f con f(x) o f(x,a) dependiendo del present(a) en el cuerpo de la subrutina. Pero cambiar la subrutina no es elegante.

Answer 1

En Fortran 2008 puede pasar funciones internas como argumentos y gfortran lo admite.

 subroutine calling()

   a0 = ...
   call treatfunction(wrapper, input=myinput, output=myoutput)
 contains
   real function wrapper(x)
     real, intent(in) :: x
     wrapper = fun(x,a0)
   end function
 end subroutine

Por cierto, me mantendría alejado de lo external , es malo, use bloques de interfaz.