/[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.25 by jmc, Thu Dec 15 21:08:59 2005 UTC revision 1.26 by heimbach, Wed Jun 7 01:55:12 2006 UTC
# Line 168  C     ... add difference of meridional f Line 168  C     ... add difference of meridional f
168             DO i=1-Olx,sNx+Olx             DO i=1-Olx,sNx+Olx
169  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
170  C     at the southern face  C     at the southern face
171              hFacCtmp=max(hFacC(I,J,K-1,bi,bj)-Half,0. _d 0)              hFacCtmp=max( _hFacC(I,J,K-1,bi,bj)-Half,0. _d 0 )
172       &           +   min(hFacC(I,J,K  ,bi,bj),Half)       &           +   min( _hFacC(I,J,K  ,bi,bj)     ,Half )
173              IF (hFacCtmp .GT. 0.) THEN              IF (hFacCtmp .GT. 0.) THEN
174               recip_hFacCtmp = 1./hFacCtmp               recip_hFacCtmp = 1./hFacCtmp
175              ELSE              ELSE
# Line 202  C Flux on Southern face Line 202  C Flux on Southern face
202            DO I=iMin,iMax            DO I=iMin,iMax
203  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
204  C     at the southern face  C     at the southern face
205             hFacStmp=max(hFacS(I,J,K-1,bi,bj)-Half,0. _d 0)             hFacStmp=max(_hFacS(I,J,K-1,bi,bj)-Half,0. _d 0)
206       &         +    min(hFacS(I,J,K  ,bi,bj),Half)       &         +    min(_hFacS(I,J,K  ,bi,bj),Half)
207             tmp_VbarZ=Half*(             tmp_VbarZ=Half*(
208       &          _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)
209       &         +_hFacS(I,J, K ,bi,bj)*vVel( I ,J, K ,bi,bj))       &         +_hFacS(I,J, K ,bi,bj)*vVel( I ,J, K ,bi,bj))
# Line 234  C Flux on Western face Line 234  C Flux on Western face
234            DO I=iMin,iMax+1            DO I=iMin,iMax+1
235  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
236  C     at the western face  C     at the western face
237             hFacWtmp=max(hFacW(I,J,K-1,bi,bj)-Half,0. _d 0)             hFacWtmp=max(_hFacW(I,J,K-1,bi,bj)-Half,0. _d 0)
238       &         +    min(hFacW(I,J,K  ,bi,bj),Half)       &         +    min(_hFacW(I,J,K  ,bi,bj),Half)
239             tmp_UbarZ=Half*(             tmp_UbarZ=Half*(
240       &         _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)
241       &        +_hFacW(I,J, K ,bi,bj)*uVel( I ,J, K ,bi,bj))       &        +_hFacW(I,J, K ,bi,bj)*uVel( I ,J, K ,bi,bj))
Legend:
Removed from v.1.25  
changed lines
  Added in v.1.26
  ViewVC Help
Powered by ViewVC 1.1.22