Skip to content

Commit

Permalink
Fortran-Lib: remove NR references
Browse files Browse the repository at this point in the history
  • Loading branch information
@MuellerSeb
MuellerSeb committed Jan 7, 2020
1 parent bc054bb commit f20a1e7
Show file tree
Hide file tree
Showing 6 changed files with 892 additions and 3,006 deletions.
128 changes: 85 additions & 43 deletions pre-proc/process_soilgrid/mo_kind.f90
View file
Original file line number Diff line number Diff line change
@@ -1,51 +1,93 @@
MODULE mo_kind
!> \file mo_kind.f90

!> \brief Define number representations

!> \details This module declares the desired ranges and precisions of the number representations,
!> such as single precision or double precision, 32-bit or 46-bit integer, etc.
!> It confirms mostly with the nrtype module of Numerical Recipes in f90.

!> \authors Juliane Mai, Matthias Cuntz, Nov 2011
!> \date 2011-2014
!> \copyright GNU Lesser General Public License http://www.gnu.org/licenses/

! Number model from which the SELECTED_REAL_KIND are requested:
! 4 byte REAL 8 byte REAL
! IEEE: precision = 6 precision = 15
! exponent = 37 exponent = 307
! CRAY: - precision = 13
! exponent = 2465

! Written Juliane Mai, Matthias Cuntz, Nov 2011
! Modified Matthias Cuntz, Nov 2011 - private/public
! - documentation
! - removed tab characters
! Matthias Cuntz, May 2014 - iso_fortran_env and iso_c_binding
! S. Mueller, Dec 2019 - remove NR specific types

! This module declares the desired ranges and precisions of the number representations,
! such as single precision or double precision, 32-bit or 46-bit integer, etc.
! It confirms mostly with the nrtype module of Numerical Recipes in f90.
! License
! -------
! This file is part of the UFZ Fortran library.

! Number model from which the SELECTED_REAL_KIND are requested:
! 4 byte REAL 8 byte REAL
! IEEE: precision = 6 precision = 15
! exponent = 37 exponent = 307
! CRAY: - precision = 13
! exponent = 2465
! The UFZ Fortran library is free software: you can redistribute it and/or modify
! it under the terms of the GNU Lesser General Public License as published by
! the Free Software Foundation, either version 3 of the License, or
! (at your option) any later version.

! The UFZ Fortran library is distributed in the hope that it will be useful,
! but WITHOUT ANY WARRANTY; without even the implied warranty of
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
! GNU Lesser General Public License for more details.

! You should have received a copy of the GNU Lesser General Public License
! along with the UFZ Fortran library (cf. gpl.txt and lgpl.txt).
! If not, see <http://www.gnu.org/licenses/>.

! Copyright 2011-2014 Matthias Cuntz, Juliane Mai

MODULE mo_kind

! Written Juliane Mai, Nov 2011
! Modified Matthias Cuntz, Nov 2011 - private/public
! - documentation
! - removed tab characters
! Does not work with compilers intel v11 and sun v12.2
! use, intrinsic :: iso_fortran_env, only: &
! int8!, int16, int32, int64, real32, real64
use, intrinsic :: iso_c_binding, only: &
c_short, c_int, c_long_long, c_float, c_double, c_float_complex, c_double_complex, c_bool

IMPLICIT NONE

PRIVATE

PUBLIC :: i1, i2, i4, i8, sp, dp, spc, dpc, lgt, sprs2_sp, sprs2_dp

INTEGER, PARAMETER :: i8 = SELECTED_INT_KIND(18) ! 8 Byte Integer
INTEGER, PARAMETER :: i4 = SELECTED_INT_KIND(9) ! 4 Byte Integer
INTEGER, PARAMETER :: i2 = SELECTED_INT_KIND(4) ! 2 Byte Integer
INTEGER, PARAMETER :: i1 = SELECTED_INT_KIND(2) ! 1 Byte Integer
INTEGER, PARAMETER :: sp = SELECTED_REAL_KIND(6,37) ! Single Precision Real
INTEGER, PARAMETER :: dp = SELECTED_REAL_KIND(15,307) ! Double Precision Real
! Note: Intrinsic KIND is an extension to the Fortran 90 standard.
! Matthias C thinks that spc=sp and dpc=dp would be enough.
INTEGER, PARAMETER :: spc = KIND((1.0_sp,1.0_sp)) ! Single Precision Complex
INTEGER, PARAMETER :: dpc = KIND((1.0_dp,1.0_dp)) ! Double Precision Complex
INTEGER, PARAMETER :: lgt = KIND(.true.) ! Logical

! Numerical recipes types for sparse arrays
TYPE sprs2_sp
INTEGER(I4) :: n, len
REAL(SP), DIMENSION(:), POINTER :: val
INTEGER(I4), DIMENSION(:), POINTER :: irow
INTEGER(I4), DIMENSION(:), POINTER :: jcol
END TYPE sprs2_sp
TYPE sprs2_dp
INTEGER(I4) :: n, len
REAL(DP), DIMENSION(:), POINTER :: val
INTEGER(I4), DIMENSION(:), POINTER :: irow
INTEGER(I4), DIMENSION(:), POINTER :: jcol
END TYPE sprs2_dp
!> 1 Byte Integer Kind
INTEGER, PARAMETER :: i1 = SELECTED_INT_KIND(2)
! INTEGER, PARAMETER :: i1 = int8 ! c_word does not exist; should be c_bool, probably; but see above
!> 2 Byte Integer Kind
! INTEGER, PARAMETER :: i2 = SELECTED_INT_KIND(4)
! INTEGER, PARAMETER :: i2 = int16
INTEGER, PARAMETER :: i2 = c_short
!> 4 Byte Integer Kind
! INTEGER, PARAMETER :: i4 = SELECTED_INT_KIND(9)
! INTEGER, PARAMETER :: i4 = int32
INTEGER, PARAMETER :: i4 = c_int
!> 8 Byte Integer Kind
! INTEGER, PARAMETER :: i8 = SELECTED_INT_KIND(18)
! INTEGER, PARAMETER :: i8 = int64
INTEGER, PARAMETER :: i8 = c_long_long
!> Single Precision Real Kind
! INTEGER, PARAMETER :: sp = SELECTED_REAL_KIND(6,37)
! INTEGER, PARAMETER :: sp = real32
INTEGER, PARAMETER :: sp = c_float
!> Double Precision Real Kind
! INTEGER, PARAMETER :: dp = SELECTED_REAL_KIND(15,307)
! INTEGER, PARAMETER :: dp = real64
INTEGER, PARAMETER :: dp = c_double
!> Single Precision Complex Kind
! INTEGER, PARAMETER :: spc = KIND((1.0_sp,1.0_sp))
! INTEGER, PARAMETER :: spc = sp
INTEGER, PARAMETER :: spc = c_float_complex
!> Double Precision Complex Kind
! INTEGER, PARAMETER :: dpc = KIND((1.0_dp,1.0_dp))
! INTEGER, PARAMETER :: dpc = dp
INTEGER, PARAMETER :: dpc = c_double_complex
!> Logical Kind
INTEGER, PARAMETER :: lgt = KIND(.true.)
!INTEGER, PARAMETER :: lgt = c_bool
! Types have to be in a public section for doxygen

END MODULE mo_kind
Loading

0 comments on commit f20a1e7

Please sign in to comment.