/[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.12 by adcroft, Tue Mar 14 17:47:26 2000 UTC revision 1.13 by adcroft, Thu Jun 29 18:29:15 2000 UTC
# Line 73  C--   Save previous solution. Line 73  C--   Save previous solution.
73  C--   Add source term arising from w=d/dt (p_s + p_nh)  C--   Add source term arising from w=d/dt (p_s + p_nh)
74        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
75         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
76          K=1  #ifdef ALLOW_NONHYDROSTATIC
77          DO j=1,sNy          DO j=1,sNy
78           DO i=1,sNx           DO i=1,sNx
79            cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj)            cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj)
80       &       +freeSurfFac*_rA(I,J,bi,bj)*horiVertRatio*(       &       +freeSurfFac*_rA(I,J,bi,bj)*horiVertRatio*(
81       &          -cg2d_x(I,J,bi,bj)       &          -cg2d_x(I,J,bi,bj)
82  #ifdef ALLOW_NONHYDROSTATIC       &          -cg3d_x(I,J,1,bi,bj)
      &          -cg3d_x(I,J,K,bi,bj)  
 #endif  
83       &        )/deltaTMom/deltaTMom       &        )/deltaTMom/deltaTMom
84              cg3d_b(i,j,1,bi,bj) = cg3d_b(i,j,1,bi,bj)
85         &      +freeSurfFac*_rA(I,J,bi,bj)*horiVertRatio*(
86         &         -cg2d_x(I,J,bi,bj)
87         &         -cg3d_x(I,J,1,bi,bj)
88         &       )/deltaTMom/deltaTMom
89           ENDDO           ENDDO
90          ENDDO          ENDDO
91  #ifdef ALLOW_NONHYDROSTATIC  #else
         K=1  
92          DO j=1,sNy          DO j=1,sNy
93           DO i=1,sNx           DO i=1,sNx
94            cg3d_b(i,j,k,bi,bj) = cg3d_b(i,j,k,bi,bj)            cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj)
95       &      +freeSurfFac*_rA(I,J,bi,bj)*horiVertRatio*(       &       +freeSurfFac*_rA(I,J,bi,bj)*horiVertRatio*(
96       &         -cg2d_x(I,J,bi,bj)       &          -cg2d_x(I,J,bi,bj)
97       &         -cg3d_x(I,J,K,bi,bj)       &        )/deltaTMom/deltaTMom
      &       )/deltaTMom/deltaTMom  
98           ENDDO           ENDDO
99          ENDDO          ENDDO
100  #endif  #endif
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22