186 |
STOP 'S/R CALC_OCE_MXLAYER: invalid method' |
STOP 'S/R CALC_OCE_MXLAYER: invalid method' |
187 |
ENDIF |
ENDIF |
188 |
|
|
189 |
|
IF ( hMixSmooth .GT. 0. _d 0 ) THEN |
190 |
|
tmpFac = (1. _d 0 - hMixSmooth ) / 4. _d 0 |
191 |
|
DO j=1-Oly+1,sNy+Oly-1 |
192 |
|
DO i=1-Olx+1,sNx+Olx-1 |
193 |
|
rhoLoc(i,j)=(hMixSmooth * hMixLayer(i,j,bi,bj) + |
194 |
|
& tmpFac * ( hMixLayer(i-1,j,bi,bj) + |
195 |
|
& hMixLayer(i+1,j,bi,bj) + |
196 |
|
& hMixLayer(i,j-1,bi,bj) + |
197 |
|
& hMixLayer(i,j+1,bi,bj) ) |
198 |
|
& ) |
199 |
|
& /(hMixSmooth + |
200 |
|
& tmpFac * ( maskC(i-1,j,1,bi,bj) + |
201 |
|
& maskC(i+1,j,1,bi,bj) + |
202 |
|
& maskC(i,j-1,1,bi,bj) + |
203 |
|
& maskC(i,j+1,1,bi,bj) ) |
204 |
|
& ) * maskC(i,j,1,bi,bj) |
205 |
|
ENDDO |
206 |
|
ENDDO |
207 |
|
DO j=1-Oly+1,sNy+Oly-1 |
208 |
|
DO i=1-Olx+1,sNx+Olx-1 |
209 |
|
hMixLayer(i,j,bi,bj) = rhoLoc(i,j) |
210 |
|
ENDDO |
211 |
|
ENDDO |
212 |
|
ENDIF |
213 |
|
|
214 |
#ifdef ALLOW_DIAGNOSTICS |
#ifdef ALLOW_DIAGNOSTICS |
215 |
IF ( useDiagnostics ) THEN |
IF ( useDiagnostics ) THEN |
216 |
CALL DIAGNOSTICS_FILL( hMixLayer, 'MXLDEPTH', |
CALL DIAGNOSTICS_FILL( hMixLayer, 'MXLDEPTH', |