/[MITgcm]/MITgcm/model/src/external_forcing.F
ViewVC logotype

Diff of /MITgcm/model/src/external_forcing.F

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

revision 1.25 by heimbach, Fri May 14 21:08:28 2004 UTC revision 1.26 by jmc, Sun Jul 18 01:04:23 2004 UTC
# Line 79  C     Add windstress momentum impulse in Line 79  C     Add windstress momentum impulse in
79         DO j=jMin,jMax         DO j=jMin,jMax
80          DO i=iMin,iMax          DO i=iMin,iMax
81           gU(i,j,kLev,bi,bj) = gU(i,j,kLev,bi,bj)           gU(i,j,kLev,bi,bj) = gU(i,j,kLev,bi,bj)
82       &   +foFacMom*surfaceTendencyU(i,j,bi,bj)       &   +foFacMom*surfaceForcingU(i,j,bi,bj)
83       &   *_maskW(i,j,kLev,bi,bj)       &   *recip_drF(kLev)*recip_hFacW(i,j,kLev,bi,bj)
84          ENDDO          ENDDO
85         ENDDO         ENDDO
86        ENDIF        ENDIF
# Line 167  C     Add windstress momentum impulse in Line 167  C     Add windstress momentum impulse in
167         DO j=jMin,jMax         DO j=jMin,jMax
168          DO i=iMin,iMax          DO i=iMin,iMax
169           gV(i,j,kLev,bi,bj) = gV(i,j,kLev,bi,bj)           gV(i,j,kLev,bi,bj) = gV(i,j,kLev,bi,bj)
170       &   +foFacMom*surfaceTendencyV(i,j,bi,bj)       &   +foFacMom*surfaceForcingV(i,j,bi,bj)
171       &   *_maskS(i,j,kLev,bi,bj)       &   *recip_drF(kLev)*recip_hFacS(i,j,kLev,bi,bj)
172          ENDDO          ENDDO
173         ENDDO         ENDDO
174        ENDIF        ENDIF
# Line 264  C     Add heat in top-layer Line 264  C     Add heat in top-layer
264         DO j=jMin,jMax         DO j=jMin,jMax
265          DO i=iMin,iMax          DO i=iMin,iMax
266           gT(i,j,kLev,bi,bj)=gT(i,j,kLev,bi,bj)           gT(i,j,kLev,bi,bj)=gT(i,j,kLev,bi,bj)
267       &     +maskC(i,j,kLev,bi,bj)*surfaceTendencyT(i,j,bi,bj)       &     +surfaceForcingT(i,j,bi,bj)
268         &     *recip_drF(kLev)*recip_hFacC(i,j,kLev,bi,bj)
269          ENDDO          ENDDO
270         ENDDO         ENDDO
271        ENDIF        ENDIF
# Line 288  C Penetrating SW radiation Line 289  C Penetrating SW radiation
289       &   -Qsw(i,j,bi,bj)*(swfracb(1)*maskC(i,j,klev,bi,bj)       &   -Qsw(i,j,bi,bj)*(swfracb(1)*maskC(i,j,klev,bi,bj)
290       &                   -swfracb(2)*maskC(i,j,kp1, bi,bj))       &                   -swfracb(2)*maskC(i,j,kp1, bi,bj))
291       &    *recip_Cp*recip_rhoConst*recip_drF(klev)       &    *recip_Cp*recip_rhoConst*recip_drF(klev)
292    c    &                            *recip_hFacC(i,j,kLev,bi,bj)
293         ENDDO         ENDDO
294        ENDDO        ENDDO
295  #endif  #endif
# Line 377  C     Add fresh-water in top-layer Line 379  C     Add fresh-water in top-layer
379         DO j=jMin,jMax         DO j=jMin,jMax
380          DO i=iMin,iMax          DO i=iMin,iMax
381           gS(i,j,kLev,bi,bj)=gS(i,j,kLev,bi,bj)           gS(i,j,kLev,bi,bj)=gS(i,j,kLev,bi,bj)
382       &   +maskC(i,j,kLev,bi,bj)*surfaceTendencyS(i,j,bi,bj)       &     +surfaceForcingS(i,j,bi,bj)
383         &     *recip_drF(kLev)*recip_hFacC(i,j,kLev,bi,bj)
384          ENDDO          ENDDO
385         ENDDO         ENDDO
386        ENDIF        ENDIF

Legend:
Removed from v.1.25  
changed lines
  Added in v.1.26

  ViewVC Help
Powered by ViewVC 1.1.22