pure function funSqrt2(b) result(ret) real(dp), intent(in) :: b real(dp) :: ret unused_var(b) ret = Root2 end function funSqrt2