kappa_val.F90

Go to the documentation of this file.

!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
!
!  Function :  k a p p a _ v a l
!
!> @file
!!
!! Heat conductivity of ice:
!! Linear interpolation of tabulated values in KAPPA(.).
!!
!! @section Copyright
!!
!! Copyright 2009-2013 Ralf Greve
!!
!! @section License
!!
!! This file is part of SICOPOLIS.
!!
!! SICOPOLIS is free software: you can redistribute it and/or modify
!! it under the terms of the GNU General Public License as published by
!! the Free Software Foundation, either version 3 of the License, or
!! (at your option) any later version.
!!
!! SICOPOLIS 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 General Public License for more details.
!!
!! You should have received a copy of the GNU General Public License
!! along with SICOPOLIS.  If not, see <http://www.gnu.org/licenses/>.
!<
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

!-------------------------------------------------------------------------------
!> Heat conductivity of ice:
!! Linear interpolation of tabulated values in KAPPA(.).
!<------------------------------------------------------------------------------
function kappa_val(temp_val)

use sico_types
use sico_variables

implicit none
real(dp)             :: kappa_val
real(dp), intent(in) :: temp_val

real(dp) :: kappa_ice

!-------- Heat conductivity of pure ice --------

#if defined(FRAC_DUST)
kappa_ice = kappa(floor(temp_val)) &
           +(kappa(floor(temp_val)+1)-kappa(floor(temp_val))) &
            *(temp_val-real(floor(temp_val),dp))
#else
kappa_val = kappa(floor(temp_val)) &
           +(kappa(floor(temp_val)+1)-kappa(floor(temp_val))) &
            *(temp_val-real(floor(temp_val),dp))
#endif

!-------- If dust is present (polar caps of Mars):
!         Heat conductivity of ice-dust mixture --------

#if defined(FRAC_DUST)
kappa_val = (1.0_dp-frac_dust)*kappa_ice + frac_dust*kappa_c
#endif

end function kappa_val
!