& streamice_vert_shear_uz, streamice_vert_shear_vz
_RL streamice_vert_shear_uz (Nr)
streamice_vert_shear_uz (m) = streamice_taubx(i,j,bi,bj) /
& 0.25 * streamice_vert_shear_uz(m)**2 +
& streamice_vert_shear_uz (m) * streamice_delsigma (m) *