/[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.17 by jmc, Tue Mar 6 16:57:10 2001 UTC revision 1.18 by jmc, Thu Mar 8 20:40:57 2001 UTC
# Line 49  C--   Save previous solution & Initialis Line 49  C--   Save previous solution & Initialis
49  #ifdef INCLUDE_CD_CODE  #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  #endif
52            cg2d_x(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
55       &     + freeSurfFac*_rA(i,j,bi,bj)*horiVertRatio*       &     + freeSurfFac*_rA(i,j,bi,bj)*
56       &       EmPmR(I,J,bi,bj)/deltaTMom       &       EmPmR(I,J,bi,bj)/deltaTMom
57  #endif  #endif
58           ENDDO           ENDDO
# Line 87  C--   Add source term arising from w=d/d Line 87  C--   Add source term arising from w=d/d
87          DO j=1,sNy          DO j=1,sNy
88           DO i=1,sNx           DO i=1,sNx
89            cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj)            cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj)
90       &       +freeSurfFac*_rA(I,J,bi,bj)*horiVertRatio*(       &      -freeSurfFac*_rA(i,j,bi,bj)/deltaTMom/deltaTMom
91       &          -cg2d_x(I,J,bi,bj)       &        *( etaN(i,j,bi,bj)
92       &          -cg3d_x(I,J,1,bi,bj)       &          +cg3d_x(i,j,1,bi,bj)*horiVertRatio/gravity )
      &        )/deltaTMom/deltaTMom  
93            cg3d_b(i,j,1,bi,bj) = cg3d_b(i,j,1,bi,bj)            cg3d_b(i,j,1,bi,bj) = cg3d_b(i,j,1,bi,bj)
94       &      +freeSurfFac*_rA(I,J,bi,bj)*horiVertRatio*(       &      -freeSurfFac*_rA(i,j,bi,bj)/deltaTMom/deltaTMom
95       &         -cg2d_x(I,J,bi,bj)       &        *( etaN(i,j,bi,bj)
96       &         -cg3d_x(I,J,1,bi,bj)       &          +cg3d_x(i,j,1,bi,bj)*horiVertRatio/gravity )
97       &       )/deltaTMom/deltaTMom  C-jmc
98    c    &      -freeSurfFac*_rA(i,j,bi,bj)*recip_Bo(i,j,bi,bj)
99    c    &        *( cg2d_x(i,j,bi,bj) + cg3d_x(i,j,1,bi,bj) )
100    c    &        /deltaTMom/deltaTMom
101    C-jmc
102           ENDDO           ENDDO
103          ENDDO          ENDDO
104  #else  #else
105          DO j=1,sNy          DO j=1,sNy
106           DO i=1,sNx           DO i=1,sNx
107            cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj)            cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj)
108       &       +freeSurfFac*_rA(I,J,bi,bj)*horiVertRatio*(       &      -freeSurfFac*_rA(i,j,bi,bj)/deltaTMom/deltaTMom
109       &          -cg2d_x(I,J,bi,bj)       &        * etaN(i,j,bi,bj)
      &        )/deltaTMom/deltaTMom  
110           ENDDO           ENDDO
111          ENDDO          ENDDO
112  #endif  #endif
# Line 152  C--   Transfert the 2D-solution to "etaN Line 154  C--   Transfert the 2D-solution to "etaN
154         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
155          DO j=1-OLy,sNy+OLy          DO j=1-OLy,sNy+OLy
156           DO i=1-OLx,sNx+OLx           DO i=1-OLx,sNx+OLx
157            etaN(i,j,bi,bj) = cg2d_x(i,j,bi,bj)            etaN(i,j,bi,bj) = recip_Bo(i,j,bi,bj)*cg2d_x(i,j,bi,bj)
158           ENDDO           ENDDO
159          ENDDO          ENDDO
160         ENDDO         ENDDO
# Line 167  C     see CG3D.h for the interface to th Line 169  C     see CG3D.h for the interface to th
169          DO bi=myBxLo(myThid),myBxHi(myThid)          DO bi=myBxLo(myThid),myBxHi(myThid)
170           DO j=1,sNy+1           DO j=1,sNy+1
171            DO i=1,sNx+1            DO i=1,sNx+1
172             uf(i,j)=-gBaro*_recip_dxC(i,j,bi,bj)*             uf(i,j)=-_recip_dxC(i,j,bi,bj)*
173       &         (cg2d_x(i,j,bi,bj)-cg2d_x(i-1,j,bi,bj))       &         (cg2d_x(i,j,bi,bj)-cg2d_x(i-1,j,bi,bj))
174             vf(i,j)=-gBaro*_recip_dyC(i,j,bi,bj)*             vf(i,j)=-_recip_dyC(i,j,bi,bj)*
175       &         (cg2d_x(i,j,bi,bj)-cg2d_x(i,j-1,bi,bj))       &         (cg2d_x(i,j,bi,bj)-cg2d_x(i,j-1,bi,bj))
176            ENDDO            ENDDO
177           ENDDO           ENDDO
# Line 207  C Western boundary Line 209  C Western boundary
209       &       -dRF(K)*dYG( i ,j,bi,bj)*hFacW( i ,j,k,bi,bj)*uf( i ,j)       &       -dRF(K)*dYG( i ,j,bi,bj)*hFacW( i ,j,k,bi,bj)*uf( i ,j)
210       &       +dRF(K)*dXG(i,j+1,bi,bj)*hFacS(i,j+1,k,bi,bj)*vf(i,j+1)       &       +dRF(K)*dXG(i,j+1,bi,bj)*hFacS(i,j+1,k,bi,bj)*vf(i,j+1)
211       &       -dRF(K)*dXG(i, j ,bi,bj)*hFacS(i, j ,k,bi,bj)*vf(i, j )       &       -dRF(K)*dXG(i, j ,bi,bj)*hFacS(i, j ,k,bi,bj)*vf(i, j )
212       &       +(       &       +( freeSurfFac*etaN(i,j,bi,bj)/deltaTMom
213       &         -wVel(i,j,k+1,bi,bj)       &          -wVel(i,j,k+1,bi,bj)
214       &        )*_rA(i,j,bi,bj)/deltaTmom       &        )*_rA(i,j,bi,bj)/deltaTmom
      &       +freeSurfFac*_rA(I,J,bi,bj)*horiVertRatio*(  
      &          +cg2d_x(I,J,bi,bj)  
      &        )/deltaTMom/deltaTMom  
215            ENDDO            ENDDO
216           ENDDO           ENDDO
217           DO K=2,Nr-1           DO K=2,Nr-1
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22