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

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

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

revision 1.16 by heimbach, Tue Sep 11 01:32:02 2012 UTC revision 1.17 by gforget, Tue Sep 18 21:29:36 2012 UTC
# Line 94  C errCode   :: > 0 if singular matrix Line 94  C errCode   :: > 0 if singular matrix
94        _RL df (1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RL df (1-Olx:sNx+Olx,1-Oly:sNy+Oly)
95        _RL af (1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RL af (1-Olx:sNx+Olx,1-Oly:sNy+Oly)
96        _RL div(1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RL div(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
97          _RL flx(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
98  #endif  #endif
99  CEOP  CEOP
100    
# Line 297  C--   Set diagnostic suffix for the curr Line 298  C--   Set diagnostic suffix for the curr
298          IF ( diagDif .OR. diagAdv ) THEN          IF ( diagDif .OR. diagAdv ) THEN
299           DO j=1-OLy,sNy+OLy           DO j=1-OLy,sNy+OLy
300            DO i=1-OLx,sNx+OLx            DO i=1-OLx,sNx+OLx
301              af(i,j) = 0. _d 0              flx(i,j) = 0. _d 0
302            ENDDO            ENDDO
303           ENDDO           ENDDO
304           DO k= Nr,1,-1           DO k= Nr,1,-1
# Line 396  C--    use current hFac (consistent with Line 397  C--    use current hFac (consistent with
397             ENDIF             ENDIF
398             DO j=jMin,jMax             DO j=jMin,jMax
399              DO i=iMin,iMax              DO i=iMin,iMax
400                af(i,j) = af(i,j)                flx(i,j) = flx(i,j)
401  #ifdef ALLOW_AUTODIFF_OPENAD  #ifdef ALLOW_AUTODIFF_OPENAD
402       &                - rkSign*div(i,j)*rA(i,j,bi,bj)%v/deltaTLev(k)       &                - rkSign*div(i,j)*rA(i,j,bi,bj)%v/deltaTLev(k)
403  #else  #else
404       &                - rkSign*div(i,j)*rA(i,j,bi,bj)/deltaTLev(k)       &                - rkSign*div(i,j)*rA(i,j,bi,bj)/deltaTLev(k)
405  #endif  #endif
406       &                - df(i,j)                af(i,j) = flx(i,j) - df(i,j)
407              ENDDO              ENDDO
408             ENDDO             ENDDO
409             diagName = 'ADVr'//diagSufx             diagName = 'ADVr'//diagSufx
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22