235 |
ENDDO |
ENDDO |
236 |
ENDDO |
ENDDO |
237 |
|
|
238 |
C Dummy values for the edges |
C Dummy values for the edges. This does not affect the results |
239 |
C This avoids weirdness in gmredi_calc_eigs |
C but avoids problems when solving for the eigenvalues. |
240 |
i=1-Olx |
i=1-Olx |
241 |
DO j=1-Oly,sNy+Oly |
DO j=1-Oly,sNy+Oly |
242 |
kLow_U(i,j) = kLow_C(i,j) |
kLow_U(i,j) = 0 |
243 |
ENDDO |
ENDDO |
244 |
j=1-Oly |
j=1-Oly |
245 |
DO i=1-Olx,sNx+Olx |
DO i=1-Olx,sNx+Olx |
246 |
kLow_V(i,j) = kLow_C(i,j) |
kLow_V(i,j) = 0 |
247 |
ENDDO |
ENDDO |
248 |
|
|
249 |
g_reciprho_sq = (gravity*recip_rhoConst)**2 |
g_reciprho_sq = (gravity*recip_rhoConst)**2 |