mod_itm.f90

Go to the documentation of this file.

module mod_itm

  use itm_types

  implicit none

  interface exists_entry
    module procedure exists_entry_real, exists_entry_integer
  end interface

  real(r8), parameter :: itm_double = -9.0e40_r8
  integer(itm_i4), parameter :: itm_integer = -999999999

  contains

!----------------------------------------------------------------------------
  function exists_entry_real(v) result(v_exists)
!----------------------------------------------------------------------------
! tests whether real data base entry v exists
!----------------------------------------------------------------------------

    real(r8), intent(in) :: v
    logical :: v_exists

    if (v > (1._r8 - epsilon(v)) * itm_double) then
      v_exists = .true.
    else
      v_exists = .false.
    end if

  end function exists_entry_real

!----------------------------------------------------------------------------
  function exists_entry_integer(v) result(v_exists)
!----------------------------------------------------------------------------
! tests whether integer data base entry v exists
!----------------------------------------------------------------------------

    integer(itm_i4), intent(in) :: v
    logical :: v_exists

    if (v /= itm_integer) then
      v_exists = .true.
    else
      v_exists = .false.
    end if

  end function exists_entry_integer

end module mod_itm