7 !! Computation of temperature for an ice-free column.
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 temperature for an ice-free column.
34 !<------------------------------------------------------------------------------
41 integer(i4b) :: i, j, kc, kt, kr
42 real(dp) :: atr1, alb1, ctr1, clb1
43 real(dp) :: lgs_a0(0:kcmax+ktmax+krmax+imax+jmax), &
44 lgs_a1(0:kcmax+ktmax+krmax+imax+jmax), &
45 lgs_a2(0:kcmax+ktmax+krmax+imax+jmax), &
46 lgs_x(0:kcmax+ktmax+krmax+imax+jmax), &
47 lgs_b(0:kcmax+ktmax+krmax+imax+jmax)
52 clb1 = alb1*q_geo(j,i)
66 lgs_a1(kr) = 1.0_dp + 2.0_dp*ctr1
68 lgs_b(kr) = temp_r(kr,j,i)
78 lgs_b(kr) = 2.0_dp*clb1
86 lgs_b(kr) = temp_s(j,i)
90 call
tri_sle(lgs_a0, lgs_a1, lgs_a2, lgs_x, lgs_b, krmax)
95 temp_r_neu(kr,j,i) = lgs_x(kr)
101 omega_t_neu(kt,j,i) = 0.0_dp
102 age_t_neu(kt,j,i) = 0.0_dp
108 temp_c_neu(kc,j,i) = temp_s(j,i)
109 age_c_neu(kc,j,i) = 0.0_dp