78 |
DO j=1,sNy |
DO j=1,sNy |
79 |
DO i=1,sNx |
DO i=1,sNx |
80 |
IF ( hOceMxL(i,j,bi,bj).NE.0. _d 0 ) THEN |
IF ( hOceMxL(i,j,bi,bj).NE.0. _d 0 ) THEN |
81 |
|
sOceMxL(i,j,bi,bj) = |
82 |
|
& (sOceMxL(i,j,bi,bj) + dtFac*sMxL_default)*dtFacR |
83 |
|
oceTfreez = -mu_Tf*sOceMxL(i,j,bi,bj) |
84 |
locTemp = ( aim_sWght0*aim_sst0(i,j,bi,bj) |
locTemp = ( aim_sWght0*aim_sst0(i,j,bi,bj) |
85 |
& + aim_sWght1*aim_sst1(i,j,bi,bj) |
& + aim_sWght1*aim_sst1(i,j,bi,bj) |
86 |
& ) - celsius2K |
& ) - celsius2K |
87 |
locTemp = MAX( locTemp , oceTfreez ) |
locTemp = MAX( locTemp , oceTfreez ) |
88 |
tOceMxL(i,j,bi,bj) = |
tOceMxL(i,j,bi,bj) = |
89 |
& (tOceMxL(i,j,bi,bj) + dtFac*locTemp)*dtFacR |
& (tOceMxL(i,j,bi,bj) + dtFac*locTemp)*dtFacR |
|
sOceMxL(i,j,bi,bj) = |
|
|
& (sOceMxL(i,j,bi,bj) + dtFac*sMxL_default)*dtFacR |
|
90 |
ENDIF |
ENDIF |
91 |
ENDDO |
ENDDO |
92 |
ENDDO |
ENDDO |