221 |
enddo |
enddo |
222 |
enddo |
enddo |
223 |
endif |
endif |
224 |
|
if (iudiag2.ne.0) then |
225 |
|
do j=jm1,jm2 |
226 |
|
do i=im1,im2 |
227 |
|
qdiag(i,j,iudiag2+l-1,bi,bj) = qdiag(i,j,iudiag2+l-1,bi,bj) + |
228 |
|
. lwdt(i,j,l,bi,bj) * pk(i,j,l,bi,bj)*pinv(i,j)*86400 |
229 |
|
enddo |
230 |
|
enddo |
231 |
|
endif |
232 |
|
if (iudiag3.ne.0) then |
233 |
|
do j=jm1,jm2 |
234 |
|
do i=im1,im2 |
235 |
|
qdiag(i,j,iudiag3+l-1,bi,bj) = qdiag(i,j,iudiag3+l-1,bi,bj) + |
236 |
|
. dlwdtg (i,j,L,bi,bj)*(tgz(i,j,bi,bj)-tg0(i,j,bi,bj)) |
237 |
|
. * pk(i,j,l,bi,bj)*pinv(i,j)*86400 |
238 |
|
enddo |
239 |
|
enddo |
240 |
|
endif |
241 |
|
|
242 |
c Longwave Heating Clear-Sky (deg/day) |
c Longwave Heating Clear-Sky (deg/day) |
243 |
c ------------------------------------ |
c ------------------------------------ |
346 |
ntmpu = ntmpu + 1 |
ntmpu = ntmpu + 1 |
347 |
ntke = ntke + 1 |
ntke = ntke + 1 |
348 |
nsphu = nsphu + 1 |
nsphu = nsphu + 1 |
349 |
|
|
350 |
|
nudiag3 = nudiag3 + 1 |
351 |
|
nudiag2 = nudiag2 + 1 |
352 |
endif |
endif |
353 |
|
|
354 |
C ********************************************************************** |
C ********************************************************************** |