| 46 |
_RL qdyn2(idim1:idim2,jdim1:jdim2,Lmplume,Nsx,Nsy) |
_RL qdyn2(idim1:idim2,jdim1:jdim2,Lmplume,Nsx,Nsy) |
| 47 |
|
|
| 48 |
integer i,j,L,iplume |
integer i,j,L,iplume |
| 49 |
_RL qplumeav(i1,j2,Lmplume) |
_RL qplumeav(i2,j2,Lmplume) |
| 50 |
_RL sqrtarg |
_RL sqrtarg |
| 51 |
|
|
| 52 |
C First step - compute the average of qplume over Nxplume |
C First step - compute the average of qplume over Nxplume |
| 88 |
else |
else |
| 89 |
qdyn1(i,j,L,bi,bj) = qplumeav(i,j,L) |
qdyn1(i,j,L,bi,bj) = qplumeav(i,j,L) |
| 90 |
qdyn2(i,j,L,bi,bj) = 0. |
qdyn2(i,j,L,bi,bj) = 0. |
| 91 |
|
endif |
| 92 |
enddo |
enddo |
| 93 |
enddo |
enddo |
| 94 |
enddo |
enddo |