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

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

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

revision 1.17 by cnh, Wed Sep 26 18:09:13 2001 UTC revision 1.18 by jmc, Thu Apr 17 13:40:06 2003 UTC
# Line 64  C     pulls div(U) + dh/dt back toward z Line 64  C     pulls div(U) + dh/dt back toward z
64  C     Fully Implicit treatment of the Barotropic Flow Divergence  C     Fully Implicit treatment of the Barotropic Flow Divergence
65          DO j=1,sNy          DO j=1,sNy
66           DO i=1,sNx+1           DO i=1,sNx+1
67            pf(i,j) = xA(i,j)*gUNm1(i,j,k,bi,bj) / deltaTmom            pf(i,j) = xA(i,j)*gU(i,j,k,bi,bj) / deltaTmom
68           ENDDO           ENDDO
69          ENDDO          ENDDO
70        ELSEIF (exactConserv) THEN        ELSEIF (exactConserv) THEN
# Line 73  C     Implicit treatment of the Barotrop Line 73  C     Implicit treatment of the Barotrop
73          DO j=1,sNy          DO j=1,sNy
74           DO i=1,sNx+1           DO i=1,sNx+1
75            pf(i,j) = implicDiv2Dflow            pf(i,j) = implicDiv2Dflow
76       &             *xA(i,j)*gUNm1(i,j,k,bi,bj) / deltaTmom       &             *xA(i,j)*gU(i,j,k,bi,bj) / deltaTmom
77           ENDDO           ENDDO
78          ENDDO          ENDDO
79        ELSE        ELSE
# Line 93  C        ENDIF Line 93  C        ENDIF
93  C#endif  C#endif
94          DO j=1,sNy          DO j=1,sNy
95           DO i=1,sNx+1           DO i=1,sNx+1
96            pf(i,j) = ( implicDiv2Dflow * gUNm1(i,j,k,bi,bj)            pf(i,j) = ( implicDiv2Dflow * gU(i,j,k,bi,bj)
97       &          + (1.-implicDiv2Dflow) * uVel(i,j,k,bi,bj)       &          + (1.-implicDiv2Dflow) * uVel(i,j,k,bi,bj)
98       &               ) * xA(i,j) / deltaTmom       &               ) * xA(i,j) / deltaTmom
99           ENDDO           ENDDO
# Line 120  C#endif Line 120  C#endif
120  C     Fully Implicit treatment of the Barotropic Flow Divergence  C     Fully Implicit treatment of the Barotropic Flow Divergence
121          DO j=1,sNy+1          DO j=1,sNy+1
122           DO i=1,sNx           DO i=1,sNx
123            pf(i,j) = yA(i,j)*gVNm1(i,j,k,bi,bj) / deltatmom            pf(i,j) = yA(i,j)*gV(i,j,k,bi,bj) / deltatmom
124           ENDDO           ENDDO
125          ENDDO          ENDDO
126        ELSEIF (exactConserv) THEN        ELSEIF (exactConserv) THEN
# Line 129  C     Implicit treatment of the Barotrop Line 129  C     Implicit treatment of the Barotrop
129          DO j=1,sNy+1          DO j=1,sNy+1
130           DO i=1,sNx           DO i=1,sNx
131            pf(i,j) = implicDiv2Dflow            pf(i,j) = implicDiv2Dflow
132       &             *yA(i,j)*gVNm1(i,j,k,bi,bj) / deltatmom       &             *yA(i,j)*gV(i,j,k,bi,bj) / deltatmom
133           ENDDO           ENDDO
134          ENDDO          ENDDO
135        ELSE        ELSE
136  C     Explicit+Implicit part of the Barotropic Flow Divergence  C     Explicit+Implicit part of the Barotropic Flow Divergence
137          DO j=1,sNy+1          DO j=1,sNy+1
138           DO i=1,sNx           DO i=1,sNx
139            pf(i,j) = ( implicDiv2Dflow * gVNm1(i,j,k,bi,bj)            pf(i,j) = ( implicDiv2Dflow * gV(i,j,k,bi,bj)
140       &          + (1.-implicDiv2Dflow) * vVel(i,j,k,bi,bj)       &          + (1.-implicDiv2Dflow) * vVel(i,j,k,bi,bj)
141       &               ) * yA(i,j) / deltaTmom       &               ) * yA(i,j) / deltaTmom
142           ENDDO           ENDDO
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22