include sfr.h routine to calculing star formation rate