& (1-streamice_density_firn*i_rhow) *
& streamice_density**2*i_rhow*
& (1-streamice_density*i_rhow) *
_RL i_rhow
i_rhow = 1./streamice_density_ocean_avg
_RL face_factor, grd_below_sl, i_rhow
_RL unconf_stress, i_rhow