add T and tau cut parameters in the code and fixed GNUmakefile
add function to calculate thermal photon source emission function at a given spatial point (tau, x, y, eta_s) with photon momentum (q, phi_q, rap_q)