/[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.9 by adcroft, Tue May 18 17:42:01 1999 UTC revision 1.10 by adcroft, Tue Mar 14 17:47:25 2000 UTC
# Line 50  C     pulls div(U) + dh/dt back toward z Line 50  C     pulls div(U) + dh/dt back toward z
50    
51        IF ( k .EQ. Nr ) THEN        IF ( k .EQ. Nr ) THEN
52  C      Initialise source term on first pass  C      Initialise source term on first pass
53         DO j=jMin,jMax         DO j=1,sNy
54          DO i=iMin,iMax          DO i=1,sNx
55    C Note: The source term containing cg2d_x and cg3d_x (at k=1)
56    C       has been moved to solve_for_pressure.F for convenience.
57           cg2d_b(i,j,bi,bj) =           cg2d_b(i,j,bi,bj) =
58       &     freeSurfFac*_rA(i,j,bi,bj)*horiVertRatio*(       &     freeSurfFac*_rA(i,j,bi,bj)*horiVertRatio*( 0.
      &      -cg2d_x(I,J,bi,bj)/deltaTMom/deltaTMom  
59  #ifdef USE_NATURAL_BCS  #ifdef USE_NATURAL_BCS
60       &      +EmPmR(I,J,bi,bj)/deltaTMom       &      +EmPmR(I,J,bi,bj)/deltaTMom
61  #endif  #endif
# Line 63  C      Initialise source term on first p Line 64  C      Initialise source term on first p
64         ENDDO         ENDDO
65        ENDIF        ENDIF
66    
67        DO j=jMin,jMax        DO j=1,sNy
68         DO i=iMin,iMax         DO i=1,sNx+1
69          pf(i,j) = xA(i,j)*gUNm1(i,j,k,bi,bj) / deltaTmom          pf(i,j) = xA(i,j)*gUNm1(i,j,k,bi,bj) / deltaTmom
70         ENDDO         ENDDO
71        ENDDO        ENDDO
# Line 85  C      Initialise source term on first p Line 86  C      Initialise source term on first p
86        ENDIF        ENDIF
87  #endif  #endif
88    
89        DO j=jMin,jMax        DO j=1,sNy+1
90         DO i=iMin,iMax         DO i=1,sNx
91          pf(i,j) = yA(i,j)*gVNm1(i,j,k,bi,bj) / deltatmom          pf(i,j) = yA(i,j)*gVNm1(i,j,k,bi,bj) / deltatmom
92         ENDDO         ENDDO
93        ENDDO        ENDDO
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22