isopycK = GM_isopycK
Kgm_tmp = isopycK*GM_isoFac2d(i,j,bi,bj)
Kwz(i,j,k,bi,bj)= ( isopycK*GM_isoFac2d(i,j,bi,bj)
_RL Kgm_tmp, isopycK, bolus_K