7 !! Rate factor for cold ice:
8 !! Linear interpolation of tabulated values in RF(.).
12 !! Copyright 2009-2013 Ralf Greve
16 !! This file is part of SICOPOLIS.
18 !! SICOPOLIS is free software: you can redistribute it and/or modify
19 !! it under the terms of the GNU General Public License as published by
20 !! the Free Software Foundation, either version 3 of the License, or
21 !! (at your option) any later version.
23 !! SICOPOLIS is distributed in the hope that it will be useful,
24 !! but WITHOUT ANY WARRANTY; without even the implied warranty of
25 !! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
26 !! GNU General Public License for more details.
28 !! You should have received a copy of the GNU General Public License
29 !! along with SICOPOLIS. If not, see <http://www.gnu.org/licenses/>.
31 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
34 !> Rate factor for cold ice:
35 !! Linear interpolation of tabulated values in RF(.).
36 !<------------------------------------------------------------------------------
44 real(dp),
intent(in) :: temp_val, temp_m_val
46 real(dp) :: temp_h_val
48 temp_h_val = temp_val-temp_m_val
50 ratefac = rf(floor(temp_h_val)) &
51 +(rf(floor(temp_h_val)+1)-rf(floor(temp_h_val))) &
52 *(temp_h_val-
real(floor(temp_h_val),dp))