From e2ea3cdff6729c7ad3e467e7812bf28c15f3a797 Mon Sep 17 00:00:00 2001
From: Martin	Diehl <m.diehl@mpie.de>
Date: Tue, 2 Oct 2012 16:53:03 +0000
Subject: [PATCH] again some changes regarding creation of NaN, ieee-conforming
 way is not in the complilers yet.

---
 code/prec.f90 | 25 +++++++++++--------------
 1 file changed, 11 insertions(+), 14 deletions(-)

diff --git a/code/prec.f90 b/code/prec.f90
index abbaecfe0..8101df841 100644
--- a/code/prec.f90
+++ b/code/prec.f90
@@ -25,13 +25,7 @@
 !> @author Martin Diehl, Max-Planck-Institut für Eisenforschung GmbH
 !> @brief setting precision for real and int type depending on makros "FLOAT" and "INT"
 !--------------------------------------------------------------------------------------------------
-#ifdef __INTEL_COMPILER
-#if __INTEL_COMPILER<1200
-#define LEGACY_COMPILER
-#endif
-#endif
 
-!--------------------------------------------------------------------------------------------------
 module prec
 
  implicit none
@@ -42,16 +36,18 @@ module prec
  SPECTRAL SOLVER DOES NOT SUPPORT SINGLE PRECISION, STOPING COMPILATION
 #endif
  integer,     parameter, public :: pReal = 4                                                        !< floating point single precition (was selected_real_kind(6,37), number with 6 significant digits, up to 1e+-37)
-#ifdef LEGACY_COMPILER
+#ifdef __INTEL_COMPILER
  real(pReal), parameter, public :: DAMASK_NaN = Z'7F800001'                                         !< quiet NaN for single precision (from http://www.hpc.unimelb.edu.au/doc/f90lrm/dfum_035.html, copy can be found in documentation/Code/Fortran)
-#else
+#endif
+#ifdef __GFORTRAN__
  real(pReal), parameter, public :: DAMASK_NaN = real(Z'7F800001', pReal)                            !< quiet NaN for single precision (from http://www.hpc.unimelb.edu.au/doc/f90lrm/dfum_035.html, copy can be found in documentation/Code/Fortran)
 #endif
 #elif (FLOAT==8)
  integer,     parameter, public :: pReal = 8                                                        !< floating point double precision (was selected_real_kind(15,300), number with 15 significant digits, up to 1e+-300)
-#ifdef LEGACY_COMPILER
+#ifdef __INTEL_COMPILER
  real(pReal), parameter, public :: DAMASK_NaN = Z'7FF8000000000000'                                 !< quiet NaN for double precision (from http://www.hpc.unimelb.edu.au/doc/f90lrm/dfum_035.html, copy can be found in documentation/Code/Fortran)
-#else
+#endif
+#ifdef __GFORTRAN__
  real(pReal), parameter, public :: DAMASK_NaN = real(Z'7FF8000000000000', pReal)                    !< quiet NaN for double precision (from http://www.hpc.unimelb.edu.au/doc/f90lrm/dfum_035.html, copy can be found in documentation/Code/Fortran)
 #endif
 #else
@@ -88,9 +84,6 @@ subroutine prec_init
  implicit none
 
 !$OMP CRITICAL (write2out)
-#ifndef LEGACY_COMPILER
- open (6, encoding='UTF-8')  
-#endif
 
  write(6,*)
  write(6,*) '<<<+-  prec init  -+>>>'
@@ -101,8 +94,12 @@ subroutine prec_init
  write(6,'(a,i3)')    ' Bytes for pLongInt: ',pLongInt
  write(6,'(a,e10.3)') ' NaN:           ',     DAMASK_NaN
  write(6,'(a,l3)')    ' NaN /= NaN:         ',DAMASK_NaN/=DAMASK_NaN
- write(6,*)
  if (DAMASK_NaN == DAMASK_NaN) call quit(9000)
+#ifdef Spectral
+ open(6, encoding='UTF-8')                                                                          ! modern fortran compilers (gfortran >4.4, ifort >11 support it)
+ write(6,*) 'using UTF-8 coded output'
+#endif
+ write(6,*)
 !$OMP END CRITICAL (write2out)
 
 end subroutine prec_init