_RL oceTs, oceV2s, snowPr
snowPr = snowP(i,j)
fresh = fresh + snowPr
qleft = qleft - snowPr*Lfresh
hs = hs + dt*snowPr/rhos
fresh = fresh + (prcAtm(i,j)-snowPr)