7 !! Computation of the melting temperatures.
11 !! Copyright 2009-2013 Ralf Greve
15 !! This file is part of SICOPOLIS.
17 !! SICOPOLIS is free software: you can redistribute it and/or modify
18 !! it under the terms of the GNU General Public License as published by
19 !! the Free Software Foundation, either version 3 of the License, or
20 !! (at your option) any later version.
22 !! SICOPOLIS is distributed in the hope that it will be useful,
23 !! but WITHOUT ANY WARRANTY; without even the implied warranty of
24 !! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
25 !! GNU General Public License for more details.
27 !! You should have received a copy of the GNU General Public License
28 !! along with SICOPOLIS. If not, see <http://www.gnu.org/licenses/>.
30 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
33 !> Computation of the melting temperatures.
34 !<------------------------------------------------------------------------------
41 integer(i4b) :: i, j, kc, kt
42 real(dp) :: atm1(0:kcmax), atm2(0:ktmax)
47 atm1(kc) = beta*(1.0_dp-eaz_c_quotient(kc))
51 atm2(kt) = 1.0_dp-zeta_t(kt)
60 temp_t_m(kt,j,i) = -beta*(h_c(j,i)+h_t(j,i)*atm2(kt))
64 temp_c_m(kc,j,i) = -atm1(kc)*h_c(j,i)