10 subroutine tridi1(l,n,cl,cm,cu,r1,au,a1)
12 use machine ,
only : kind_phys
14 integer,
parameter :: one = 1.0_kind_phys
16 real(kind=kind_phys) fk
18 real(kind=kind_phys) cl(l,2:n),cm(l,n),cu(l,n-1),r1(l,n), &
28 fk = one / (cm(i,k)-cl(i,k)*au(i,k-1))
30 a1(i,k) = fk*(r1(i,k)-cl(i,k)*a1(i,k-1))
34 fk = one / (cm(i,n)-cl(i,n)*au(i,n-1))
35 a1(i,n) = fk*(r1(i,n)-cl(i,n)*a1(i,n-1))
39 a1(i,k) = a1(i,k)-au(i,k)*a1(i,k+1)
49 subroutine tridi2(l,n,cl,cm,cu,r1,r2,au,a1,a2)
51 use machine ,
only : kind_phys
53 integer,
parameter :: one = 1.0_kind_phys
55 real(kind=kind_phys) fk
57 real(kind=kind_phys) cl(l,2:n),cm(l,n),cu(l,n-1),r1(l,n),r2(l,n), &
58 & au(l,n-1),a1(l,n),a2(l,n)
68 fk = one / (cm(i,k)-cl(i,k)*au(i,k-1))
70 a1(i,k) = fk*(r1(i,k)-cl(i,k)*a1(i,k-1))
71 a2(i,k) = fk*(r2(i,k)-cl(i,k)*a2(i,k-1))
75 fk = one / (cm(i,n)-cl(i,n)*au(i,n-1))
76 a1(i,n) = fk*(r1(i,n)-cl(i,n)*a1(i,n-1))
77 a2(i,n) = fk*(r2(i,n)-cl(i,n)*a2(i,n-1))
81 a1(i,k) = a1(i,k)-au(i,k)*a1(i,k+1)
82 a2(i,k) = a2(i,k)-au(i,k)*a2(i,k+1)
94 subroutine tridin(l,n,nt,cl,cm,cu,r1,r2,au,a1,a2)
96 use machine ,
only : kind_phys
98 integer,
parameter :: one = 1.0_kind_phys
99 integer is,k,kk,n,nt,l,i
100 real(kind=kind_phys) fk(l)
102 real(kind=kind_phys) cl(l,2:n), cm(l,n), cu(l,n-1), &
103 & r1(l,n), r2(l,n*nt), &
104 & au(l,n-1), a1(l,n), a2(l,n*nt), &
108 fk(i) = one / cm(i,1)
109 au(i,1) = fk(i)*cu(i,1)
110 a1(i,1) = fk(i)*r1(i,1)
115 a2(i,1+is) = fk(i) * r2(i,1+is)
120 fkk(i,k) = one / (cm(i,k)-cl(i,k)*au(i,k-1))
121 au(i,k) = fkk(i,k)*cu(i,k)
122 a1(i,k) = fkk(i,k)*(r1(i,k)-cl(i,k)*a1(i,k-1))
129 a2(i,k+is) = fkk(i,k)*(r2(i,k+is)-cl(i,k)*a2(i,k+is-1))
134 fk(i) = one / (cm(i,n)-cl(i,n)*au(i,n-1))
135 a1(i,n) = fk(i)*(r1(i,n)-cl(i,n)*a1(i,n-1))
140 a2(i,n+is) = fk(i)*(r2(i,n+is)-cl(i,n)*a2(i,n+is-1))
145 a1(i,k) = a1(i,k) - au(i,k)*a1(i,k+1)
152 a2(i,k+is) = a2(i,k+is) - au(i,k)*a2(i,k+is+1)
163 subroutine tridit(l,n,nt,cl,cm,cu,rt,au,at)
166 use machine ,
only : kind_phys
168 integer,
parameter :: one = 1.0_kind_phys
169 integer is,k,kk,n,nt,l,i
170 real(kind=kind_phys) fk(l)
172 real(kind=kind_phys) cl(l,2:n), cm(l,n), cu(l,n-1), &
174 & au(l,n-1), at(l,n*nt), &
178 fk(i) = one / cm(i,1)
179 au(i,1) = fk(i)*cu(i,1)
184 at(i,1+is) = fk(i) * rt(i,1+is)
189 fkk(i,k) = one / (cm(i,k)-cl(i,k)*au(i,k-1))
190 au(i,k) = fkk(i,k)*cu(i,k)
197 at(i,k+is) = fkk(i,k)*(rt(i,k+is)-cl(i,k)*at(i,k+is-1))
202 fk(i) = one / (cm(i,n)-cl(i,n)*au(i,n-1))
207 at(i,n+is) = fk(i)*(rt(i,n+is)-cl(i,n)*at(i,n+is-1))
214 at(i,k+is) = at(i,k+is) - au(i,k)*at(i,k+is+1)
subroutine tridi1(l, n, cl, cm, cu, r1, au, a1)
Routine to solve the tridiagonal system to calculate temperature and moisture at .
subroutine tridit(l, n, nt, cl, cm, cu, rt, au, at)
This subroutine solves tridiagonal problem for TKE.
subroutine tridin(l, n, nt, cl, cm, cu, r1, r2, au, a1, a2)
Routine to solve the tridiagonal system to calculate u- and v-momentum at .
subroutine tridi2(l, n, cl, cm, cu, r1, r2, au, a1, a2)
This subroutine ..