Fortran Wiki
set_exponent