37 |
|
|
38 |
#ifdef ALLOW_DIAGNOSTICS |
#ifdef ALLOW_DIAGNOSTICS |
39 |
#include "SIZE.h" |
#include "SIZE.h" |
40 |
#include "diagnostics_SIZE.h" |
#include "DIAGNOSTICS_SIZE.h" |
41 |
#include "diagnostics.h" |
#include "DIAGNOSTICS.h" |
42 |
#endif |
#endif |
43 |
|
|
44 |
integer myid,im,jm,lm,bi,bj |
integer myid,im,jm,lm,bi,bj |
71 |
enddo |
enddo |
72 |
enddo |
enddo |
73 |
|
|
|
do L=1,lm |
|
|
do j=1,jm |
|
|
do i=1,im |
|
|
dz(i,j,L) = cpog * (pekap(i,j,L+1)-pekap(i,j,L)) * t(i,j,L) * |
|
|
. (1.+virtcon*s(i,j,L)) |
|
|
enddo |
|
|
enddo |
|
|
enddo |
|
|
|
|
74 |
do j=1,jm |
do j=1,jm |
75 |
do i=1,im |
do i=1,im |
76 |
z(i,j,lm) = zsurf(i,j) + 0.5 * dz(i,j,lm) |
z(i,j,lm) = zsurf(i,j) + 0.5 * dz(i,j,lm) |
88 |
do L=1,lm |
do L=1,lm |
89 |
do j=1,jm |
do j=1,jm |
90 |
do i=1,im |
do i=1,im |
91 |
rf(i,j,L) = dampcoef*(1+tanh((z(i,j,L)-80000.)/5000.))/86400. |
rf(i,j,L) = dampcoef*(1+tanh((z(i,j,L)-50000.)/5000.))/86400. |
92 |
rfu(i,j,L) = - rf(i,j,L) * u(i,j,L) |
rfu(i,j,L) = - rf(i,j,L) * u(i,j,L) |
93 |
rfv(i,j,L) = - rf(i,j,L) * v(i,j,L) |
rfv(i,j,L) = - rf(i,j,L) * v(i,j,L) |
94 |
rft(i,j,L) = -(u(i,j,L)*rfu(i,j,L) + v(i,j,L)*rfv(i,j,L) )*cpinv |
rft(i,j,L) = -(u(i,j,L)*rfu(i,j,L) + v(i,j,L)*rfv(i,j,L) )*cpinv |