Use format e25.16 for output of doubles to avoid warnings with ifort.