& wmean_hflux, wmean_sflux, wmean_tau,
wmean_sflux = 0. _d 0
& + (fctilemm(bi,bj)/wmean_sflux/nyearsrec)**2
& wmean_sflux, wsflux,
_RL wmean_sflux