module my_stuff implicit none double precision rawSE end module my_stuff ********* program test_rawSE use my_stuff implicit none write (*, *) rawSE (2.0d0) end program test_rawSE ******** double precision FUNCTION rawSE(x) double precision x rawSE = 2.0d0 * x return end FUNCTION