CALL WRITE_0D_R8( GM_background_K*GM_skewflx, INDEX_NONE,
_RL GM_skewflx
& GM_skewflx, GM_advect
GM_skewflx = 0.
GM_skewflx = 1.
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