/[MITgcm]/MITgcm/pkg/generic_advdiff/gad_calc_rhs.F
ViewVC logotype

Diff of /MITgcm/pkg/generic_advdiff/gad_calc_rhs.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.57 by jmc, Mon Jun 6 15:44:00 2011 UTC revision 1.58 by jmc, Thu Dec 1 14:16:30 2011 UTC
# Line 161  C--   Set diagnostic suffix for the curr Line 161  C--   Set diagnostic suffix for the curr
161        advFac  = 0. _d 0        advFac  = 0. _d 0
162        IF (calcAdvection) advFac = 1. _d 0        IF (calcAdvection) advFac = 1. _d 0
163        rAdvFac = rkSign*advFac        rAdvFac = rkSign*advFac
164        IF (implicitAdvection) rAdvFac = 0. _d 0        IF (implicitAdvection) rAdvFac = rkSign
165    
166        DO j=1-OLy,sNy+OLy        DO j=1-OLy,sNy+OLy
167         DO i=1-OLx,sNx+OLx         DO i=1-OLx,sNx+OLx
# Line 724  C     Anelastic: scale vertical fluxes b Line 724  C     Anelastic: scale vertical fluxes b
724       &   *( (fZon(i+1,j)-fZon(i,j))       &   *( (fZon(i+1,j)-fZon(i,j))
725       &     +(fMer(i,j+1)-fMer(i,j))       &     +(fMer(i,j+1)-fMer(i,j))
726       &     +(fVerT(i,j,kDown)-fVerT(i,j,kUp))*rkSign       &     +(fVerT(i,j,kDown)-fVerT(i,j,kUp))*rkSign
727       &     -localT(i,j)*( (uTrans(i+1,j)-uTrans(i,j))       &     -localT(i,j)*( (uTrans(i+1,j)-uTrans(i,j))*advFac
728       &                   +(vTrans(i,j+1)-vTrans(i,j))       &                   +(vTrans(i,j+1)-vTrans(i,j))*advFac
729       &                   +(rTransKp1(i,j)-rTrans(i,j))*rAdvFac       &                   +(rTransKp1(i,j)-rTrans(i,j))*rAdvFac
730       &                  )*advFac       &                  )
731       &    )       &    )
732         ENDDO         ENDDO
733        ENDDO        ENDDO
Legend:
Removed from v.1.57  
changed lines
  Added in v.1.58
  ViewVC Help
Powered by ViewVC 1.1.22