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