| 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 |
Parent Directory
| 
Revision Log
| 
Revision Graph
| 
Patch