/[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.2 by cnh, Sat Aug 22 17:51:07 1998 UTC revision 1.12 by jmc, Tue Feb 20 15:06:21 2001 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #include "CPP_EEOPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6  C     /==========================================================\  C     /==========================================================\
7  C     | S/R CALC_DIV_GHAT                                        |  C     | S/R CALC_DIV_GHAT                                        |
8  C     | o Form the right hand-side of the surface pressure eqn.  |  C     | o Form the right hand-side of the surface pressure eqn.  |
9    C     |==========================================================|
10  C     \==========================================================/  C     \==========================================================/
11        SUBROUTINE CALC_DIV_GHAT(        SUBROUTINE CALC_DIV_GHAT(
12       I        bi,bj,iMin,iMax,jMin,jMax,       I        bi,bj,iMin,iMax,jMin,jMax,
# Line 22  C     == Global variables == Line 24  C     == Global variables ==
24  #include "PARAMS.h"  #include "PARAMS.h"
25  #include "GRID.h"  #include "GRID.h"
26  #include "CG2D.h"  #include "CG2D.h"
27    #ifdef ALLOW_NONHYDROSTATIC
28    #include "CG3D.h"
29    #endif
30    
31  C     == Routine arguments ==  C     == Routine arguments ==
 C     pH - Hydrostatic pressure  
32  C     bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation  C     bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation
33  C                                      results will be set.  C                                      results will be set.
34  C     kUp, kDown, kM1                - Index for upper and lower layers.  C     k                              - Index of layer.
35    C     xA, yA                         - Cell face areas
36  C     myThid - Instance number for this innvocation of CALC_MOM_RHS  C     myThid - Instance number for this innvocation of CALC_MOM_RHS
37        INTEGER bi,bj,iMin,iMax,jMin,jMax        INTEGER bi,bj,iMin,iMax,jMin,jMax
38        INTEGER K        INTEGER K
# Line 46  C     pulls div(U) + dh/dt back toward z Line 51  C     pulls div(U) + dh/dt back toward z
51    
52        IF ( k .EQ. Nr ) THEN        IF ( k .EQ. Nr ) THEN
53  C      Initialise source term on first pass  C      Initialise source term on first pass
54         DO j=jMin,jMax         DO j=1,sNy
55          DO i=iMin,iMax          DO i=1,sNx
56  C        cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj)  C Note: The source term containing cg2d_x and cg3d_x (at k=1)
57  C    &    -freeSurfFac*_rA(i,j,bi,bj)*  C       has been moved to solve_for_pressure.F for convenience.
58  C    &     cg2d_x(I  ,J  ,bi,bj)/deltaTMom/deltaTMom           cg2d_b(i,j,bi,bj) = 0.
59           cg2d_b(i,j,bi,bj) =  #ifdef USE_NATURAL_BCS
60       &    -freeSurfFac*_rA(i,j,bi,bj)*       &     + freeSurfFac*_rA(i,j,bi,bj)*horiVertRatio*
61       &     cg2d_x(I  ,J  ,bi,bj)/deltaTMom/deltaTMom       &       EmPmR(I,J,bi,bj)/deltaTMom
62    #endif
63          ENDDO          ENDDO
64         ENDDO         ENDDO
65        ENDIF        ENDIF
66    
67        DO j=jMin,jMax        IF (implicDiv2Dflow.EQ.1.) then
68         DO i=iMin,iMax  C     Fully Implicit treatment of the Barotropic Flow Divergence
69          pf(i,j) = xA(i,j)*gUNm1(i,j,k,bi,bj) / deltaTmom          DO j=1,sNy
70         ENDDO           DO i=1,sNx+1
71        ENDDO            pf(i,j) = xA(i,j)*gUNm1(i,j,k,bi,bj) / deltaTmom
72             ENDDO
73            ENDDO
74          ELSE
75    C     Explicit+Implicit part of the Barotropic Flow Divergence
76    C      => Filtering of uVel,vVel is necessary
77    #ifdef ALLOW_ZONAL_FILT
78            IF (zonal_filt_lat.LT.90.) THEN
79              CALL ZONAL_FILTER(
80         &      uVel, hFacW, 1-1, sNy+1, k, k, bi, bj, 1, myThid)
81              CALL ZONAL_FILTER(
82         &      vVel, hFacS, 1-1, sNy+1, k, k, bi, bj, 2, myThid)
83            ENDIF
84    #endif
85            DO j=1,sNy
86             DO i=1,sNx+1
87              pf(i,j) = ( implicDiv2Dflow * gUNm1(i,j,k,bi,bj)
88         &          + (1.-implicDiv2Dflow) * uVel(i,j,k,bi,bj)
89         &               ) * xA(i,j) / deltaTmom
90             ENDDO
91            ENDDO
92          ENDIF
93        DO j=1,sNy        DO j=1,sNy
94         DO i=1,sNx         DO i=1,sNx
95          cg2d_b(i,j,bi,bj) =  cg2d_b(i,j,bi,bj) +          cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj) +
96       &   pf(i+1,j)-pf(i,j)       &   pf(i+1,j)-pf(i,j)
97         ENDDO         ENDDO
98        ENDDO        ENDDO
99    
100        DO j=jMin,jMax  #ifdef ALLOW_NONHYDROSTATIC
101         DO i=iMin,iMax        IF (nonHydrostatic) THEN
102          pf(i,j) = yA(i,j)*gVNm1(i,j,k,bi,bj) / deltatmom         DO j=1,sNy
103            DO i=1,sNx
104             cg3d_b(i,j,k,bi,bj) = pf(i+1,j)-pf(i,j)
105            ENDDO
106         ENDDO         ENDDO
107        ENDDO        ENDIF
108    #endif
109    
110          IF (implicDiv2Dflow.EQ.1.) then
111    C     Fully Implicit treatment of the Barotropic Flow Divergence
112            DO j=1,sNy+1
113             DO i=1,sNx
114              pf(i,j) = yA(i,j)*gVNm1(i,j,k,bi,bj) / deltatmom
115             ENDDO
116            ENDDO
117          ELSE
118    C     Explicit+Implicit part of the Barotropic Flow Divergence
119            DO j=1,sNy+1
120             DO i=1,sNx
121              pf(i,j) = ( implicDiv2Dflow * gVNm1(i,j,k,bi,bj)
122         &          + (1.-implicDiv2Dflow) * vVel(i,j,k,bi,bj)
123         &               ) * yA(i,j) / deltaTmom
124             ENDDO
125            ENDDO
126          ENDIF
127        DO j=1,sNy        DO j=1,sNy
128         DO i=1,sNx         DO i=1,sNx
129          cg2d_b(i,j,bi,bj) =  cg2d_b(i,j,bi,bj) +          cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj) +
130       &   pf(i,j+1)-pf(i,j)       &   pf(i,j+1)-pf(i,j)
131         ENDDO         ENDDO
132        ENDDO        ENDDO
133    
134    #ifdef ALLOW_NONHYDROSTATIC
135          IF (nonHydrostatic) THEN
136           DO j=1,sNy
137            DO i=1,sNx
138             cg3d_b(i,j,k,bi,bj) = cg3d_b(i,j,k,bi,bj) +
139         &    pf(i,j+1)-pf(i,j)
140            ENDDO
141           ENDDO
142          ENDIF
143    #endif
144    
145        RETURN        RETURN
146        END        END
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22