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

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

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

revision 1.25 by adcroft, Fri Jun 29 17:14:49 2001 UTC revision 1.26 by jmc, Wed Sep 19 13:58:08 2001 UTC
# Line 46  C--   Save previous solution & Initialis Line 46  C--   Save previous solution & Initialis
46         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
47          DO j=1-OLy,sNy+OLy          DO j=1-OLy,sNy+OLy
48           DO i=1-OLx,sNx+OLx           DO i=1-OLx,sNx+OLx
49    #ifdef INCLUDE_CD_CODE
50            etaNm1(i,j,bi,bj) = etaN(i,j,bi,bj)            etaNm1(i,j,bi,bj) = etaN(i,j,bi,bj)
51    #endif
52            cg2d_x(i,j,bi,bj) = Bo_surf(i,j,bi,bj)*etaN(i,j,bi,bj)            cg2d_x(i,j,bi,bj) = Bo_surf(i,j,bi,bj)*etaN(i,j,bi,bj)
53            cg2d_b(i,j,bi,bj) = 0.            cg2d_b(i,j,bi,bj) = 0.
54  #ifdef USE_NATURAL_BCS  #ifdef USE_NATURAL_BCS
# Line 92  C--   Add source term arising from w=d/d Line 94  C--   Add source term arising from w=d/d
94       &      -freeSurfFac*_rA(i,j,bi,bj)/deltaTMom/deltaTMom       &      -freeSurfFac*_rA(i,j,bi,bj)/deltaTMom/deltaTMom
95       &        *( etaN(i,j,bi,bj)       &        *( etaN(i,j,bi,bj)
96       &          +phi_nh(i,j,1,bi,bj)*horiVertRatio/gravity )       &          +phi_nh(i,j,1,bi,bj)*horiVertRatio/gravity )
 C-jmc  
 c    &      -freeSurfFac*_rA(i,j,bi,bj)*recip_Bo(i,j,bi,bj)  
 c    &        *( cg2d_x(i,j,bi,bj) + phi_nh(i,j,1,bi,bj) )  
 c    &        /deltaTMom/deltaTMom  
 C-jmc  
97           ENDDO           ENDDO
98          ENDDO          ENDDO
99  #else  #else
100          DO j=1,sNy          IF ( exactConserv ) THEN
101           DO i=1,sNx  c       IF (nonlinFreeSurf.GT.0) THEN
102            cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj)           DO j=1,sNy
103       &      -freeSurfFac*_rA(i,j,bi,bj)/deltaTMom/deltaTMom            DO i=1,sNx
104       &        * etaN(i,j,bi,bj)             cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj)
105         &       -freeSurfFac*_rA(i,j,bi,bj)/deltaTMom/deltaTMom
106         &         * etaH(i,j,bi,bj)
107              ENDDO
108           ENDDO           ENDDO
109          ENDDO          ELSE
110             DO j=1,sNy
111              DO i=1,sNx
112               cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj)
113         &       -freeSurfFac*_rA(i,j,bi,bj)/deltaTMom/deltaTMom
114         &         * etaN(i,j,bi,bj)
115              ENDDO
116             ENDDO
117            ENDIF
118  #endif  #endif
119    
120  #ifdef ALLOW_OBCS  #ifdef ALLOW_OBCS
Legend:
Removed from v.1.25  
changed lines
  Added in v.1.26
  ViewVC Help
Powered by ViewVC 1.1.22