pure function funBplus1(b) result(ret) real(dp), intent(in) :: b real(dp) :: ret ret = b + 1.0_dp end function funBplus1