509 |
C anelastic: all transports & advect. fluxes are scaled by rhoFac |
C anelastic: all transports & advect. fluxes are scaled by rhoFac |
510 |
DO j=jMin,jMax |
DO j=jMin,jMax |
511 |
DO i=iMin,iMax |
DO i=iMin,iMax |
512 |
|
C to recover old (before 2009/11/30) results (since flxAdvUp(k=2) was zero) |
513 |
|
c IF (k.EQ.2) flxAdvUp(i,j) = 0. |
514 |
gW(i,j,k,bi,bj) = |
gW(i,j,k,bi,bj) = |
515 |
& -( ( flx_EW(i+1,j)-flx_EW(i,j) ) |
& -( ( flx_EW(i+1,j)-flx_EW(i,j) ) |
516 |
& + ( flx_NS(i,j+1)-flx_NS(i,j) ) |
& + ( flx_NS(i,j+1)-flx_NS(i,j) ) |