& streamice_sigma_coord, streamice_delsigma,
_RL streamice_sigma_coord (Nr)
streamice_sigma_coord (m) = 1./Nr * (m-0.5)
& * streamice_sigma_coord(m)
& streamice_sigma_coord(m) * streamice_delsigma(m) /