grd_below_sl = 1. _d 0
grd_below_sl = 0. _d 0
& streamice_density_ocean_avg * grd_below_sl *
& streamice_density_ocean_avg*grd_below_sl *
_RL grd_below_sl
_RL face_factor, grd_below_sl, i_rhow