183 |
|
|
184 |
c----------------------------------------------------------------------- |
c----------------------------------------------------------------------- |
185 |
c combine interior and boundary layer coefficients and nonlocal term |
c combine interior and boundary layer coefficients and nonlocal term |
186 |
|
c !!!NOTE!!! In shallow (2-level) regions and for shallow mixed layers |
187 |
|
c (< 1 level), diffusivity blmc can become negative. The max's below |
188 |
|
c are a hack until this problem is properly diagnosed and fixed. |
189 |
c----------------------------------------------------------------------- |
c----------------------------------------------------------------------- |
|
|
|
190 |
do k = 1, Nr |
do k = 1, Nr |
191 |
do i = 1, imt |
do i = 1, imt |
192 |
if (k .lt. kbl(i)) then |
if (k .lt. kbl(i)) then |
193 |
do md = 1, mdiff |
diffus(i,k,1) = max ( blmc(i,k,1), viscAr ) |
194 |
diffus(i,k,md) = blmc(i,k,md) |
diffus(i,k,2) = max ( blmc(i,k,2), diffKrS ) |
195 |
end do |
diffus(i,k,3) = max ( blmc(i,k,3), diffKrT ) |
196 |
else |
else |
197 |
ghat(i,k) = 0. |
ghat(i,k) = 0. |
198 |
endif |
endif |