Get the current precision for real numbers
function get_precision() result(precision) !! Get the current precision for real numbers integer(default_int) :: precision precision = dp_precision end function get_precision