156 |
c ustar (nx,ny) - surface friction velocity (m/s) |
c ustar (nx,ny) - surface friction velocity (m/s) |
157 |
c bo (nx,ny) - surface turbulent buoyancy forcing (m^2/s^3) |
c bo (nx,ny) - surface turbulent buoyancy forcing (m^2/s^3) |
158 |
c bosol (nx,ny) - surface radiative buoyancy forcing (m^2/s^3) |
c bosol (nx,ny) - surface radiative buoyancy forcing (m^2/s^3) |
159 |
c boplume(nx,ny,Nr) - surface haline buoyancy forcing (m^2/s^3) |
c boplume(nx,ny,Nrp1) - surface haline buoyancy forcing (m^2/s^3) |
160 |
c shsq (nx,ny,Nr) - local velocity shear squared |
c shsq (nx,ny,Nr) - local velocity shear squared |
161 |
c at interfaces for ri_iwmix (m^2/s^2) |
c at interfaces for ri_iwmix (m^2/s^2) |
162 |
c dVsq (nx,ny,Nr) - velocity shear re surface squared |
c dVsq (nx,ny,Nr) - velocity shear re surface squared |
185 |
#ifdef ALLOW_SALT_PLUME |
#ifdef ALLOW_SALT_PLUME |
186 |
_RL temp1 ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr ) |
_RL temp1 ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr ) |
187 |
_RL temp2 ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr ) |
_RL temp2 ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr ) |
188 |
_RL boplume ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr ) |
_RL boplume ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, 0:Nr ) |
189 |
#endif /* ALLOW_SALT_PLUME */ |
#endif /* ALLOW_SALT_PLUME */ |
190 |
_RL shsq ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr ) |
_RL shsq ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr ) |
191 |
_RL dVsq ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr ) |
_RL dVsq ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr ) |
404 |
temp1(i,j,1) = saltPlumeFlux(i,j,bi,bj) |
temp1(i,j,1) = saltPlumeFlux(i,j,bi,bj) |
405 |
temp2(i,j,1) = 0. _d 0 |
temp2(i,j,1) = 0. _d 0 |
406 |
DO k=2,Nr |
DO k=2,Nr |
407 |
temp1(i,j,k) = 0. _ 0 |
temp1(i,j,k) = 0. _d 0 |
408 |
temp2(i,j,k) = 0. _ 0 |
temp2(i,j,k) = 0. _d 0 |
409 |
ENDDO |
ENDDO |
410 |
#else /* def SALT_PLUME_VOLUME */ |
#else /* def SALT_PLUME_VOLUME */ |
411 |
DO k=1,Nr |
DO k=1,Nr |