The *fft are binary files with the Fourier transform of the propagator displayed in the following way PARAMETER(L=24,LT=48) parameter(Ll=6,l0=12) parameter(nmom=(2*Ll+1)**3*(2*L0+1)) !total number of momenta parameter(norb3=2100,norb4=210,maxorb=48) complex*16 prop(0:11,0:11,-Ll:Ll,-Ll:Ll,-Ll:Ll,-L0:L0) c sink source nz ny nx nt remember that the inner loops are at the left. The 0:11 of course represent the spin-color in the standard way s=0,c=0:2,s=1,c=0:2 etc In other words the inner loop is the sp-col at the sink, the the sp_col at the source, etc. The outer loop is on nt The momenta are 2 pi nx/L, 2 pi ny/L, 2 pi nz/L, pi(2nt-1)/L T