22 |
_RL vphy(im2,jm2,Nrphys,Nbi,Nbj) |
_RL vphy(im2,jm2,Nrphys,Nbi,Nbj) |
23 |
_RL thphy(im2,jm2,Nrphys,Nbi,Nbj) |
_RL thphy(im2,jm2,Nrphys,Nbi,Nbj) |
24 |
_RL sphy(im2,jm2,Nrphys,Nbi,Nbj) |
_RL sphy(im2,jm2,Nrphys,Nbi,Nbj) |
25 |
_RL qq(im2,jm2,Nrphys),pk(im2,jm2,Nrphys,Nbi,Nbj) |
_RL qq(im2,jm2,Nrphys,Nbi,Nbj),pk(im2,jm2,Nrphys,Nbi,Nbj) |
26 |
_RL dp(im2,jm2,Nrphys,Nbi,Nbj) |
_RL dp(im2,jm2,Nrphys,Nbi,Nbj) |
27 |
_RL radswt(im2,jm2,Nbi,Nbj),radswg(im2,jm2,Nbi,Nbj) |
_RL radswt(im2,jm2,Nbi,Nbj),radswg(im2,jm2,Nbi,Nbj) |
28 |
_RL swgclr(im2,jm2,Nbi,Nbj),osr(im2,jm2,Nbi,Nbj) |
_RL swgclr(im2,jm2,Nbi,Nbj),osr(im2,jm2,Nbi,Nbj) |
297 |
if( itke.ne.0 ) then |
if( itke.ne.0 ) then |
298 |
do j=jm1,jm2 |
do j=jm1,jm2 |
299 |
do i=im1,im2 |
do i=im1,im2 |
300 |
qdiag(i,j,itke+L-1,bi,bj) = qdiag(i,j,itke+L-1,bi,bj) + qq(i,j,L) |
qdiag(i,j,itke+L-1,bi,bj) = qdiag(i,j,itke+L-1,bi,bj) |
301 |
|
& + qq(i,j,L,bi,bj) |
302 |
enddo |
enddo |
303 |
enddo |
enddo |
304 |
endif |
endif |