/[MITgcm]/MITgcm/model/src/impldiff.F
ViewVC logotype

Diff of /MITgcm/model/src/impldiff.F

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

revision 1.12 by heimbach, Mon Nov 13 16:30:32 2000 UTC revision 1.13 by adcroft, Wed Nov 29 22:38:31 2000 UTC
# Line 6  C     /================================= Line 6  C     /=================================
6  C     | S/R IMPLDIFF                                             |  C     | S/R IMPLDIFF                                             |
7  C     | o Solve implicit diffusion equation for vertical         |  C     | o Solve implicit diffusion equation for vertical         |
8  C     |   diffusivity.                                           |  C     |   diffusivity.                                           |
9    C     | o Recoded from 2d intermediate fields to 3d to reduce    |
10    C     |   TAMC storage                                           |
11    C     | o Fixed missing masks for fields a(), c()                |
12  C     \==========================================================/  C     \==========================================================/
13        SUBROUTINE IMPLDIFF( bi, bj, iMin, iMax, jMin, jMax,        SUBROUTINE IMPLDIFF( bi, bj, iMin, iMax, jMin, jMax,
14       I                     deltaTX,KappaRX,recip_hFac,       I                     deltaTX,KappaRX,recip_hFac,
# Line 67  C--   Old aUpper Line 70  C--   Old aUpper
70          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
71            c(i,j,k) = -deltaTX*recip_hFac(i,j,k,bi,bj)*recip_drF(k)            c(i,j,k) = -deltaTX*recip_hFac(i,j,k,bi,bj)*recip_drF(k)
72       &               *KappaRX(i,j,k+1)*recip_drC(k+1)       &               *KappaRX(i,j,k+1)*recip_drC(k+1)
73              IF (recip_hFac(i,j,k+1,bi,bj).EQ.0.) c(i,j,k)=0.
74          ENDDO          ENDDO
75         ENDDO         ENDDO
76        ENDDO        ENDDO
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22