/[MITgcm]/MITgcm/verification/tidal_basin_2d/code/external_forcing.F
ViewVC logotype

Diff of /MITgcm/verification/tidal_basin_2d/code/external_forcing.F

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

revision 1.3 by adcroft, Thu May 15 15:19:02 2003 UTC revision 1.4 by jmc, Sun Jul 18 01:23:05 2004 UTC
# Line 62  C     Add windstress momentum impulse in Line 62  C     Add windstress momentum impulse in
62         DO j=jMin,jMax         DO j=jMin,jMax
63          DO i=iMin,iMax          DO i=iMin,iMax
64           gU(i,j,kLev,bi,bj) = gU(i,j,kLev,bi,bj)           gU(i,j,kLev,bi,bj) = gU(i,j,kLev,bi,bj)
65       &   +foFacMom*surfaceTendencyU(i,j,bi,bj)       &   +foFacMom*surfaceForcingU(i,j,bi,bj)
66       &   *_maskW(i,j,kLev,bi,bj)       &   *recip_drF(kLev)*recip_hFacW(i,j,kLev,bi,bj)
67          ENDDO          ENDDO
68         ENDDO         ENDDO
69        ENDIF        ENDIF
# Line 163  C     Add windstress momentum impulse in Line 163  C     Add windstress momentum impulse in
163         DO j=jMin,jMax         DO j=jMin,jMax
164          DO i=iMin,iMax          DO i=iMin,iMax
165           gV(i,j,kLev,bi,bj) = gV(i,j,kLev,bi,bj)           gV(i,j,kLev,bi,bj) = gV(i,j,kLev,bi,bj)
166       &   +foFacMom*surfaceTendencyV(i,j,bi,bj)       &   +foFacMom*surfaceForcingV(i,j,bi,bj)
167       &   *_maskS(i,j,kLev,bi,bj)       &   *recip_drF(kLev)*recip_hFacS(i,j,kLev,bi,bj)
168          ENDDO          ENDDO
169         ENDDO         ENDDO
170        ENDIF        ENDIF
# Line 243  C     Add heat in top-layer Line 243  C     Add heat in top-layer
243         DO j=jMin,jMax         DO j=jMin,jMax
244          DO i=iMin,iMax          DO i=iMin,iMax
245           gT(i,j,kLev,bi,bj)=gT(i,j,kLev,bi,bj)           gT(i,j,kLev,bi,bj)=gT(i,j,kLev,bi,bj)
246       &     +maskC(i,j,kLev,bi,bj)*surfaceTendencyT(i,j,bi,bj)       &     +surfaceForcingT(i,j,bi,bj)
247         &     *recip_drF(kLev)*recip_hFacC(i,j,kLev,bi,bj)
248          ENDDO          ENDDO
249         ENDDO         ENDDO
250        ENDIF        ENDIF
# Line 334  C     Add fresh-water in top-layer Line 335  C     Add fresh-water in top-layer
335         DO j=jMin,jMax         DO j=jMin,jMax
336          DO i=iMin,iMax          DO i=iMin,iMax
337           gS(i,j,kLev,bi,bj)=gS(i,j,kLev,bi,bj)           gS(i,j,kLev,bi,bj)=gS(i,j,kLev,bi,bj)
338       &   +maskC(i,j,kLev,bi,bj)*surfaceTendencyS(i,j,bi,bj)       &     +surfaceForcingS(i,j,bi,bj)
339         &     *recip_drF(kLev)*recip_hFacC(i,j,kLev,bi,bj)
340          ENDDO          ENDDO
341         ENDDO         ENDDO
342        ENDIF        ENDIF
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22