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)