In Evol/GAUGE/meas2.f
zpol=zpol/(NV*NC) -> zpol=zpol/(NV*NC/NT) at the forth last line