/[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.20 by jmc, Sat Jul 30 23:38:54 2005 UTC revision 1.21 by jmc, Tue Aug 16 22:42:49 2005 UTC
# Line 61  C     flx_Dn       :: Temp. used for fVo Line 61  C     flx_Dn       :: Temp. used for fVo
61        _RL    fZon(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL    fZon(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
62        _RL    fMer(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL    fMer(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
63        _RL    del2w(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL    del2w(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
64  C     I,J,K - Loop counters  C     i,j,k - Loop counters
65        INTEGER i,j,k, kP1, kUp        INTEGER i,j,k, kP1
66        _RL  wOverride        _RL  wOverride
67        _RS  hFacWtmp        _RS  hFacWtmp
68        _RS  hFacStmp        _RS  hFacStmp
# Line 102  C        slipSideFac = 0. _d 0 Line 102  C        slipSideFac = 0. _d 0
102          DO K=1,Nr          DO K=1,Nr
103           DO j=1-OLy,sNy+OLy           DO j=1-OLy,sNy+OLy
104            DO i=1-OLx,sNx+OLx            DO i=1-OLx,sNx+OLx
105             gWNM1(i,j,k,bi,bj) = gW(i,j,k,bi,bj)             gwNm1(i,j,k,bi,bj) = gW(i,j,k,bi,bj)
106             gW(i,j,k,bi,bj) = 0.             gW(i,j,k,bi,bj) = 0.
107            ENDDO            ENDDO
108           ENDDO           ENDDO
# Line 182  C     First compute the fraction of open Line 182  C     First compute the fraction of open
182  C     at the southern face  C     at the southern face
183              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)
184       &           +   min(hFacC(I,J,K  ,bi,bj),Half)       &           +   min(hFacC(I,J,K  ,bi,bj),Half)
             recip_hFacCtmp = 0. _d 0  
185              IF (hFacCtmp .GT. 0.) THEN              IF (hFacCtmp .GT. 0.) THEN
186               recip_hFacCtmp = 1./hFacCtmp               recip_hFacCtmp = 1./hFacCtmp
187              ELSE              ELSE
# Line 201  CML ************************************ Line 200  CML ************************************
200  CML   No-slip Boundary conditions for bi-harmonic dissipation  CML   No-slip Boundary conditions for bi-harmonic dissipation
201  CML   need to be implemented here!  CML   need to be implemented here!
202  CML ************************************************************  CML ************************************************************
203             ELSE
204    C-    Initialize del2w to zero:
205              DO j=1-Oly,sNy+Oly
206               DO i=1-Olx,sNx+Olx
207                del2w(i,j) = 0. _d 0
208               ENDDO
209              ENDDO
210           ENDIF           ENDIF
211    
212  C Flux on Southern face  C Flux on Southern face
# Line 242  C     First compute the fraction of open Line 248  C     First compute the fraction of open
248  C     at the western face  C     at the western face
249             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)
250       &         +    min(hFacW(I,J,K  ,bi,bj),Half)       &         +    min(hFacW(I,J,K  ,bi,bj),Half)
251                   tmp_UbarZ=Half*(             tmp_UbarZ=Half*(
252       &         _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)
253       &        +_hFacW(I,J, K ,bi,bj)*uVel( I ,J, K ,bi,bj))       &        +_hFacW(I,J, K ,bi,bj)*uVel( I ,J, K ,bi,bj))
254             Flx_EW(I,J,bi,bj)=             Flx_EW(I,J,bi,bj)=
# Line 299  caja           and an hFacUS (above V po Line 305  caja           and an hFacUS (above V po
305         ENDDO         ENDDO
306        ENDDO        ENDDO
307    
308        
309        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
310         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
311          DO K=2,Nr          DO K=2,Nr
# Line 307  caja           and an hFacUS (above V po Line 313  caja           and an hFacUS (above V po
313            DO i=iMin,iMax            DO i=iMin,iMax
314             wVel(i,j,k,bi,bj) = wVel(i,j,k,bi,bj)             wVel(i,j,k,bi,bj) = wVel(i,j,k,bi,bj)
315       &     +deltatMom*nh_Am2*( ab15*gW(i,j,k,bi,bj)       &     +deltatMom*nh_Am2*( ab15*gW(i,j,k,bi,bj)
316       &                 +ab05*gWNM1(i,j,k,bi,bj) )       &                 +ab05*gwNm1(i,j,k,bi,bj) )
317             IF (hFacC(I,J,K,bi,bj).EQ.0.) wVel(i,j,k,bi,bj)=0.             IF (hFacC(I,J,K,bi,bj).EQ.0.) wVel(i,j,k,bi,bj)=0.
318            ENDDO            ENDDO
319           ENDDO           ENDDO
Legend:
Removed from v.1.20  
changed lines
  Added in v.1.21
  ViewVC Help
Powered by ViewVC 1.1.22