_RL recip_dhSnowLin
recip_dhSnowLin = 1. _d 0 / dhSnowLin
recip_dhSnowLin = 0. _d 0
frsnow = hs*recip_dhSnowLin/iceMask(i,j)