_RL exf_outscal_hflux
& exf_outscal_hflux,
Qnet(i,j,bi,bj) = exf_outscal_hflux*hflux(i,j,bi,bj)
& exf_outscal_hflux * ( hflux_exfremo_intercept +
& exf_outscal_hflux, exf_outscal_ustress, exf_outscal_vstress,
exf_outscal_hflux = 1. _d 0