Fortran - subprogramas en fortran

saben que tengo un problemas con funciones o subprogramas, alguien me puede mostar un codigo con una funcion que sume 2 numeros, por fa. no logro que me corran las funciones.

Espero que este programita te sirva

Program Sumar2n
implicit none

real*8 :: n1, n2, resultado

n1=10.0
n2=10.0

resultado = Sumar(n1,n2)

write(*,*) resultado

contains

function sumar(A,B)
real*8 :: sumar, A, B
sumar=A+B
end function sumar

end Program Sumar2n

Si es un procedimiento tienes que usar subroutine en ves de function y este no te regresa un resultado, lo puedes usar como por ejemplo para llenar una matriz o algo asi y lo mandas llamar como

call procedimiento(parametros)