_RL GM_skewflx
& GM_skewflx, GM_advect
Kgm_tmp = GM_isopycK + GM_skewflx*GM_background_K
& + VisbeckK(i,j,bi,bj)*(1. _d 0 + GM_skewflx)
& ( GM_isopycK - GM_skewflx*GM_background_K
GM_skewflx = 0.
GM_skewflx = 1.
CALL WRITE_0D_R8( GM_background_K*GM_skewflx, INDEX_NONE,