replace custom is_finite, is_nan, is_special by ieee_arithmetic
http://fortranwiki.org/fortran/show/ieee_arithmetic is the current standard and all compilers do by now have it, so much of the content in mo_utils.f90
can either go directly reference ieee_arithmetic without compiler-dependent pre-processors.