325 |
C- express the Tensor in term of Diffusivity (= m**2 / s ) |
C- express the Tensor in term of Diffusivity (= m**2 / s ) |
326 |
DO j=1-Oly+1,sNy+Oly-1 |
DO j=1-Oly+1,sNy+Oly-1 |
327 |
DO i=1-Olx+1,sNx+Olx-1 |
DO i=1-Olx+1,sNx+Olx-1 |
328 |
|
#ifdef ALLOW_KAPGM_CONTROL |
329 |
|
Kgm_tmp = GM_isopycK + GM_skewflx*kapgm(i,j,k,bi,bj) |
330 |
|
#else |
331 |
Kgm_tmp = GM_isopycK + GM_skewflx*GM_background_K |
Kgm_tmp = GM_isopycK + GM_skewflx*GM_background_K |
332 |
|
#endif |
333 |
#ifdef GM_VISBECK_VARIABLE_K |
#ifdef GM_VISBECK_VARIABLE_K |
334 |
& + VisbeckK(i,j,bi,bj)*(1. _d 0 + GM_skewflx) |
& + VisbeckK(i,j,bi,bj)*(1. _d 0 + GM_skewflx) |
335 |
#endif |
#endif |
418 |
DO j=1-Oly+1,sNy+Oly-1 |
DO j=1-Oly+1,sNy+Oly-1 |
419 |
DO i=1-Olx+1,sNx+Olx-1 |
DO i=1-Olx+1,sNx+Olx-1 |
420 |
Kuz(i,j,k,bi,bj) = |
Kuz(i,j,k,bi,bj) = |
421 |
|
#ifdef ALLOW_KAPGM_CONTROL |
422 |
|
& ( GM_isopycK - GM_skewflx*kapgm(i,j,k,bi,bj) |
423 |
|
#else |
424 |
& ( GM_isopycK - GM_skewflx*GM_background_K |
& ( GM_isopycK - GM_skewflx*GM_background_K |
425 |
|
#endif |
426 |
#ifdef GM_VISBECK_VARIABLE_K |
#ifdef GM_VISBECK_VARIABLE_K |
427 |
& +op5*(VisbeckK(i,j,bi,bj)+VisbeckK(i-1,j,bi,bj))*GM_advect |
& +op5*(VisbeckK(i,j,bi,bj)+VisbeckK(i-1,j,bi,bj))*GM_advect |
428 |
#endif |
#endif |
541 |
DO j=1-Oly+1,sNy+Oly-1 |
DO j=1-Oly+1,sNy+Oly-1 |
542 |
DO i=1-Olx+1,sNx+Olx-1 |
DO i=1-Olx+1,sNx+Olx-1 |
543 |
Kvz(i,j,k,bi,bj) = |
Kvz(i,j,k,bi,bj) = |
544 |
|
#ifdef ALLOW_KAPGM_CONTROL |
545 |
|
& ( GM_isopycK - GM_skewflx*kapgm(i,j,k,bi,bj) |
546 |
|
#else |
547 |
& ( GM_isopycK - GM_skewflx*GM_background_K |
& ( GM_isopycK - GM_skewflx*GM_background_K |
548 |
|
#endif |
549 |
#ifdef GM_VISBECK_VARIABLE_K |
#ifdef GM_VISBECK_VARIABLE_K |
550 |
& +op5*(VisbeckK(i,j,bi,bj)+VisbeckK(i,j-1,bi,bj))*GM_advect |
& +op5*(VisbeckK(i,j,bi,bj)+VisbeckK(i,j-1,bi,bj))*GM_advect |
551 |
#endif |
#endif |