76 |
do j = jdim1,jdim2 |
do j = jdim1,jdim2 |
77 |
do i = idim1,idim2 |
do i = idim1,idim2 |
78 |
do L = 1,Lmplume |
do L = 1,Lmplume |
79 |
sqrtarg = (qplumeav(i,j,L)*qplumeav(i,j,L)) / |
if(vref(i,j,L,bi,bj).ne.0.) then |
80 |
|
sqrtarg = (qplumeav(i,j,L)*qplumeav(i,j,L)) / |
81 |
. ( ( (uref(i,j,L,bi,bj)*uref(i,j,L,bi,bj)) / |
. ( ( (uref(i,j,L,bi,bj)*uref(i,j,L,bi,bj)) / |
82 |
. (vref(i,j,L,bi,bj)*vref(i,j,L,bi,bj)) ) + 1. ) |
. (vref(i,j,L,bi,bj)*vref(i,j,L,bi,bj)) ) + 1. ) |
83 |
qdyn2(i,j,L,bi,bj) = sqrt(sqrtarg) |
qdyn2(i,j,L,bi,bj) = sqrt(sqrtarg) |
84 |
qdyn1(i,j,L,bi,bj) = qdyn2(i,j,L,bi,bj) * |
qdyn1(i,j,L,bi,bj) = qdyn2(i,j,L,bi,bj) * |
85 |
. (uref(i,j,L,bi,bj)/vref(i,j,L,bi,bj)) |
. (uref(i,j,L,bi,bj)/vref(i,j,L,bi,bj)) |
86 |
|
else |
87 |
|
qdyn1(i,j,L,bi,bj) = qplumeav(i,j,L) |
88 |
|
qdyn2(i,j,L,bi,bj) = 0. |
89 |
enddo |
enddo |
90 |
enddo |
enddo |
91 |
enddo |
enddo |