/[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.26 by jmc, Mon Nov 30 16:26:48 2009 UTC revision 1.27 by jmc, Fri Nov 9 22:37:05 2012 UTC
# Line 29  C     == Global variables == Line 29  C     == Global variables ==
29  #include "PARAMS.h"  #include "PARAMS.h"
30  #include "GRID.h"  #include "GRID.h"
31  #include "DYNVARS.h"  #include "DYNVARS.h"
32    #ifdef ALLOW_ADDFLUID
33    # include "FFIELDS.h"
34    #endif
35    
36  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
37  C     == Routine arguments ==  C     == Routine arguments ==
# Line 77  C     pulls div(U) + dh/dt back toward z Line 80  C     pulls div(U) + dh/dt back toward z
80  C     Fully Implicit treatment of the Barotropic Flow Divergence  C     Fully Implicit treatment of the Barotropic Flow Divergence
81          DO j=1,sNy          DO j=1,sNy
82           DO i=1,sNx+1           DO i=1,sNx+1
83            pf(i,j) = xA(i,j)*gU(i,j,k,bi,bj) / deltaTmom            pf(i,j) = xA(i,j)*gU(i,j,k,bi,bj) / deltaTMom
84           ENDDO           ENDDO
85          ENDDO          ENDDO
86        ELSEIF (exactConserv) THEN        ELSEIF (exactConserv) THEN
# Line 86  C     Implicit treatment of the Barotrop Line 89  C     Implicit treatment of the Barotrop
89          DO j=1,sNy          DO j=1,sNy
90           DO i=1,sNx+1           DO i=1,sNx+1
91            pf(i,j) = implicDiv2Dflow            pf(i,j) = implicDiv2Dflow
92       &             *xA(i,j)*gU(i,j,k,bi,bj) / deltaTmom       &             *xA(i,j)*gU(i,j,k,bi,bj) / deltaTMom
93           ENDDO           ENDDO
94          ENDDO          ENDDO
95        ELSE        ELSE
# Line 111  c#endif Line 114  c#endif
114           DO i=1,sNx+1           DO i=1,sNx+1
115            pf(i,j) = ( implicDiv2Dflow * gU(i,j,k,bi,bj)            pf(i,j) = ( implicDiv2Dflow * gU(i,j,k,bi,bj)
116       &     + (1. _d 0-implicDiv2Dflow)* uVel(i,j,k,bi,bj)       &     + (1. _d 0-implicDiv2Dflow)* uVel(i,j,k,bi,bj)
117       &               ) * xA(i,j) / deltaTmom       &               ) * xA(i,j) / deltaTMom
118           ENDDO           ENDDO
119          ENDDO          ENDDO
120        ENDIF        ENDIF
# Line 154  C     Explicit+Implicit part of the Baro Line 157  C     Explicit+Implicit part of the Baro
157           DO i=1,sNx           DO i=1,sNx
158            pf(i,j) = ( implicDiv2Dflow * gV(i,j,k,bi,bj)            pf(i,j) = ( implicDiv2Dflow * gV(i,j,k,bi,bj)
159       &     + (1. _d 0-implicDiv2Dflow)* vVel(i,j,k,bi,bj)       &     + (1. _d 0-implicDiv2Dflow)* vVel(i,j,k,bi,bj)
160       &               ) * yA(i,j) / deltaTmom       &               ) * yA(i,j) / deltaTMom
161           ENDDO           ENDDO
162          ENDDO          ENDDO
163        ENDIF        ENDIF
# Line 181  C     Explicit+Implicit part of the Baro Line 184  C     Explicit+Implicit part of the Baro
184          DO j=1,sNy          DO j=1,sNy
185           DO i=1,sNx           DO i=1,sNx
186            cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj)            cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj)
187       &         - addMass(i,j,k,bi,bj)*mass2rUnit/deltaTmom       &         - addMass(i,j,k,bi,bj)*mass2rUnit/deltaTMom
188           ENDDO           ENDDO
189          ENDDO          ENDDO
190  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
# Line 189  C     Explicit+Implicit part of the Baro Line 192  C     Explicit+Implicit part of the Baro
192          DO j=1,sNy          DO j=1,sNy
193           DO i=1,sNx           DO i=1,sNx
194            cg3d_b(i,j,k,bi,bj) = cg3d_b(i,j,k,bi,bj)            cg3d_b(i,j,k,bi,bj) = cg3d_b(i,j,k,bi,bj)
195       &         - addMass(i,j,k,bi,bj)*mass2rUnit/deltaTmom       &         - addMass(i,j,k,bi,bj)*mass2rUnit/deltaTMom
196           ENDDO           ENDDO
197          ENDDO          ENDDO
198         ENDIF         ENDIF
Legend:
Removed from v.1.26  
changed lines
  Added in v.1.27
  ViewVC Help
Powered by ViewVC 1.1.22