added gamma_source.f