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

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

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

revision 1.15 by mlosch, Tue May 25 07:48:22 2004 UTC revision 1.16 by edhill, Fri Jun 25 13:09:40 2004 UTC
# Line 128  C Flux on Southern face Line 128  C Flux on Southern face
128            DO I=iMin,iMax            DO I=iMin,iMax
129  C     First compute the fraction of open water for the w-control volume  C     First compute the fraction of open water for the w-control volume
130  C     at the southern face  C     at the southern face
131             hFacStmp=max(hFacS(I,J,K-1,bi,bj)-Half,0)             hFacStmp=max(hFacS(I,J,K-1,bi,bj)-Half,0. _d 0)
132       &         +    min(hFacS(I,J,K  ,bi,bj),Half)       &         +    min(hFacS(I,J,K  ,bi,bj),Half)
133             tmp_VbarZ=Half*(             tmp_VbarZ=Half*(
134       &          _hFacS(I,J,K-1,bi,bj)*vVel( I ,J,K-1,bi,bj)       &          _hFacS(I,J,K-1,bi,bj)*vVel( I ,J,K-1,bi,bj)
# Line 152  C Flux on Western face Line 152  C Flux on Western face
152            DO I=iMin,iMax+1            DO I=iMin,iMax+1
153  C     First compute the fraction of open water for the w-control volume  C     First compute the fraction of open water for the w-control volume
154  C     at the western face  C     at the western face
155             hFacWtmp=max(hFacW(I,J,K-1,bi,bj)-Half,0)             hFacWtmp=max(hFacW(I,J,K-1,bi,bj)-Half,0. _d 0)
156       &         +    min(hFacW(I,J,K  ,bi,bj),Half)       &         +    min(hFacW(I,J,K  ,bi,bj),Half)
157                   tmp_UbarZ=Half*(                   tmp_UbarZ=Half*(
158       &         _hFacW(I,J,K-1,bi,bj)*uVel( I ,J,K-1,bi,bj)       &         _hFacW(I,J,K-1,bi,bj)*uVel( I ,J,K-1,bi,bj)
Legend:
Removed from v.1.15  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22