12 |
|
|
13 |
integer myid,im1,im2,jm1,jm2,Nrphys,Nbi,Nbj,bi,bj,ntracer |
integer myid,im1,im2,jm1,jm2,Nrphys,Nbi,Nbj,bi,bj,ntracer |
14 |
_RL p(im2,jm2,Nbi,Nbj) |
_RL p(im2,jm2,Nbi,Nbj) |
15 |
_RL uphy(im2,jm2,Nrphys,Nbi,Nbj) |
_RL uphy(im2,jm2,Nrphys) |
16 |
_RL vphy(im2,jm2,Nrphys,Nbi,Nbj) |
_RL vphy(im2,jm2,Nrphys) |
17 |
_RL thphy(im2,jm2,Nrphys,Nbi,Nbj) |
_RL thphy(im2,jm2,Nrphys) |
18 |
_RL sphy(im2,jm2,Nrphys,Nbi,Nbj) |
_RL sphy(im2,jm2,Nrphys) |
19 |
_RL qq(im2,jm2,Nrphys,Nbi,Nbj),pk(im2,jm2,Nrphys,Nbi,Nbj) |
_RL qq(im2,jm2,Nrphys,Nbi,Nbj),pk(im2,jm2,Nrphys,Nbi,Nbj) |
20 |
_RL dp(im2,jm2,Nrphys,Nbi,Nbj) |
_RL dp(im2,jm2,Nrphys,Nbi,Nbj) |
21 |
_RL radswt(im2,jm2,Nbi,Nbj),radswg(im2,jm2,Nbi,Nbj) |
_RL radswt(im2,jm2,Nbi,Nbj),radswg(im2,jm2,Nbi,Nbj) |
275 |
if(diagnostics_is_on('UWND ',myid) ) then |
if(diagnostics_is_on('UWND ',myid) ) then |
276 |
do j=jm1,jm2 |
do j=jm1,jm2 |
277 |
do i=im1,im2 |
do i=im1,im2 |
278 |
tmpdiag(i,j) = uphy(i,j,L,bi,bj) |
tmpdiag(i,j) = uphy(i,j,L) |
279 |
enddo |
enddo |
280 |
enddo |
enddo |
281 |
call diagnostics_fill(tmpdiag,'UWND ',L,1,3,bi,bj,myid) |
call diagnostics_fill(tmpdiag,'UWND ',L,1,3,bi,bj,myid) |
286 |
if(diagnostics_is_on('VWND ',myid) ) then |
if(diagnostics_is_on('VWND ',myid) ) then |
287 |
do j=jm1,jm2 |
do j=jm1,jm2 |
288 |
do i=im1,im2 |
do i=im1,im2 |
289 |
tmpdiag(i,j) = vphy(i,j,L,bi,bj) |
tmpdiag(i,j) = vphy(i,j,L) |
290 |
enddo |
enddo |
291 |
enddo |
enddo |
292 |
call diagnostics_fill(tmpdiag,'VWND ',L,1,3,bi,bj,myid) |
call diagnostics_fill(tmpdiag,'VWND ',L,1,3,bi,bj,myid) |
297 |
if(diagnostics_is_on('TMPU ',myid) ) then |
if(diagnostics_is_on('TMPU ',myid) ) then |
298 |
do j=jm1,jm2 |
do j=jm1,jm2 |
299 |
do i=im1,im2 |
do i=im1,im2 |
300 |
tmpdiag(i,j) = thphy(i,j,L,bi,bj)*pk(i,j,L,bi,bj) |
tmpdiag(i,j) = thphy(i,j,L)*pk(i,j,L,bi,bj) |
301 |
enddo |
enddo |
302 |
enddo |
enddo |
303 |
call diagnostics_fill(tmpdiag,'TMPU ',L,1,3,bi,bj,myid) |
call diagnostics_fill(tmpdiag,'TMPU ',L,1,3,bi,bj,myid) |
319 |
if(diagnostics_is_on('SPHU ',myid) ) then |
if(diagnostics_is_on('SPHU ',myid) ) then |
320 |
do j=jm1,jm2 |
do j=jm1,jm2 |
321 |
do i=im1,im2 |
do i=im1,im2 |
322 |
tmpdiag(i,j) = sphy(i,j,L,bi,bj) * 1000. |
tmpdiag(i,j) = sphy(i,j,L) * 1000. |
323 |
enddo |
enddo |
324 |
enddo |
enddo |
325 |
call diagnostics_fill(tmpdiag,'SPHU ',L,1,3,bi,bj,myid) |
call diagnostics_fill(tmpdiag,'SPHU ',L,1,3,bi,bj,myid) |
440 |
do j=jm1,jm2 |
do j=jm1,jm2 |
441 |
do i=im1,im2 |
do i=im1,im2 |
442 |
qbar(i,j) = qbar(i,j) + |
qbar(i,j) = qbar(i,j) + |
443 |
. sphy(i,j,L,bi,bj)*dp(i,j,L,bi,bj) |
. sphy(i,j,L)*dp(i,j,L,bi,bj) |
444 |
enddo |
enddo |
445 |
enddo |
enddo |
446 |
enddo |
enddo |