43 integer(i4b),
intent(in) :: i, j
44 real(dp),
intent(in) :: atr1, alb1
46 integer(i4b) :: kc, kt, kr
47 real(dp) :: ctr1, clb1
48 real(dp) :: lgs_a0(0:kcmax+ktmax+krmax+imax+jmax), &
49 lgs_a1(0:kcmax+ktmax+krmax+imax+jmax), &
50 lgs_a2(0:kcmax+ktmax+krmax+imax+jmax), &
51 lgs_x(0:kcmax+ktmax+krmax+imax+jmax), &
52 lgs_b(0:kcmax+ktmax+krmax+imax+jmax)
57 clb1 = alb1*q_geo(j,i)
71 lgs_a1(kr) = 1.0_dp + 2.0_dp*ctr1
73 lgs_b(kr) = temp_r(kr,j,i)
83 lgs_b(kr) = 2.0_dp*clb1
91 lgs_b(kr) = temp_s(j,i)
95 call
tri_sle(lgs_a0, lgs_a1, lgs_a2, lgs_x, lgs_b, krmax)
100 temp_r_neu(kr,j,i) = lgs_x(kr)
106 omega_t_neu(kt,j,i) = 0.0_dp
107 age_t_neu(kt,j,i) = 0.0_dp
113 temp_c_neu(kc,j,i) = temp_s(j,i)
114 age_c_neu(kc,j,i) = 0.0_dp