18 |
|
|
19 |
integer myThid,im1,im2,jm1,jm2,Nrphys,Nbi,Nbj,bi,bj,ntracer |
integer myThid,im1,im2,jm1,jm2,Nrphys,Nbi,Nbj,bi,bj,ntracer |
20 |
_RL p(im2,jm2,Nbi,Nbj) |
_RL p(im2,jm2,Nbi,Nbj) |
21 |
_RL uphy(im2,jm2,Nrphys,Nbi,Nbj),vphy(im2,jm2,Nrphys,Nbi,Nbj) |
_RL uphy(im2,jm2,Nrphys,Nbi,Nbj) |
22 |
_RL thphy(im2,jm2,Nrphys,Nbi,Nbj),sphy(im2,jm2,Nrphys,Nbi,Nbj) |
_RL vphy(im2,jm2,Nrphys,Nbi,Nbj) |
23 |
|
_RL thphy(im2,jm2,Nrphys,Nbi,Nbj) |
24 |
|
_RL sphy(im2,jm2,Nrphys,Nbi,Nbj) |
25 |
_RL qq(im2,jm2,Nrphys),pk(im2,jm2,Nrphys,Nbi,Nbj) |
_RL qq(im2,jm2,Nrphys),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) |
30 |
_RL dst4(im2,jm2,Nbi,Nbj),tgz(im2,jm2,Nbi,Nbj) |
_RL dst4(im2,jm2,Nbi,Nbj),tgz(im2,jm2,Nbi,Nbj) |
31 |
_RL tg0(im2,jm2,Nbi,Nbj),radlwg(im2,jm2,Nbi,Nbj) |
_RL tg0(im2,jm2,Nbi,Nbj),radlwg(im2,jm2,Nbi,Nbj) |
32 |
_RL lwgclr(im2,jm2,Nbi,Nbj) |
_RL lwgclr(im2,jm2,Nbi,Nbj) |
33 |
_RL turbu(im2,jm2,Nrphys,Nbi,Nbj),turbv(im2,jm2,Nrphys,Nbi,Nbj) |
_RL turbu(im2,jm2,Nrphys,Nbi,Nbj) |
34 |
|
_RL turbv(im2,jm2,Nrphys,Nbi,Nbj) |
35 |
_RL turbt(im2,jm2,Nrphys,Nbi,Nbj) |
_RL turbt(im2,jm2,Nrphys,Nbi,Nbj) |
36 |
_RL turbq(im2,jm2,Nrphys,ntracer,Nbi,Nbj) |
_RL turbq(im2,jm2,Nrphys,ntracer,Nbi,Nbj) |
37 |
_RL moistu(im2,jm2,Nrphys,Nbi,Nbj),moistv(im2,jm2,Nrphys,Nbi,Nbj) |
_RL moistu(im2,jm2,Nrphys,Nbi,Nbj) |
38 |
|
_RL moistv(im2,jm2,Nrphys,Nbi,Nbj) |
39 |
_RL moistt(im2,jm2,Nrphys,Nbi,Nbj) |
_RL moistt(im2,jm2,Nrphys,Nbi,Nbj) |
40 |
_RL moistq(im2,jm2,Nrphys,ntracer,Nbi,Nbj) |
_RL moistq(im2,jm2,Nrphys,ntracer,Nbi,Nbj) |
41 |
_RL lwdt(im2,jm2,Nrphys,Nbi,Nbj),swdt(im2,jm2,Nrphys,Nbi,Nbj) |
_RL lwdt(im2,jm2,Nrphys,Nbi,Nbj) |
42 |
|
_RL swdt(im2,jm2,Nrphys,Nbi,Nbj) |
43 |
_RL lwdtclr(im2,jm2,Nrphys,Nbi,Nbj) |
_RL lwdtclr(im2,jm2,Nrphys,Nbi,Nbj) |
44 |
_RL swdtclr(im2,jm2,Nrphys,Nbi,Nbj) |
_RL swdtclr(im2,jm2,Nrphys,Nbi,Nbj) |
45 |
_RL dlwdtg(im2,jm2,Nrphys,Nbi,Nbj) |
_RL dlwdtg(im2,jm2,Nrphys,Nbi,Nbj) |
49 |
|
|
50 |
C ********************************************************************** |
C ********************************************************************** |
51 |
|
|
52 |
|
#ifdef ALLOW_DIAGNOSTICS |
53 |
do j=jm1,jm2 |
do j=jm1,jm2 |
54 |
do i=im1,im2 |
do i=im1,im2 |
55 |
pinv(i,j) = 1.0 / p(i,j,bi,bj) |
pinv(i,j) = 1.0 / p(i,j,bi,bj) |
431 |
nvdtradlw = nvdtradlw + 1 |
nvdtradlw = nvdtradlw + 1 |
432 |
nvdtradsw = nvdtradsw + 1 |
nvdtradsw = nvdtradsw + 1 |
433 |
|
|
434 |
|
#endif |
435 |
return |
return |
436 |
end |
end |