&     streamice_sigma_coord, streamice_delsigma,

      _RL streamice_delsigma (Nr)

        streamice_delsigma (m) = 1./Nr

     &       H_streamice(i,j,bi,bj) * streamice_delsigma (m) *

     &       streamice_sigma_coord(m) * streamice_delsigma(m) /

     &       streamice_vert_shear_uz (m) * streamice_delsigma (m) *

     &       streamice_vert_shear_vz (m) * streamice_delsigma (m) *

     &       streamice_delsigma(m)*(omega_temp(m)+omega_temp(m+1))*.5

     &       streamice_delsigma(m)*(u_shear(m)+u_shear(m+1))*.5

     &       streamice_delsigma(m)*(v_shear(m)+v_shear(m+1))*.5