/[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.32 by jmc, Sat Jun 15 03:18:07 2002 UTC revision 1.33 by adcroft, Fri Jun 21 18:36:06 2002 UTC
# Line 112  C--   Add source term arising from w=d/d Line 112  C--   Add source term arising from w=d/d
112           DO j=1,sNy           DO j=1,sNy
113            DO i=1,sNx            DO i=1,sNx
114             cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj)             cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj)
115       &       -freeSurfFac*_rA(i,j,bi,bj)/deltaTMom/deltaTMom       &       -freeSurfFac*_rA(i,j,bi,bj)/deltaTMom/deltaTfreesurf
116       &         *( etaN(i,j,bi,bj)       &         *( etaN(i,j,bi,bj)
117       &           +phi_nh(i,j,1,bi,bj)*horiVertRatio/gravity )       &           +phi_nh(i,j,1,bi,bj)*horiVertRatio/gravity )
118             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)
119       &       -freeSurfFac*_rA(i,j,bi,bj)/deltaTMom/deltaTMom       &       -freeSurfFac*_rA(i,j,bi,bj)/deltaTMom/deltaTfreesurf
120       &         *( etaN(i,j,bi,bj)       &         *( etaN(i,j,bi,bj)
121       &           +phi_nh(i,j,1,bi,bj)*horiVertRatio/gravity )       &           +phi_nh(i,j,1,bi,bj)*horiVertRatio/gravity )
122            ENDDO            ENDDO
# Line 128  C--   Add source term arising from w=d/d Line 128  C--   Add source term arising from w=d/d
128           DO j=1,sNy           DO j=1,sNy
129            DO i=1,sNx            DO i=1,sNx
130             cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj)             cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj)
131       &       -freeSurfFac*_rA(i,j,bi,bj)/deltaTMom/deltaTMom       &       -freeSurfFac*_rA(i,j,bi,bj)/deltaTMom/deltaTfreesurf
132       &         * etaH(i,j,bi,bj)       &         * etaH(i,j,bi,bj)
133            ENDDO            ENDDO
134           ENDDO           ENDDO
# Line 136  C--   Add source term arising from w=d/d Line 136  C--   Add source term arising from w=d/d
136           DO j=1,sNy           DO j=1,sNy
137            DO i=1,sNx            DO i=1,sNx
138             cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj)             cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj)
139       &       -freeSurfFac*_rA(i,j,bi,bj)/deltaTMom/deltaTMom       &       -freeSurfFac*_rA(i,j,bi,bj)/deltaTMom/deltaTfreesurf
140       &         * etaN(i,j,bi,bj)       &         * etaN(i,j,bi,bj)
141            ENDDO            ENDDO
142           ENDDO           ENDDO
Legend:
Removed from v.1.32  
changed lines
  Added in v.1.33
  ViewVC Help
Powered by ViewVC 1.1.22