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

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

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

revision 1.2 by jmc, Wed Jun 20 23:06:51 2007 UTC revision 1.3 by jmc, Mon Aug 11 22:25:52 2008 UTC
# Line 99  c       hMixCriteria  = -0.8 _d 0 Line 99  c       hMixCriteria  = -0.8 _d 0
99             rhoKm1(i,j) = rhoSurf(i,j)             rhoKm1(i,j) = rhoSurf(i,j)
100             rhoMxL(i,j) = rhoSurf(i,j)             rhoMxL(i,j) = rhoSurf(i,j)
101       &                 + MAX( rhoMxL(i,j)*hMixCriteria, dRhoSmall )       &                 + MAX( rhoMxL(i,j)*hMixCriteria, dRhoSmall )
102             hMixLayer(i,j,bi,bj) = rF(1)-R_low(I,J,bi,bj)             hMixLayer(i,j,bi,bj) = rF(1)-R_low(i,j,bi,bj)
103           ENDDO           ENDDO
104          ENDDO          ENDDO
105          DO k = 2,Nr          DO k = 2,Nr
106  C-    potential density (reference level = surface level)  C-    potential density (reference level = surface level)
107           CALL FIND_RHO(           CALL FIND_RHO_2D(
108       I             bi, bj, 1-OLx, sNx+OLx, 1-OLy, sNy+OLy, K, 1,       I        1-OLx, sNx+OLx, 1-OLy, sNy+OLy, 1,
109       I             theta, salt,       I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
110       O             rhoLoc, myThid )       O        rhoLoc,
111         I        k, bi, bj, myThid )
112    
113           DO j=1-Oly,sNy+Oly           DO j=1-Oly,sNy+Oly
114            DO i=1-Olx,sNx+Olx            DO i=1-Olx,sNx+Olx
# Line 135  C     where the local stratification exc Line 136  C     where the local stratification exc
136  C     (from surface down to here) by factor hMixCriteria  C     (from surface down to here) by factor hMixCriteria
137    
138  c       hMixCriteria  = 2. _d 0  c       hMixCriteria  = 2. _d 0
139  C-Note: dRhoSmall is hard coded for now but should become run-time parameter  C-Note: dRhoSmall is hard coded for now but should become run-time parameter
140          dRhoSmall = 1. _d -2          dRhoSmall = 1. _d -2
141          DO j=1-Oly,sNy+Oly          DO j=1-Oly,sNy+Oly
142           DO i=1-Olx,sNx+Olx           DO i=1-Olx,sNx+Olx
# Line 150  C-Note: dRhoSmall is hard coded for now Line 151  C-Note: dRhoSmall is hard coded for now
151          ENDDO          ENDDO
152          DO k = 2,Nr-1          DO k = 2,Nr-1
153  C-    potential density (reference level = surface level)  C-    potential density (reference level = surface level)
154           CALL FIND_RHO(           CALL FIND_RHO_2D(
155       I             bi, bj, 1-OLx, sNx+OLx, 1-OLy, sNy+OLy, K, 1,       I        1-OLx, sNx+OLx, 1-OLy, sNy+OLy, 1,
156       I             theta, salt,       I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
157       O             rhoLoc, myThid )       O        rhoLoc,
158         I        k, bi, bj, myThid )
159    
160           DO j=1-Oly,sNy+Oly           DO j=1-Oly,sNy+Oly
161            DO i=1-Olx,sNx+Olx            DO i=1-Olx,sNx+Olx
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22