double precision Q(-nf:nf),tau(-nf:nf) common/ewcharge/Q,tau save /ewcharge/