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 |
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 |
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 |