!
! Type_Kinds
!
! Module to hold specification kinds for variable declaration, as well as 
! associated descriptors.
!
! CREATION HISTORY:
!       Written by:     Paul van Delst, CIMSS/SSEC 12-Jun-2000
!                       paul.vandelst@ssec.wisc.edu
!

MODULE Type_Kinds

  ! ---------------------------
  ! Disable all implicit typing
  ! ---------------------------
  IMPLICIT NONE


  ! ------------
  ! Visibilities
  ! ------------
  ! Everything is private by default
  PRIVATE
  ! The integer types
  PUBLIC :: Byte   , n_Bytes_Byte 
  PUBLIC :: Short  , n_Bytes_Short
  PUBLIC :: Long   , n_Bytes_Long 
  PUBLIC :: LLong  , n_Bytes_LLong
  PUBLIC :: IP_Kind, n_Bytes_IP_Kind  ! Default integer set by IIP
  PUBLIC :: IP     , n_Bytes_IP       ! Aliases for IP_Kind
  ! The floating point types
  PUBLIC :: Single , n_Bytes_Single
  PUBLIC :: Double , n_Bytes_Double
  PUBLIC :: Quad   , n_Bytes_Quad  
  PUBLIC :: FP_Kind, n_Bytes_FP_Kind  ! Default integer set by IFP
  PUBLIC :: FP_Kind4,n_Bytes_FP_Kind4  ! Default integer set by IFP
  PUBLIC :: FP     , n_Bytes_FP       ! Aliases for FP_Kind
  PUBLIC :: FP4     , n_Bytes_FP4       ! Aliases for FP_Kind4


  ! -------------------------------------------------------------------
  ! THE DEFAULT INTEGER INDEX. Change the value of IIP for the required
  ! integer kind. The following chart details the correspondence:
  !
  !    IIP        INTEGER(IP)
  !  ==============================
  !     1        Byte 
  !     2       Short (2 bytes)
  !     3        Long (4 bytes)
  !     4       LLong (8 bytes)  **IF AVAILABLE, Long OTHERWISE**
  !
  ! -------------------------------------------------------------------
  INTEGER, PARAMETER :: IIP = 3  ! 1=Byte, 2=Short, 3=Long, 4=LLong


  ! -------------------------------------------------------------------
  ! THE DEFAULT FLOATING POINT INDEX. Change the value of IFP for the
  ! required floating point kind. The following chart details the
  ! correspondence:
  !
  !    IFP          REAL(FP)
  !  ==============================
  !     1       Single (4  bytes)
  !     2       Double (8  bytes)
  !     3       Quad   (16 bytes)  **IF AVAILABLE, Double OTHERWISE**
  !
  ! -------------------------------------------------------------------
  INTEGER, PARAMETER :: IFP = 2  ! 1=Single, 2=Double, 3=Quad
!tks
  INTEGER, PARAMETER :: IFP4 = 1  ! 1=Single, 2=Double, 3=Quad


  ! -------------------
  ! Integer definitions
  ! -------------------
  ! Integer types
  INTEGER, PARAMETER :: Byte    = SELECTED_INT_KIND(1)   ! Byte  integer
!tks Byte is reserved word in new ifort
  INTEGER, PARAMETER :: Byte2    = SELECTED_INT_KIND(1)   ! Byte  integer
!tks
  INTEGER, PARAMETER :: Short   = SELECTED_INT_KIND(4)   ! Short integer
  INTEGER, PARAMETER :: Long    = SELECTED_INT_KIND(8)   ! Long  integer
  INTEGER, PARAMETER :: LLong_t = SELECTED_INT_KIND(16)  ! LLong integer
  INTEGER, PARAMETER :: LLong   = ((( 1 + SIGN(1,LLong_t)) / 2) * LLong_t) + &
                                  ((( 1 - SIGN(1,LLong_t)) / 2) * Long   )
  ! Expected 8-bit byte sizes of the integer kinds
  INTEGER, PARAMETER :: n_Bytes_Byte  = 1
  INTEGER, PARAMETER :: n_Bytes_Short = 2
  INTEGER, PARAMETER :: n_Bytes_Long  = 4
  INTEGER, PARAMETER :: n_Bytes_LLong = 8
  ! Define arrays for default definition
  INTEGER, PARAMETER :: N_IP = 4
! INTEGER, PARAMETER, DIMENSION(N_IP) :: IP_KIND_TYPES = (/ Byte,  &
  INTEGER, PARAMETER, DIMENSION(N_IP) :: IP_KIND_TYPES = (/ Byte2,  &
                                                            Short, &
                                                            Long,  &
                                                            LLong  /) 
  INTEGER, PARAMETER, DIMENSION(N_IP) :: IP_BYTE_SIZES = (/ n_Bytes_Byte,  &
                                                            n_Bytes_Short, &
                                                            n_Bytes_Long,  &
                                                            n_Bytes_LLong  /)
  ! Default values
  INTEGER, PARAMETER :: IP_Kind        =IP_KIND_TYPES(IIP)
  INTEGER, PARAMETER :: n_Bytes_IP_Kind=IP_BYTE_SIZES(IIP)
  INTEGER, PARAMETER :: IP        =IP_Kind
  INTEGER, PARAMETER :: n_Bytes_IP=n_Bytes_IP_Kind


  ! --------------------------
  ! Floating point definitions
  ! --------------------------
  ! Floating point types
  INTEGER, PARAMETER :: Single = SELECTED_REAL_KIND(6)  ! Single precision
  INTEGER, PARAMETER :: Double = SELECTED_REAL_KIND(15) ! Double precision
  INTEGER, PARAMETER :: Quad_t = SELECTED_REAL_KIND(20) ! Quad precision
  INTEGER, PARAMETER :: Quad   = ((( 1 + SIGN(1,Quad_t)) / 2) * Quad_t) + &
                                 ((( 1 - SIGN(1,Quad_t)) / 2) * Double)
  ! Expected 8-bit byte sizes of the floating point kinds
  INTEGER, PARAMETER :: n_Bytes_Single = 4
  INTEGER, PARAMETER :: n_Bytes_Double = 8
  INTEGER, PARAMETER :: n_Bytes_Quad   = 16
  ! Define arrays for default definition
  INTEGER, PARAMETER :: N_FP = 3
  INTEGER, PARAMETER, DIMENSION(N_FP) :: FP_KIND_TYPES = (/ Single, &
                                                            Double, &
                                                            Quad    /) 
  INTEGER, PARAMETER, DIMENSION(N_FP) :: FP_BYTE_SIZES = (/ n_Bytes_Single, &
                                                            n_Bytes_Double, &
                                                            n_Bytes_Quad    /)


  ! Default values
  INTEGER, PARAMETER :: FP_Kind         = FP_KIND_TYPES(IFP)
  INTEGER, PARAMETER :: n_Bytes_FP_Kind = FP_BYTE_SIZES(IFP)
  INTEGER, PARAMETER :: FP        =FP_Kind
  INTEGER, PARAMETER :: n_Bytes_FP=n_Bytes_FP_Kind
!tks
  INTEGER, PARAMETER :: FP_Kind4         = FP_KIND_TYPES(IFP4)
  INTEGER, PARAMETER :: n_Bytes_FP_Kind4 = FP_BYTE_SIZES(IFP4)
  INTEGER, PARAMETER :: FP4       =FP_Kind4
  INTEGER, PARAMETER :: n_Bytes_FP4=n_Bytes_FP_Kind4
  
END MODULE Type_Kinds