FUNCTION gaunt_ff,nu,T

;nu in Hz
nu40=nu/1.e9/40   ;in units of 40 GHz
T4=T/1.e4

mean_gff=4.4*(T4)^0.21*(nu40)^(-0.14)

RETURN,mean_gff

END