actually calculate the ice season

calc_ice_season(yfile, threshval = 15)