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 |