125 |
c ---------------------------------------------- |
c ---------------------------------------------- |
126 |
do j=1,jm |
do j=1,jm |
127 |
do i=1,im |
do i=1,im |
128 |
phis_std(i,j) = min( 400.0 _d 0, sqrt( max(0.0 _d 0,phis_var(i,j)) )/grav ) |
phis_std(i,j) = min( 400.0 _d 0, sqrt( max(0.0 _d 0, |
129 |
|
$ phis_var(i,j)) )/grav ) |
130 |
enddo |
enddo |
131 |
enddo |
enddo |
132 |
|
|
184 |
do L = 1,Lm |
do L = 1,Lm |
185 |
do j = 1,jm |
do j = 1,jm |
186 |
do i = 1,im |
do i = 1,im |
187 |
dragu(i,j,L) = sign( min(0.006 _d 0,abs(dragu(i,j,L))),dragu(i,j,L) ) |
dragu(i,j,L) = sign( min(0.006 _d 0,abs(dragu(i,j,L))), dragu(i |
188 |
dragv(i,j,L) = sign( min(0.006 _d 0,abs(dragv(i,j,L))),dragv(i,j,L) ) |
$ ,j,L) ) |
189 |
|
dragv(i,j,L) = sign( min(0.006 _d 0,abs(dragv(i,j,L))), dragv(i |
190 |
|
$ ,j,L) ) |
191 |
dragt(i,j,L) = -( uz(i,j,L)*dragu(i,j,L)+vz(i,j,L)*dragv(i,j,L) ) |
dragt(i,j,L) = -( uz(i,j,L)*dragu(i,j,L)+vz(i,j,L)*dragv(i,j,L) ) |
192 |
. *cpinv |
. *cpinv |
193 |
dudt(i,j,L) = dudt(i,j,L) + dragu(i,j,L) |
dudt(i,j,L) = dudt(i,j,L) + dragu(i,j,L) |