calc_fd function

Calculate sap flux density time series